Next Article in Journal
Recent Advancements in Cyclodextrin-Based Adsorbents for the Removal of Hazardous Pollutants from Waters
Next Article in Special Issue
Structure and Diffusion of Ionic PDMS Melts
Previous Article in Journal
The Role of Roller Rotation Pattern in the Spreading Process of Polymer/Short-Fiber Composite Powder in Selective Laser Sintering
Previous Article in Special Issue
Kinetic, Isotherm, and Equilibrium Investigation of Cr(VI) Ion Adsorption on Amine-Functionalized Porous Silica Beads
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mesoscale Modeling of Agglomeration of Molecular Bottlebrushes: Focus on Conformations and Clustering Criteria

Department of Materials Science and Engineering, Clemson University, Clemson, SC 29634, USA
*
Author to whom correspondence should be addressed.
Current address: Prescience Insilico Pvt. Ltd., Bengaluru 560037, Karnataka, India.
Current address: Oakwood Labs, Oakwood, OH 44146, USA.
Polymers 2022, 14(12), 2339; https://doi.org/10.3390/polym14122339
Submission received: 3 May 2022 / Revised: 3 June 2022 / Accepted: 6 June 2022 / Published: 9 June 2022
(This article belongs to the Special Issue Molecular Simulation and Modeling of Polymers)

Abstract

:
Using dissipative particle dynamics, we characterize dynamics of aggregation of molecular bottlebrushes in solvents of various qualities by tracking the number of clusters, the size of the largest cluster, and an average aggregation number. We focus on a low volume fraction of bottlebrushes in a range of solvents and probe three different cutoff criteria to identify bottlebrushes belonging to the same cluster. We demonstrate that the cutoff criteria which depend on both the coordination number and the length of the side chain allows one to correlate the agglomeration status with the structural characteristics of bottlebrushes in solvents of various qualities. We characterize conformational changes of the bottlebrush within the agglomerates with respect to those of an isolated bottlebrush in the same solvents. The characterization of bottlebrush conformations within the agglomerates is an important step in understanding the relationship between the bottlebrush architecture and material properties. An analysis of three distinct cutoff criteria to identify bottlebrushes belonging to the same cluster introduces a framework to identify both short-lived transient and long-lived agglomerates; the same approach could be further extended to characterize agglomerates of various macromolecules with complex architectures beyond the specific bottlebrush architecture considered herein.

Graphical Abstract

1. Introduction

Molecular bottlebrushes (MBBs) are branched macromolecules with the side chains densely grafted to the linear backbone [1,2,3,4,5,6]. The structural characteristics of the bottlebrushes can be tailored during their synthesis by regulating their grafting density, degrees of polymerization of the side chains and the backbone, and chemical composition [6,7,8,9]. Due to their unique properties, MBBs can find their use in various applications, from electronic and photonic materials [10,11,12] to strain-adaptive stiffening materials [2,7] and drug delivery applications [13,14]. High grafting density leads to the pronounced effects of steric repulsions, which results in the stretched conformations of a backbone and side chains compared with the linear chain conformations. The extent of stretching, as well as equilibrium conformation of the MBB, strongly depends on its architecture, chemical nature of the moieties within the backbone and side chains, solvent quality, and concentration of MBBs within the solvent [1,6,15,16]. Relatively high backbone tension can develop depending on bottlebrush conformation [17,18]; this tension can be further amplified upon bottlebrush adsorption onto the interface [19,20,21].
For a chosen architecture and chemical nature of a bottlebrush, its conformation can be tuned either by regulating the solvent quality or by regulating the concentration of bottlebrushes in solution. The phase separation and agglomeration processes in blends and solutions incorporating MBBs—as well as migration of bottlebrushes to the interfaces and surfaces in multi-phase systems—introduce further complexity, so that the resulting behavior depends on the interplay between the contributions from the components with various affinities [22,23,24,25]. The self-assembly process can also be tailored by external stimuli, such as temperature, pH, ionic strength, and external fields [2,26,27,28,29]. For example, the conformation of the bottlebrush with side chains exhibiting lower critical solution temperature (LCST) behavior depends strongly on the architectural properties of the bottlebrush as well as on temperature [30,31,32]. At low temperatures, the side chains are hydrophilic, hence the bottlebrush attains cylindrical morphology provided that the backbones are significantly longer than the side chains [30]. An increase in temperature results in a pronounced decrease in the bottlebrush radius of gyration and in cylinder-to-globule transition at sufficiently high temperatures [30,32]. Furthermore, a two-step conformational change with an increase in temperature was reported for polymers with molecular architecture intermediate between a star-like and a cylindrical bottlebrush: an increase in temperature first resulted in stretching and only then in collapse of the macromolecules at sufficiently high temperatures [33]. The chemical nature of the end groups along with the length of the side chain also significantly affects phase behavior and interfacial properties of the bottlebrushes [31,34]. Varying the architecture of the stimuli-responsive side chains by choosing either random or block copolymers further allows one to tailor bottlebrush conformation and aggregation behavior [35]. Recent experiments showed that the bottlebrushes within aggregates formed at high temperatures (above LCST) are either loosely packed if the side chains are random copolymers or strongly interpenetrate if the side chains are block copolymers [35].
Of particular interest is using a relatively low fraction of specifically designed bottlebrushes to controllably tailor bulk and surface properties of various systems. Bottlebrushes can be effectively used as compatibilizers to control morphologies and characteristic domain sizes of immiscible polymer blends; specifically, MBBs were shown to encapsulate one of the components of the polymer blend, thereby suppressing coarsening at elevated temperatures [36,37]. Bottlebrushes were also shown to encapsulate proteins such as lysozyme, preventing its denaturation at high temperatures [38]. Tailoring architecture and composition of MBBs allows one to form micelles with controllable sizes [14,39] and imparts the capability to carry drugs with varied hydrophobicity [40,41].
In our recent study [42], oleophobic perfluoropolyether (PFPE) methacrylate molecular bottlebrush (PFM) was synthesized and used as an additive to commodity thermoplastics (nylon 6 and poly(methyl methacrylate)). Our studies showed that an addition of only 1% to 5% of the PFM bottlebrush additives imparts high water and oil repellency to the host thermoplastic matrix. While similar oleophobic properties can be achieved via utilizing long-chain perfluoroalkyl substances (LCPFAS), these materials have currently been phased out due to their negative environmental impact [43,44,45]. PFM bottlebrush additives introduce a feasible and potentially safer alternative to LCPFAS due to the lower toxicity and high oxidative and thermal stability of perfluoropolyethers [46,47,48,49]. We showed that the observed high oil repellency of the surfaces is due to the migration of PFM bottlebrush additives to the thermoplastics surfaces due to the mismatch in the affinity between the densely grafted side chains and a host polymer [42]. Specifically, we showed that first the bottlebrushes aggregate into the clusters with various number of brushes depending on the affinity between the moieties comprising the backbone/side chain and the host matrix. Upon approaching the interface, the bottlebrush clusters spread over the interface exposing highly oleophobic side chains [42].
The characterization of bottlebrush conformations is an important step in understanding the relationship between the bottlebrush architecture and material properties [50]. A number of techniques is available for the direct visualization of the conformation of MBBs in experiments, including atomic force microscopy [7,51], transmission electron microscopy, small-angle neutron scattering [52], and single-molecule localization microscopy [53]. However, due to the limited resolution of available experimental techniques [53,54], it remains challenging to distinguish the bottlebrushes within clusters and characterize their conformations. Computational modeling allows one to understand and predict bottlebrush conformations under various conditions [5,22,23,51,52,55,56,57,58,59,60,61,62,63,64]; current modeling efforts focusing on characterizing self-assembly of bottlebrushes are surveyed in a recent review by Mohmmadi et al. [65].
Of particular relevance to the current work are recent studies characterizing conformations of MBBs in solvents of various quality. Recent integrated computational and experimental study quantified the effects of the backbone and side chain lengths on the conformations of bottlebrushes in a good solvent for a range of architectural parameters of bottlebrushes [66]. Specifically, the variation in the backbone molecular weight for the fixed length of the side chain was shown to significantly affect the resulting MBB conformation, ranging from a star-like structure for a short backbone, to an extended rod, and then to a coiled chain with an increase in backbone length; importantly, these changes were reflected in intrinsic viscosity measured experimentally [66]. Furthermore, the decrease in solvent quality from good to poor solvent for a fixed bottlebrush architecture was shown to induce transitions from the coil regime to the dense collapsed state [23,55,57,60]. This transition is dependent upon the architectural parameters, such as the length of the side chain and grafting density [55,57]. In addition, the conformation of the bottlebrush also reflects variations in their concentration with scaling relations affected by the screening of excluded volume interactions [58,59].
For multiple bottlebrushes with various affinities between the backbone and solvent, the bottlebrushes agglomerate into clusters of various shapes—including spherical, cylindrical, or necklace-type structures—depending on the side chain grafting densities and chemical nature of the monomers [67]. For bottlebrush block copolymer with two-component side chains, the phase separation behavior in solutions was quantified with respect to the arrangement and relative fraction of the two components [22,56,61]. Recently, Gumerov and Potemkin [23] utilized dissipative particle dynamics (DPD) simulations to characterize conformations of a single MBB with responsive inner blocks of the side chains, as well as self-assembly of these MBBs, in solvents of various quality. They demonstrated that MBBs form continuous or branched aggregates or cylindrical micelles depending upon solvent quality and grafting density; an increase in grafting density was shown to change the aggregate shape and required a lower solvent quality to induce an aggregation of the MBBs [23]. Chang et al. [68] demonstrated a transition of micellar aggregates into unilamellar vesicles and then to “onion-like” multilayered vesicles with the decrease in the affinities between the side chain and solvent. This morphology was observed for the relatively short MBB molecules [68,69] and hyperbranched multi-arm copolymers [70,71].
Herein, we focus on characterizing the effects of the solvent quality on the average size of agglomerates in the system with low volume fraction of bottlebrushes in solvents of various quality with respect to the side chains. We characterize the changes in the conformation of the bottlebrush within the agglomerate with respect to that in the limit of infinite dilution (single bottlebrush). In what follows, we first introduce our modeling approach and the choice of model parameters; this choice is informed by our prior experimental results [42] providing experimentally relevant architectural characteristics of the bottlebrush. We then characterize radius of gyration, shape anisotropy, and an acylindricity of a single bottlebrush in solvents of various quality and compare these characteristics with the respective values calculated during the agglomeration process. We track the agglomeration process by tracking the time evolution of cluster distribution, the size of the largest cluster, and an average aggregation number. We probe three different cutoff criteria to identify whether two bottlebrushes belong to the same cluster. The results below show that the cutoff criterion which depends on both the coordination number and the length of the side chain allows one to correlate the agglomeration status with the structural characteristics of the bottlebrushes in solvents of various qualities.

2. Methods

We use the dissipative particle dynamics (DPD) approach [72,73,74,75] to model the time evolution and self-assembly in bottlebrush-solvent systems. DPD is a computationally efficient mesoscale approach that has been utilized to model a broad variety of polymer systems [76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92,93,94], including studies of the effects of solvent quality on structural characteristics [95,96,97,98] and self-assembly in various polymer systems [23,79,89]. Herein, we briefly introduce the DPD approach, the details of this approach can be found in [72,73,74]. The beads in DPD represent groups of atoms with the motion of these beads governed by the Newton’s equations of motion [74]. There are three contributions to the pairwise additive force F i j exerted on a bead i by a bead j : conservative ( F i j C ), dissipative ( F i j D ) , and random ( F i j R ) force. The total force acting on a bead i, F t o t = j i F i j , is calculated via summation over all forces acting on this bead from the beads located within the characteristic cutoff distance r c ; the value of r c effectively defines the intrinsic length scale of the model. We chose the conservative force corresponding to the soft repulsion between the beads, which is the most common choice of the conservative force in DPD, as [74]
F i j C = { a i j ( 1 r i j r c ) e i j 0   ( r i j < r c ) ( r i j r c )
where a i j is the repulsion coefficient that defines the magnitude of maximum repulsion between the beads i and j , r i j = | r i j | is the distance between the centers of masses of these beads, r i j = r i r j , and e i j = r i j r i j . The repulsion coefficient between the same types of beads, a i i is often chosen based on the degree of coarse-graining [74,99], while the repulsion parameter between the beads iand j representing chemically distinct species can be chosen based on the affinities between the species these beads, represented as [74]
a i j = a i i + 3.27 χ i j
where χ i j is the Flory–Huggins interaction parameter between the respective species; the above expression holds for the bead number density ρ = 3 , which is the common choice on the number density in DPD.
The dissipative and random forces are coupled to satisfy the fluctuation-dissipation theorem and are taken as [74] F i j D = γ ω 2 ( r i j ) ( e i j · v i j ) e i j , and F i j R = σ ω ( r i j ) ζ i j · t 1 / 2 e i j , respectively, where γ and σ are the strengths of the dissipative and random forces, v i j = v i v j is the relative bead velocity, Δ t is the simulation time step, and ζ i j   is a symmetric Gaussian random variable with unit variance and zero mean. We use the most common choice of the weight function [74,99] ω ( r i j ) = 1 r i j r c ; other choices of the weight function are also permitted [100] within the DPD framework.
The bonded beads are connected by harmonic bonds with the interaction potential U b o n d = K b o n d 2 ( r i j r b ) 2 , where K b o n d is a spring constant and r b is an equilibrium bond length. To minimize topological crossings of bonded polymer chains, a known limitation of DPD, we utilize a modified segmental repulsive potential (mSRP) formulation [101]. Within the mSRP formulation, an additional repulsive force
F i j S R P = b ( 1 d i j / d c S ) e i j S
is introduced between the centers of the bonds if the distance between these centers, d i j = | d i j | is below the cutoff distance [101] d c S ; herein, e i j S = d i j d i j .
The reference parameters in our DPD simulations are chosen as following [74,99]: the beads number density is ρ = 3 , the cutoff distance is r c = 1 , the strengths of the dissipative and random forces are γ = 4.5 and σ = 3.0, respectively, and the temperature and bead mass are set to unity in reduced DPD units. For the bonded interactions, we set [87] K b = 10 3   and r b = 0.7 ; notably, high spring constant values were used in a number of recent DPD studies [102,103,104]. The parameters defining mSRP interactions are set as [101] b = 80 and d c = 0.8 . These mSRP parameters were shown to effectively minimize topology violations in the original publication [101] and in the subsequent studies [105]. For the repulsion parameter between the same type of beads, we choose a i i = 78 in reduced DPD units, this value is calculated to reproduce water compressibility and corresponds to a single bead representing three water molecules [99]. This choice of the degree of coarse-graining allows one to relate the reduced DPD length scale, r c = 1 , to a dimensional value of r c     0.65   nm , and the reduced unit of time, τ , to the dimensional value of τ   88   p s [99]. The choice of the repulsion parameters between the beads representing chemically distinct moieties ( a i j ,   where   i j )   is discussed further below. The parameters above and all the simulation values below are provided in reduced DPD units [74,99]. The equations of motion are integrated using the velocity-Verlet algorithm with the time step Δt = 0.01 τ . The simulations were carried out using the Large Scale Atomic Molecular Massively Parallel Simulator (LAMMPS) package [106,107] with mSRP code [101]. The simulation snapshots are visualized using the Visual Molecular Dynamics (VMD) software [108].
The bottlebrush architecture in the reference case is shown in Figure 1a, where the beads representing the backbone are shown in green and the beads within the side chains are shown in blue. The number of backbone and side chain beads are denoted as Nbb = 62 and Nsc = 6, respectively, and the grafting density is set to z = 1. The grafting density equal to one and the ratio between the number of beads within the backbone and within the side chain is chosen based on the experimental data [42]. We model the cases where the backbone beads have high affinity with the solvent beads (abb-sol = 78) and low affinity with the side chain beads (abb-sc = 89). The choice of these parameters based on the measured interfacial tensions for the examples of molecular bottlebrushes synthesized in the concurrent experiments is provided in our recent study [42]. In what follows, we vary the affinities between side chain and solvents (asc-sol) within the relatively wide range.
The box size in all the simulations below involving multiple bottlebrushes is chosen as 60 × 60 × 60. The periodic boundary conditions are applied in all directions. The fraction of bottlebrush beads is set at 5% of total beads within the system. For the reference architecture shown in Figure 1a, the bottlebrush consists of N T = 422 beads. nMBB = 77 bottlebrushes are dispersed within the simulation box of reference size (Figure 1b). In the simulations below, we also vary the bottlebrush architecture keeping the grafting density fixed but varying the degree of polymerization of the side chain (Nsc = 10 and Nsc = 20); in these series of simulations nMBB is varied correspondingly to keep the fraction of bottlebrushes fixed at 5%. Initially, we set all the repulsion parameters equal to aij = 78 between all the beads and equilibrate the system for 2 × 106 time steps to ensure that at the initial time for the production runs the bottlebrushes are well dispersed within the box as shown in Figure 1b. For the production runs, we use 3 × 106 time steps. In all the studies involving multiple bottlebrushes, four independent simulations were analyzed. In all the simulations involving single bottlebrushes, the simulation box size is chosen as 30 × 30 × 30 and the last 1200 frames are used for the analysis.
To characterize the conformations of the bottlebrushes dependent upon solvent quality, we calculate the following characteristics for each bottlebrush at a chosen simulation time step t: the squared radius of gyration, R g 2 , shape anisotropy, κ 2 , and acylindricity, c . The squared radius of gyration of each bottlebrush within the system is calculated as a first invariant of a diagonalized gyration tensor S ^ as
R g 2 = I 1 = t r   [ S ^ ]   = λ 1 + λ 2 + λ 3 ,  
where λ 3 λ 2 λ 1 are the eigenvalues of the gyration tensor with the components S x y = 1 N T 2 Σ i = 1 N T Σ j > i N T ( x i x j ) ( y i y j ) ; herein the summation is taken over all beads within the bottlebrush, N T . A relative shape anisotropy of a bottlebrush is defined via first ( I 1 ) and second invariants of gyration tensor, I 2 = λ 1 λ 2 + λ 2 λ 3 + λ 3 λ 1 , as [109]
κ 2 = 1 3 I 2 I 1 2
Shape anisotropy is used to characterize shapes of objects with complex geometries such as polymers with various architectures [109,110,111,112]. The relative shape anisotropy ranges between κ 2 = 0 for an object with high symmetry and κ 2 = 1 for the points located on a line [109]. For a planar object [109] (with λ 1 = 0   and λ 2 = λ 3 ) κ 2 = 0.25 . For linear polymer chains, the shape anisotropy is calculated as κ 2 0.43 in a good solvent and κ 2 0.39   in the theta solvent [113,114,115].
Finally, we also calculate acylindricity
c = λ 2 λ 1
and the ratio between the largest and smallest eigenvalues for each bottlebrush λ 3 / λ 1 .
To analyze all the above structural characteristics upon equilibration, the data are averaged within the last 6 × 106 time steps (using the last 1200 frames) for the case of a single bottlebrush. For the cases of multiple bottlebrushes, the data are averaged over all the bottlebrushes, nMBB, (nMBB = 77 in the reference case) within 60 frames (last 3 × 105 time steps) in four independent simulation sets.

3. Results and Discussions

3.1. Characterizing Conformation of a Single Bottlebrush Dependent upon Solvent Quality

We first characterized equilibrium conformation of a single bottlebrush with a chosen architecture dependent upon solvent quality. Specifically, we calculated the radius of gyration, relative shape anisotropy, eigenvalue ratio, and acylindricity of a single bottlebrush (Equations (4)–(6)) at late times upon reaching an equilibrated conformation in solvents of various qualities. These structural characteristics are then compared with the same structural characteristics of bottlebrushes within the agglomerates formed during the self-assembly of multiple bottlebrushes (Section 3.2 below) to identify the effect of agglomeration on these characteristics. We vary the repulsion parameter asc-sol between the side chain and solvent beads to tailor the side chain affinity to the solvent. The side chain–solvent repulsion parameters are chosen in independent series of simulations as asc-sol = {78,80,82,83,84,85,87,89}, where asc-sol = 78 represents a high affinity corresponding to a good solvent and asc-sol = 89 indicates a low affinity corresponding to a poor solvent. Figure 2a shows representative late time snapshots of a bottlebrush in solvents of different qualities, solvent beads are not shown for clarity of representation. Upon the decrease in the affinity of the side chain and solvent beads, the transition from the random coil to the globule conformation is clearly observed by comparing bottlebrush conformations I–IV in Figure 2a.
The collapse of the bottlebrush upon decreasing the side chains’ solubility in the solvent for various lengths of the side chain Nsc = 1, 6, 10, 20 is evident from the decrease in the average radius of gyration upon equilibration for all the Nsc considered (Figure 2b). This collapse is observed at conditions close to the theta condition for the linear chains interacting with solvent with asc-sol = 83 (see Figure S1 of SI). While for the shortest side chain (Nsc = 1) this decrease is gradual, the effect becomes more pronounced with an increase in the side chain length, Nsc. This effect can be attributed to both, an increase in the relative fraction of solvophobic beads at higher Nsc, and more pronounced backbone stretching at this high grafting density for the longer side chains. The effect of the backbone stretching is also apparent in Figure S1 of SI, where the average end-to-end distance of the entire bottlebrush, <Re>, and an average radius of gyration of the backbone are shown to exhibit the same trends with varying solvent quality. Similar behavior as shown in Figure 2b was previously reported in a coarse-grained MD study [57] of brushes with solvophobic backbones (recall that the backbone in our study is solvophilic), which showed a coil-to-globule transition close to the theta condition for the linear chain upon reducing side chain affinity with more pronounced effect for the longer side chains.
To gain further insights into the conformation of a bottlebrush in solvents of various quality with respect to the side chains, we next characterize the relative shape anisotropy, κ 2 ,   acylindricity c , and the ratio between the largest and smallest eigenvalues, λ 3 / λ 1 , for the same series of simulations (Figure 3). For all the cases considered, the values of κ 2 decrease with the decrease in the side chain affinity to the solvent (Figure 3a). When the side chains are in a good solvent (asc-sol = 78), for the bottlebrushes with Nsc = 1, 6, and 10 the calculated average values of the shape anisotropy κ 2 are close to 0.43, the values previously reported for linear chains in a good solvent [113,114] (see the representative snapshot in Figure 3d(I)). For the longest side chain length considered, Nsc = 20, the average value of κ 2 is somewhat lower than that for all the shorter side chains. This difference can be attributed to the relatively low ratio between the backbone and the side chain degrees of polymerization (see the snapshot in Figure 3d(V)). Correspondingly, the average value of κ 2 in this case is somewhat closer to that for the star polymer [116]. For the bottlebrushes in poor solvent with respect to the side chains (asc-sol = 89), the shape anisotropy κ 2 is close to zero for Nsc = 6, 10, and 20, indicating formation of highly symmetric globules for all these cases (Figure 3d(II,VI)). For the shortest side chain length of Nsc = 1, the average κ 2 is notably higher in poor solvent than that for all the other architectures, but it still remains relatively low ( < κ 2 >   0.1). The latter indicates collapse of the bottlebrush into a structure with relatively low symmetry (see the snapshot in Figure 3d(IV)).
Same trends are further confirmed by comparing an average ratio between the largest and smallest eigenvalues, < λ 3 / λ 1 > (Figure 3b), and an average acylindricity normalized on the smallest eigenvalue, c / λ 1 (Figure 3c). The normalization on the λ 1 in the latter plot is needed for more direct comparison between the structural characteristics of macromolecules with highly disparate dimensions; the acylindricity c , as well as λ 2 / λ 1 ratios are provided in Figure S2. These plots confirm that the bottlebrushes remain highly globular for Nsc = 6, 10, and 20 at high asc-sol while for Nsc = 1 both < λ 3 / λ 1 > and c / λ 1 remain relatively high even at the highest asc-sol considered. For the solvents with high affinity to the side chains, the values c / λ 1 decrease with an increase in Nsc, indicating the lowest acylindricity normalized per smallest linear size in the system for the macromolecules with the longest side chains. We also note large error bars in the calculated values in Figure 3a–c for all the cases except the cases corresponding to nearly globular conformations. These large error bars indicate a broad distribution of macromolecular conformations upon equilibration and are on the same scale as observed in the prior studies of shape anisotropy of various polymers [114,117].

3.2. Characterizing Bottlebrush Conformations in Solutions

In the next series of simulations, we focus on the dynamics of the bottlebrushes during the self-assembly in solvents of various qualities. Initially, the bottlebrushes are well dispersed and equilibrated as detailed in the Methods section. Representative simulation snapshots of time evolution of 5% solution of bottlebrushes (nMBB = 77) in a good solvent (asc-sol = 78, Case I), in solvents of intermediate quality (asc-sol = 83 and asc-sol = 84, Cases II and III, respectively), and in a poor solvent (asc-sol = 89, case IV) are shown in Figure 4. The bottlebrushes with the side chains with high affinity to the solvent (asc-sol = 78, Figure 4a) remain well dispersed throughout the entire simulation run. For the intermediate scenarios, the bottlebrushes begin to form clusters but remain mainly dispersed for asc-sol = 83 (Figure 4b), and are mainly aggregated for asc-sol = 84 (Figure 4c). In the poor solvents (Figure 4d), the bottlebrushes first form small clusters relatively quickly (see snapshot at t = 105) and then these smaller clusters coalesce forming larger close-to-spherical aggregates (remaining snapshots in Figure 4d).
To characterize the self-assembly of bottlebrushes into aggregates under various conditions quantitatively, we begin with characterizing clusters formation with time for all the above cases. Since our scenarios encompass a range of systems from well-dispersed to agglomerated bottlebrushes, a clear protocol defining a cluster that is suitable for a range of scenarios needs to be identified. The intuitively simple Stillinger criterion [118], which identifies two monomers as belonging to the same cluster if they are separated by a distance smaller than a chosen cutoff distance, is successfully used in a number of studies including studies involving molecular modeling of the self-assembly of bottlebrushes [61,119]. Depending on the system considered, the Stillinger criterion may however strongly overestimate the size of small molecular clusters [120,121,122], counting passing molecules within the given cutoff distance. In the later scenarios, the ten Wolde–Frenkel cluster definition [123], which in addition to the distance criterion accounts for a prescribed number of neighbors, can be utilized [120,121,122]. In DPD simulations, the cutoff distance, which determines whether two beads belong to the same cluster is often chosen [23,124,125] as the model cutoff distance rc defined above. Hence, counting two bottlebrushes as part of the same cluster provided that there is at least a single contact (nc = 1) within the cutoff distance rc between the pair of beads belonging to each bottlebrush is an equivalent of the Stillinger criterion for our system. However, it is reasonable to anticipate that using a single contact would result in strong overestimation of the cluster sizes in the cases when the bottlebrushes are well dispersed or even partially dispersed due to counting short-lived (or passing by) contacts. Hence, we introduce two more criteria for identifying whether two bottlebrushes belong to the same cluster and compare the utility and robustness of all three criteria.
Since all interactions in DPD are considered only within the model cutoff distance rc, an extension of the clustering criteria based on the single contact defined above needs to incorporate a specific required number of neighbors within this distance. The minimum required number of neighbors in turn can be defined based on the coordination number calculated up to the cutoff distance rc. The coordination number of a single bead within the entire system calculated using radial distribution function up to the cutoff distance rc is 12.52. Hence, an additional criterion to identify whether two bottlebrushes belong to the same cluster can be introduced by requiring that at least one bead within one bottlebrush is effectively fully surrounded by the beads of another bottlebrush. This criterion corresponds to the minimum number of contacts required to consider both bottlebrushes belonging to the same cluster to be based on the coordination number, namely nc = 13 contacts between the beads belonging to two different bottlebrushes. Finally, we define another, an even more strict criterion that effectively requires all the beads within the single side chain of one bottlebrush to be fully surrounded by the beads of another bottlebrush for these two brushes to belong to the same cluster. This number is based on the coordination number of a single bead multiplied by the degree of the polymerization of the side chain, nc = 13 Nsc (or nc = 78 in the reference case scenario for Nsc = 6).
Using all three criteria introduced above (i.e., requiring at least nc = 1, 13, or 13 Nsc contacts between the beads of two bottlebrushes for them to belong to the same cluster), we first tracked the distribution of clusters during the time evolution in Cases I–IV. The average size of the largest cluster (the number of the bottlebrushes that belong to this cluster) and the distribution of clusters in selected cases are shown in Figure 5 and Figure S3 respectively. For the well-dispersed Case I (snapshots in Figure 4a), using the single contact cutoff criterion (nc = 1) leads to extensive variations in the size of the largest cluster (from 5 to 46 bottlebrushes, black dotted line in Figure 5a). Hence, as anticipated, this criterion results in significant overestimation of the size of the largest cluster due to counting passing-by brushes as short-lived clusters. With nc = 13 (red dash-dotted line in Figure 5a) and nc = 78 (blue solid line), the calculated sizes of the largest cluster are comparable. While some variations in the largest cluster size are observed with nc = 13 (in some cases counting a small number of transient agglomerates of 2–3 bottlebrushes), for nc = 78 the largest cluster consists of a single bottlebrush during the entire simulation run with no transient agglomerates being accounted for (Figure 5a, black solid line). Similar trends are observed for Case II (Figure 5b); however, the size of the largest cluster correspondingly increases to approximately 25 bottlebrushes on average for nc = 13 (red dash-dotted line). It should be noted that the size of the largest cluster strongly varies with time from 2 to 40 bottlebrushes indicating that the short-lived contacts between the neighboring macromolecules are accounted for while this criterion is used. With nc = 78, the largest cluster size for Case II is around 3 (blue solid line in Figure 5b) indicating a transient nature of visually larger agglomerates (intermediate to late time snapshots in Figure 4b). Similar trends distinguishing various criteria are observed at early times in Case III (Figure 5c), while at late times all the criteria give the same average size of the largest clusters. Finally, the size of the largest cluster increases in a stepwise manner in the Case IV (Figure 5d); notably, all three criteria result in the same largest cluster size. The stepwise largest cluster size increase clearly indicates the non-transient nature of the agglomeration process: namely, if a given bottlebrush joined the largest cluster, it remains a member of this cluster for the entire simulation run. In other words, no contribution to the size of the largest cluster from the transient short-lived aggregates is observed in Case IV.
The distribution of all cluster sizes with time for the time instances corresponding to the snapshots in Figure 4a,d (Case I and Case IV) are shown in Figure S3. These plots further illustrate that while with nc = 1, a large number of transient short-lived agglomerates are accounted for in the good solvent (Figure S3a); with nc = 78, none of these transient agglomerates are accounted for (Figure S3b). On the contrary, for the poor solvent (Case IV) all criteria give essentially identical clusters distribution (Figure S3c,d).
Next, we calculate an average number of clusters and the aggregation number defined as [23,56]
N a g g =   n i P i /   n i ,  
where Pi is the size of the cluster i (number of bottlebrushes within this cluster) and ni is the number of clusters of this size; the summation is taken over all the clusters. The time evolution data are summarized in Figure 6 for Cases I–IV (in black, red, blue, and green, respectively) and for all the above cutoff criteria (nc = 1, 13, and 78 are shown by circles, squares, and solid lines of the respective color as marked in the legend). The data points are averaged over four independent simulations, the error bars represent standard deviations. We first focus on the solid curves for all four cases (nc = 78). The dispersed bottlebrushes (Case I) remain well dispersed in all the runs, solid black curves correspond to the number of bottlebrushes, nMBB = 77, (Figure 6a) and respectively N a g g = 1   (Figure 6b). In other words, using the criteria nc = 78 identifies each bottlebrush as an independent cluster. However, if one uses a single contact cutoff (nc = 1, black circles in Figure 6), the number of clusters is significantly lower and varies extensively with time due to the formation of transient short-lived agglomerates as discussed above for a single simulation run.
Upon agglomeration (Cases II–IV), the bottlebrushes form larger clusters hence the average number of clusters decreases with time (solid red, blue, and green curves in Figure 6a, correspondingly). For Case II, a number of clusters calculated using nc = 78 remains just somewhat lower than the total number of bottlebrushes nMBB (Figure 6a, solid red curve). Using cutoff criteria nc = 1 or 13 and following the time evolution for N a g g for Case II (red symbols in Figure 6b) further indicates formation of primarily short-lived agglomerates in this scenario.
The agglomeration process is distinctly different in Cases III and IV. The process is notably faster for Case IV (green curve, solid line for nc = 78) than that for Case III (in blue, solid line for nc = 78) since the driving force for the agglomeration due to the differences in affinities of the side chains and solvent is the strongest in Case IV. Furthermore, an average aggregation number—which reflects an average size of the cluster—changes approximately stepwise for the average   N a g g (and strictly stepwise for each individual run) for the Case IV, confirming the non-reversibility of the agglomeration process noted above (i.e., a particular bottlebrush joining a cluster belongs to this cluster for the remainder of the simulation run). On the contrary, large absolute values of N a g g and extensive variations in these values with time are observed for Case III (blue solid line in Figure 6b), indicating that even using nc = 78 criteria accounts for the notable contribution of short-lived relatively large agglomerates since the transient agglomeration is an inherent feature of the aggregation process for this solvent of intermediate quality.
To identify the applicability and the utility of the chosen criteria, we next analyze an average number of clusters (Figure 7a) and an average aggregation number (Figure 7b) for all three cutoff numbers of contacts as a function of solvent quality, asc-sol. Dotted, dash-dotted, and solid curves correspond to nc = 1, 13, and 78, respectively. These plots clearly show that for the relatively low solvent quality with respect to the side chains, asc-sol  85 , all the above cutoff criteria lead to identical results.
For the relatively good solvent quality (78 asc-sol 82), using nc = 1 results in a significant contribution of transient agglomerates into a cluster count, while nc = 13 already allows one to disregard the majority of short-lived transient clusters, with only small differences in measured values of N a g g   observed while using nc = 13 or nc = 78 (Figure 7b). Recall that, when using nc = 13 contacts as a cutoff number of contacts for the two bottlebrushes to belong to a single cluster ensures, effectively at least one bead from a chosen bottlebrush is fully surrounded by the beads of another bottlebrush; while using nc =78 contacts ensures that effectively all the beads within the single side chain of one bottlebrush are fully surrounded by the beads of another bottlebrush. Using the criterion nc =78 results in a number of clusters equal to the number of bottlebrushes nMBB in the system (Figure 7a); hence, no passing by bottlebrushes is accounted for if this criterion is used.
For the solvents of intermediate quality (82  asc-sol  85), transient short-lived clusters largely contribute to the dynamics of the self-assembly. Within the transient (short-lived) cluster, a bottlebrush that is at a given time instant belongs to this cluster may migrate away from it and become a member of another cluster. As discussed above, this transient agglomeration is an inherent feature of the aggregation process for the solvents of intermediate quality. Hence, characterizing the agglomeration process using various cutoff criteria provides further insights into the relative contributions of transient aggregation and formation of the long-lived clusters.
For the same series of the simulations as introduced above (Cases I–IV, Figure 4), we next characterize conformations of bottlebrushes during the assembly process as a function of time by calculating the radius of gyration Rg and the shape anisotropy κ 2 (Figure 8). In Case I (black curves), both Rg and κ 2 remain approximately constant during all the simulation runs with the values close to those for a single bottlebrush in a good solvent with respect to the side chains. Both these values notably decrease at a relatively early time for Case II (red curves) and maintain these lower values (with respect to Case I) for the remainder of the simulation runs. The equilibrated values of Rg and κ 2 in Case II remain close to—but somewhat lower than—the respective values for a single bottlebrush (see direct comparison below).
Unlike the monotonous decrease and saturation of Rg and κ 2 in Case II, a shallow dip in both these characteristics is observed in Case III (in blue) and a more pronounced dip is observed in Case IV (in green) at relatively early times prior to the saturation of both characteristics at their equilibrium values. The sharp decrease in both Rg and κ 2 in Case IV clearly reflects formation of small nearly globular structures containing only a few bottlebrushes (as evident in Figure 4d at early times). Notably, κ 2 drops to 0.1 at early times (in green Figure 8b), while this value is low it is notably higher than κ 2 corresponding to globular structures formed by a single bottlebrush (Figure 3a). An increase in κ 2   with time up to 0.35 and a corresponding increase in Rg indicate the conformational changes during the agglomeration process. While the long-lived clusters are formed at late times in both Cases III and IV, our results indicate differences in aggregation processes in these two scenarios. In Case III, a number of short-lived agglomerates are formed and nearly immediately broken-up (as can also be seen from clusters characterization in Figure 6a,b). In other words, relatively higher affinity between the solvent and the side chains in Case III compared to that in Case IV promotes short-lived clusters and hinders formation of long-lived clusters. Furthermore, upon equilibration κ 2 saturates at approximately the same value in Cases III and IV (Figure 8b), indicating similar shapes of the bottlebrushes within the clusters. The average values of Rg however remain higher upon equilibration in Case III than that in Case IV, indicating that the bottlebrushes are somewhat more stretched in Case III.

3.3. Characterizing Equilibrium Conformations of Bottlebrushes upon Agglomeration

We now focus on directly quantifying the effects of the agglomeration of bottlebrushes in various solvents on their equilibrium conformations. We first compare equilibrium radius of gyration of bottlebrushes in solutions containing 5% of bottlebrushes (reference case architecture with Nsc = 6 shown in Figure 1a) with that for a single bottlebrush (Figure 9a). As anticipated from the above discussion of the time evolution of structural properties, the value of Rg remains the same as that for a single bottlebrush when the bottlebrushes are well-dispersed (asc-sol = 78, Case I in Figure 4a). With the gradual increase in asc-sol (a decrease in solvent quality with respect to the side chains), Rg decreases similarly to the behavior observed for a single bottlebrush. This decrease is however significantly less pronounced for multiple bottlebrushes (black curve in Figure 9a). The radius of gyration for the multiple bottlebrushes is notably larger than that for a single bottlebrush even for the intermediate solvent quality, for which a large number of transient agglomerates is formed (for example, asc-sol = 83 in Case II discussed above). Furthermore, an average radius of gyration Rg is found to be significantly higher than that for a single bottlebrush in all the cases where the self-assembly into the stable clusters at late times is observed (asc-sol = [84:89], where asc-sol = 84 and asc-sol =89 are the Cases III and IV analyzed above).
It is instructive to characterize the distribution of the values of Rg in the case of multiple bottlebrushes. For the dispersed case (asc-sol = 78), the distribution of Rg is symmetric (Figure 9c) with the peak at the average value reported in Figure 9a. For the case of the poor solvent (asc-sol = 89), the peak is clearly shifted to the lower value (Figure 9d) and is notably lower than the average Rg (Figure 9a). The minimum value in this distribution Rg     3.0 corresponds to that observed for a single bottlebrush in the same solvent (i.e., globular configuration of a single bottlebrush), while the highest value at the end of the long tail (Rg    6.2) corresponds to the value for the same chain in a coiled conformation. This asymmetric distribution of the frequencies of occurrences of various values of Rg clearly illustrates broad distribution of chain conformations for the chains belonging to bottlebrush aggregates. Specifically, while a large number of bottlebrushes attain conformations close to globular (but remain partially stretched), a notable fraction also remains in close-to-coiled conformations.
Further details on the configurations of the bottlebrushes in these systems are revealed via analysis of the shape descriptors. The average shape anisotropy,   κ 2 , attains significantly higher values for all the cases corresponding to the formation of relatively large agglomerates than that for a single bottlebrush, as can be seen by comparing red and black curves in Figure 10a for asc-sol  83 . These higher values of κ 2 indicate that the shape of a bottlebrush belonging to an agglomerate, on average, remains somewhat closer to that for the polymer chains in a theta-solvent (recall that κ 2 0.39   for   the   liner   chain   in the theta solvent [113,114,115]) than to the chains in globular conformation. For asc-sol  85 , a single bottlebrush clearly attains globular conformation (low κ 2 , red curve in Figure 10a). Consistent with the distribution of the probability of occurrences of Rg (Figure 9d), the probability of occurrences of κ 2 is also broadly distributed and shifted to the left (Figure 10c) for the case of poor solvent with respect to the side chains, while the distribution of κ 2 remains significantly more symmetric, though also broadly distributed for asc-sol = 78 (Figure 10b). The broad distribution of κ 2 results in the large error bars in Figure 10a. Finally, backbone stretching for the bottlebrushes belonging to the agglomerates is also evident from the relatively high average values of < λ 3 / λ 1 > at asc-sol  85 compared to the same ratio for the single bottlebrushes (compare black and red curves in Figure 10d). Large error bars in Figure 10d at high asc-sol for the case of multiple bottlebrushes correspond to the broad distribution of their conformations as discussed above (Figure 9d and Figure 10c). To summarize, relating the clustering results using three criteria (Figure 7) and structural characteristics of the bottlebrushes in solutions (with respect to the same characteristics of the single bottlebrush, Figure 9 and Figure 10) show that using nc = 78 allows one to relate bottlebrush conformation in solution for a wide range of solvent qualities to the corresponding agglomeration status. These results also identify a range of solvent qualities (relatively poor solvents) for which any of the criteria may be used to characterize clustering.

3.4. Effect of the Side Chain Length on the Agglomeration and Bottlebrush Conformation within Agglomerate

To examine the effects of the bottlebrush architecture on the agglomeration process, multiple bottlebrushes constituting 5% of total volume fractions are placed in a good solvent (Case I, asc-sol = 78), solvents of intermediate quality (asc-sol = 83, asc-sol = 84 for Cases II and III, respectively), and poor solvent (asc-sol = 89, Case IV). In these series of simulations, the longer side chains are chosen with respect to that in the reference case, specifically Nsc = 10 and Nsc =20. To keep the total number of bottlebrush beads at 5% of the total number of beads, the numbers of bottlebrushes are reduced with respect to that in the reference case to nMBB = 49 for Nsc = 10 (NT = 662) and nMBB = 26 for Nsc = 20 (NT = 1262). The representative time evolution simulation snapshots are shown in Figure S4 and Figure S5, respectively.
We track the time evolution of the number of clusters, aggregation number, radius of gyration Rg, and shape anisotropy κ 2 to characterize the agglomeration process using the same procedures as introduced above for the reference case architecture (Nsc = 6). The trends discussed for the reference bottlebrush architecture are robust and hold for the longer side chains (Figure S6 and Figure S7, respectively). Note that the dip in κ 2 at early times for Case IV is more pronounced for the longer side chains (green curves in Figure S7b,d) than that for the reference bottlebrush architecture. Furthermore, the average shape anisotropy at equilibrium for Case IV (Nsc = 10 and 20) and Case III (Nsc =10) is also notably lower than that observed in the respective cases for Nsc = 6 (Figure S7b,d). This can be attributed to the fact that the clusters containing a smaller number of bottlebrushes are formed if bottlebrushes with longer side chains are used (both at early times and at equilibrium), hence there is effectively smaller relative influence of the neighboring bottlebrushes on the structural characteristics of a chosen bottlebrush.
We next compare the agglomerates at equilibrium by calculating the number of clusters and an average aggregation number for all side chain lengths considered (Nsc = 6, 10, and 20). The number of clusters at equilibrium normalized on the number of bottlebrushes, nMBB, as a function of the side chain affinities using all three cutoff criteria is provided in Figure 11a; circles, squares, and solid lines of the same color correspond to nc = 1, 13, and nc = 13 Nsc, respectively. This plot clearly shows that effectively requiring at least one side chain of one bottlebrush to be surrounded by the beads of another bottlebrush for the two brushes to belong to the same cluster (i.e., using the cut off criterion nc = 13 Nsc) for the cases where the bottlebrushes are relatively well dispersed (asc-sol    82 ) results in the number of clusters equal to the number of bottlebrushes. However, using the cutoff criteria nc = 1 or 13 for the same side chain affinities to the solvent results in a significantly smaller number of clusters for all Nsc considered (black, red, and blue symbols and solid lines in Figure 11a for Nsc = 6, 10, and 20, respectively). This behavior reflects the relative contribution of formation and breaking of transient short-lived agglomerates as discussed above for the reference case scenario, the same behavior holds for the longer side chains. Note that, for Nsc = 6, the same data as discussed above (Figure 7) are used; herein however the data are normalized on nMBB for clarity of comparison between different side chain lengths.
For the simulations of bottlebrush conformations in relatively poor solvents (asc-sol  > 84 ), all three criteria defined above result in an identical number of clusters upon equilibration (Figure 11a). These results show that the number of clusters in equilibrium normalized on nMBB decreases with an increase in Nsc in poor solvent. Note that an average size of the cluster in equilibrium expressed via the number of bottlebrush beads—as well as the number of clusters in equilibrium without the normalization on nMBB—remain approximately the same for all Nsc considered (Figure S8) for the solvents with asc-sol     84 ; this is due to the fact that the same volume fraction (5%) of bottlebrushes with different architectures is used. Finally, the aggregation number—as defined in Equation (7) (expressed in the number of bottlebrushes) with the cutoff criterion nc = 13 Nsc—is notably lower for higher Nsc (for asc-sol  > 84 ) due to the larger number of beads comprising a single bottlebrush. For the solvents with relatively high affinities to the side chains (asc-sol  < 82 ) , the lowest cutoff value (nc  = 1 , dotted curves) results in highest N a g g for the longest side chain, indicating a more pronounced contribution of short-lived agglomerates, while nc = 13 Nsc results in N a g g = 1 (Figure 11b). For the solvents of intermediate quality (82   asc-sol 84), the transient clusters are more pronounced for the longer side chains, as can be seen in shifting the transient regions to the left (to the lower asc-sol) in both plots with an increase in Nsc (Figure 11).
Finally, to quantify the effect of clustering on the structural properties of the bottlebrushes with longer side chains, we track structural characteristics of bottlebrushes in solvents of different qualities and compare these characteristics with that for a single bottlebrush. For the solvents with relatively high affinities to the side chains (asc-sol  < 82 ) , an average radius of gyration R g , average shape anisotropy κ 2 , a ratio of eigenvalues < λ 3 / λ 1 > , and an acylindricity normalized by the smallest eigenvalue c / λ 1 (Figure 12) all remain the same (within the error bars) for the multiple bottlebrushes (black curves, solid for Nsc  = 10 and dashed for Nsc  = 20 ) as those characteristics for the single bottlebrush (red curves, solid for Nsc = 10 and dashed for Nsc = 20 ). These results clearly indicate that, despite the increased contribution of formation of transient short-lived agglomerates with an increase in asc-sol within this range, this transient agglomeration does not affect the structural characteristics of the bottlebrushes. A comparison of structural characteristics with the characterization of the agglomeration at equilibrium (Figure 11) clearly indicates that using the cutoff value of nc = 13 Nsc allows one to predict the expected trends in bottlebrush conformation based on the calculated average aggregation number. Specifically, this cutoff value primarily captures long-lived clusters, hence it allows one to predict that the structural characteristics of bottlebrushes in solutions will remain the same as those of single bottlebrushes if N a g g = 1 . While using lower cut off values (nc  = 1   or 13) allows one to capture an existence of transient agglomerates at the same range of solvent qualities, these criteria are not sufficient to differentiate short-lived or transient and long-lived clusters.
On the contrary, in relatively poor solvents (asc-sol > 84 ) , the average radius of gyration R g , average shape anisotropy κ 2 , a ratio of eigenvalues < λ 3 / λ 1 > , and an acylindricity normalized by the smallest eigenvalue c / λ 1 (Figure 12) differ significantly for the multiple bottlebrushes (black curves, solid for Nsc = 10 and dashed for) Nsc = 20 and the single bottlebrush (red curves, solid for Nsc = 10 and dashed for Nsc = 20). These trends are the same as observed in the reference case (Nsc = 6) and correspond to the relative stretching of the bottlebrushes within the long-lived clusters. The average κ 2 is lowest for Nsc = 20 reflecting lower relative influence of the neighboring bottlebrushes on the structural characteristics of a chosen bottlebrush within the clusters since these clusters on average encompass a smaller number of bottlebrushes. (Recall that an average size of the cluster in equilibrium expressed via the number of bottlebrush beads remain approximately the same for all Nsc considered for the solvents with asc-sol    84 since the same volume fraction (5%) of bottlebrushes with different architectures is used.) Our results show that, in relatively poor solvents, all three criteria for the critical number of contacts can be used to describe the clustering process and predict an effect of clustering on the structural characteristics. For the solvents of intermediate quality (82 asc-sol 84 ), the transient clusters are more pronounced for the longer side chains as discussed above. Using the highest cutoff nc = 13 Nsc allows one to primarily focus on the effects of long-lived clusters, while using smaller cutoff allows one to also account for contribution of transient clustering processes. To summarize, these studies show that the cutoff nc = 13 Nsc can be used to predict the anticipated changes in structural characteristics of bottlebrushes upon agglomeration with respect to the same characteristics for the single bottlebrushes for all the solvent qualities considered herein. Finally, we replot the data shown in Figure 11 and Figure 12 parametrically for the average shape anisotropy κ 2 as a function of the aggregation number, N a g g (Figure S9). Herein, we only consider N a g g   calculated using highest critical number of contacts, nc = 13 Nsc. A consistent trend showing a decrease in κ 2 with the decrease in N a g g   is seen in Figure S9 for all the relatively poor solvent cases.

4. Conclusions

We characterized agglomeration of molecular bottlebrushes by tracking the time evolution of the number of clusters, the size of the largest cluster, and an average aggregation number. We systematically analyzed three different cutoff criteria to identify whether two bottlebrushes belong to the same cluster. This analysis allows one to introduce a clear protocol to identify both short-lived (or transient) and long-lived agglomerates. Our results demonstrate that the cutoff criterion which depends on both the coordination number and the length of the side chain allows one to largely disregard transient or short-lived clusters. Furthermore, this criterion allows one to correlate the agglomeration status with the structural characteristics of bottlebrushes. We characterized conformational changes of the bottlebrush within the agglomerates with respect to that of an isolated bottlebrush in the same solvent. Specifically, we tracked the radius of gyration, shape anisotropy, acylindricity, and the ratios of eigenvalues of the gyration tensor. These studies allowed us to clearly identify a range of solvent qualities at which conformational changes remain negligible; we also showed that only transient short-lived agglomerates are formed in these scenarios. Furthermore, we identified a range of affinities at which distinct differences in conformations of bottlebrushes within the agglomerates and that of an isolated bottlebrush in the same solvent are observed. Specifically, our results demonstrated a distinct increase in both the radius of gyration and the shape anisotropy for the bottlebrush within the agglomerate, which in turn indicates stretching of a bottlebrush within the agglomerates. The characterization of bottlebrush conformations within the agglomerates is an important step in characterizing the relationship between the bottlebrush architecture and material properties. A comparison between the three distinct cutoff criteria to identify whether two bottlebrushes belong to the same cluster introduces a protocol to identify both short-lived transient and long-lived agglomerates and could be further adapted to characterize the agglomeration process of macromolecules with various complex architectures.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/polym14122339/s1, Figure S1: End-to-end distance (<Re>) and Radius of gyration (<Rg>) of the backbone of a bottlebrush; Figure S2: Acylindricity c and ratio of eigenvalues λ 2 / λ 1 of a bottlebrush; Figure S3: Time evolution of distribution of clusters in case I and case IV; Figure S4: Representative snapshots of time evolution of nMBB = 49 bottlebrushes with Nsc = 10 in solvents of various qualities; Figure S5: Representative snapshots of time evolution of nMBB = 26 bottlebrushes with Nsc = 20 in solvents of various qualities; Figure S6: Time evolution of number of clusters and aggregation number ( N a g g ) for the side chain length Nsc = 10/20; Figure S7: Time evolution of radius of gyration and shape anisotropy for sidechain length Nsc = 10/20; Figure S8: Number of clusters and an average size of cluster in number of beads as a function of side chain affinity for Nsc = 10 and 20; Figure S9: Shape anisotropy κ2 as a function of average aggregation number and snapshots of bottlebrush cluster.

Author Contributions

Conceptualization, O.K., S.T., C.K.C. and I.L.; Methodology, S.T. and C.K.C.; Validation, S.T., C.K.C. and O.K.; Formal analysis, S.T., C.K.C., M.G. and O.K.; Investigation, S.T., C.K.C., I.L. and O.K.; Resources, O.K.; Data curation, S.T. and C.K.C.; Writing—original draft preparation, S.T.; Writing—review and editing, S.T., C.K.C., I.L. and O.K.; Visualization, S.T.; Supervision, O.K.; Funding acquisition, I.L. and O.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Science Foundation EPSCoR Program under NSF Award no. OIA-1655740. M.G. contribution was supported by National Science Foundation under Award no. 1950557 (REU site). Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect those of the National Science Foundation.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Clemson University is acknowledged for generous allotment of compute time on The Palmetto Cluster.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sheiko, S.S.; Sumerlin, B.S.; Matyjaszewski, K. Cylindrical molecular brushes: Synthesis, characterization, and properties. Prog. Polym. Sci. 2008, 33, 759–785. [Google Scholar] [CrossRef]
  2. Lee, H.I.; Pietrasik, J.; Sheiko, S.S.; Matyjaszewski, K. Stimuli-responsive molecular brushes. Prog. Polym. Sci. 2010, 35, 24–44. [Google Scholar] [CrossRef]
  3. Zhang, N.; Huber, S.; Schulz, A.; Luxenhofer, R.; Jordan, R. Cylindrical Molecular Brushes of Poly(2-oxazoline)s from 2-Isopropenyl-2-oxazoline. Macromolecules 2009, 42, 2215–2221. [Google Scholar] [CrossRef]
  4. Rzayev, J. Molecular Bottlebrushes: New Opportunities in Nanomaterials Fabrication. ACS Macro Lett. 2012, 1, 1146–1149. [Google Scholar] [CrossRef]
  5. Tu, S.; Choudhury, C.K.; Luzinov, I.; Kuksenok, O. Recent advances towards applications of molecular bottlebrushes and their conjugates. Curr. Opin. Solid State Mater. Sci. 2019, 23, 50–61. [Google Scholar] [CrossRef]
  6. Xie, G.; Martinez, M.R.; Olszewski, M.; Sheiko, S.S.; Matyjaszewski, K. Molecular bottlebrushes as novel materials. Biomacromolecules 2018, 20, 27–54. [Google Scholar] [CrossRef]
  7. Daniel, W.F.M.; Burdyńska, J.; Vatankhah-Varnoosfaderani, M.; Matyjaszewski, K.; Paturej, J.; Rubinstein, M.; Dobrynin, A.V.; Sheiko, S.S. Solvent-free, supersoft and superelastic bottlebrush melts and networks. Nat. Mater. 2016, 15, 183–189. [Google Scholar] [CrossRef]
  8. Paturej, J.; Sheiko, S.S.; Panyukov, S.; Rubinstein, M. Molecular structure of bottlebrush polymers in melts. Sci. Adv. 2016, 2, e1601478. [Google Scholar] [CrossRef] [Green Version]
  9. Pan, T.; Dutta, S.; Kamble, Y.; Patel, B.B.; Wade, M.A.; Rogers, S.A.; Diao, Y.; Guironnet, D.; Sing, C.E. Materials Design of Highly Branched Bottlebrush Polymers at the Intersection of Modeling, Synthesis, Processing, and Characterization. Chem. Mater. 2022, 34, 1990–2024. [Google Scholar] [CrossRef]
  10. Motoyanagi, J.; Ishikawa, T.; Minoda, M. Stimuli-responsive brush-shaped conjugated polymers with pendant well-defined poly(vinyl ether)s. J. Polym. Sci. Part A Polym. Chem. 2016, 54, 3318–3325. [Google Scholar] [CrossRef]
  11. Liberman-Martin, A.L.; Chu, C.K.; Grubbs, R.H. Application of Bottlebrush Block Copolymers as Photonic Crystals. Macromol. Rapid Commun. 2017, 38, 1700058. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Chae, C.-G.; Yu, Y.-G.; Seo, H.-B.; Kim, M.-J.; Grubbs, R.H.; Lee, J.-S. Experimental Formulation of Photonic Crystal Properties for Hierarchically Self-Assembled POSS–Bottlebrush Block Copolymers. Macromolecules 2018, 51, 3458–3466. [Google Scholar] [CrossRef] [Green Version]
  13. Müllner, M. Molecular Polymer Brushes in Nanomedicine. Macromol. Chem. Phys. 2016, 217, 2209–2222. [Google Scholar] [CrossRef]
  14. Fenyves, R.; Schmutz, M.; Horner, I.J.; Bright, F.V.; Rzayev, J. Aqueous Self-Assembly of Giant Bottlebrush Block Copolymer Surfactants as Shape-Tunable Building Blocks. J. Am. Chem. Soc. 2014, 136, 7762–7770. [Google Scholar] [CrossRef] [PubMed]
  15. López-Barrón, C.R.; Vargas-Lara, F.; Kang, S. Single-Chain Conformation of Poly (α-olefins) in Dilute Solutions at the Crossover between Linear and Bottlebrush Architectures. Macromolecules 2021, 54, 6854–6866. [Google Scholar] [CrossRef]
  16. Rathgeber, S.; Pakula, T.; Wilk, A.; Matyjaszewski, K.; Beers, K.L. On the shape of bottle-brush macromolecules: Systematic variation of architectural parameters. J. Chem. Phys. 2005, 122, 124904. [Google Scholar] [CrossRef] [Green Version]
  17. Sheiko, S.S.; Sun, F.C.; Randall, A.; Shirvanyants, D.; Rubinstein, M.; Lee, H.-i.; Matyjaszewski, K. Adsorption-induced scission of carbon–carbon bonds. Nature 2006, 440, 191–194. [Google Scholar] [CrossRef]
  18. Panyukov, S.V.; Sheiko, S.S.; Rubinstein, M. Amplification of tension in branched macromolecules. Phys. Rev. Lett. 2009, 102, 148301. [Google Scholar] [CrossRef] [Green Version]
  19. Panyukov, S.; Zhulina, E.B.; Sheiko, S.S.; Randall, G.C.; Brock, J.; Rubinstein, M. Tension Amplification in Molecular Brushes in Solutions and on Substrates. J. Phys. Chem. B 2009, 113, 3750–3768. [Google Scholar] [CrossRef] [Green Version]
  20. Kutnyanszky, E.; Hempenius, M.A.; Vancso, G.J. Polymer bottlebrushes with a redox responsive backbone feel the heat: Synthesis and characterization of dual responsive poly(ferrocenylsilane)s with PNIPAM side chains. Polym. Chem. 2014, 5, 771–783. [Google Scholar] [CrossRef]
  21. Paturej, J.; Kuban, L.; Milchev, A.; Vilgis, T.A. Tension enhancement in branched macromolecules upon adhesion on a solid substrate. EPL Europhys. Lett. 2012, 97, 58003. [Google Scholar] [CrossRef] [Green Version]
  22. Wessels, M.G.; Jayaraman, A. Molecular dynamics simulation study of linear, bottlebrush, and star-like amphiphilic block polymer assembly in solution. Soft Matter 2019, 15, 3987–3998. [Google Scholar] [CrossRef] [PubMed]
  23. Gumerov, R.A.; Potemkin, I.I. Computer simulations of comb-like macromolecules with responsive diblock copolymer side chains. Colloid Polym. Sci. 2021, 299, 407–418. [Google Scholar] [CrossRef]
  24. Li, H.; Liu, H.; Nie, T.; Chen, Y.; Wang, Z.; Huang, H.; Liu, L.; Chen, Y. Molecular bottlebrush as a unimolecular vehicle with tunable shape for photothermal cancer therapy. Biomaterials 2018, 178, 620–629. [Google Scholar] [CrossRef] [PubMed]
  25. Jiang, L.; Nykypanchuk, D.; Pastore, V.J.; Rzayev, J. Morphological Behavior of Compositionally Gradient Polystyrene–Polylactide Bottlebrush Copolymers. Macromolecules 2019, 52, 8217–8226. [Google Scholar] [CrossRef]
  26. Li, X.; Wang, B.; Liu, Q.-J.; Zhao, R.; Song, D.-P.; Li, Y. Supersoft Elastic Bottlebrush Microspheres with Stimuli-Responsive Color-Changing Properties in Brine. Langmuir 2021, 37, 6744–6753. [Google Scholar] [CrossRef] [PubMed]
  27. Magurudeniya, H.D.; Ringstrand, B.S.; Seifert, S.; Firestone, M.A. Reversible hierarchical structure induced by solvation and temperature modulation in an ionic liquid-based random bottlebrush copolymer. Polym. Chem. 2018, 9, 5200–5214. [Google Scholar] [CrossRef]
  28. Kent, E.W.; Henn, D.M.; Zhao, B. Shape-changing linear molecular bottlebrushes with dually pH-and thermo-responsive diblock copolymer side chains. Polym. Chem. 2018, 9, 5133–5144. [Google Scholar] [CrossRef]
  29. Yao, K.; Chen, Y.; Zhang, J.; Bunyard, C.; Tang, C. Cationic Salt-Responsive Bottle-Brush Polymers. Macromol. Rapid Commun. 2013, 34, 645–651. [Google Scholar] [CrossRef]
  30. Li, C.M.; Gunari, N.; Fischer, K.; Janshoff, A.; Schmidt, M. New perspectives for the design of molecular actuators: Thermally induced collapse of single macromolecules from cylindrical brushes to spheres. Angew. Chem. Int. Edit. 2004, 43, 1101–1104. [Google Scholar] [CrossRef]
  31. Li, X.Y.; ShamsiJazeyi, H.; Pesek, S.L.; Agrawal, A.; Hammouda, B.; Verduzco, R. Thermoresponsive PNIPAAM bottlebrush polymers with tailored side-chain length and end-group structure. Soft Matter 2014, 10, 2008–2015. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Bejagam, K.K.; Singh, S.K.; Ahn, R.; Deshmukh, S.A. Unraveling the Conformations of Backbone and Side Chains in Thermosensitive Bottlebrush Polymers. Macromolecules 2019, 52, 9398–9408. [Google Scholar] [CrossRef]
  33. Kang, J.J.; Shehu, K.; Sachse, C.; Jung, F.A.; Ko, C.H.; Barnsley, L.C.; Jordan, R.; Papadakis, C.M. A molecular brush with thermoresponsive poly(2-ethyl-2-oxazoline) side chains: A structural investigation. Colloid Polym. Sci. 2021, 299, 193–203. [Google Scholar] [CrossRef]
  34. Ohnsorg, M.L.; Prendergast, P.C.; Robinson, L.L.; Bockman, M.R.; Bates, F.S.; Reineke, T.M. Bottlebrush Polymer Excipients Enhance Drug Solubility: Influence of End-Group Hydrophilicity and Thermoresponsiveness. ACS Macro Lett. 2021, 10, 375–381. [Google Scholar] [CrossRef] [PubMed]
  35. Kang, J.J.; Jung, F.A.; Ko, C.H.; Shehu, K.; Barnsley, L.C.; Kohler, F.; Dietz, H.; Zhao, J.P.; Pispas, S.; Papadakis, C.M. Thermoresponsive Molecular Brushes with Propylene Oxide/Ethylene Oxide Copolymer Side Chains in Aqueous Solution. Macromolecules 2020, 53, 4068–4081. [Google Scholar] [CrossRef]
  36. Mah, A.H.; Afzali, P.; Qi, L.; Pesek, S.; Verduzco, R.; Stein, G.E. Bottlebrush copolymer additives for immiscible polymer blends. Macromolecules 2018, 51, 5665–5675. [Google Scholar] [CrossRef] [Green Version]
  37. Wang, Y.; Zhang, S.; Jiang, H.; Lin, Y.; Xing, H.; Tang, T. Insight on compatibilization of LLDPE/PS blends from morphology, interfacial state, mechanical properties and melt properties: Comb-like copolymer vs diblock copolymer. Polymer 2021, 218, 123540. [Google Scholar] [CrossRef]
  38. Yadavalli, N.S.; Borodinov, N.; Choudhury, C.K.; Quiñones-Ruiz, T.; Laradji, A.M.; Tu, S.; Lednev, I.K.; Kuksenok, O.; Luzinov, I.; Minko, S. Thermal stabilization of enzymes with molecular brushes. ACS Catal. 2017, 7, 8675–8684. [Google Scholar] [CrossRef]
  39. Xu, Z.; Liu, S.; Liu, H.; Yang, C.; Kang, Y.; Wang, M. Unimolecular micelles of amphiphilic cyclodextrin-core star-like block copolymers for anticancer drug delivery. Chem. Commun. 2015, 51, 15768–15771. [Google Scholar] [CrossRef]
  40. Unsal, H.; Onbulak, S.; Calik, F.; Er-Rafik, M.; Schmutz, M.; Sanyal, A.; Rzayev, J. Interplay between molecular packing, drug loading, and core cross-linking in bottlebrush copolymer micelles. Macromolecules 2017, 50, 1342–1352. [Google Scholar] [CrossRef]
  41. Shi, X.; Zhang, Y.; Tian, Y.; Xu, S.; Ren, E.; Bai, S.; Chen, X.; Chu, C.; Xu, Z.; Liu, G. Multi-Responsive Bottlebrush-Like Unimolecules Self-Assembled Nano-Riceball for Synergistic Sono-Chemotherapy. Small Methods 2021, 5, 2000416. [Google Scholar] [CrossRef]
  42. Wei, L.; Caliskan, T.D.; Tu, S.; Choudhury, C.K.; Kuksenok, O.; Luzinov, I. Highly Oil-Repellent Thermoplastic Boundaries via Surface Delivery of CF3 Groups by Molecular Bottlebrush Additives. ACS Appl. Mater. Interfaces 2020, 12, 38626–38637. [Google Scholar] [CrossRef] [PubMed]
  43. Conder, J.M.; Hoke, R.A.; Wolf, W.d.; Russell, M.H.; Buck, R.C. Are PFCAs Bioaccumulative? A Critical Review and Comparison with Regulatory Criteria and Persistent Lipophilic Compounds. Environ. Sci. Technol. 2008, 42, 995–1003. [Google Scholar] [CrossRef] [PubMed]
  44. US Environmental Protection Agency: Long-chain Perfluorinated Chemicals (PFCs) Action Plan. Available online: https://www.epa.gov/assessing-and-managing-chemicals-under-tsca/long-chain-perfluorinated-chemicals-pfcs-action-plan (accessed on 28 April 2020).
  45. Available online: https://www.nationalgeographic.com/science/2020/01/pfas-contamination-safe-drinking-water-study/ (accessed on 28 April 2020).
  46. Buck, R.C.; Franklin, J.; Berger, U.; Conder, J.M.; Cousins, I.T.; de Voogt, P.; Jensen, A.A.; Kannan, K.; Mabury, S.A.; van Leeuwen, S.P.J. Perfluoroalkyl and polyfluoroalkyl substances in the environment: Terminology, classification, and origins. Integr. Environ. Asses. 2011, 7, 513–541. [Google Scholar] [CrossRef] [PubMed]
  47. Camaiti, M.; Brizi, L.; Bortolotti, V.; Papacchini, A.; Salvini, A.; Fantazzini, P. An Environmental Friendly Fluorinated Oligoamide for Producing Nonwetting Coatings with High Performance on Porous Surfaces. ACS Appl. Mater. Interfaces 2017, 9, 37279–37288. [Google Scholar] [CrossRef] [PubMed]
  48. Toselli, M.; Messori, M.; Bongiovanni, R.; Malucelli, G.; Priola, A.; Pilati, F.; Tonelli, C. Poly(ϵ-caprolactone)-poly(fluoroalkylene oxide)-poly(ϵ-caprolactone) block copolymers. 2. Thermal and surface properties. Polymer 2001, 42, 1771–1779. [Google Scholar] [CrossRef]
  49. Lopez, G.; Ameduri, B.; Habas, J.-P. A Versatile Strategy to Synthesize Perfluoropolyether-Based Thermoplastic Fluoropolymers by Alkyne-Azide Step-Growth Polymerization. Macromol. Rapid Commun. 2016, 37, 711–717. [Google Scholar] [CrossRef]
  50. Abbasi, M.; Faust, L.; Wilhelm, M. Comb and bottlebrush polymers with superior rheological and mechanical properties. Adv. Mater. 2019, 31, 1806484. [Google Scholar] [CrossRef]
  51. Zhulina, E.; Borisov, O.; Pryamitsyn, V.A.; Birshtein, T. Coil-globule type transitions in polymers. 1. Collapse of layers of grafted polymer chains. Macromolecules 1991, 24, 140–149. [Google Scholar] [CrossRef]
  52. Pesek, S.L.; Xiang, Q.; Hammouda, B.; Verduzco, R. Small-angle neutron scattering analysis of bottlebrush backbone and side chain flexibility. J. Polym. Sci. Part B Polym. Phys. 2017, 55, 104–111. [Google Scholar] [CrossRef]
  53. Chan, J.M.; Kordon, A.C.; Zhang, R.; Wang, M. Direct visualization of bottlebrush polymer conformations in the solid state. Proc. Natl. Acad. Sci. USA 2021, 118, e2109534118. [Google Scholar] [CrossRef] [PubMed]
  54. Horkay, F.; Chremos, A.; Douglas, J.F.; Jones, R.; Lou, J.; Xia, Y. Comparative experimental and computational study of synthetic and natural bottlebrush polyelectrolyte solutions. J. Chem. Phys. 2021, 155, 074901. [Google Scholar] [CrossRef] [PubMed]
  55. Fytas, N.G.; Theodorakis, P.E. Molecular dynamics simulations of single-component bottle-brush polymers with flexible backbones under poor solvent conditions. J. Phys. Condens. Matter 2013, 25, 285105. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Lyubimov, I.; Wessels, M.G.; Jayaraman, A. Molecular dynamics simulation and PRISM theory study of assembly in solutions of amphiphilic bottlebrush block copolymers. Macromolecules 2018, 51, 7586–7599. [Google Scholar] [CrossRef]
  57. Pannuzzo, M.; Tilton, R.D.; Deserno, M. Responsive behavior of a branched-chain polymer network: A molecular dynamics study. Soft Matter 2018, 14, 6485–6495. [Google Scholar] [CrossRef]
  58. Paturej, J.; Kreer, T. Hierarchical excluded volume screening in solutions of bottlebrush polymers. Soft Matter 2017, 13, 8534–8541. [Google Scholar] [CrossRef] [Green Version]
  59. Sunday, D.F.; Chremos, A.; Martin, T.B.; Chang, A.B.; Burns, A.B.; Grubbs, R.H. Concentration Dependence of the Size and Symmetry of a Bottlebrush Polymer in a Good Solvent. Macromolecules 2020, 53, 7132–7140. [Google Scholar] [CrossRef]
  60. Theodorakis, P.; Paul, W.; Binder, K. Microphase separation in bottlebrush polymers under poor-solvent conditions. EPL Europhys. Lett. 2009, 88, 63002. [Google Scholar] [CrossRef] [Green Version]
  61. Theodorakis, P.; Paul, W.; Binder, K. Analysis of the cluster formation in two-component cylindrical bottle-brush polymers under poor solvent conditions. A simulation study. Eur. Phys. J. E 2011, 34, 52. [Google Scholar] [CrossRef]
  62. Verduzco, R.; Li, X.; Pesek, S.L.; Stein, G.E. Structure, function, self-assembly, and applications of bottlebrush copolymers. Chem. Soc. Rev. 2015, 44, 2405–2420. [Google Scholar] [CrossRef] [Green Version]
  63. Binder, K.; Hsu, H.-P.; Paul, W. Understanding the stiffness of macromolecules: From linear chains to bottle-brushes. Eur. Phys. J. Spec. Top. 2016, 225, 1663–1671. [Google Scholar] [CrossRef]
  64. Lebedeva, I.O.; Zhulina, E.B.; Borisov, O.V. Self-Assembly of Bottlebrush Block Copolymers in Selective Solvent: Micellar Structures. Polymers 2021, 13, 1351. [Google Scholar] [CrossRef] [PubMed]
  65. Mohammadi, E.; Joshi, S.Y.; Deshmukh, S.A. A review of computational studies of bottlebrush polymers. Comput. Mater. Sci. 2021, 199, 110720. [Google Scholar] [CrossRef]
  66. Dutta, S.; Wade, M.A.; Walsh, D.J.; Guironnet, D.; Rogers, S.A.; Sing, C.E. Dilute solution structure of bottlebrush polymers. Soft Matter 2019, 15, 2928–2941. [Google Scholar] [CrossRef] [Green Version]
  67. Košovan, P.; Kuldová, J.; Limpouchová, Z.; Procházka, K.; Zhulina, E.B.; Borisov, O.V. Amphiphilic graft copolymers in selective solvents: Molecular dynamics simulations and scaling theory. Macromolecules 2009, 42, 6748–6760. [Google Scholar] [CrossRef]
  68. Chang, H.-Y.; Lin, Y.-L.; Sheng, Y.-J.; Tsao, H.-K. Structural Characteristics and Fusion Pathways of Onion-Like Multilayered Polymersome Formed by Amphiphilic Comb-Like Graft Copolymers. Macromolecules 2013, 46, 5644–5656. [Google Scholar] [CrossRef]
  69. Luo, Z.; Li, Y.; Wang, B.; Jiang, J. pH-sensitive vesicles formed by amphiphilic grafted copolymers with tunable membrane permeability for drug loading/release: A multiscale simulation study. Macromolecules 2016, 49, 6084–6094. [Google Scholar] [CrossRef]
  70. Tan, H.; Yu, C.; Lu, Z.; Zhou, Y.; Yan, D. A dissipative particle dynamics simulation study on phase diagrams for the self-assembly of amphiphilic hyperbranched multiarm copolymers in various solvents. Soft Matter 2017, 13, 6178–6188. [Google Scholar] [CrossRef]
  71. Tan, H.; Wang, W.; Yu, C.; Zhou, Y.; Lu, Z.; Yan, D. Dissipative particle dynamics simulation study on self-assembly of amphiphilic hyperbranched multiarm copolymers with different degrees of branching. Soft Matter 2015, 11, 8460–8470. [Google Scholar] [CrossRef]
  72. Hoogerbrugge, P.J.; Koelman, J. Simulating Microscopic Hydrodynamic Phenomena with Dissipative Particle Dynamics. Europhys. Lett. 1992, 19, 155–160. [Google Scholar] [CrossRef]
  73. Español, P.; Warren, P. Statistical Mechanics of Dissipative Particle Dynamics. Europhys. Lett. EPL 1995, 30, 191–196. [Google Scholar] [CrossRef] [Green Version]
  74. Groot, R.D.; Warren, P.B. Dissipative particle dynamics: Bridging the gap between atomistic and mesoscopic simulation. J. Chem. Phys. 1997, 107, 4423–4435. [Google Scholar] [CrossRef]
  75. Yong, X. Hydrodynamic Interactions and Entanglements of Polymer Solutions in Many-Body Dissipative Particle Dynamics. Polymers 2016, 8, 426. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Santo, K.P.; Neimark, A.V. Effects of metal-polymer complexation on structure and transport properties of metal-substituted polyelectrolyte membranes. J. Colloid Interface Sci. 2021, 602, 654–668. [Google Scholar] [CrossRef] [PubMed]
  77. Feng, Y.H.; Zhang, X.P.; Zhao, Z.Q.; Guo, X.D. Dissipative Particle Dynamics Aided Design of Drug Delivery Systems: A Review. Mol. Pharm. 2020, 17, 1778–1799. [Google Scholar] [CrossRef] [PubMed]
  78. Wang, X.; Santo, K.P.; Neimark, A.V. Modeling Gas–Liquid Interfaces by Dissipative Particle Dynamics: Adsorption and Surface Tension of Cetyl Trimethyl Ammonium Bromide at the Air–Water Interface. Langmuir 2020, 36, 14686–14698. [Google Scholar] [CrossRef]
  79. Shin, J.M.; Kim, Y.; Yun, H.; Yi, G.-R.; Kim, B.J. Morphological Evolution of Block Copolymer Particles: Effect of Solvent Evaporation Rate on Particle Shape and Morphology. ACS Nano 2017, 11, 2133–2142. [Google Scholar] [CrossRef]
  80. Salib, I.; Yong, X.; Crabb, E.J.; Moellers, N.M.; McFarlin, G.T.; Kuksenok, O.; Balazs, A.C. Harnessing Fluid-Driven Vesicles to Pick Up and Drop Off Janus Particles. ACS Nano 2013, 7, 1224–1238. [Google Scholar] [CrossRef]
  81. Dutt, M.; Kuksenok, O.; Nayhouse, M.J.; Little, S.R.; Balazs, A.C. Modeling the Self-Assembly of Lipids and Nanotubes in Solution: Forming Vesicles and Bicelles with Transmembrane Nanotube Channels. ACS Nano 2011, 5, 4769–4782. [Google Scholar] [CrossRef]
  82. Gumerov, R.A.; Filippov, S.A.; Richtering, W.; Pich, A.; Potemkin, I.I. Amphiphilic microgels adsorbed at oil-water interfaces as mixers of two immiscible liquids. Soft Matter 2019, 15, 3978–3986. [Google Scholar] [CrossRef]
  83. Vishnyakov, A.; Lee, M.T.; Neimark, A.V. Prediction of the Critical Micelle Concentration of Nonionic Surfactants by Dissipative Particle Dynamics Simulations. J. Phys. Chem. Lett. 2013, 4, 797–802. [Google Scholar] [CrossRef] [PubMed]
  84. Mills, Z.G.; Mao, W.B.; Alexeev, A. Mesoscale modeling: Solving complex flows in biology and biotechnology. Trends Biotechnol. 2013, 31, 426–434. [Google Scholar] [CrossRef] [PubMed]
  85. Chen, S.S.; Yong, X. Dissipative particle dynamics modeling of hydrogel swelling by osmotic ensemble method. J. Chem. Phys. 2018, 149, 094904. [Google Scholar] [CrossRef] [PubMed]
  86. Singh, A.; Kuksenok, O.; Johnson, J.A.; Balazs, A.C. Tailoring the structure of polymer networks with iniferter-mediated photo-growth. Polym. Chem. 2016, 7, 2955–2964. [Google Scholar] [CrossRef] [Green Version]
  87. Choudhury, C.K.; Palkar, V.; Kuksenok, O. Computational Design of Nanostructured Soft Interfaces: Focus on Shape Changes and Spreading of Cubic Nanogels. Langmuir 2020, 36, 7109–7123. [Google Scholar] [CrossRef]
  88. Chen, S.S.; Yong, X. Elastocapillary interactions of thermoresponsive microgels across the volume phase transition temperatures. J. Colloid Interface Sci. 2021, 584, 275–280. [Google Scholar] [CrossRef]
  89. Nikolov, S.V.; Fernandez-Nieves, A.; Alexeev, A. Behavior and mechanics of dense microgel suspensions. Proc. Natl. Acad. Sci. USA 2020, 117, 27096–27103. [Google Scholar] [CrossRef]
  90. Anakhov, M.V.; Gumerov, R.A.; Potemkin, I.I. Stimuli-responsive aqueous microgels: Properties and applications. Mendeleev Commun. 2020, 30, 555–562. [Google Scholar] [CrossRef]
  91. Espanol, P.; Warren, P.B. Perspective: Dissipative Particle Dynamics. J. Chem. Phys 2017, 146, 150901. [Google Scholar] [CrossRef]
  92. Santo, K.P.; Neimark, A.V. Dissipative particle dynamics simulations in colloid and interface science: A review. Adv. Colloid Interface Sci. 2021, 298, 102545. [Google Scholar] [CrossRef]
  93. Teng, C.-Y.; Sheng, Y.-J.; Tsao, H.-K. Boundary-induced segregation in nanoscale thin films of athermal polymer blends. Soft Matter 2016, 12, 4603–4610. [Google Scholar] [CrossRef] [PubMed]
  94. Wang, Z.; Wang, X.; Ji, Y.; Qiang, X.; He, L.; Li, S. Bottlebrush block polymers in solutions: Self-assembled microstructures and interactions with lipid membranes. Polymer 2018, 140, 304–314. [Google Scholar] [CrossRef]
  95. Nardai, M.M.; Zifferer, G. Simulation of dilute solutions of linear and star-branched polymers by dissipative particle dynamics. J. Chem. Phys. 2009, 131, 124903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  96. Guo, J.; Liang, H.; Wang, Z.-G. Coil-to-globule transition by dissipative particle dynamics simulation. J. Chem. Phys. 2011, 134, 244904. [Google Scholar] [CrossRef]
  97. Liu, Y.; Dai, Y.; Xu, X. Dynamic Feature of Incipient Polymer Collapse below the Theta Point. J. Phys. Chem. B 2017, 121, 9469–9475. [Google Scholar] [CrossRef]
  98. Xu, X.; Gao, X. Polymer Knot in Solution near the θ Point. J. Phys. Chem. B 2020, 124, 2723–2729. [Google Scholar] [CrossRef]
  99. Groot, R.D.; Rabone, K.L. Mesoscopic Simulation of Cell Membrane Damage, Morphology Change and Rupture by Nonionic Surfactants. Biophys. J. 2001, 81, 725–736. [Google Scholar] [CrossRef] [Green Version]
  100. Nikunen, P.; Karttunen, M.; Vattulainen, I. How would you integrate the equations of motion in dissipative particle dynamics simulations? Comput. Phys. Commun. 2003, 153, 407–423. [Google Scholar] [CrossRef] [Green Version]
  101. Sirk, T.W.; Slizoberg, Y.R.; Brennan, J.K.; Lisal, M.; Andzelm, J.W. An enhanced entangled polymer model for dissipative particle dynamics. J. Chem. Phys. 2012, 136, 134903. [Google Scholar] [CrossRef]
  102. Vishnyakov, A.; Mao, R.; Lee, M.-T.; Neimark, A.V. Coarse-grained model of nanoscale segregation, water diffusion, and proton transport in Nafion membranes. J. Chem. Phys. 2018, 148, 024108. [Google Scholar] [CrossRef]
  103. Bray, D.J.; Anderson, R.L.; Warren, P.B.; Lewtas, K. Wax Formation in Linear and Branched Alkanes with Dissipative Particle Dynamics. J. Chem. Theory Comput. 2020, 16, 7109–7122. [Google Scholar] [CrossRef] [PubMed]
  104. Choudhury, C.K.; Kuksenok, O. Native-Based Dissipative Particle Dynamics Approach for α-Helical Folding. J. Phys. Chem. B 2020, 124, 11379–11386. [Google Scholar] [CrossRef] [PubMed]
  105. Palkar, V.; Kuksenok, O. Controlling Degradation and Erosion of Polymer Networks: Insights from Mesoscale Modeling. J. Phys. Chem. B 2021, 126, 336–346. [Google Scholar] [CrossRef] [PubMed]
  106. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  107. LAMMPS Molecular Dynamics Simulator. Available online: http://lammps.sandia.gov (accessed on 29 October 2020).
  108. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. Model. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  109. Theodorou, D.N.; Suter, U.W. Shape of Unperturbed Linear-Polymers—Polypropylene. Macromolecules 1985, 18, 1206–1214. [Google Scholar] [CrossRef]
  110. Moreno, N.; Nunes, S.P.; Peinemann, K.-V.; Calo, V.M. Topology and Shape Control for Assemblies of Block Copolymer Blends in Solution. Macromolecules 2015, 48, 8036–8044. [Google Scholar] [CrossRef] [Green Version]
  111. Bosko, J.T.; Todd, B.D.; Sadus, R.J. Analysis of the shape of dendrimers under shear. J. Chem. Phys. 2006, 124, 044910. [Google Scholar] [CrossRef]
  112. Lee, H.; Baker, J.R.; Larson, R.G. Molecular Dynamics Studies of the Size, Shape, and Internal Structure of 0% and 90% Acetylated Fifth-Generation Polyamidoamine Dendrimers in Water and Methanol. J. Phys. Chem. B 2006, 110, 4014–4019. [Google Scholar] [CrossRef] [Green Version]
  113. Jagodzinski, O.; Eisenriegler, E.; Kremer, K. Universal shape properties of open and closed polymer chains: Renormalization group analysis and Monte Carlo experiments. J. Phys. I 1992, 2, 2243–2279. [Google Scholar] [CrossRef]
  114. Narros, A.; Moreno, A.J.; Likos, C.N. Effects of Knots on Ring Polymers in Solvents of Varying Quality. Macromolecules 2013, 46, 3654–3668. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  115. Steinhauser, M.O. A molecular dynamics study on universal properties of polymer chains in different solvent qualities. Part I. A review of linear chain properties. J. Chem. Phys. 2005, 122, 094901. [Google Scholar] [CrossRef] [PubMed]
  116. Larin, S.; Lyulin, S.; Lyulin, A.; Darinskii, A. Computer Simulations of Interpolyelectrolyte Complexes Formed by Star-like Polymers and Linear Polyelectrolytes. Macromol. Symp. 2009, 278, 40–47. [Google Scholar] [CrossRef]
  117. Morphew, D.; Chakrabarti, D. Hierarchical self-assembly of colloidal magnetic particles into reconfigurable spherical structures. Nanoscale 2015, 7, 8343–8350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  118. Stillinger, F.H., Jr. Rigorous Basis of the Frenkel-Band Theory of Association Equilibrium. J. Chem. Phys. 1963, 38, 1486–1494. [Google Scholar] [CrossRef] [Green Version]
  119. Theodorakis, P.E.; Paul, W.; Binder, K. Pearl-necklace structures of molecular brushes with rigid backbone under poor solvent conditions: A simulation study. J. Chem. Phys. 2010, 133, 104901. [Google Scholar] [CrossRef]
  120. Oh, K.J.; Zeng, X.C. A small-system ensemble Monte Carlo simulation of supersaturated vapor: Evaluation of barrier to nucleation. J. Chem. Phys. 2000, 112, 294–300. [Google Scholar] [CrossRef]
  121. Wedekind, J.; Strey, R.; Reguera, D. New method to analyze simulations of activated processes. J. Chem. Phys. 2007, 126, 134103. [Google Scholar] [CrossRef]
  122. Wedekind, J.; Reguera, D. What is the best definition of a liquid cluster at the molecular scale? J. Chem. Phys. 2007, 127, 154516. [Google Scholar] [CrossRef]
  123. Wolde, P.R.T.; Frenkel, D. Computer simulation study of gas–Liquid nucleation in a Lennard-Jones system. J. Chem. Phys. 1998, 109, 9901–9918. [Google Scholar] [CrossRef]
  124. Johnston, M.A.; Swope, W.C.; Jordan, K.E.; Warren, P.B.; Noro, M.G.; Bray, D.J.; Anderson, R.L. Toward a Standard Protocol for Micelle Simulation. J. Phys. Chem. B 2016, 120, 6337–6351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  125. Lee, M.-T.; Vishnyakov, A.; Neimark, A.V. Calculations of Critical Micelle Concentration by Dissipative Particle Dynamics Simulations: The Role of Chain Rigidity. J. Phys. Chem. B 2013, 117, 10304–10310. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (a) Bottlebrush architecture, backbone (bb), and side chain (sc) beads are shown in green and blue, Nbb and Nsc denote numbers of beads within the backbone and side chain, respectively. (b) Initial distribution of nMBB = 77 bottlebrushes equilibrated within the good solvent; solvent beads are shown in purple using reduced bead size.
Figure 1. (a) Bottlebrush architecture, backbone (bb), and side chain (sc) beads are shown in green and blue, Nbb and Nsc denote numbers of beads within the backbone and side chain, respectively. (b) Initial distribution of nMBB = 77 bottlebrushes equilibrated within the good solvent; solvent beads are shown in purple using reduced bead size.
Polymers 14 02339 g001
Figure 2. (a) Late time snapshots of a bottlebrush conformations (Nsc = 6) in solvents of various quality: asc-sol = 78 in (I), asc-sol = 83 in (II), asc-sol = 85 in (III), and asc-sol = 89 in (IV). (b) Radius of gyration of the bottlebrush upon equilibration in solvents of various quality from the good solvent (asc-sol = 78) to the poor solvent (asc-sol = 89); black, red, blue, and green curves correspond to Nsc = 1, 6, 10, and 20, respectively. The data are averaged over last 6 × 106 steps (1200 frames), error bars represent standard deviations.
Figure 2. (a) Late time snapshots of a bottlebrush conformations (Nsc = 6) in solvents of various quality: asc-sol = 78 in (I), asc-sol = 83 in (II), asc-sol = 85 in (III), and asc-sol = 89 in (IV). (b) Radius of gyration of the bottlebrush upon equilibration in solvents of various quality from the good solvent (asc-sol = 78) to the poor solvent (asc-sol = 89); black, red, blue, and green curves correspond to Nsc = 1, 6, 10, and 20, respectively. The data are averaged over last 6 × 106 steps (1200 frames), error bars represent standard deviations.
Polymers 14 02339 g002
Figure 3. (a) Shape anisotropy κ2, (b) ratio of eigenvalues λ 3 / λ 1 , and (c) acylindricity c normalized by the eigenvalue λ 1 . The data are averaged over last 6 × 106 steps (1200 frames), error bars represent standard deviations. Black, red, blue, and green curves correspond to Nsc = 1, 6, 10, and 20, respectively. (d) Late time snapshots for bottlebrushes with Nsc = 1 in good and poor solvents (asc-sol = 78 in (I), asc-sol = 89 in (IV), respectively), bottlebrushes with Nsc = 10 in good and poor solvents (asc-sol = 78 in (II), asc-sol = 89 in (V), respectively), and bottlebrushes with Nsc = 20 in good and poor solvents (asc-sol = 78 in (IV) and asc-sol = 89 in (VI), respectively).
Figure 3. (a) Shape anisotropy κ2, (b) ratio of eigenvalues λ 3 / λ 1 , and (c) acylindricity c normalized by the eigenvalue λ 1 . The data are averaged over last 6 × 106 steps (1200 frames), error bars represent standard deviations. Black, red, blue, and green curves correspond to Nsc = 1, 6, 10, and 20, respectively. (d) Late time snapshots for bottlebrushes with Nsc = 1 in good and poor solvents (asc-sol = 78 in (I), asc-sol = 89 in (IV), respectively), bottlebrushes with Nsc = 10 in good and poor solvents (asc-sol = 78 in (II), asc-sol = 89 in (V), respectively), and bottlebrushes with Nsc = 20 in good and poor solvents (asc-sol = 78 in (IV) and asc-sol = 89 in (VI), respectively).
Polymers 14 02339 g003
Figure 4. Representative snapshots of time evolution of nMBB = 77 bottlebrushes with Nsc = 6 in solvents of various qualities: (a) asc-sol = 78 (Case I), (b) asc-sol = 83 (Case II), (c) asc-sol = 84 (Case III), and (d) asc-sol = 89 (Case IV). The snapshots correspond to the time steps provided below each column, the solvent is not shown for clarity of representation.
Figure 4. Representative snapshots of time evolution of nMBB = 77 bottlebrushes with Nsc = 6 in solvents of various qualities: (a) asc-sol = 78 (Case I), (b) asc-sol = 83 (Case II), (c) asc-sol = 84 (Case III), and (d) asc-sol = 89 (Case IV). The snapshots correspond to the time steps provided below each column, the solvent is not shown for clarity of representation.
Polymers 14 02339 g004
Figure 5. Time evolution of the size of the largest cluster size (number of bottlebrushes) in various solvents: (a) asc-sol = 78 (Case I), (b) asc-sol = 83 (Case II), (c) asc-sol = 84 Case (III), and (d) asc-sol = 89 (Case IV). Black dotted, red dash-dotted, and blue solid lines of the same color represent the number of contacts used to identify clusters (nc = 1, 13, and 78, respectively).
Figure 5. Time evolution of the size of the largest cluster size (number of bottlebrushes) in various solvents: (a) asc-sol = 78 (Case I), (b) asc-sol = 83 (Case II), (c) asc-sol = 84 Case (III), and (d) asc-sol = 89 (Case IV). Black dotted, red dash-dotted, and blue solid lines of the same color represent the number of contacts used to identify clusters (nc = 1, 13, and 78, respectively).
Polymers 14 02339 g005
Figure 6. Time evolution of (a) number of clusters, and (b) aggregation number N a g g for asc-sol = 78 (Case I, in black), asc-sol = 83 (Case II, in red), asc-sol = 84 Case (III, in blue), and asc-sol = 89 (Case IV, in green). Circles, squares, and solid lines of the same color represent the number of contacts used to identify clusters (nc = 1, 13, and 78, respectively). The data are averaged over nMBB = 77 bottlebrushes in four independent simulation sets for each time step, the error bars represent standard deviation.
Figure 6. Time evolution of (a) number of clusters, and (b) aggregation number N a g g for asc-sol = 78 (Case I, in black), asc-sol = 83 (Case II, in red), asc-sol = 84 Case (III, in blue), and asc-sol = 89 (Case IV, in green). Circles, squares, and solid lines of the same color represent the number of contacts used to identify clusters (nc = 1, 13, and 78, respectively). The data are averaged over nMBB = 77 bottlebrushes in four independent simulation sets for each time step, the error bars represent standard deviation.
Polymers 14 02339 g006
Figure 7. Late time average of (a) number of clusters and (b) aggregation number for sidechain length Nsc = 6. Dotted, dash-dotted, and solid lines of the same color represent the number of contacts used to identify clusters (nc = 1, 13, and 78, respectively). The data are averaged over nMBB = 77 bottlebrushes in four independent simulation sets for each time step, the error bars represent standard deviation.
Figure 7. Late time average of (a) number of clusters and (b) aggregation number for sidechain length Nsc = 6. Dotted, dash-dotted, and solid lines of the same color represent the number of contacts used to identify clusters (nc = 1, 13, and 78, respectively). The data are averaged over nMBB = 77 bottlebrushes in four independent simulation sets for each time step, the error bars represent standard deviation.
Polymers 14 02339 g007
Figure 8. Time evolution of (a) radius of gyration and (b) shape anisotropy. At each time step, the data are averaged over all nMBB = 77 bottlebrushes in four independent simulation sets, error bars represent standard deviation. Black, red, blue, and green curves correspond to affinities asc-sol =78, 83, 84, and 89, respectively.
Figure 8. Time evolution of (a) radius of gyration and (b) shape anisotropy. At each time step, the data are averaged over all nMBB = 77 bottlebrushes in four independent simulation sets, error bars represent standard deviation. Black, red, blue, and green curves correspond to affinities asc-sol =78, 83, 84, and 89, respectively.
Polymers 14 02339 g008
Figure 9. (a) Radius of gyration (<Rg>) of bottlebrushes upon equilibration as a function of sidechains affinity for nMBB = 77 (black curve) and for the single bottlebrush (red curve). The values are averaged over all nMBB = 77 bottlebrushes in four independent simulation sets within last 3 × 105 timesteps (60 frames for each set), error bars represent standard deviation. (b) Representative snapshot highlighting a single bottlebrush that is a part of the larger cluster (asc-sol = 89, t = 105). (c,d) Distribution of Rg in good solvent asc-sol = 78 (in (c)) and in poor solvent asc-sol = 89 (in (d)) for nMBB = 77 bottlebrushes in four sets within the last 3 × 105 timesteps.
Figure 9. (a) Radius of gyration (<Rg>) of bottlebrushes upon equilibration as a function of sidechains affinity for nMBB = 77 (black curve) and for the single bottlebrush (red curve). The values are averaged over all nMBB = 77 bottlebrushes in four independent simulation sets within last 3 × 105 timesteps (60 frames for each set), error bars represent standard deviation. (b) Representative snapshot highlighting a single bottlebrush that is a part of the larger cluster (asc-sol = 89, t = 105). (c,d) Distribution of Rg in good solvent asc-sol = 78 (in (c)) and in poor solvent asc-sol = 89 (in (d)) for nMBB = 77 bottlebrushes in four sets within the last 3 × 105 timesteps.
Polymers 14 02339 g009
Figure 10. (a) Shape anisotropy κ2 upon equilibration as a function of side chains affinity for nMBB = 77 (black curve) and for the single bottlebrush (red curve). (b,c) Distribution of the κ2 in good solvent asc-sol = 78 (in (b)) and in poor solvent asc-sol = 89 (in (c)) for 77 bottlebrushes in four sets with the last 3 × 105 time steps. (d) Ratio of eigenvalues λ31 upon equilibration as a function of side chains affinity for nMBB = 77 (black curve) and for the single bottlebrush (red curve). In (a,d), the data are averaged over all 77 bottlebrushes in four sets within the last 3 × 105 time steps (60 frames for each set), error bars represent standard deviation. In (b,c), data used to plot distributions includes relevant characteristics of all 77 bottlebrushes in four sets within the last 3 × 105 timesteps.
Figure 10. (a) Shape anisotropy κ2 upon equilibration as a function of side chains affinity for nMBB = 77 (black curve) and for the single bottlebrush (red curve). (b,c) Distribution of the κ2 in good solvent asc-sol = 78 (in (b)) and in poor solvent asc-sol = 89 (in (c)) for 77 bottlebrushes in four sets with the last 3 × 105 time steps. (d) Ratio of eigenvalues λ31 upon equilibration as a function of side chains affinity for nMBB = 77 (black curve) and for the single bottlebrush (red curve). In (a,d), the data are averaged over all 77 bottlebrushes in four sets within the last 3 × 105 time steps (60 frames for each set), error bars represent standard deviation. In (b,c), data used to plot distributions includes relevant characteristics of all 77 bottlebrushes in four sets within the last 3 × 105 timesteps.
Polymers 14 02339 g010
Figure 11. (a) Number of clusters scaled by total number of bottlebrushes and (b) average aggregation number as a function of solvent quality. Black, red, and blue curves correspond to the side chain length Nsc = 6, 10, and 20. Dotted, dash-dotted, and solid lines of the same color represent the minimum cutoff value used to identify clusters (1, 13, and 13 multiplied by Nsc, respectively).
Figure 11. (a) Number of clusters scaled by total number of bottlebrushes and (b) average aggregation number as a function of solvent quality. Black, red, and blue curves correspond to the side chain length Nsc = 6, 10, and 20. Dotted, dash-dotted, and solid lines of the same color represent the minimum cutoff value used to identify clusters (1, 13, and 13 multiplied by Nsc, respectively).
Polymers 14 02339 g011
Figure 12. (a) Radius of gyration Rg, (b) shape anisotropy κ2, (c) ratio of eigenvalues λ31, and (d) acylindricity normalized on the smallest eigenvalue, c/λ1, of the bottlebrushes as a function of side chain affinities for multiple bottlebrushes (black curves) and single bottlebrushes (red curves). The solid and dash-dotted lines correspond to the sidechain length Nsc = 10 and 20, correspondingly. For the multiple bottlebrushes (black curves), the data are averaged over all nMBB bottlebrushes within 60 frames (last 3 × 105 time steps) in four independent simulation runs, error bars represent standard deviations.
Figure 12. (a) Radius of gyration Rg, (b) shape anisotropy κ2, (c) ratio of eigenvalues λ31, and (d) acylindricity normalized on the smallest eigenvalue, c/λ1, of the bottlebrushes as a function of side chain affinities for multiple bottlebrushes (black curves) and single bottlebrushes (red curves). The solid and dash-dotted lines correspond to the sidechain length Nsc = 10 and 20, correspondingly. For the multiple bottlebrushes (black curves), the data are averaged over all nMBB bottlebrushes within 60 frames (last 3 × 105 time steps) in four independent simulation runs, error bars represent standard deviations.
Polymers 14 02339 g012
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tu, S.; Choudhury, C.K.; Giltner, M.; Luzinov, I.; Kuksenok, O. Mesoscale Modeling of Agglomeration of Molecular Bottlebrushes: Focus on Conformations and Clustering Criteria. Polymers 2022, 14, 2339. https://doi.org/10.3390/polym14122339

AMA Style

Tu S, Choudhury CK, Giltner M, Luzinov I, Kuksenok O. Mesoscale Modeling of Agglomeration of Molecular Bottlebrushes: Focus on Conformations and Clustering Criteria. Polymers. 2022; 14(12):2339. https://doi.org/10.3390/polym14122339

Chicago/Turabian Style

Tu, Sidong, Chandan K. Choudhury, Michaela Giltner, Igor Luzinov, and Olga Kuksenok. 2022. "Mesoscale Modeling of Agglomeration of Molecular Bottlebrushes: Focus on Conformations and Clustering Criteria" Polymers 14, no. 12: 2339. https://doi.org/10.3390/polym14122339

APA Style

Tu, S., Choudhury, C. K., Giltner, M., Luzinov, I., & Kuksenok, O. (2022). Mesoscale Modeling of Agglomeration of Molecular Bottlebrushes: Focus on Conformations and Clustering Criteria. Polymers, 14(12), 2339. https://doi.org/10.3390/polym14122339

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