## Exploiting Scaling Laws for Designing Polymeric Bottle Brushes: a Theoretical Coarse Graining for Homopolymeric Branched Polymers

Pietro Corsi, 1 Elia Roma, 1 Tecla Gasperi, 1 Fabio Bruni, 1 and Barbara Capone 1, ∗

1 Dipartimento di Scienze, Universit` degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy a (Dated: August 9, 2019)

Bottle brushes are polymeric macromolecules made of a linear polymeric backbone grafted with side chains. The choice of the grafting density s g , of the length ns of the grafted side chains and their chemical nature, fully determines the properties of each macromolecule, such as its elasticity and its folding behaviour. Typically experimental bottle brushes are systems made of tens of thousands of monomeric units, rendering a computational approach extremely expensive, especially in the case of bottle brushes solutions. A proper coarse graining description of these macromolecules, thus appears essential. We here present a theoretical approach able to develop a general, transferable and analytical multi-scale coarse graining of homopolymeric bottle brush polymers under good solvent condition. Starting from scaling theories, each macromolecule is mapped onto a chain of tethered star polymers, whose effective potential is known from scaling predictions, computational and experimental validations and can be expressed as a function of the number of arms f , and on the length na of each arm. Stars are then tethered to one another and the effective potentials between them is shown to only depend on the key parameters of the original bottle brush polymer ( s g , ns ). The generalised form of the effective potential is then used to reproduce properties of the macromolecules obtained both with scaling theories and with simulations. The general form of the effective potentials derived in the current study, allows to describe theoretically and computationally properties of homopolymeric bottle brush polymers for all grafting densities and all lengths of both backbone and grafted arms, opening the path for a manifold of applications.

PACS numbers: 82.70.-y, 81.16.Fg, 05.10.-a, 83.80.Uv

## INTRODUCTION

Scaling theories are a powerful tool allowing to map properties of extremely complex systems onto simple relations linking a selection of microscopic details to a subset of general features. Polymers, due to their fractal nature, are an excellent example of systems that can extensively be described by means of scaling theories for the most diversified architectures and chemical compositions [1]. The average size (or radius of gyration Rg ) of a single polymeric macromolecule can be, for example, mapped onto a power law linking the number of monomers n to a critical exponent n that grasps the monomersolvent interaction [2]: Rg ∼ n n . Scaling laws also apply to finite density systems either made of free macromolecules in solution or characterised by complex geometrical constraints. In this context, general properties of polymeric brushes can be described by means of simple scaling considerations [1, 3, 4] as well as through mean-field self-consistent calculations [57]. For polymeric systems grafted on a planar surface, scaling predictions allow to estimate the dependence of average height of a brush in terms of the grafting density s g = N A / - expressed as the average number N of polymers grafted per unit area A - the length n of the grafted polymers, and of the interaction between the monomers and solvent [1, 3]:

$$h \sim n \sigma _ { g } ^ { \frac { 1 - \nu } { 2 \nu } }. \quad \ \ \ \ ( 1 ) \quad \ \ c c$$

The results obtained for the planar case can then be generalised to brushes with spherical symmetry [8], that find their macromolecular correspondent in star-like polymers [9-11]. Similarly to how the average height of a planar brush was computed, it is possible to estimate the scaling law of the av- erage radius of gyration of such macromolecules in terms of the number na of monomers belonging to each of the f arms grafted on a central core :

$$R _ { s } \sim n _ { a } ^ { \nu } f ^ { 1 / 5 },$$

where for simplicity we assume the radius of gyration of the star to scale with the 1 / 5 exponent [8], neglecting finite size corrections to the scaling exponent linked to anisotropy that can arise in stars with a small number of arms.

Further works extended the theories of planar and spherical brushes to different geometries and to more complex grafting units [12, 14, 15]. The present manuscript focuses on the properties of the so-called bottle brush polymers, i.e. macromolecules made of a polymeric backbone of nb monomers, grafted with linear chains each one made of ns monomeric units. The grafting density is expressed as s g = 1 / ( bm ) , where bm is the unit of length between each grafting site expressed as the product of the number of m monomers between the diverse grafting sites, and the equilibrium distance b between bonded neighbouring monomers. In the following part of this work b is set to one, and all the simulations performed in good solvent condition, hence n ∼ 0 588. .

As side-chains fluctuate, the average monomer distribution around the backbone and the average height of the grafted homopolymeric chains can be predicted by assuming a tubular conformation for the macromolecule [16, 17].

The tube diameter Dc has been shown to scale as [12]:

$$D _ { c } = 2 R _ { c } \sim n _ { s } ^ { \nu } \tau ^ { 1 / 5 } \left ( \sigma _ { g } n _ { s } \right ) ^ { \nu / 5 },$$

where the solvent quality is taken into account both in the Flory exponent n and the reduced temperature t =( T -q ) / T ,

where T is the temperature of the system and q is the q -temperature e.g. the transition temperature between good solvent and bad solvent, at which polymers behave as random walks. As we do not change nor the temperature nor the solvent quality throughout this work, t is constant. By neglecting all contributions arising from constant terms, we obtain the following scaling law for the cylinder diameter: Dc ∼ n s n ( s gns ) n / 5 .

The radius of gyration Rbb of the whole macromolecule can be written as [12]:

$$R _ { b b } \sim N ^ { \nu } \tau ^ { 1 / 5 } \left ( \sigma _ { g } n _ { s } \right ) ^ { - 2 \nu / 5 } \sim n _ { b } ^ { \nu } \left ( 1 + \sigma _ { g } n _ { s } \right ) ^ { \nu } \left ( \sigma _ { g } n _ { s } \right ) ^ { - 2 \nu / 5 }, \ ( 4 ) \quad \text{with $s$} \quad.$$

where in the last term we ignore all constant terms, and write explicitly the total number of monomers N = nb ( 1 + s gns ) .

The large number of microscopic units constituting each bottle-brush renders an extensive computational study of these macromolecules prohibitive from a computational point of view. On the other hand, despite the progresses obtained in describing the system by means of scaling laws, detailed analyses either at finite density, or when particular phenomena have to be investigated, are still difficult to access. An effective description invokes a reduction of the degrees of freedom to be analysed, making of paramount importance the development of an efficient and quantitative multi-scale strategy.

This work is structured in two distinct sub-sequent sections; we first use full monomer Molecular Dynamics (MD) simulations to extract the key information and the set of effective potentials that we will need to perform the proposed coarse graining (CG). We then use the newly derived potentials to reproduce and compare the results obtained with CG to the full monomeric MD results, and to scaling predictions.

## RESULTS AND DISCUSSION

## The full monomer representation: deriving scaling laws and effective potentials

All calculations, both in the full monomer case as well as for the coarse grained representation, are done by means of MD simulations, performed with LAMMPS [19] in the NVT ensemble. In the full monomer representation, monomers in the simulations interact through a purely repulsive truncated and shifted Lennard-Jones potential (i.e, good solvent conditions),

$$V _ { L J } ( r ) = \begin{cases} 4 \epsilon \left [ \left ( \frac { \sigma } { r } \right ) ^ { 1 2 } - \left ( \frac { \sigma } { r } \right ) ^ { 6 } \right ] + \epsilon & \text{for $r\leq 2^{1/6} \sigma$} \\ 0 & \text{for $r>2^{1/6} \sigma$} \end{cases} \quad ( 5 ) \quad \quad 3.$$

with a cutoff set at 2 1 / 6 s .

Bonded monomers are tethered by means of a Finite Extensible Nonlinear Elastic (FENE) potential [13]:

$$\text{$n$ good} \\ \text{$as$ ran} _ { \text{ure nor} } \quad V _ { \text{FENE} } ( r ) = - 0. 5 k R _ { 0 } ^ { 2 } \ln \left [ 1 - \left ( \frac { r } { R _ { 0 } } \right ) ^ { 2 } \right ] + 4 \varepsilon \left [ \left ( \frac { \varsigma } { r } \right ) ^ { 1 2 } - \left ( \frac { \varsigma } { r } \right ) ^ { 6 } \right ] + \varepsilon \\ \text{nt. $By} \\ \text{mc we} \quad \text{where $l=2\infty/\varepsilon^2$ and $D$}. = 1 \, \varepsilon$$

where k = 30 e / s 2 and R 0 = 1 5 . s .

In the following part of this work, we use the reduced parameters T (temperature), e (interaction energy), and s = 1 (monomer diameter) .

The full monomeric simulations have been performed on bottle brushes made of a backbone of nb = 300 monomers, with side chains of ns monomers in the range ns ∈ [ 50 300 , ] for the various analysed cases, with an increment of d ns = 25 monomers between one case and another, and grafting densities s g = 0 2 0 25 0 33 0 5. . , . , . , .

Simulations are performed on macromolecules made of a total number of monomers in the 10 3 -10 5 range for at least 10 8 MDsteps.

Wefirst compute, within the full monomeric representation, general properties of the single bottle brush, such as the average cylinder radius Rc and the average radius of gyration Rbb for a wide set of ( s g , ns ) combinations, as shown in Figure 1 (a) and (b).

Rc can be extrapolated as the average distance of all monomers belonging to the grafted arms from the backbone. The probability distribution P r ( ⊥ ) of the monomer height, computed as the distance of all monomers along the orthogonal direction to the backbone, is shown in panel (c) of Figure 1. By rescaling the P r ( ⊥ ) computed for all of the ( s g , ns ) combinations, with the corresponding cylinder radius Rc ( s g , ns ) , we obtain the perfect rescaling of all curves onto a general mastercurve.

Early works on scaling behaviours of bottle brush polymers, [14, 15, 18] suggested that the length of the rigid portion of a bottle brush macromolecule is of the order of its cylinder diameter, or thickness. This renders possible to define within each macromolecule spherical symmetric elements - or super blobs -characterised by having a constant rigidity on a scale of the order of the cylinder diameter Dc . This work aims at developing a consistent and general coarse grained representation for homopolymeric bottle brushes in good solvent condition. The idea that will be exploited, is that each macromolecule can be seen as a chain of star-like super blobs , hence the coarse graining procedure will require three major steps:

- 1. identification of the number x of super blobs that will be used to represent the macromolecule;
- 2. mapping of each super blob onto an equivalent star polymer;
- 3. derivation of a general form for the set of effective potentials acting between the coarse grained beads.

We sketch the coarse graining scheme in Figure 2; panel (a) highlights the division of the bottle brush into super blobs of radius R x , while panel (b) shows the mapping of each super

/5

/

)

(

)

(

Figure 1: All plots are obtained for bottle-brushes made of nb = 300 monomers per backbone, with grafting densities s g = 0 2 0 25 0 33 . , . , . and 0 5 and arms of length . ns ∈ [ 50 300 , ] with d ns = 25 monomers per chain. Top left panels shows scaling laws for Rc as in (3), top right shows the scaling of the radius of gyration Rbb as in (4) and the bottom left panel shows the distribution P r ( ⊥ ) of monomers around the backbone, rescaled with the scaling law in equation (3). Snapshots in the top left panel show a bottle brush made by a backbone of nb = 300, with grafting density s g = 0 5 and arms made . of ns = 300 monomers, and the same brush enclosed in a cylinder of radius Rc ( s g = 0 5 . , ns = 300 . ) The snapshots in the bottom left figure highlight - in green - one of the arms of a bottle brush made of nb = 300 backbone monomers, with grafting density s g = 0 2 and . arms made of ns = 125 monomers. Side and top view show the way in which arms distribute around the backbone

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01183v2-figure-545ebb4a-4970-4728-854f-6355a2967859.jpg)

blob into a star polymer of radius Rs = R x made of a number f eq of arms each one with na = ns monomers.

=

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01183v2-figure-325a3c28-dd6c-4a2e-8546-c260014b6024.jpg)

=

Figure 2: Schematic representation of the coarse graining procedure. Panel (a) sketches the super blobbing of the bottle brush macromolecules into x spheric-symmetrical blobs. Each blob contains a number n x of backbone monomers as defined in equation (7) and higlighted in red. Panel (b) is a zoom on one of the x super blobs that will become basic units of the coarse grained representations. Each super blob is represented as a star polymer of radius Rs = R x , made of f = f eq arm of na = ns monomers. Blobs interact with one each other by means of the effective star-star potential (9).

g

## 1. Counting the number of super blobs

To identify the number of blobs, we exploit the spherical symmetry that each bottle brush presents on sub-portions of the macromolecule of lengths comparable to the cylinder radius. We therefore divide the brush in a number x of identical sub-bottle brush segments, each one made of n x backbone monomers, such that the radius of gyration R x of the segment is identical to the radius of the cylinder.

By imposing R x = Rc in (4) we obtain the number of n x of backbone monomers belonging to each of the of the x subsegments of the bottle brush

$$n _ { \xi } \sim \frac { R _ { c } ^ { 1 / \vee } } { \left ( 1 + \sigma _ { g } n _ { s } \right ) \left ( \sigma _ { g } n _ { s } \right ) ^ { - 2 / 5 } }. \quad \quad \quad \left ( 7 \right )$$

The number x = nb / n x of super blobs per bottle brush is then obtained as the integer ratio between the number of monomers per backbone and the number of monomers per blob.

## 2. Identify the equivalent star

We now want to identify the equivalent star on which we will map each super blob. In particular we want to derive the number f eq of equivalent arms, made of ns monomers each, needed so that the radius of gyration of the equivalent star polymer is identical to the radius of gyration of a bottle brush made of n x backbone monomers, with arms of ns monomers and a grafting density s g .

Equating (2) and (3) we obtain Rs = R x = Rc . It can be easily derived that the equivalent star has a number of arms:

$$f _ { \text{eq} } \sim \left ( R _ { c } / n _ { s } ^ { \nu } \right ) ^ { 5 } \sim \left [ \left ( n _ { s } ^ { \nu } \left ( \sigma _ { g } n _ { s } \right ) ^ { \nu / 5 } \right ) / n _ { s } ^ { \nu } \right ] ^ { 5 } \sim \left ( \sigma _ { g } n _ { s } \right ) ^ { \nu } \quad ( 8 )$$

each made of ns monomers, as sketched in Figure 2.

## 3. Deriving the effective potential

The first step into the proposed coarse graining procedure is to assume that each super blob interacts with other super blobs by means of a star-like effective pair potential Vs [10, 20]:

$$\text{:red}. \quad & \text{y means or a state-like executive pain poementual } v _ { s } \{ \mathbb { U }, \mathbb { U }. \, \mathbb { U }. \\ \text{macro-} \\ \text{ins a} \\ \mathbb { U } \, \text{as} \\ \mathbb { U } \, \text{each} \\ \mathbb { U } \, \text{each} \\ \mathbb { U } \, \text{each} \\ \mathbb { U } \, \text{each} \\ \mathbb { U } \, \text{each} \\ \mathbb { U } \, \text{for} \\ \mathbb { U } \, \text{$\mathbb { U }_{c}=\left(4/3\right)R_{s}=\left(4/3\right)R_{c}$ is the corona diameter of the}$$

where s c =( 4 3 / ) Rs =( 4 / 3 ) Rc is the corona diameter of the star polymer.

To estimate the tethering potential linking neighbouring coarse grained beads, we compute the distribution of distances P r i ( cm -r i + 1 cm ) between the centres of masses of the neighbouring super blobs, e.g. between the centres of masses of consecutive sub-parts of the bottle brush made of n x backbone beads. The logarithmic inversion of the probability distribution leads to the total potential acting between two neighbouring beads. Such an interaction will be modelled as the superposition of a repulsive potential and a tethering term F th ( r ) :

$$V _ { \xi } ( r ) = V _ { s } ( r ) + \Phi _ { \text{th} } ( r ), \text{ \quad \ \ } ( 1 0 )$$

where Vs ( r ) is the effective potential between two equivalent stars as in equation (9), and F th a tethering term that has the form:

$$\Phi _ { \text{th} } ( r ) = k ( r - r _ { 0 } ) ^ { 2 } \,, \quad \quad \quad ( 1 1 )$$

where k = 5 / 18 s 2 g ns and r 0 = s c .

The k and r 0 values have been determined by analysing all of the ( s g , ns ) combinations that we defined earlier on in this section. For clarity we show in panel (a) of Figure 3 only one amongst the possible sets of effective potentials that can be obtained for all the ( s g , ns ) combinations; in particular we present the potentials that are used to map a bottle brush with s g = 0 5 and . ns = 300 onto a chain of stars of f eq = 19 equivalent arms. The inset shows the comparison between the theoretical predictions for the total potential acting between two tethered beads, and the raw data obtained by logarithmically inverting the probability distributions for four ( s g , ns ) combinations as specified in the caption. Panel (b) of Figure 3 shows the fitting result for the elastic constant k obtained with all of the extracted potentials.

Remarkably, the set of effective potentials obtained in (10) only depend on a few details of the microscopic system, and can be described through a functional form. The sets of effective interactions obtained are thus general: all systems with the most diverse microscopic characteristics can be reproduced by the set of potentials just introduced, where (9) is the potential acting between non bonded units, while bonded units interact via the effective potential in (10).

## Testing the coarse graining procedure: bottle brushes as linear homopolymeric chains made of tethered star polymers

To test the coarse graining procedure, we performed a series of computational analyses by representing classes of generality of bottle brushes as linear homopolymeric chains. The monomeric units of each coarse grained molecule are coarsegrained star polymers, that interact with one each other by means of the purely repulsive potential (9), and are tethered by means of the effective interaction (11).

The mapping between the underlying full monomer representation and the coarse grained one, lies in the choice of the microscopic parameters ( s g , ns ) that univocally determine the bottle brush cylinder radius Rc . The latter becomes the unit

Figure 3: Panel (a): the effective potential acting between two neighbouring super blobs (black solid curve) split into the equivalent star polymer contribution (red dashed curve) and the tethering therm F th ( r ) (blue dash-dotted curve). The inset of the figure shows raw data vs theoretical fitting of the total tethering potentials between neighbouring subparts of the bottle brushes for s g = 0 2 . , ns = 100, (orange), s g = 0 25 . , ns = 150 (red), s g = 0 33 . , ns = 125 (magenta), s g = 0 5 . , ns = 50 (blue) Panel (b) shows the fitting of all of the spring constants k leading to the theoretical form in equation (11).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01183v2-figure-8ecee2bb-d34d-4d8e-92bc-a026df232826.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01183v2-figure-0ebcd80d-194c-4cac-b91e-2abd78644aa0.jpg)

length for the monomeric beads of the coarse grained chain Rs = Rc = R x . This reflects on the f eq value, thus on the particular set of effective interactions to be used between all super blobs. As for the mapping of the backbone length, since the number n x of backbone monomers per blob is fixed by the choice of Rc in (7), a variation of the number x of super blobs implies a change in the total microscopic backbone monomers per molecule nb = x n x . The more the blobs used, the longer the corresponding backbone.

We performed MD simulations, and computed the radius of gyration for coarse grained homopolymeric chains made of x = 5 10 15 20 25 30 50 100 200 effective monomers, for , , , , , , , , systems that correspond - on the microscopic level - to all

of the combinations of s g = 0 2 0 25 0 33 0 5 1 and . , . , . , . , ns = 50 100 200 400 600 800. , , , , ,

To check the consistency between the two representations, we compute the radius of gyration for all of the aforementioned systems in the coarse grained representation and map it to the monomer resolved bottle brush one (4). To perform the mapping, the first step is to compute from the x coarse grained beads the underlying monomer resolved number of backbone monomers:

$$n _ { b } = \xi n _ { \xi } \sim \xi \frac { R _ { c } ^ { 1 / \nu } } { \left ( 1 + \sigma _ { g } n _ { s } \right ) \left ( \sigma _ { g } n _ { s } \right ) ^ { - 2 / 5 } }. \quad \ \ \stackrel { \sim } { \asymp }$$

The radius of gyration of linear homopolymeric chain made of n beads of radius Rc is know to scale as [2]

$$R _ { g } \sim R _ { c } n ^ { \nu }. \quad \ \ \ ( 1 3 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

To compare the radii of gyration obtained with the two models, we thus need to compare (13) and (4) for all of the measured systems:

$$R _ { b b } \sim n _ { b } ^ { \vee } ( 1 + \sigma _ { g } n _ { s } ) ^ { \vee } ( \sigma _ { g } n _ { s } ) ^ { - 2 \vee / 5 } \sim R _ { g }. \quad ( 1 4 ) \quad \text{of $\sigma_{g}$} \\ \sigma _ { \cdot \, = }$$

Wethus obtain a general relation between the two representations:

$$R _ { g } / \left ( 1 + \sigma _ { g } n _ { s } \right ) ^ { \nu } \left ( \sigma _ { g } n _ { s } \right ) ^ { - 2 \nu / 5 } \sim n _ { b } ^ { \nu }, \quad \ \ ( 1 5 ) \quad \text{and} \, \mathbb { 1 }$$

where Rg in the left hand of equation (15) is computed in the coarse grained representation.

We show in Figure 4 all of the data obtained with the here introduced coarse grained method, falling on the predicted mastercurve (15), together with a line highlighting the correct dependence of Rg / ( 1 -s ns ) n ( s ns ) -2 n / 5 on nb . Three shapshots are also drawn: the top one shows a full monomer represented bottle brush, the middle coloured one is the full monomeric mapping - in terms of star polymers - corresponding to the original bottle brush, and the bottom one emphasises how each star is represented by means of its isotropic effective potential, shown as a sphere. The perfect agreement between all of the simulated systems, demonstrates the effectiveness of the coarse graining process. To assess the predictive power of the procedure, we extended simulations to the s g = 1 grafting density case, that had not been included in the derivation of the effective potential, reaching a reasonable agreement even for such a stiff molecule.

## CONCLUSIONS

Summarising, we introduced an accurate and extremely simple coarse grained methodology for homopolymeric bottle brushes based on scaling arguments. We derived a set of general effective potentials that can be used to represent the complex macromolecular systems as well as their functional

Figure 4: The plot shows the radius of gyration obtained by means of a coarse grained representation for the underlying bottle brush. The quantities reported on the y -axis has been rescaled according to equation (4), while the quantity plotted on x -axis has been rescaled using equation (12) to map results obtained by means of a homopolymeric representation on the predictions obtained with the scaling law of equation (4). Circles correspond to a microscopic grafting density of s g = 0 2, squares are . s g = 0 25, diamond are . s g = 0 33, stars are . s g = 0 5, triangles are . s g = 1, while the corresponding ns are colour coded in the legend of the graph (black for ns =50, red for 100, green for 200, blue for 400, magenta for 600 and orange for 800. The two snapshots show a bottle brush with nb = 300 , s g = 1 and ns = 225 and the full monomeric view of a bottle brush made of a series of star polymers. a is a fitting constant.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01183v2-figure-dacb577b-6ac6-4750-a6b7-a8f1ff982ef3.jpg)

dependence on a few microscopic details, and lay the basis for a further exploration of physical properties featured in complex macromolecules in diversified environments. We show that, within our coarse graining framework, macromolecules made of tens of thousands of monomers can be mapped onto a chain of a few effective beads, retaining the link between the coarse grained representation and the underlying microscopic one in the functional form of the effective potentials used. We tested the multiscale methodology by reproducing properties of the full monomeric molecules by means of our novel coarse grained representation. Notwithstanding the coarse grained molecules do not even retain the geometry of the original macromolecules, the results obtained with the multiscale approach have a quantitative agreement with both the full monomeric computational results, and with the predictions obtained with scaling laws, as shown in Figure 4. We in fact show that it is possible to map a full monomeric system characterised by a cylindrical geometry into a chain of spherical effective beads. The successful achievement of our quantitative coarse graining makes possible to access prohibitive regions, such as finite density analyses, that the large number of monomeric units involved would not allow to investigate with an monomeric description (each macromolecule is made of ∼ 10 3 -10 5 monomers). Consequently, the initial theoretical prediction of all interactions between the titled

macromolecule and its environment should give a facile access to the design and synthesis of bottle brushes suitable for a plethora of widespread uses. For example, the opportunity of foreseeing how to selectively tune the key interactions between each homopolymer and various bioactive compounds by simply changing either pH or temperature [21, 22], makes bottle brushes an extremely effective and efficient tool in the field of innovative drug delivery systems [23].

On a different perspective, the possibility of controlling the elastic properties of the macromolecules by tuning the ( s g , ns ) combination, also renders finite density solutions of bottle brushes of extreme interest due to their peculiar rheological properties [24].

The present methodology will allow for numerous applications. It can, for instance, be easily generalised to cases of polydisperse bottle brushes, where a different local grafting density s g and/or a different local length ns of the grafted arms will lead to a change in the elastic properties of the macromolecule. The coarse graining will allow to map such a case onto a chain of stars characterised by different radii of gyration, each represented by its corresponding effective potential, following the generalisation of the effective potential, as shown for example in reference [25]. The possibility to use this coarse graining representation to study the adsorbption properties of the bottle brushes, and to compare them to those of star polymer [11] is currently being investigated by the authors of this paper. Additionally, in order to experimentally validate our theoretical findings, the syntheses of various star polymers, as well as the corresponding bottle brushes, are ongoing in our laboratories and will be reported in due course.

## ACKNOWLEDGMENTS

BC acknowledges funding from the Marie Curie Individual Fellowship, project ID 751255. ER acknowledge the financial support from Erasmus+ Unipharma Graduates 2017/2018 and the hospitality of University of Natural Resources and Life Sciences Vienna. The Grant of Excellence Departments, MIUR-Italy (ARTICOLO 1, COMMI 314 - 337 LEGGE 232/2016) is gratefully acknowledged.

∗ Electronic address: barbara.capone@uniroma3.it

- [1] P. G. de Gennes, Scaling Concepts in Polymer Physics (Cornell University Press, 1980).
- [2] P. Flory, Principles of Polymer Chemistry (Cornell University Press, 1953).
- [3] S. Alexander, J. Phys. France 38 , 983 (1977).
- [4] I. Coluzza, B. Capone, and J.-P. Hansen, Soft Matter 7 , 5255 (2011).
- [5] S. T. Milner, T. A. Witten, and M. E. Cates, Macromolecules 21 , 2610 (1988).
- [6] C. M. Wijmans, J. M. H. M. Scheutjens, and Y. B. Zhulina, Macromolecules 25 , 2657 (1992).
- [7] R. R. Netz and M. Schick, Macromolecules 31 , 5105 (1998).
- [8] M. Daoud and J. P. Cotton, J. Phys. France 43 , 531 (1982).
- [9] M. Watzlawek, C. N. Likos, and H. L¨ owen, Phys. Rev. Lett. 82 , 5289 (1999).
- [10] C. N. Likos, H. L¨ owen, M. Watzlawek, B. Abbas, O. Jucknischke, J. Allgaier, and D. Richter, Phys. Rev. Lett. 80 , 4450 (1998).
- [11] D. Marzi, C. N. Likos, and B. Capone, The Journal of Chemical Physics 137 , 014902 (2012).
- [12] T. M. Birshtein, O. V. Borisov, Y. B. Zhulina, A. R. Khokhlov, and T. A. Yurasova, Polymer Science U.S.S.R. 29 , 1293 (1987), ISSN 0032-3950.
- [13] K. Kremer, G. S. Grest, J Chem. Phys. 94 , 4103 (1990)
- [14] Y. B. Zhulina, Polymer Science U.S.S.R. 26 , 885 (1984), ISSN 0032-3950.
- [15] Y. B. Zhulina and T. M. Birshtein, Polymer Science U.S.S.R. 27 , 570 (1985), ISSN 0032-3950.
- [16] H.-P. Hsu, W. Paul, and K. Binder, Macromolecular Theory and Simulations 20 , 510 (2011).
- [17] H.-P. Hsu, W. Paul, and K. Binder, Polymer Science Series C 55 , 39 (2013), ISSN 1555-614X.
- [18] H.-P. Hsu, W. Paul, and K. Binder, Europhysics Letters 92 , 28003 (2010).
- [19] S. Plimpton, Journal of Computational Physics 117 , 1 (1995), ISSN 0021-9991.
- [20] C. N. Likos, Physics Reports 348 , 267 (2001), ISSN 03701573.
- [21] X.-Y. Tu, C. Meng, X.-L. Zhang, M.-G. Jin, X.-S. Zhang, X.Z. Zhao, Y.-F. Wang, L.-W. Ma, B.-Y. Wang, M.-Z. Liu, et al., Macromolecular Bioscience 18 , 1800022 (2018).
- [22] X.-Y. Tu, C. Meng, Y.-F. Wang, L.-W. Ma, B.-Y. Wang, J.-L. He, P.-H. Ni, X.-L. Ji, M.-Z. Liu, and H. Wei, Macromolecular Rapid Communications 39 , 1870014 (2018).
- [23] R. Verduzco, X. Li, S. L. Pesek, and G. E. Stein, Chem. Soc. Rev. 44 , 2405 (2015).
- [24] W. F. M. Daniel, J. Burdynska, M. Vatankhah-Varnoosfaderani, K. Matyjaszewski, J. Paturej, M. Rubinstein, A. V. Dobrynin, and S. S. Sheiko, Nature Materials 15 , 183 (2016).
- [25] C. Mayer, F. Sciortino, C. N. Likos, P. Tartaglia, H. L¨wen, and o E. Zaccarelli, Macromolecules 42 , 423 (2009).
P. Corsi E. Roma T. Gasperi F. Bruni B. Capone
Exploiting Scaling Laws for Polymeric Bottle Brushes: a Theoretical Coarse-Graining for Homopolymeric Branched Polymers
Bottle brushes are polymeric macromolecules made of a linear polymeric
backbone grafted with side chains. The choice of the grafting density
{\sigma}g, the length ns the grafted side chains and their chemical nature
fully determines the properties of each macromolecule, such as its elasticity
and its folding behaviour. Typically, experimental bottle brushes are systems
made of tens of thousands of monomeric units, rendering a computational
approach extremely expensive, especially in the case of bottle brush solutions.
A proper coarse graining description of these macromolecules thus appears
essential. We present here a theoretical approach able to develop a general,
transferable and analytical multi-scale coarse graining of homopolymeric bottle
brush polymers under good solvent conditions. Starting from scaling theories,
each macromolecule is mapped onto a chain of tethered star polymers, whose
effective potential is known from scaling predictions, computational and
experimental validations and can be expressed as a function of the number of
arms f, and the length na of each arm. Stars are then tethered to one another
and the effective potential between them is shown to only depend on the key
parameters of the original bottle brush polymer ({\sigma}g, ns). The
generalised form of the effective potential is then used to reproduce
properties of the macromolecules obtained both with scaling theories and with
simulations. The general form of the effective potentials derived in the
current study allows a theoretical and computational description of the
properties of homopolymeric bottle brush polymers for all grafting densities
and all lengths of both backbone and grafted arms, opening the path for a
manifold of applications.