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)
, the DP of the main chain
, DPs of side chains
, and number of monomers in a spacer separating two neighboring grafting points
. The number of side chains emanating from each branching point is
. The ratio
(number of side chains per monomer unit of the main chain) quantifies grafting density. Total number of side chains in the block A is
. The set of architectural parameters for the block
A will be abbreviated as
. Block
B with total DP
is characterized by the corresponding set of parameters
. Below we focus on comb-shaped copolymers with densely grafted side chains,
, and, following standard nomenclature, term them as molecular bottlebrushes. Evidently,
or
correspond to linear respective blocks. The DPs of the blocks are expressed through the sets of architectural parameters as
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 () 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
Here, D and R are the corona thickness and radius of the core, respectively ( is overall radius of micelle), or are the distances from the core–corona interface, is the thermal energy, and .
The topological coefficients
depend on the set of topological parameters
of each respective block. Below, instead of the topological coefficients
, we use the topological ratios
which are equal to unity for linear blocks (“bare backbones”),
at
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
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
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
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
for
or
. For that, Equation (
4) can be presented as
or
Obviously, is limited from above as .
Alternatively, we can analyze parameters of micelles as a function of grafting density
for
or
, so that
or
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,
which comprises contribution of the corona,
, and the core,
, domains, as well as the free energy of the core–corona interface,
.
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
as
where
is concentration (volume fraction) profile of monomer units of block
A in the corona.
Under good solvent conditions,
where
v is dimensionless (normalized by
) second virial coefficient, and, therefore
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
where
,
and
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
is defined as
where signs “+” and “−” refer to the corona (
) and core (
) domains, respectively
The condition of uniform packing of blocks
B in the core imposes the relation between the core–corona
interface area
s per block copolymer and the radius
R of the core domain
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
where
is the volume fraction of monomer units
B in the core.
Conservation of the number of monomer units of the
A blocks,
, provides (with account of Equation (
11)) the relationship between
and
D
where
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
where
is the flux of elastic tension per unit area in the corresponding domains at distance
z from
interface.
It can be presented as
where the upper limits of the integral are equal to
D or
R for the corona and core domains, respectively, and
is given by Equation (
11), while
.
After performing integration, one gets
where
is the free energy (per molecule) in a planar brush with area
s per
A block and Equation (
3) for
was used.
The elastic free energy of block
B in the core with radius
R can be expressed using Equations (
12) and (
19) as
The excess free energy of the core–corona interface (per block copolymer chain) can be expressed as
where
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
As follows from Equation (
24), the surface area
s per block copolymer chain is expressed as a function of dimensionless ratio
The free energy of the core–corona interface is given by
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
and
By minimizing total free energy (per chain) as a function of
x
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
where
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,
.
4. Starlike and Crew-Cut Micelles: Asymptotic Results
Following standard nomenclature, we distinguish starlike,
, and crew-cut,
, 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,
and the free energy of the core–corona interface,
. 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
.
The asymptotic power law dependencies of structural and thermodynamic properties of starlike and crew-cut micelles on the set of architectural parameters of A-block and are presented below. They differ for the cases of bottlebrush, , and miktostar-like, , blocks A.
4.1. Starlike Micelles,
Aggregation number:
and core radius
.
4.2. Crew-Cut Micelles,
Aggregation number:
and core radius
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 of the main chain of the soluble block is kept constant and the number and the DP of the side chains is increasing; in the latter case D is an increasing function of .
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 . When the DP of the main chain of the soluble block is short enough (miktostar regime), all these properties become a power law function of the number 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 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 leads to an increase in 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.