Disk-shaped bicelles in block copolymer/homopolymer blends

Disk-shaped bicelles in block copolymer/homopolymer blends

M. J. Greenall Institute of Mathematics, Physics and Computer Science, Physical Sciences Building, Aberystwyth University, Aberystwyth SY23 3BZ, United Kingdom. School of Mathematics and Physics, University of Lincoln, Brayford Pool, Lincoln LN6 7TS, United Kingdom. mjgreenall@physics.org

Mixtures of micelle-forming and lamella-forming amphiphiles in solution can form disk-shaped bilayers, sometimes referred to as bicelles. Using self-consistent field theory (SCFT), we investigate the structure and stability of these aggregates in a blend of two species of PS-PDMS diblock with PDMS homopolymer at . We find that the center of each disk is mainly composed of lamella-forming diblocks, while its thicker rim is mostly formed of micelle-forming diblocks. However, this segregation is not perfect, and the concentration of micelle formers is of the order of on the flat central surface of the bicelle. We also find that the addition of micelle former to the mixture of lamella former and homopolymer is necessary for disk-like bicelles to be stable. Specifically, the free energy density of the disk has a minimum as a function of the disk radius when both micelle- and lamella-forming diblocks are present, indicating that the bicelles have a preferred, finite radius. However, it decays monotonically when only lamella former is present, indicating that the bicelle structure is always unstable with respect to further aggregation in these systems. Finally, we identify a concentration range where the bicelle is predicted to have a lower free energy density than the simple cylindrical and lamellar aggregates, and so might be thermodynamically stable.

I Introduction

Mixtures of lamella-forming and micelle-forming amphiphiles can form structures in solution that are not seen when only one type of amphiphile is present Jain and Bates (2004); Safran et al. (1990); Zidovska et al. (2009). These structures, which include disks, ribbons, and perforated lamellae Dürr et al. (2013); Katsaras et al. (2005), have been studied in particular depth in lipid-detergent mixtures, where they are known collectively as bicelles Dürr et al. (2012). Bicelles are widely used in biophysical experiments as model membranes Dürr et al. (2013) for the solubilization of proteins Sanders and Landis (1995); Seddon et al. (2004), and have several advantages over alternative membrane structures such as multilayers and micelles. In particular, they avoid the curvature and strain effects that can occur in micelles Matsumori and Murata (2010); Prosser et al. (2006), and, unlike oriented multilayers, are usually optically tranparent Sanders and Prosser (1998). They contain smaller amounts of detergent, which can be harmful to membrane proteins, than mixed micelles Seddon et al. (2004), and, under certain circumstances Loudet et al. (2010); Diller et al. (2009); Macdonald et al. (2013); Dürr et al. (2013); Ram and Prestegard (1988), can be aligned in a magnetic field. Which of the various bicellar morphologies forms depends on a variety of factors Sanders et al. (2004), with the disk morphology appearing at lower temperatures Ramos et al. (2011) and lower lamella-former concentrations Dürr et al. (2013), and perforated lamellae forming at higher temperatures Ramos et al. (2011) and higher concentrations of lamella former Dürr et al. (2013).

Much progress has also been made in the modeling of bicelles, using techniques ranging from simple geometric models Vold and Prosser (1996) to microscopic simulations Vácha and Frenkel (2014). Recently, a theory based on the chemical potentials of the amphiphiles in differently shaped aggregates has been used to model light-scattering data on disk formation in mixed surfactant systems Anachkov et al. (2014). Disk-shaped bilayers and ribbon-like structures have been found in coarse-grained molecular dynamics simulations de Joannis et al. (2006); Noguchi (2013); Lee and Pastor (2011); Vácha and Frenkel (2014), and similar techniques have been used to investigate the transformation of small vesicles into disks Shinoda et al. (2011); Noguchi (2012) and vice versa Marrink and Mark (2003). However, the complexity of lipid-detergent systems, and the simplifying assumptions that must be made when modeling them, mean that detailed comparison of theoretical models with experimental results can prove difficult.

Fortunately, it can sometimes be possible to perform more detailed comparisons of theory with experiment for the self-assembly of micelles and bilayers in systems that are composed of a diblock copolymer and a homopolymer, rather than a lipid or detergent dissolved in water Kinning et al. (1991); Roe (1986); Leibler et al. (1983); Mayes and de la Cruz (1988). In these systems, scaling theories Leibler et al. (1983); Roe (1986) and self-consistent field theory Greenall et al. (2009a) give a good description of the size of the micelles, and self-consistent field theory can also predict the shape of the aggregates that will be formed for a given set of polymer parameters Greenall et al. (2009b).

In the current paper, we take a similar approach to the study of disk-shaped bicelles, and use self-consistent field theory Edwards (1965) to investigate the structures formed in a blend of two types of poly(styrene)-poly(dimethylsiloxane) (PS-PDMS) copolymer with poly(dimethylsiloxane) (PDMS) homopolymer. We have a number of reasons for focusing on this system. First, the two polymers segregate strongly Cochran et al. (2003); Eitouni and Balsara (2007), with the result that the system will mimic the lower temperatures where the disk-shaped bicelle morphology is most likely to be seen. Second, the large difference in electron absorption and scattering of PS and PDMS means that the microstructures formed in this system can be studied directly by electron microscopy, without the need for staining Saam and Fearon (1971). The third reason is the possibility of modifying the properties of the PDMS. The two polymers have very different glass transition temperatures Chow et al. (2006); Krause et al. (1982), and the PS blocks will be become glassy when the system is cooled to room temperature, while the PDMS will remain viscoelastic, yielding a dispersion of hard disks. Such systems Dykes et al. (2010) can display a range of behavior, including shear-induced phase separation Brown and Rennie (2001) and the formation of networks of platelets Nicolai and Cocard (2001); Patil et al. (2008).

The paper is organized as follows. In the following section, we describe our PDMS/PS system. Next, we give a brief introduction to the technique to be used, self-consistent field theory (SCFT). We then present and discuss our results, and give our conclusions in the final section.

Ii Details of the system

As discussed in the Introduction, blends of PS and PDMS are strongly segregated, and the Flory parameter between the two species is high. It is given as a function of temperature (measured in K) by Cochran et al. (2003); Eitouni and Balsara (2007) . This expression for is defined with respect to a reference volume of , and is valid for temperatures from 165 to 225C. To keep our SCFT algorithm numerically stable, we carry out our calculations at the upper end of this range (225C), where the polymers will be slightly more weakly segregated and the interfaces less sharp. The lamella-forming species is chosen to have a PDMS block of molar mass 6000g/mol and a PS block of 13000g/mol, and the molar mass of the PDMS homopolymer is set to 4500g/mol. The homopolymer has been chosen to be shorter than the PDMS block of the copolymer to ensure thermodynamic stability Mayes and Delacruz (1988); Leibler et al. (1983), and the ratio of the block lengths of the copolymer has been set so that this species will have a strong tendency to form lamellae Kinning et al. (1988); Greenall et al. (2009b). A polymer chain with ‘hydrophobic’ and ‘hydrophilic’ components in a ratio of roughly two to one, as is the case here, also provides a simple model of the long-chain DMPC molecules that form the body of bicelles in lipid-detergent mixtures Dürr et al. (2013); Vácha and Frenkel (2014). To ensure that it has a clear preference for forming micelles, the other species has a long PDMS block, with molar mass 24000g/mol, and a PS block of 6000g/mol. The overall weight fraction of copolymer, including both species, is set to . Using this relatively low concentration has two advantages. First, it will allow us to neglect interactions between the aggregates when modeling the thermodynamics of the system. Second, it will avoid interference between the aggregate and the boundary of the calculation box. We do not risk dropping below the critical micelle concentration by using a copolymer concentration of this magnitude, because the system under consideration is very strongly segregatedKinning et al. (1991). We also need to know the mean-square end-to-end distance, , of each polymer species as a function of its molar mass, . For polymers in a melt, this is given by an expression of the form , where the ratio is close to constant for a given polymer Brandrup and Immergut (1989). For PS Brandrup and Immergut (1989), we have that . The dimensions of PDMS Konishi et al. (1991) are rather similar, and .

Certain of the parameters above cannot be used directly as input to an SCFT calculation, and need to be converted into the appropriate forms. First, the molar masses listed above are converted to molar volumes Kinning et al. (1991) using the specific volumes (in ) of PS and PDMS. Empirical expressions for the specific volumes of PS Richardson and Savill (1977) and PDMS Shih and Flory (1972) as a function of temperature are taken from the literature. The molar volumes are then used to calculate the volume fraction of PDMS in each type of diblock and the volume ratios of the various species. The overall weight fraction of diblocks is converted to a volume fraction by a similar procedure, and the volume of a single molecule of each species in is calculated by dividing the appropriate molar volume by , a numerical constant that incorporates Avogadro’s number and the conversion from to . Finally, the number, , of repeat units in a polymer chain can be calculated by dividing the volume of the molecule by the reference volume, .

Iii Self-consistent field theory

Self-consistent field theory (SCFT) Edwards (1965) is a mean-field model that has been used with success to calculate the form and free energy of equilibrium Maniadis et al. (2007); Drolet and Fredrickson (1999); Matsen (2006) and metastable Duque (2003); Katsov et al. (2004) structures in systems composed of homopolymers Werner et al. (1999), copolymers Müller and Gompper (2002); Wang et al. (2010) and mixtures of these Denesyuk and Gompper (2006). SCFT has several features that make it especially suited to the current problem. First, as stated above, it has been shown to give a good description of the shape and size of micelles in blends of a single species of block copolymer with a homopolymer Greenall et al. (2009a, b), a system closely related to the current one. Second, it is faster than simulation techniques such as Monte Carlo, but can yield comparably accurate predictions of micelle size and shape Cavallo et al. (2006); Wijmans and Linse (1995); Leermakers and Scheutjens (1990). This will allow bicelles of a range of sizes to be studied in a reasonable period of time. Finally, it does not require any assumptions to be made regarding the segregation of copolymers of different architecture within an aggregate, meaning that any such effects that we observe are a natural prediction of the theory and have not been added by hand.

We now give a short overview of the application of SCFT to our system of two copolymers and a homopolymer, and refer the reader to reviews Matsen (2006); Fredrickson (2006); Schmid (1998) for in-depth presentations of the theory and to earlier papers Greenall and Gompper (2011); Greenall et al. (2009b, a) for a detailed description of our calculations. SCFT is a coarse-grained theory, and individual molecules are modeled as random walks in space Schmid (1998). An ensemble of many such molecules is considered, and the intermolecular forces are modeled by introducing contact potentials between the molecules and assuming that the blend is incompressible Matsen (2006). The Flory parameter discussed above is used to specify the strength of the repulsion between the two chemical species. In order to reduce the computational difficulty of the problem, fluctuations are neglected; that is, a mean-field approximation is made Matsen (2006). In the case of long molecules, this approximation is quantitatively accurate Cavallo et al. (2006); Fredrickson (2006); Matsen (2006).

SCFT can be used to perform calculations in different thermodynamical ensembles Duque (2003); Zhou and Shi (2011). In this paper, we perform all calculations in the canonical ensemble, which corresponds to keeping the amounts of copolymer and homopolymer in the simulation box fixed. This approach makes it easier for us to access more complex aggregates, such as bicelles. Such structures are more difficult to find in ensembles where the concentrations of the various species are able to change, and sometimes need to be stabilized by applying geometric constraints to the density profile Katsov et al. (2004).

Applying the mean-field approximationMatsen (2006) to our system, we find that the SCFT expression for the free energy of a system of copolymer species 1 and 2 in homopolymer is given by


where the are the mean volume fractions of the various components, with or PDMS2 for the poly(dimethyl siloxane) components of species and , or PS2 for the poly(styrene) components and for the homopolymer solvent, and the are the local volume fractions. is the total volume, is the number of repeat units in species 1, and is the SCFT free energy of a homogeneous system of the same composition. The architectures of the individual molecules enter through the single-chain partition functions , which are calculated from the propagators and Matsen (2006). These latter quantities satisfy diffusion equations with a field term that incorporates the polymer interactions. One field is associated with the PDMS segments, and one with the PS segments. This means that, to calculate the copolymer partition functions, and , the diffusion equation must be solved with the field appropriate to each of the two blocks of the copolymer Fredrickson (2006); Matsen (2006). In addition, the difference between the expressions for the mean-square end-to-end distances of PDMS and PS Konishi et al. (1991); Brandrup and Immergut (1989) must be taken into account when calculating the prefactor of the term in each block. The polymer density profiles are computed from integrals over the propagators Matsen (2006); Fredrickson (2006), with the volume fractions of PS and PDMS in each copolymer species entering via the limits of integration.

We perform all our calculations in cylindrical polar coordinates. Since we mainly focus on disk-shaped aggregates, we reduce the problem to a two-dimensional one by assuming that the system has rotational invariance about the -axis, and carry out all our calculations in a cylindrical box. We impose reflecting boundary conditions at all edges of the box, and the center of the bicelle lies at the origin of the coordinate system. The vertical height of this box is set to (meaning that its effective height is ), and its radius is varied according to a procedure that will be described later. The diffusion equations are solved using a finite difference method Press et al. (1992) with a step size of . The curve parameter that specifies the distance along the polymer backboneFredrickson (2006) runs from to , and its step size is set to for the copolymers and for the homopolymer.

The derivation of the mean-field free energy also yields a set of simultaneous equations relating the fields and to the densities . To calculate the density profiles for a given set of mean volume fractions , we make an initial guess for the fields that has the approximate form of the structure we wish to study, and then solve the diffusion equations to calculate the propagators and density profiles corresponding to these fields. The new are then substituted into the simultaneous equations to compute updated fields Matsen (2004), which are then used in turn to compute new by solving the diffusion equation as described above. For the algorithm to remain stable, the iteration needs to be damped, and, instead of using the updated values of directly to calculate the , we use the linear combination where . This process is repeated until convergence is achieved. For smaller systems, this simple mixing is sufficient. However, in larger calculation boxes, it can stall after an initial period of convergence. When performing calculations on these systems, we follow Thompson and co-workers Thompson et al. (2004) in passing the fields generated by the simple mixing iterations to an Anderson mixing algorithm Anderson (1965). This algorithm has greater flexibility in the iteration steps it can take, since it stores a history of previous values of the , and calculates the next estimate for the fields by adding these together in a linear combination Ng (1974). The pure Anderson mixing procedure is not stable in our case, and has to be damped Schmid (1997); Stasiak and Matsen (2011). We find that a history of 50–100 previous values of the and damping parameter of 0.1–0.2 yield good results.

We now need to address the issue of how to relate the thermodynamics of a single bicelle to those of a larger system containing many aggregates. To do this, we adapt a procedure that has been developed to study simple micelles and bilayers Greenall et al. (2009a, b); Greenall and Marques (2013); Liaw et al. (2015). First, we compute the free-energy density of a cylindrical box containing a single disk-shaped aggregate. Since the copolymer concentration is low, the aggregate is surrounded by a large volume of homopolymer, and the shape of the aggregate is not influenced by contact with the boundaries of the system. The radius of the simulation box is then varied, keeping the total volume fractions of both types of copolymer constant. Changing the size of the box in this way causes the bicelle to grow in the radial direction. The free-energy density of the system is calculated for each box radius. As is the case when this method is applied to micelles, there will be a minimum in the free energy as a function of the radiusLiaw et al. (2015) if the bicelle is stable. This solution of the SCFT equations corresponds to the optimum size of the bicelle. If, on the other hand, the bicelle is unstable with respect to further aggregation, the free energy density will decay monotonically as the box radius is increased, meaning that the bicelle can always move to a more energetically favorable state by growing in the radial direction, eventually forming an extended bilayer.

This approach is designed to mimic the behavior of a larger system (of fixed volume and fixed copolymer volume fraction) containing many bicelles. The reason for this is that this larger system minimizes its total free energy by changing the number of aggregates and hence the volume (‘box size’) occupied by each. Minimizing the free energy density in this way locates the bicelle that is the most energetically favorable and therefore the most likely to be observed in a sample containing many aggregates.

This minimum of the free-energy density with respect to the box size corresponds to the absolute free energy minimum of a solution of aggregates, and a point on the curve corresponds to a monodisperse solution of aggregates of a given size. In an earlier publication Greenall and Marques (2013), we showed how to use these curves to take into account simple fluctuations around the free energy minimum and so calculate the width, , of the size distribution of aggregates in a system containing only one type of amphiphile. This was achieved by relating the curvature of to the second derivative, , of the free energy of an aggregate containing molecules, which was then used Puvvada and Blankschtein (1990) to calculate via . We began by writing down an expression for the free-energy density of a monodisperse system of aggregates, each containing copolymers. This was given by


where is the number density of copolymers, is the number density of aggregates, is the free energy of a copolymer in solution, is the free energy of an aggregate of copolymers, is the volume of a single copolymer, and first term arises from the entropy of the free copolymers in solution. We then noted that a single SCFT calculation finds the local minimum of the free energy density for an aggregate in a box of volume . In the process, it determines the optimum number of molecules in the aggregate for this box size and so corresponds to minimizing with respect to at a given aggregate number density . Varying then yields the curve , from which we can read off . Remembering that this derivative is evaluated along the line where , we found that


where we have written , , and converted the number density to the volume fraction to express in terms of quantities that are either input to or results of our SCFT calculations.

We now extend this calculation to a system containing two amphiphile species. As in the calculation above, we suppose that there is only one type of aggregate in the system. We make the further simplifying assumption that the number fractions of the two species in the aggregate are the same as the overall number fractions. This turns out to hold accurately in our numerical results, and allows us to write the free-energy density of the two-species system as


where the number densities of species 1 and 2 are and respectively (so that ), and are the free energies of single copolymers of species 1 and 2 in solution, and is the volume of a single copolymer of species 2. This assumption also allows us to write the number of copolymers in the aggregate as , where


The calculation then proceeds as before Greenall and Marques (2013), and we find that


where, as before, we have converted number densities to volume fractions. The extra factor of in equation 6 arises from the normalization of the free energy in equation 1.

To see if the bicelle is likely to form, we also need to compare its free energy density with those of other aggregates. If we can establish that, for a given concentration range, the bicelle has a lower free-energy density than the simple spherical, cylindrical and lamellar aggregates, there is a clear possibility that it will form in this region. To calculate the free energy density of spherical micelles, we simply continue our calculations to smaller values of the calculation box radius until the disk-like bicelle shrinks to a sphere. It turns out that, for the systems considered here, where the parameter is large and a high proportion of thepolymers have a strong preference for forming lamellae, spherical micelles have much higher free-energy densities than the other aggregates. In consequence, we omit them from our results. The free energies of the optimum cylinder and lamella are found by a similar method to that described above Greenall et al. (2009b). These structures are assumed to be of infinite extent, so calculating their optimum free energies using SCFT is a one-dimensional problem, with the box size being varied in the - and -directions respectivelyGreenall et al. (2009b). For convenience, and for consistency with our earlier calculations, we perform these calculations using the same two-dimensional algorithm as before, but with the calculation box made very thin in the direction in which the density profiles do not vary.

Iv Results and discussion

We begin this section by demonstrating that the disk-like bicelle structure is a solution to the SCFT equations. Next, we investigate the dependence of the free energy density of the disk-like bicelle on its radius, to determine in which systems the bicelle has a preferred size. Finally, we compare the free energy of the bicelle to those of the simple cylindrical and lamellar structures for a range of concentrations of the two copolymer species, to find the range of parameters for which it might form in an experiment.

iv.1 The disk-shaped bicelle morphology

In Figure 1, we show a ray-traced plot of the surface of a disk-like bicelle obtained as a solution to the SCFT equations in the PS/PDMS system described above. by weight of diblocks are sphere formers, and are lamella formers. The surface is defined as the locus of the points where the local solvent volume fraction , and has a biconcave disk shape reminiscent of a red blood cell.

Figure 1: Ray-traced plot of the surface of the bicelle.

We now investigate the distribution of the two species, and their hydrophilic and hydrophobic components, within the bicelle. In Figure 2a, we show radial cuts through the volume fraction profiles of the various blocks that make up the bicelle. The division of the bicelle into a well defined “hydrophobic” PS core and “hydrophilic” PDMS corona is clearly visible, with the interface between the two regions being located at around . A marked difference in behavior between the sphere-forming and lamella-forming species is also seen. Specifically, the sphere formers are concentrated at the edge of the bicelle, with the volume fraction profile of their hydrophobic components being sharply peaked just before the core-corona interface, and their relatively long hydrophilic components stretching out into the solvent. This structure, in which the species with a preference for forming curved membranes segregates to the rim of the bicelle, is often sketched in the literature Sanders and Prosser (1998), and our results show that it can be reproduced in explicit calculations.

The curves in Figure 2b are cuts through the same volume fraction profiles in the -direction. Again, clear core and corona regions can be seen. The difference between the two plots lies in the fact that the peak in the hydrophobic profile of the sphere-forming species just before the core-corona interface is markedly less pronounced in the -direction, confirming the clear preference of the sphere formers for the bicelle rim. However, the presence of this peak, which attains a maximum height of 0.12, shows that the segregation of the sphere formers to the rim is not perfect, with an appreciable concentration of this species remaining near , particularly on the flat surface of the bicelle. Again, this result agrees with the pictorial model of bicelles often suggested in the biophysics literature Sanders and Prosser (1998).

Figure 2: Cuts through the volume fraction profiles of the various blocks making up the bicelle: sphere former PS blocks (thick solid line); sphere former PDMS blocks (thick dashed line); lamella former PS blocks (thin solid line); lamella former PDMS blocks (thin dashed line). The volume fraction profile of the PDMS homopolymer “solvent” is omitted for clarity. Panel (a) shows cuts in the -direction, and panel (b) shows cuts in the -direction. Note that neither panel shows the full range of the calculation box.

iv.2 The preferred radius of the bicelle

Having shown that the bicelle exists as a solution to SCFT and investigated the distribution of the two copolymer species within it, we now study the dependence of its free energy density on its radius, with the aim of finding whether it has a preferred size. The size of the bicelle is varied by changing the calculation box radius as described above, and we plot its free energy density against the radius of its core. We focus on the core radius because this is the measure of micelle size that is most easily measured in experiment Kinning et al. (1988). Since the interface between the core and the corona is sharp, all reasonable definitions of the core radius will yield similar values, and we define it here as the value of at which the core and corona densities are equal, so that .

Figure 3a shows the results of this calculation for the system studied above, where by weight of all amphiphiles are sphere formers and are lamella formers. The curve of the free energy density shows a clear minimum as a function of the bicelle radius, showing that the aggregate has a preferred size.

Figure 3: Plots of the free-energy density against bicelle core radius for (a) a mixed system in which by weight of all amphiphiles are sphere formers and are lamella formers and (b) a pure system in which all amphiphiles are lamella formers. The narrow range of the axes arises because the system is dilute and the free energy density is calculated for the entire calculation box, with the result that the differences in free energy density between one size of bicelle and another appear small. However, a more detailed calculation, in which the curvature of the free energy per chain in the bicelle is extracted from the above curves, reveals that the size selection in this system is significant, with a relative polydispersity, , of approximately in the aggregation number, corresponding to a relative polydispersity of around in the radius.

We can now use Equation 6 to estimate the relative polydispersity in the aggregation number of the micelle, . This is found to be approximately , corresponding Greenall and Marques (2013) to a relative polydispersity of around in the bicelle radius. This demonstrates that significant size selection takes place in this system, with a clear preferred size for the bicelles.

iv.3 Concentration dependence

We now compare the free energy density of the bicelle with those of the cylindrical and lamellar morphologies over a range of concentrations. Since the bicelle is a hybrid structure that contains elements of both the lamella and the cylinder, we expect Vácha and Frenkel (2014) that it will form at compositions around that at which the free energy densities of these two structures are the same. We therefore begin by locating this composition, and focus our attention on its vicinity. Figure 4 shows the results of these calculations. At higher weight fractions of sphere-forming amphiphile, the cylindrical micelle has the lowest free-energy density and the lamella has the highest. As the amount of sphere former is reduced, the free-energy density of the bicelle becomes the lowest of the three aggregates: it drops below that of the cylinder, while remaining below that of the lamella.

Figure 4: Free energy densities of the lamella (squares), bicelle (triangles) and cylinder (circles), normalised with respect to that of the bicelle and plotted against the percentage of amphiphiles that are sphere formers.

As the weight fraction of sphere former is decreased further, the free-energy densities of the bicelle and lamella become closer, although they do not cross within the range of compositions that we are able to study. It is possible that the two quantities approach each other asymptotically as the amount of sphere former tends to zero. This is consistent with the growth in the bicelle radius seen as the amount of sphere former is reduced (Figure 5). As the bicelle becomes very large, it approaches the lamella in shape, and the free-energy densities of the two structures can also be expected to become very close.

Figure 5: The bicelle core radius plotted against the percentage of amphiphiles that are sphere formers.

V Conclusions

Using self-consistent field theory, we have studied the structure and stability of disk-like bicelles in a blend of lamella-forming and sphere-forming PS-PDMS diblocks with PDMS homopolymer. We have found that these structures have a characteristic biconcave disk shape, like red blood cells. Furthermore, we have shown that the two amphiphile species are unequally distributed within the aggregate, with the center being mainly composed of lamella formers, and the rim being mostly formed of micelle formers. This picture of bicelles is often sketched in the biophysics literature Sanders and Prosser (1998), and is reproduced here in calculations on a well defined model.

We also find that the presence of micelle former is necessary for disk-like bicelles to be stable, and that they will be unstable with respect to further aggregation when only lamella former is present. Finally, we locate a region of parameter space range where the bicelle is predicted to have a lower free energy density than the competing cylindrical and lamellar aggregates.

There are are a number of ways in which our calculations could be extended. First, we used the rather high temperature of in order to keep the parameter at a level that our existing numerical methods were able to deal with. Further refinements of the Anderson mixing method and the algorithm used to solve the diffusion equation might allow the calculations to be extended to lower temperatures, where the concentration range at which the bicelle has a lower free-energy density than the other structures should be broader. Second, perhaps using these extended methods, we could look for other polymer blends in which the bicelle is predicted to be stable and in which its presence could have an effect on the mechanical properties. Finally, a similar SCFT approach could be applied to lipid systems Kik et al. (2010), where bicelles were originally observed. This could be used to investigate the degree of detergent penetration into the center of the bicelle, an important issue for the biophysical experiments that use bicelles as model membranes.


  • Jain and Bates (2004) S. Jain and F. S. Bates, Macromolecules 37, 1511 (2004).
  • Safran et al. (1990) S. A. Safran, P. Pincus, and D. Andelman, Science 248, 354 (1990).
  • Zidovska et al. (2009) A. Zidovska, K. K. Ewert, J. Quispe, B. Carragher, C. S. Potter, and C. R. Safinya, Langmuir 25, 2979 (2009).
  • Dürr et al. (2013) U. H. Dürr, R. Soong, and A. Ramamoorthy, Prog. Nucl. Mag. Res. Sp. 69, 1 (2013).
  • Katsaras et al. (2005) J. Katsaras, T. A. Harroun, J. Pencer, and M.-P. Nieh, Naturwissenschaften 92, 355 (2005).
  • Dürr et al. (2012) U. H. Dürr, M. Gildenberg, and A. Ramamoorthy, Chem. Rev. 112, 6054 (2012).
  • Sanders and Landis (1995) C. R. Sanders and G. C. Landis, Biochemistry 34, 4030 (1995).
  • Seddon et al. (2004) A. M. Seddon, P. Curnow, and P. J. Booth, Biochim. Biophys. Acta 1666, 105 (2004).
  • Matsumori and Murata (2010) N. Matsumori and M. Murata, Nat. Prod. Rep. 27, 1480 (2010).
  • Prosser et al. (2006) R. S. Prosser, F. Evanics, J. L. Kitevski, and M. S. Al-Abdul-Wahid, Biochemistry 45, 8453 (2006).
  • Sanders and Prosser (1998) C. R. Sanders and R. S. Prosser, Structure 6, 1227 (1998).
  • Loudet et al. (2010) C. Loudet, A. Diller, A. Grélard, R. Oda, and E. J. Dufourc, Prog. Lipid Res. 49, 289 (2010).
  • Diller et al. (2009) A. Diller, C. Loudet, F. Aussenac, G. Raffard, S. Fournier, M. Laguerre, A. Grélard, S. J. Opella, F. M. Marassi, and E. J. Dufourc, Biochimie 91, 744 (2009).
  • Macdonald et al. (2013) P. M. Macdonald, Q. Saleem, A. Lai, and H. H. Morales, Chem. Phys. Lipids 166, 31 (2013).
  • Ram and Prestegard (1988) P. Ram and J. H. Prestegard, Biochim. Biophys. Acta 940, 289 (1988).
  • Sanders et al. (2004) C. Sanders, A. K. Hoffmann, D. N. Gray, M. H. Keyes, and C. D. Ellis, ChemBioChem 5, 423 (2004).
  • Ramos et al. (2011) J. Ramos, A. Imaz, J. Callejas-Fernández, L. Barbosa-Barros, J. Estelrich, M. Quesada-Pérez, and J. Forcada, Soft Matter 7, 5067 (2011).
  • Vold and Prosser (1996) R. R. Vold and R. S. Prosser, J. Magn. Reson. Ser. B 113, 267 (1996).
  • Vácha and Frenkel (2014) R. Vácha and D. Frenkel, Langmuir 30, 4229 (2014).
  • Anachkov et al. (2014) S. E. Anachkov, P. A. Kralchevsky, K. D. Danov, G. S. Georgieva, and K. P. Ananthapadmanabhan, J. Colloid Interf. Sci. 416, 258 (2014).
  • de Joannis et al. (2006) J. de Joannis, F. Y. Jiang, and J. T. Kindt, Langmuir 22, 998 (2006).
  • Noguchi (2013) H. Noguchi, J. Chem. Phys. 138, 024907 (2013).
  • Lee and Pastor (2011) H. Lee and R. W. Pastor, J. Phys. Chem. B 115, 7830 (2011).
  • Shinoda et al. (2011) W. Shinoda, T. Nakamura, and S. O. Nielsen, Soft Matter 7, 9012 (2011).
  • Noguchi (2012) H. Noguchi, Soft Matter 8, 8926 (2012).
  • Marrink and Mark (2003) S. J. Marrink and A. E. Mark, J. Am. Chem. Soc. 125, 15233 (2003).
  • Kinning et al. (1991) D. J. Kinning, E. L. Thomas, and L. J. Fetters, Macromolecules 24, 3893 (1991).
  • Roe (1986) R.-J. Roe, Macromolecules 19, 728 (1986).
  • Leibler et al. (1983) L. Leibler, H. Orland, and J. C. Wheeler, J. Chem. Phys. 79, 3550 (1983).
  • Mayes and de la Cruz (1988) A. M. Mayes and M. O. de la Cruz, Macromolecules 21, 2543 (1988).
  • Greenall et al. (2009a) M. J. Greenall, D. M. A. Buzza, and T. C. B. McLeish, Macromolecules 42, 5873 (2009a).
  • Greenall et al. (2009b) M. J. Greenall, D. M. A. Buzza, and T. C. B. McLeish, J. Chem. Phys. 131, 034904 (2009b).
  • Edwards (1965) S. F. Edwards, Proc. Phys. Soc. 85, 613 (1965).
  • Cochran et al. (2003) E. W. Cochran, D. C. Morse, and F. S. Bates, Macromolecules 36, 782 (2003).
  • Eitouni and Balsara (2007) H. B. Eitouni and N. P. Balsara, Physical Properties of Polymers Handbook (Springer-Verlag, New York, 2007), chap. 19, 2nd ed.
  • Saam and Fearon (1971) J. C. Saam and F. W. G. Fearon, Ind. Eng. Chem. Prod. Res. Develop. 10, 10 (1971).
  • Chow et al. (2006) W. W. Y. Chow, K. F. Lei, G. Shi, W. J. Li, and Q. Huang, Smart Mater. Struct. 15, S112 (2006).
  • Krause et al. (1982) S. Krause, Z.-H. Lu, and M. Iskandar, Macromolecules 15, 1076 (1982).
  • Dykes et al. (2010) L. M. C. Dykes, J. M. Torkelson, W. R. Burghardt, and R. Krishnamoorti, Polymer 51, 4916 (2010).
  • Brown and Rennie (2001) A. B. D. Brown and A. R. Rennie, Chem. Eng. Sci. 56, 2999 (2001).
  • Nicolai and Cocard (2001) T. Nicolai and S. Cocard, Eur. Phys. J. E 5, 221 (2001).
  • Patil et al. (2008) S. P. Patil, R. Mathew, T. G. Ajithkumar, P. R. Rajamohanan, T. S. Mahesh, and G. Kumaraswamy, J. Phys. Chem. B 112, 4536 (2008).
  • Mayes and Delacruz (1988) A. M. Mayes and M. O. Delacruz, Macromolecules 21, 2543 (1988).
  • Kinning et al. (1988) D. J. Kinning, K. I. Winey, and E. L. Thomas, Macromolecules 21, 3502 (1988).
  • Brandrup and Immergut (1989) J. Brandrup and E. H. Immergut, eds., Polymer Handbook (Wiley, New York, 1989).
  • Konishi et al. (1991) T. Konishi, T. Yoshizaki, and H. Yamakawa, Macromolecules 24, 5614 (1991).
  • Richardson and Savill (1977) M. J. Richardson and N. G. Savill, Polymer 18, 3 (1977).
  • Shih and Flory (1972) H. Shih and P. J. Flory, Macromolecules 5, 758 (1972).
  • Maniadis et al. (2007) P. Maniadis, T. Lookman, E. M. Kober, and K. O. Rasmussen, Phys. Rev. Lett. 99, 048302 (2007).
  • Drolet and Fredrickson (1999) F. Drolet and G. H. Fredrickson, Phys. Rev. Lett. 83, 4317 (1999).
  • Matsen (2006) M. W. Matsen, in Soft Matter, edited by G. Gompper and M. Schick (Wiley-VCH, Weinheim, 2006), chap. 2.
  • Duque (2003) D. Duque, J. Chem. Phys. 119, 5701 (2003).
  • Katsov et al. (2004) K. Katsov, M. Müller, and M. Schick, Biophys. J. 87, 3277 (2004).
  • Werner et al. (1999) A. Werner, M. Müller, F. Schmid, and K. Binder, J. Chem. Phys. 110, 1221 (1999).
  • Müller and Gompper (2002) M. Müller and G. Gompper, Phys. Rev. E 66, 041805 (2002).
  • Wang et al. (2010) J. F. Wang, K. K. Guo, L. J. An, M. Müller, and Z. G. Wang, Macromolecules 43, 2037 (2010).
  • Denesyuk and Gompper (2006) N. A. Denesyuk and G. Gompper, Macromolecules 39, 5497 (2006).
  • Cavallo et al. (2006) A. Cavallo, M. Müller, and K. Binder, Macromolecules 39, 9539 (2006).
  • Wijmans and Linse (1995) C. M. Wijmans and P. Linse, Langmuir 11, 3748 (1995).
  • Leermakers and Scheutjens (1990) F. A. M. Leermakers and J. M. H. M. Scheutjens, J. Colloid Interface Sci. 136, 231 (1990).
  • Fredrickson (2006) G. H. Fredrickson, The Equilibrium Theory of Inhomogeneous Polymers (Oxford University Press, Oxford, 2006).
  • Schmid (1998) F. Schmid, J. Phys.: Condens. Matter 10, 8105 (1998).
  • Greenall and Gompper (2011) M. J. Greenall and G. Gompper, Langmuir 27, 3416 (2011).
  • Zhou and Shi (2011) J. Zhou and A.-C. Shi, Macromol. Theory Simul. 20, 690 (2011).
  • Press et al. (1992) W. H. Press, B. P. Flannery, S. A. Teukolsky, and W. T. Vetterling, Numerical Recipes in C (Cambridge University Press, Cambridge, 1992), 2nd ed.
  • Matsen (2004) M. W. Matsen, J. Chem. Phys. 121, 1938 (2004).
  • Thompson et al. (2004) R. B. Thompson, K. Ø. Rasmussen, and T. Lookman, J. Chem. Phys. 120, 31 (2004).
  • Anderson (1965) D. G. Anderson, J. Assoc. Comput. Mach. 12, 547 (1965).
  • Ng (1974) K. C. Ng, J. Chem. Phys. 61, 2680 (1974).
  • Schmid (1997) F. Schmid, Phys. Rev. E 55, 5774 (1997).
  • Stasiak and Matsen (2011) P. Stasiak and M. W. Matsen, Eur. Phys. J. E 34, 110 (2011).
  • Greenall and Marques (2013) M. J. Greenall and C. M. Marques, Phys. Rev. Lett 110, 088301 (2013).
  • Liaw et al. (2015) C. Y. Liaw, K. J. Henderson, W. R. Burghardt, J. Wang, and K. R. Shull, Macromolecules 48, 173 (2015).
  • Puvvada and Blankschtein (1990) S. Puvvada and D. Blankschtein, J. Chem. Phys. 92, 3710 (1990).
  • Kik et al. (2010) R. A. Kik, F. A. M. Leermakers, and J. M. Kleijn, Phys. Rev. E 81, 021915 (2010).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description