Next Article in Journal
Chemical-Physical Behaviour of Microgels Made of Interpenetrating Polymer Networks of PNIPAM and Poly(acrylic Acid)
Next Article in Special Issue
Organosilica-Modified Multiblock Copolymers for Membrane Gas Separation
Previous Article in Journal
Preparation of Molded Fiber Products from Hydroxylated Lignin Compounded with Lewis Acid-Modified Fibers Its Analysis
Previous Article in Special Issue
Mikto-Arm Stars as Soft-Patchy Particles: From Building Blocks to Mesoscopic Structures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Self-Assembly of Bottlebrush Block Copolymers in Selective Solvent: Micellar Structures

by
Inna O. Lebedeva
1,
Ekaterina B. Zhulina
2 and
Oleg V. Borisov
1,2,*
1
Institut des Sciences Analytiques et de Physico-Chimie pour l’Environnement et les Matériaux, UMR 5254 CNRS UPPA, 64000 Pau, France
2
Institute of Macromolecular Compounds of the Russian Academy of Sciences, 190121 St. Petersburg, Russia
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(9), 1351; https://doi.org/10.3390/polym13091351
Submission received: 20 March 2021 / Revised: 15 April 2021 / Accepted: 16 April 2021 / Published: 21 April 2021
(This article belongs to the Special Issue Self-Assembly of Block Copolymers: Experiment and Modelling)

Abstract

:
Block copolymers comprising chemically different bottlebrush blocks can self-assemble in selective solvents giving rise to micellar-like solution nanostructures. The self-consistent field theoretical approach is used for predicting relation between architectural parameters of both bottlebrush blocks (polymerization degrees of the main and side chains, density of grafting of the side chains to the backbone) and structural properties of micelles as well as critical micelle concentration (CMC). As predicted by the theory, replacement of linear blocks by bottlebrush ones with the same degrees of polymerization results in a decrease in the micellar core size (in aggregation number) and extension of the corona, whereas the CMC increases. These theoretical findings are in good agreement with results of computer simulations.

1. Introduction

Diblock copolymers comprising chemically different blocks A and B are capable of self-assembly in selective solvent, which is good for blocks A, but poor for block B, thus giving rise to nanoscale micellar-like aggregates. In such aggregates, insoluble B blocks associate into solvent-free core domain decorated by solubilizing corona formed by solvophilic A blocks. Amphiphilic diblock copolymers with hydrophobic B block and hydrophilic (neutral or ionically charged) A blocks represent generic type of polymeric surfactants resembling in many aspects their low molecular weight counterparts [1].
Rational understanding of self-assembly of block copolymers in selective solvent has been achieved on the basis of existing theories [2,3,4,5,6,7] amply supported by experiments (see, e.g., reviews [8,9]). A generic feature of block copolymer assembly is coupling of inter-molecular association to conformational changes in both blocks which makes self-assembly of block copolymers more complex than that of conventional surfactants. Structure and properties of self-assembled aggregates can be efficiently controlled by DPs of blocks, as well as by tuning their solubility (so called stimuli-responsiveness) [6].
Moreover, as has been recently understood, by changing topology (mode and degree of branching) of the blocks, one can achieve new properties of the self-assembled nanostructures without affecting chemical nature of the constituent blocks. Linear-dendritic or double dendritic block copolymers (termed also as amphiphilic or Janus dendrimers) represent typical examples [10,11,12,13,14,15]. Such copolymers are most promising for biomedical applications due to large number of exposed to environment and potentially functionizable terminal groups of dendritically-branched solvophilic blocks [16].
Recently, self-assembled structures of block copolymers with bottlebrush-like blocks have attracted considerable interest. Amphiphilic diblock copolymers with bottlebrush soluble blocks can replace in many aspects copolymers with dendritic blocks because they comprise large number of functionizable end segments in multiple lateral branches. Moreover, triblock copolymers with middle bottlebrush blocks can associate in solution or in the melt state giving rise to materials (microphase segregated, mesogels) with unique mechanical properties resembling biological tissues and efficiently controllable through the generic set of architectural parameters, i.e., degrees of polymerization of the main and side chains and grafting density in the bottlebrush blocks.
Although the first computer simulations [17,18] proved the possibility of affecting the micellar structure by varying architecture of the bottlebrush blocks, theoretical understanding of systematic relations between the topological parameters set and experimentally relevant properties (aggregation number, overall dimensions) of micelles formed by bottlebrush block copolymers is still missing. Here, we aim to fill this important knowledge gap. The aim of the present paper is to develop a theory of self-assembly of diblock copolymers comprising one soluble and one insoluble bottlebrush blocks in selective solvent (linear-bottlebrush block copolymers represent a particular case).

2. Model and Method

We consider diblock copolymers comprising two chemically different comb-shaped blocks in selective solvent, Figure 1. Block A has total degree of polymerization (DP) N A , the DP of the main chain M A , DPs of side chains n A , and number of monomers in a spacer separating two neighboring grafting points m A . The number of side chains emanating from each branching point is q A . The ratio q A / m A (number of side chains per monomer unit of the main chain) quantifies grafting density. Total number of side chains in the block A is M A q A / m A . The set of architectural parameters for the block A will be abbreviated as { M A ; n A ; m A ; q A } . Block B with total DP N B is characterized by the corresponding set of parameters { M B ; n B ; m B ; q B } . Below we focus on comb-shaped copolymers with densely grafted side chains, q A , B n A , B / m A , B 1 , and, following standard nomenclature, term them as molecular bottlebrushes. Evidently, n A = 0 or n B = 0 correspond to linear respective blocks. The DPs of the blocks are expressed through the sets of architectural parameters as
N i = M i ( 1 + q i · n i m i ) i = A , B
The backbone and the side chains in each of the blocks are assumed to be chemically identical, the monomer unit size a (the same for both blocks), is on the order of the Kuhn segment length. The solvent is assumed to be poor for monomer units of blocks B and moderately good for monomer units of blocks A.
Poor solubility of blocks B drives association of block copolymers in selective solvent that leads to formation of micellar-like aggregates where blocks B constitute solvent-free core domain whereas soluble blocks A protrude in solution to form corona. Here, we focus on spherical micelles although depending on the DPs of blocks A and B, and their degrees of branching other morphologies (wormlike micelles, polymersomes) may correspond to thermodynamic equilibrium structures which will be considered elsewhere.
Poor solubility of blocks B assures narrow core–corona interface so that the corona and the core of micelle can be envisioned as convex swollen and concave dry brushes of A and B blocks tethered to the core–corona ( A / B ) interface, respectively.
We implement the strong stretching self-consistent field (SS-SCF) approximation for calculating the free energies of both corona (A) and core (B) domains. In the SS-SCF formalism, the self-consistent molecular potential acting on monomer units of A or B blocks is parabolic and presented [19,20,21] as
U i ( z ) k B T = 3 κ i 2 2 a 2 · D 2 z 2 , i = A R 2 z 2 , i = B
Here, D and R are the corona thickness and radius of the core, respectively ( R + D is overall radius of micelle), z [ 0 , D ] or z [ 0 , R ] are the distances from the core–corona interface, k B T is the thermal energy, and i = A , B .
The topological coefficients κ i depend on the set of topological parameters { M i , n i , m i , q i } of each respective block. Below, instead of the topological coefficients κ i , we use the topological ratios
η i = 2 N i κ i π
which are equal to unity for linear blocks (“bare backbones”), N i = M i at n i = 0 but are larger than unity for any branched, e.g., bottlebrush blocks. The topological ratio quantifies the relative increase in the elastic free energy per chain in the brush formed by branched macromolecules as compared to that in the brush formed by linear chains with the same degree of polymerization. The topological ratio depends on the macromolecular architectural parameters but is independent of the chain grafting density, and the character of intermolecular interactions in the brush.
The topological ratio for brushes of comblike macromolecules was calculated in ref. [22] and can be expressed analytically in the two asymptotic limits
η i ( 1 + q i n i / m i ) 1 / 2 , M i ( n i m i / q i ) 1 / 2 q i M i / m i , M i ( n i m i / q i ) 1 / 2
The first line in Equation (4) corresponds to the comblike block with long backbone and multiple relatively short side chains attached to it (bottlebrush limit, as depicted in Figure 1a), whereas the second line corresponds to the block with few long arms attached to the relatively short backbone (miktostar limit, depicted in Figure 1b). Although the topological coefficient (and the topological ratio) can be calculated for comb-shaped blocks with arbitrary set { M i , n i , m i , q i } of the parameters, the self-consistent potential exhibits the parabolic dependence on z, Equation (2), provided (i) Gaussian conformational elasticity of all linear segments of the branched blocks on all length scales and (ii) absence of “dead zone” proximal to A / B interface and depleted from the free ends of backbones of blocks. In convex spherical geometry the latter criterion can be violated but, nonetheless, Equation (2) with topological ratio given by Equation (4) can serve as a reasonable approximation.
In this paper, we focus on copolymers with comb-shaped blocks that exhibit linear (Gaussian) elasticity on all length scales. Although, due to repulsive interactions between densely grafted side chains, the backbones of comb-shaped polymers in the melt or/and in solution are elongated with respect to unperturbed Gaussian end-to-end dimensions [23], the application of the SS-SCF approach to micellar corona and core ignores pre-stretching of backbones in the individual A and B blocks. A more elaborated model taking into account renormalization of the blocks’ elasticity will be considered elsewhere.
Below we analyze how parameters of micelles (aggregation number, radius of the core, corona thickness, and the overall radius of the micelle) depend on the length of side chains n i for N i , m i = c o n s t or M i , m i = c o n s t . For that, Equation (4) can be presented as
η i ( n i ) ( 1 + q i n i / m i ) 1 / 2 , n i N i 2 / 3 ( m i / q i ) 1 / 3 N i / ( n i + m i / q i ) , n i N i 2 / 3 ( m i / q i ) 1 / 3
or
η i ( n i ) ( 1 + q i n i / m i ) 1 / 2 , n i M i 2 q i / m i q i M i / m i , n i M i 2 q i / m i
Obviously, n i is limited from above as n i N i m i .
Alternatively, we can analyze parameters of micelles as a function of grafting density m i / q i for N i , n i = c o n s t or M i , n i = c o n s t , so that
η i ( m i / q i ) ( 1 + q i n i / m i ) 1 / 2 , m i / q i N i 2 n i 3 N i / ( n i + m i / q i ) , m i / q i , N i 2 n i 3
or
η i ( m i / q i ) ( 1 + q i n i / m i ) 1 / 2 , m i / q i M i 2 n i 1 M i q i / m i , m i / q i M i 2 n i 1
The first and the second lines in Equations (5)–(8) refer to bottlebrush or miktostar shape of a block, respectively.

3. Free Energy of a Micelle

The properties of equilibrium micelles are obtained through minimization of the free energy per block copolymer molecule,
F = F c o r o n a + F c o r e + F i n t e r f a c e
which comprises contribution of the corona, F c o r o n a , and the core, F c o r e , domains, as well as the free energy of the core–corona interface, F i n t e r f a c e .
Within SS-SCF approximation, the self-consistent molecular potential acting on the monomers of corona chains is given by Equation (2), and related to the interaction free energy density f i n t , A { ϕ A ( z ) } as
U ( z ) = δ f i n t , A { ϕ A ( z ) } δ ϕ A ( z )
where ϕ A ( z ) is concentration (volume fraction) profile of monomer units of block A in the corona.
Under good solvent conditions,
f i n t , A { ϕ A ( z ) } = k B T a 3 v ϕ 2 ( z )
where v is dimensionless (normalized by a 3 ) second virial coefficient, and, therefore
ϕ A ( z ) = 3 κ A 2 4 v a 2 ( D 2 z 2 )
also exhibits a parabolic shape.
Detailed calculations of the free energies of corona and core domains in a micelle formed by copolymers with arbitrary branched architecture of blocks are presented elsewhere [15] and here we outline only the key steps.
The cumulative free energy of the corona and the core domains can be presented as
F c o r o n a + F c o r e =
0 D f i n t , A { ϕ A ( z ) } s A ( z ) d z + 0 D f e l a s t i c , A ( z ) s A ( z ) d z + 0 R f e l a s t i c , B ( z ) s B ( z ) d z
where f i n t , A { ϕ A ( z ) } , f e l a s t i c , A ( z ) and f e l a s t i c , B ( z ) are the free energy of excluded volume interactions, the elastic (conformational) free energy per unit volume in the corona, and of the core, respectively. Here s is the core surface area per block copolymer molecule and s A , B ( z ) is defined as
s i ( z ) = s · R ± z R 2
where signs “+” and “−” refer to the corona ( i = A ) and core ( i = B ) domains, respectively
The condition of uniform packing of blocks B in the core imposes the relation between the core–corona A / B interface area s per block copolymer and the radius R of the core domain
s = 3 N B a 3 R ϕ B ,
The number p of block copolymer chains in a spherical micelles is related to the radius R of the core or to the core surface area s per copolymer chain as
p = 4 π R 3 ϕ B 3 N B = 36 π N B 2 a 6 s 3 ϕ B 2
where ϕ B 1 is the volume fraction of monomer units B in the core.
Conservation of the number of monomer units of the A blocks, 0 D A ( z ) s A ( z ) d z = N A a 3 , provides (with account of Equation (11)) the relationship between s , R and D
D = D 0 ( s ) · ( 1 + 3 4 D R + 1 5 D 2 R 2 ) 1 / 3
where
D 0 ( s ) = 2 N A v a 2 s κ A 1 / 3 a ( 8 π 2 ) 1 / 3 N A v 1 / 3 ( a 2 s ) 1 / 3 η A 2 / 3 a
is the thickness of a planar brush (with zero curvature) with grafting area s per molecule.
As long as both blocks A and B exhibit Gaussian elasticity, the density of elastic free energy can be expressed as
f e l a s t i c , i ( z ) = T i ( z ) 2
where T i ( z ) is the flux of elastic tension per unit area in the corresponding domains at distance z from A / B interface.
It can be presented as
T i ( z ) k B T = 3 κ i 2 2 s i ( z ) a 3 z D , R z ϕ i ( z ) s i ( z ) d z
where the upper limits of the integral are equal to D or R for the corona and core domains, respectively, and ϕ A ( z ) is given by Equation (11), while ϕ B ( z ) c o n s t 1 .
After performing integration, one gets
F c o r o n a = F c o r o n a ( 0 ) ( s ) D D 0 ( s ) 5 · 1 + 5 6 D R + 5 21 D 2 R 2
where
F c o r o n a ( 0 ) ( s ) = k B T 9 20 κ A 4 s v a 2 D 0 ( s ) a 5 = k B T 9 π 2 / 3 10 N A v 2 / 3 ( a 2 s ) 2 / 3 η A 2 / 3
is the free energy (per molecule) in a planar brush with area s per A block and Equation (3) for η A was used.
The elastic free energy of block B in the core with radius R can be expressed using Equations (12) and (19) as
F c o r e ( s ) = 0 R f e l a s t i c , B ( z ) s B ( z ) d z = k B T 27 π 2 80 N B a 4 s 2 ϕ B 2 η B 2
The excess free energy of the core–corona interface (per block copolymer chain) can be expressed as
F i n t e r f a c e = k B T γ s
where γ k B T is the interfacial free energy per unit area, which is controlled primarily by solubility parameter of the core-forming B-blocks but is virtually independent of concentration of monomer units of A-blocks in the corona domain.
Using Equations (14), (16), and (17), the area s of the core surface per block copolymer chain can be presented as
s a 2 = π η N A 3 / 2 v 1 / 2 N B ϕ 3 / 2 · ( 3 D 2 R ) 3 / 2 ( 1 + 3 4 D R + 1 5 D 2 R 2 ) 1 / 2
As follows from Equation (24), the surface area s per block copolymer chain is expressed as a function of dimensionless ratio
x = D / R
The free energy of the core–corona interface is given by
F i n t e r f a c e ( x ) = k B T γ s a 2 = k B T γ π η N A 3 / 2 v 1 / 2 N B ϕ 3 / 2 ( 3 x 2 ) 3 / 2 ( 1 + 3 4 x + 1 5 x 2 ) 1 / 2 ,
By substituting Equation (24) for s in Equations (12) and (21) for corona and core, contributions to the free energy can be also expressed as functions of x
F c o r o n a ( x ) k B T = 3 5 v N A 2 ϕ N B ( 1 + 5 6 x + 5 21 x 2 ) x ( 1 + 3 4 x + 1 5 x 2 ) 2
and
F c o r e ( x ) k B T = N A 3 v ϕ N B 2 η B η A 2 · 1 10 x 3 ( 1 + 3 4 x + 1 5 x 2 ) 1
By minimizing total free energy (per chain) as a function of x
d [ F c o r o n a ( x ) + F i n t e r f a c e ( x ) + F c o r e ( x ) ] d x = 0
which is equivalent to the minimization with respect to s, one finds the value of s (or R) corresponding to equilibrium micelle and, using Equations (9), (26)–(28) the free energy per block copolymer chain in the equilibrium micelle. This free energy is directly related to CMC, as
k B T ln C M C = F ( s ) F 1
where
F 1 k B T γ ( N B / ϕ B ) 2 / 3
is (with the accuracy of numerical prefactor) the free energy of a unimer which is dominated by the excess free energy of the interface between single collapsed B-block and surrounding (poor) solvent. We remark that Equation (30) disregards translational entropy of micelles which is justified for sufficiently large aggregation numbers, p 1 .

4. Starlike and Crew-Cut Micelles: Asymptotic Results

Following standard nomenclature, we distinguish starlike, D R , and crew-cut, D R , micelles (Figure 2). The overall dimensions of the micelles in these limiting cases are controlled either by the extension of the corona, D, or by the core radius, R, respectively. In both limits, structural and thermodynamic properties of micelles can be derived analytically in the form of power-law dependencies by keeping only the dominant terms in the free energy, that is, the free energy of the corona, F c o r o n a and the free energy of the core–corona interface, F i n t e r f a c e . Neglecting conformational entropy of the core-forming B-blocks does not allow us to capture the effect of their architecture on the properties of micelles, which may serve as a reasonable approximation as long as the B-blocks are linear or weakly branched. As demonstrated below, the calculations with proper account of the contribution of the conformational entropy of the B-blocks, prove that increasing branching of the B-blocks leads to systematic decrease in the dimensions and aggregation number of micelles. However, this decrease cannot be described in terms of power law dependencies on η B .
The asymptotic power law dependencies of structural and thermodynamic properties of starlike and crew-cut micelles on the set { M A ; n A ; m A ; q A } of architectural parameters of A-block and N B / ϕ B are presented below. They differ for the cases of bottlebrush, M A ( n A m A / q A ) 1 / 2 , and miktostar-like, M A ( n A m A / q A ) 1 / 2 , blocks A.

4.1. Starlike Micelles, D R

Aggregation number:
p γ 15 / 11 ( N B ϕ B ) 10 / 11 v 6 / 11 ·
N A 3 / 11 ( 1 + q A n A m A ) 9 / 11 = M A 3 / 11 ( 1 + q A n A m A ) 12 / 11 , M A ( n A m A / q A ) 1 / 2 N A 21 / 11 n A 18 / 11 ( 1 + m A q A n A ) 18 / 11 = ( M A q A m A ) 21 / 11 n A 3 / 11 ( 1 + m A q A n A ) 3 / 11 , M A ( n A m A / q A ) 1 / 2
and core radius R / a N B ϕ B 1 / 3 p 1 / 3 .
Corona thickness:
D / a γ 3 / 11 ( N B ϕ B ) 2 / 11 v 1 / 11
N A 6 / 11 ( 1 + q A n A m A ) 4 / 11 = M A 6 / 11 ( 1 + q A n A m A ) 2 / 11 , M A ( n A m A / q A ) 1 / 2 N A 2 / 11 n A 8 / 11 ( 1 + m A q A n A ) 8 / 11 = ( q A M A m A ) 2 / 11 n A 6 / 11 ( 1 + m A q A n A ) 6 / 11 , M A ( n A m A / q A ) 1 / 2
The CMC:
ln C M C γ N B ϕ B 2 / 3 + γ 6 / 11 N B ϕ B 4 / 11 v 2 / 11 ·
N A 1 / 11 ( 1 + q A n A m A ) 3 / 11 = M A 1 / 11 ( 1 + q A n A m A ) 4 / 11 , M A ( n A m A / q A ) 1 / 2 N A 7 / 11 n A 6 / 11 ( 1 + m A q A n A ) 6 / 11 = ( q A M A m A ) 7 / 11 n A 1 / 11 ( 1 + m A q A n A ) 1 / 11 , M A ( n A m A / q A ) 1 / 2

4.2. Crew-Cut Micelles, D R

Aggregation number:
p γ 9 / 5 ( N B ϕ B ) 2 v 6 / 5 ·
N A 9 / 5 ( 1 + q A n A m A ) 3 / 5 = M A 9 / 5 ( 1 + q A n A m A ) 12 / 5 , M A ( n A m A / q A ) 1 / 2 N A 3 n A 6 / 5 ( 1 + m A q A n A ) 6 / 5 = ( M A q A m A ) 3 n A 9 / 5 ( 1 + m A q A n A ) 9 / 5 , M A ( n A m A / q A ) 1 / 2
and core radius R / a N B ϕ B 1 / 3 p 1 / 3
Corona thickness:
D / a γ 1 / 5 v 1 / 5 ·
N A 4 / 5 ( 1 + q A n A m A ) 2 / 5 = M A 4 / 5 ( 1 + q A n A m A ) 2 / 5 , M A ( n A m A / q A ) 1 / 2 n A 4 / 5 ( 1 + m A q A n A ) 4 / 5 , M A ( n A m A / q A ) 1 / 2
The CMC:
ln C M C γ N B ϕ B 2 / 3 + γ 2 / 5 v 2 / 5 ·
N A 3 / 5 ( 1 + q A n A m A ) 1 / 5 = M A 3 / 5 ( 1 + q A n A m A ) 4 / 5 , M A ( n A m A / q A ) 1 / 2 N A n A 2 / 5 ( 1 + m A q A n A ) 2 / 5 = q A M A m A n A 3 / 5 ( 1 + m A q A n A ) 3 / 5 , M A ( n A m A / q A ) 1 / 2

5. Discussion

5.1. Diagram of States

First of all, it is instructive to outline the ranges of stability of starlike ( s t ) and crew-cut ( c c ) micelles with corona formed by copolymers with either bottlebrush ( b b ) or miktostar ( m k s ) solvophilic blocks A. The corresponding diagrams of states are presented in ( N A , n A ) and ( M A , n A ) coordinates in Figure 3a,b, respectively. The diagrams contain regimes b b / s t , m k s / s t , b b / c c , m k s / c c of starlike or crew-cut micelles with bottlebrush and miktostar corona-forming A-blocks, respectively. The diagrams contain also the unimer regime where micellization is suppressed due to strong repulsive interactions between A-blocks. An increase in the overall DP of soluble blocks ultimately provokes stabilization of the unimer state of the block copolymers in the solution.

5.2. Structural and Thermodynamic Properties of Micelles

Below we discuss the main trends in the dependencies of the micellar properties (aggregation number, core, corona and overall dimensions, and CMC) on architectural parameters of the blocks. These trends are qualitatively captured by asymptotic Equations (32)–(37) whereas the explicit dependencies are accurately calculated at arbitrary ratio R / D according to the routine described in Section 3.
We start with analysis of the effect of varied side chains length, n A , in soluble block A at either N A , m A / q A = c o n s t or M A , m A / q A = c o n s t . An increase in n A (at N A = c o n s t or M A = c o n s t ) provokes transformation of the A-block from bottlebrush to miktostar-like and may also lead (particularly at M A = c o n s t ) to change of the shape of the micelles from crew-cut to starlike ones.
As illustrated by Figure 4, corresponding to cross-section (1) of the diagram in Figure 3a, at N A , m A / q A = c o n s t the aggregation number, core radius, and corona thickness in starlike micelles with bottlebrush A-block decrease as a function of n A . If not the total DP of the A-block, but rather DP of the main chain M A is kept constant, the aggregation number p and core radius R also decrease, whereas D and overall micellar size increase as a function of n A . In both cases ( N A = c o n s t or M A = c o n s t ) CMC is an increasing function of n A .
Figure 5 demonstrates more complex, non-monotonous dependencies of the micellar parameters on n A at N A = c o n s t , m A / q A = c o n s t corresponding to cross-section (2) of the diagram in Figure 3a. Here aggregation number, core radius and corona thickness decrease as a function of n A as long as A-block is bottlebrush-like, but increase as a function of n A when A-block acquires miktostar shape. Hence, p , R , D , and R + D pass through a minimum as a function of n A when the shape of the coronal block A changes from bottlebrush to miktostar. In Figure 5, this crossover occurs in the regime of crew-cut micelles ( R D ) , but similar trends are predicted for starlike micelles (cross-section (3) in Figure 3a). The CMC also depends non-monotonously on n A , increasing in the range of bottlebrush A-blocks and decreasing in the range of miktostar A-blocks.
In Figure 6, the micellar structural properties p , R , D , R + D are plotted as a function of q A n A / m A at N A = c o n s t for a number of selected values of q A / m A . As one can see from the figure, all the dependencies follow universal curves as long as A-blocks keep bottlebrush shape but split into a series of curves each corresponding to specific value of q A / m A at large n A when the A-block acquires miktostar shape.
The effect of DP M A of the backbone in the corona block A (at constant n A , m A / q A ) on the aggregation number and dimensions of micelles is as follows: As predicted by Equations (32), (33), (35), and (36) for both starlike and crew-cut micelles p and R do decrease as a function of M A while CMC increases. For starlike micelles, the corona thickness D and the overall micellar sise R + D increase as a function of M A in the case of bottlebrush A-block, but decrease in the case of miktoarm A-block. For crew-cut micelles, the corona thickness increases as a function of M A in the case of bottlebrush A-block, but is virtually independent of M A when A-block acquires miktostar shape. Hence, the corona thickness and the overall radius of the micelle may exhibit non-monotonic behavior as a function of M A passing through a minimum at the crossover between bottlebrush and miktostar regimes for the coronal block A, as illustrated by Figure 7 for crew-cut micelles.
The effect of branching of the core-forming block B on the structural properties of spherical micelles is less pronounced, as illustrated in Figure 8. Although asymptotic Equations (32)–(37) neglecting conformational entropy of the core-forming blocks do not predict any dependence on the branching parameters of the core-forming B-blocks, detailed calculation with the account of the contribution of the conformational entropy of the core-forming block indicate that increasing branching of block B leads to a decrease in the micelle aggregation number and an increase in CMC, which is more pronounced in the crew-cut domain.

5.3. Comparison to Molecular Dynamics Simulations

Recent coarse-grained MD simulations [18] of A B diblock self-assembly in solution demonstrated how aggregation number p of micelles depends on the architectural parameters of comb-shaped blocks. Diblock copolymers with symmetric composition A:B = 50:50 were studied at fixed total diblock DP N = 96 , equal and variable lengths of side chains in both blocks, n i = n , and equal lengths of spacers, m i = 1 , in the backbone with DP M = N / ( 1 + n ) . A set of 9 different block copolymer architectures was constructed via changing n, ranging from bottlebrush to starlike macromolecules. The average parameters of micelles, formed with inferior solvent strength for B-block, were presented as a function of branching parameter, g = ( R b r 2 ¯ ) / ( R l i n 2 ¯ ) , specified as the ratio of gyration radii, R 2 ¯ , of the branched (subscript b r ), and linear (subscript l i n ) individual chains with the same DP in a good solvent.
To relate g to molecular parameters in our model, we use the asymptotic power law dependencies for the average end-to-end distance (expressed in units of monomer length) for bottlebrush and starlike polymers with spacer length m 1 under good solvent conditions. In the mean filed approximation, the average size of macromolecule, R b r , is given by
R b r n 2 / 5 M 3 / 5 b o t t l e b r u s h p o l y m e r n 3 / 5 M 1 / 5 s t a r l i k e p o l y m e r
The first line in Equation (38) is obtained by balancing the elastic free energy of the stretched bottlebrush backbone (∼ R 2 / M ) with the free energy of repulsive binary interactions (∼ ( M n ) 2 / R 3 ). For a starlike polymer, the elastic free energy accounts for stretching of M side chains (∼ M R 2 / n ) while the contribution of binary interactions remains the same as for bottlebrush polymer (∼ ( M n ) 2 / R 3 ). The crossover between two asymptotes in Equation (38) occurs at M = M * n 1 / 2 .
For a linear chain with N monomer units,
R l i n N 3 / 5 = M 3 / 5 ( 1 + n ) 3 / 5 M 3 / 5 n 3 / 5
The theoretical branching ratio g is then specified as
g = R b r R l i n 2 n 2 / 5 b o t t l e b r u s h p o l y m e r M 4 / 5 = ( n / N ) 4 / 5 s t a r l i k e p o l y m e r
with the accuracy of prefactors on the order of unity. As it follows from Equation (40), an increase in n at fixed N leads to a non-monotonic dependence for g ( n ) , in agreement with the simulation data (Figure 2g in Ref. [18]). The crossover between asymptotes for bottlebrush and starlike polymers occurs at
n * N 2 / 3 ; g * = g ( n * ) N 4 / 15
For N = 96 , and both prefactors equal to unity in Equation (40), n * 21 and g * 0.30 .
In selective solvent, A B copolymer with equal DPs of linear blocks ( n i = 0 ) gives rise to starlike micelles in which the core of insoluble B-blocks is decorated by the corona of soluble A-blocks. As it follows from the diagram of states in Figure 3, branching of the blocks, i.e., an increase in n A at fixed value of N A , leads to the non-monotonic variation in micelle aggregation number, p. According to Equation (32) in spherical micelles formed by diblock copolymer with m A / q A 1 , and N B = N A = N / 2 , the average aggregation number varies as
p N B 10 / 11 · N A 3 / 11 ( 1 + n A ) 9 / 11 b b / s t N A 21 / 11 n A 18 / 11 m k s / s t = ( N / 2 ) 7 / 11 ( 1 + n A ) 9 / 11 b b / s t ( N / 2 ) 1 n A 18 / 11 m k s / s t
That is, aggregation number decreases from p ( n A = 0 ) = p l i n ( N / 2 ) 7 / 11 in b b / s t regime, reaches a minimum at the boundary between b b / s t and m k s / s t regimes ( n A N 2 / 3 ), and subsequently increases back to p l i n at n A = N / 2 in m k s / s t regime. By substituting g ( n ) from Equation (40) in Equation (41), one finds
p p l i n g 45 / 22 b b / s t g 45 / 22 m k s / s t
Therefore, according to Equation (42), aggregation number p in micelles formed by copolymer with fixed DP N and bottlebrush blocks is expected to decrease upon an increase in n (due to decreasing g, Equation (40)), while for block copolymer with starlike blocks, p is expected to increase with n (and g ) , making the dependence p ( g ) looplike. The loop originates at g = 1 , and turns backwards at g = g * . In the mean field approximation, the exponent in Equation (42) is the same for bottlebrush and starlike polymers, however, the omitted numerical coefficients and nonpower law dependencies could be different.
The predicted trends are in agreement with the results of MD simulations [18]. The latter demonstrated clear difference in the behavior of A B copolymers with bottlebrush (samples 1–6) and starlike (samples 7–9) blocks. In accordance with Equations (40) and (42), the data for micelle aggregation number < N a g g > formed a loop as a function of branching parameter g (Figure 3e in Ref. [18]). The crossover value of n * 21 (estimated with accuracy of prefactors on the order of unity) was between DPs N s c = 15 and N s c = 23 in starlike samples 8 and 9 (Figure 3a in Ref. [18]). The crossover value of g * 0.30 was close to g 0.3 , corresponding to reverse in < N a g g > ( g ) dependence in MD simulations. However, evaluation of exponents in the asymptotic theoretical dependencies requires more extensive simulations with longer chains, and remains a challenge at the moment.

6. Conclusions

To summarize, we have developed a theory of micellization of diblock copolymers comprising chemically different comb-shaped (bottlebrush) blocks in selective solvents. This theory enables predicting how DPs of the main and side chains and grafting densities in both soluble and insoluble blocks of the copolymers affect aggregation number and the equilibrium dimensions of self-assembled micelles and critical micellar concentration. Both limits of blocks with long main chain and multiple short side chains (bottlebrush) and short main chain and a few long side chains (miktostar) are considered. Asymptotic analytical dependencies are derived for the limiting cases of starlike and crew-cut micelles, whereas full solution is obtained for arbitrary ratios between the micellar core and corona dimensions.
We have demonstrated that the replacement of the linear soluble block by a bottlebrush one with the same degree of polymerization results in a decrease in the aggregation number and dimensions (core size, corona thickness, the overall radius) of micelles. A similar and even stronger effect on p and R is predicted when the DP M A of the main chain of the soluble block is kept constant and the number and the DP n A of the side chains is increasing; in the latter case D is an increasing function of n A .
As long as the main chain of the A-block is sufficiently long (the A-block has the bottlebrush shape), all structural properties of the micelles and CMC follow approximately power law dependencies on effective grafting density q A n A / m A . When the DP M A of the main chain of the soluble block is short enough (miktostar regime), all these properties become a power law function of the number M A q A / m A of side chains in A-block.
The described effects of architecture of the soluble blocks on the properties of micelles are in full agreement with the trends observed in Molecular Dynamics simulations in Refs. [17,18] and can be explained by enhanced repulsive interactions and larger conformational entropy penalty for stretching of the comb-shaped blocks in the micellar corona, as compared to those for linear blocks. Experimental validation of these theoretical predictions would require comprehensive study of micellization of block copolymers with systematically varied DPs of the main and side chains in the soluble blocks
Replacement of the insoluble linear block by comb-shaped ones with the same DP N B results in a weak decrease in the aggregation number with concomitant decrease in the micellar size and increase in CMC. However, these dependencies are not described by power law functions. These effects are due to larger conformational entropy penalty for radial stretching of branched insoluble block in the micellar core. On the contrary, increasing DP and grafting density of side chains in the insoluble block with constant length of the main chain M B leads to an increase in N B and, as a consequence, increasing aggregation number and decreasing CMC.
Importantly, in the present study, we used a mean-field approximation, which is justified for semi-dilute solutions of semi-flexible polymers, that is, under the conditions that Kuhn segments of the main and side chains of both blocks are larger than the monomer unit size [24,25]. A more refined scaling analysis [7] may lead to different values of power law exponents for different properties of the block copolymer micelles. Such analysis, however, is beyond the scope of the present paper.
Finally, we remark that in the present study we considered only spherical micelles. By adjusting the sets of architectural parameters of the blocks, nano-aggregates of other morphologies (i.e., wormlike micelles, polymersomes) can be obtained as thermodynamically equilibrium structures. We shall address polymorphism of nanostructures formed by diblock copolymers with bottlebrush blocks in our forthcoming publication.

Author Contributions

All the authors equally contributed to formulation of the problem, elaboration of the computational method, analytical and numerical calculations and article writing. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Russian Science Foundation, grant 20-13-00270.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This work was supported by Russian Science Foundation, grant 20-13-00270.

Conflicts of Interest

The authors declear no conflict of interests.

References

  1. Mai, Y.; Eisenberg, E. Self assembly of block copolymers. Chem. Soc. Rev. 2012, 41, 5969–5985. [Google Scholar] [CrossRef] [PubMed]
  2. De Gennes, P.-J. Macromolecules and Liquid Crystals: Reflections on Certain Lines of Research. In Solid State Physics; Academic Press: New York, NY, USA, 1978; pp. 1–17. [Google Scholar]
  3. Halperin, A. Polymeric micelles: A star model. Macromolecules 1987, 20, 2943–2946. [Google Scholar] [CrossRef]
  4. Zhulina, Y.B.; Birshtein, T.M. Conformations of block copolymer molecules in selective solvents (micellar structures). Polym. Sci. USSR 1985, 27, 570–578. [Google Scholar] [CrossRef]
  5. Zhulina, E.B.; Borisov, O.V. Theory of block copolymer micelles: Recent advances and current challenges. Macromolecules 2012, 45, 4229–4240. [Google Scholar] [CrossRef]
  6. Borisov, O.V.; Zhulina, E.B.; Leermakers, F.A.M.; Müller, A.H.E. Self-assembled structures of amphiphilic ionic block copolymers: Theory, self-consistent field modelling and experiment. Adv. Polym. Sci. 2011, 241, 57–129. [Google Scholar]
  7. Zhulina, E.B.; Adam, M.; Sheiko, S.; LaRue, I.; Rubinstein, M. Diblock copolymer micelles in a dilute solution. Macromolecules 2005, 38, 5330. [Google Scholar] [CrossRef]
  8. Schacher, F.H.; Rupar, P.A.; Manners, I. Functional block copolymers: Nanostructured materials with emerging applications. Angew. Chem. Int. Ed. 2012, 41, 5969–5985. [Google Scholar] [CrossRef]
  9. Tritschler, U.; Pearce, S.; Gwyther, J.; Whittell, G.R.; Manners, I. 50th Anniversary Perspective: Functional Nanoparticles from the Solution Self-Assembly of Block Copolymers. Macromolecules 2017, 50, 3439–3463. [Google Scholar] [CrossRef] [Green Version]
  10. Wurm, F.; Frey, H. Linear-dendritic block copolymers: The state of the art and exciting perspectives. Progress. Polym. Sci. 2011, 36, 1–52. [Google Scholar] [CrossRef]
  11. Blasco, E.; Pinol, M.; Oriol, L. Responsive Linear-Dendritic Block Copolymers. Macromol. Rapid Commun. 2014, 35, 1090–1115. [Google Scholar] [CrossRef]
  12. Garcia-Juan, H.; Nogales, A.; Blasco, E.; Martinez, J.C.; Sics, I.; Ezquerra, T.A.; Pinol, M.; Oriol, L. Self-assembly of thermo and light responsive amphiphilic linear dendritic block copolymers. Eur. Polym. J. 2016, 81, 621–633. [Google Scholar] [CrossRef] [Green Version]
  13. Mirsharghi, S.; Knudsen, K.D.; Bagherifam, S.; Niström, B.; Boas, U. Preparation and self-assembly of amphiphilic polylysine dendrons. New J. Chem. 2016, 40, 3597–3611. [Google Scholar] [CrossRef] [Green Version]
  14. Lebedeva, I.O.; Zhulina, E.B.; Borisov, O.V. Theory of Linear-Dendritic Block Copolymer Micelles. ACS Macro Lett. 2018, 7, 811–816. [Google Scholar] [CrossRef]
  15. Lebedeva, I.O.; Zhulina, E.B.; Borisov, O.V. Self-assembly of linear-dendritic and double-dendritic block copolymers: From spherical micelles to dendrimersomes. Macromolecules 2019, 52, 3655–3667. [Google Scholar] [CrossRef]
  16. Fan, X.; Zhao, Y.; Xu, W.; Li, L. Linear-Dendritic Block Copolymer for Drug and Gene Delivery. Mater. Sci. Eng. C 2016, 62, 943–959. [Google Scholar] [CrossRef]
  17. 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]
  18. 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]
  19. Pickett, G.T. Classical Path Analysis of end-Grafted Dendrimers: Dendrimer Forest. Macromolecules 2001, 34, 8784–8791. [Google Scholar] [CrossRef]
  20. Zook, T.C.; Pickett, G.T. Hollow-Core Dendrimers Revised. Phys. Rev. Lett. 2003, 90, 015502. [Google Scholar] [CrossRef] [Green Version]
  21. Zhulina, E.B.; Leermakers, F.A.M.; Borisov, O.V. Ideal mixing in multicomponent brushes of branched macromolecules. Macromolecules 2015, 48, 5614–5622. [Google Scholar] [CrossRef]
  22. Mikhailov, I.V.; Zhulina, E.B.; Borisov, O.V. Brushes and lamellar mesophases of comb-shaped (co)polymers: A self-consistent field theory. Phys. Chem. Chem. Phys. 2020, 22, 23358–23398. [Google Scholar] [CrossRef]
  23. Zhulina, E.B.; Sheiko, S.S.; Borisov, O.V. Solution and melts of barbwire bottlebrushes: Hierarchical structure and scale-dependent elasticity. Macromolecules 2019, 52, 1671–1684. [Google Scholar] [CrossRef]
  24. Birshtein, T.M. Diagram of state for a solution of semi-rigid macromolecules. Polym. Sci. USSR 1982, 24, 2416–2422. [Google Scholar] [CrossRef]
  25. Shaefer, D.W.; Joanny, J.-F.; Pincus, P. Dynamics of Semiflexible Polymers in Solution. Macromolecules 1980, 13, 1280–1289. [Google Scholar] [CrossRef]
Figure 1. Schematics of bottlebrush (a) and miktoarm star (b) diblock copolymer for a particular case of branching activity q A = q B = 1 .
Figure 1. Schematics of bottlebrush (a) and miktoarm star (b) diblock copolymer for a particular case of branching activity q A = q B = 1 .
Polymers 13 01351 g001
Figure 2. Schematics of starlike, D R , and crew-cut, D R , micelles formed by copolymers with bottlebrush blocks.
Figure 2. Schematics of starlike, D R , and crew-cut, D R , micelles formed by copolymers with bottlebrush blocks.
Polymers 13 01351 g002
Figure 3. Schematic diagram of states of micelles in ( N A , n A ) (a) and ( M A , n A ) (b) coordinates. Regions of starlike (st) and crew-cut (cc) micelles with bottlebrush (bb) or miktostar (mks) A-blocks are outlined. Equations for boundaries between different regimes are presented near the lines. Parameters of the boundaries in the diagram are C ( γ 2 ( N B / ϕ B ) 5 v 3 ) 1 / 7 ; U = γ 5 ( N B / ϕ B ) 10 / 3 v 2 . Gray area corresponds to the region of spherical micelles instability (transition to wormlike micelles and polymersomes).
Figure 3. Schematic diagram of states of micelles in ( N A , n A ) (a) and ( M A , n A ) (b) coordinates. Regions of starlike (st) and crew-cut (cc) micelles with bottlebrush (bb) or miktostar (mks) A-blocks are outlined. Equations for boundaries between different regimes are presented near the lines. Parameters of the boundaries in the diagram are C ( γ 2 ( N B / ϕ B ) 5 v 3 ) 1 / 7 ; U = γ 5 ( N B / ϕ B ) 10 / 3 v 2 . Gray area corresponds to the region of spherical micelles instability (transition to wormlike micelles and polymersomes).
Polymers 13 01351 g003
Figure 4. Aggregation number p, core radius R (dashed lines), and corona thickness D (solid lines), and overall radius, R + D , of starlike spherical micelles formed by block copolymers with N B = 200 and N A = 1000 or M A = 1000 as a function of n A (cross-section (1) in Figure 3a). Red and blue curves correspond to N A = 1000 = c o n s t (varied n A , M A ) or to M A = 1000 = c o n s t (varied n A , N A ), respectively. Other parameters are q A = 1 , m A = 2 , γ = 1 , v = 1 , ϕ B = 1 .
Figure 4. Aggregation number p, core radius R (dashed lines), and corona thickness D (solid lines), and overall radius, R + D , of starlike spherical micelles formed by block copolymers with N B = 200 and N A = 1000 or M A = 1000 as a function of n A (cross-section (1) in Figure 3a). Red and blue curves correspond to N A = 1000 = c o n s t (varied n A , M A ) or to M A = 1000 = c o n s t (varied n A , N A ), respectively. Other parameters are q A = 1 , m A = 2 , γ = 1 , v = 1 , ϕ B = 1 .
Polymers 13 01351 g004
Figure 5. Aggregation number p, core radius R (dashed lines), and corona thickness D (solid lines), and overall radius, R + D , of micelles formed by block copolymers with N B = 150 and N A = 100 as a function of n A (cross-section (2) in Figure 3a). m A is varied, as indicated in the figure. Other parameters are q A = 1 , γ = 1 , v = 1 , ϕ B = 1 .
Figure 5. Aggregation number p, core radius R (dashed lines), and corona thickness D (solid lines), and overall radius, R + D , of micelles formed by block copolymers with N B = 150 and N A = 100 as a function of n A (cross-section (2) in Figure 3a). m A is varied, as indicated in the figure. Other parameters are q A = 1 , γ = 1 , v = 1 , ϕ B = 1 .
Polymers 13 01351 g005
Figure 6. Aggregation number p, core radius R (dashed lines), and corona thickness D (solid lines), and overall radius, R + D , of micelles formed by block copolymers with N B = 150 and N A = 100 as a function of n A q A / m A . m A is varied, as indicated in the figure. Other parameters are q A = 1 , γ = 1 , v = 1 , ϕ B = 1 . The crossover between bottlebrush and miktostar shapes of the A-block is indicated by vertical dashed lines in the middle panel.
Figure 6. Aggregation number p, core radius R (dashed lines), and corona thickness D (solid lines), and overall radius, R + D , of micelles formed by block copolymers with N B = 150 and N A = 100 as a function of n A q A / m A . m A is varied, as indicated in the figure. Other parameters are q A = 1 , γ = 1 , v = 1 , ϕ B = 1 . The crossover between bottlebrush and miktostar shapes of the A-block is indicated by vertical dashed lines in the middle panel.
Polymers 13 01351 g006
Figure 7. Aggregation number p, core radius, R (red), corona thickness, D (blue) and overall radius, R + D , of spherical micelles formed by block copolymers with N B = 200 as a function of M A at n A = 100 . Other parameters are q A = 1 , γ = 1 , v = 1 , ϕ B = 1 .
Figure 7. Aggregation number p, core radius, R (red), corona thickness, D (blue) and overall radius, R + D , of spherical micelles formed by block copolymers with N B = 200 as a function of M A at n A = 100 . Other parameters are q A = 1 , γ = 1 , v = 1 , ϕ B = 1 .
Polymers 13 01351 g007
Figure 8. Aggregation number p, core radius, R (dashed lines), corona thickness, D (solid lines), and overall dimensions R + D as a function of the DP n B of side chains in the core forming block at N B = c o n s t = 1000 . m B is varied, as indicated in the figure. Other parameters are q A = q B = 1 N A = 300 , m A = 4 , n A = 10 γ = 1 , v = 1 , ϕ = 1 .
Figure 8. Aggregation number p, core radius, R (dashed lines), corona thickness, D (solid lines), and overall dimensions R + D as a function of the DP n B of side chains in the core forming block at N B = c o n s t = 1000 . m B is varied, as indicated in the figure. Other parameters are q A = q B = 1 N A = 300 , m A = 4 , n A = 10 γ = 1 , v = 1 , ϕ = 1 .
Polymers 13 01351 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lebedeva, I.O.; Zhulina, E.B.; Borisov, O.V. Self-Assembly of Bottlebrush Block Copolymers in Selective Solvent: Micellar Structures. Polymers 2021, 13, 1351. https://doi.org/10.3390/polym13091351

AMA Style

Lebedeva IO, Zhulina EB, Borisov OV. Self-Assembly of Bottlebrush Block Copolymers in Selective Solvent: Micellar Structures. Polymers. 2021; 13(9):1351. https://doi.org/10.3390/polym13091351

Chicago/Turabian Style

Lebedeva, Inna O., Ekaterina B. Zhulina, and Oleg V. Borisov. 2021. "Self-Assembly of Bottlebrush Block Copolymers in Selective Solvent: Micellar Structures" Polymers 13, no. 9: 1351. https://doi.org/10.3390/polym13091351

APA Style

Lebedeva, I. O., Zhulina, E. B., & Borisov, O. V. (2021). Self-Assembly of Bottlebrush Block Copolymers in Selective Solvent: Micellar Structures. Polymers, 13(9), 1351. https://doi.org/10.3390/polym13091351

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