Thermodynamics and Bending Energetics of Microemulsions

A comprehensive, yet simple, theoretical model for droplet microemulsions is presented. The model combines thermodynamics of self-assembly with bending elasticity theory and relates microemulsion properties, such as average droplet size, polydispersity, interfacial tension and solubilisation capacity with the three bending elasticity constants, spontaneous curvature ( H 0 ), bending rigidity ( k c ) and saddle-splay constant ( k c ). In addition, the self- association entropy constant ( k s ) explicitly determines various microemulsion properties. The average droplet size is shown to increase with increasing effective bending constant, defined as k eff ¼ 2 k c þ k c þ k s , as well as with decreasing magnitudes of H 0 . The polydis- persity decreases with increasing values of k eff , but does not at all depend on H 0 . The model predicts ultra-low interfacial tensions, the values of which decrease considerably with increasing droplet radius, in agreement with experiments. The solubilisation capacity increases as the number of droplets is increased with increasing surfactant concentration. In addition, an enhanced solubilisation effect is obtained as the size of the droplets increases with increasing surfactant concentration, as a result of self-association entropy effects. It is demonstrated that self-association entropy effects favour smaller droplet size as well as larger droplet polydispersity.


Introduction
Surfactants are amphiphilic molecules consisting of a hydrophilic head group and a hydrophobic tail. Surfactants self-assemble above a certain surfactant concentration in an aqueous solvent to form micelles or bilayers. The driving force for the self-assembly process is the tendency to minimize the interfacial contact area between water and surfactant hydrophobic tails. As a result, the tails make up an interior core of the micelles that is absent of water whereas the hydrophilic head groups are located at the interface adjacent to the aqueous solvent. Hydrophobic components that are usually insoluble in water, such as oil or fat, may be dissolved into a surfactant-water system. The oil molecules become incorporated into the hydrophobic core of the micelles to form swollen micelles. However, there is a limit how much oil that may be dissolved by micelles and above a certain amount of oil two separate liquid phases coexist, i.e. thermodynamically stable oil droplets (discrete phase) dissolved in water (continuous phase) coexisting with excess oil (cf. Figure 1). Such an oil-in-water (o/w) microemulsion phase coexisting with an excess oil phase in a two-phase system is denoted as Winsor I microemulsion. Some surfactants, on the other hand, may form reversed water-in-oil (w/o) droplets that dissolve water (discrete phase) in oil (continuous phase) in the presence of excess water (Winsor II microemulsion). Moreover, a microemulsion phase may also coexist with both excess oil and excess water in a three-phase system (Winsor III microemulsion) [1]. A Winsor III microemulsion does not usually consist of finite-sized droplets. Rather some kind of macroscopic structure is formed, either more or less planar alternating layers of oil and water separated by surfactant monolayers, or some kind of ordered or disordered bicontinuous phase with a system of separated tunnels of water and oil [2].
The microscopic structure of microemulsions (whether Winsor I, II or III), is mainly determined by the chemical structure of the surfactant. Surfactants with large hydrophilic head groups and small hydrophobic (or lipophilic if one prefers) parts tend to curve so as to form ordinary micelles or o/w microemulsion droplets whereas surfactants with a small hydrophilic part and large hydrophobic part (like phospholipids) tend to form w/o microemulsion droplets (cf. Figure 1). The molecular properties determining the curvature of a surfactant monolayer may be summed up in the quantity hydrophilic-lipophilic balance, or shorter the HLB value [3]. As a result, one expects o/w microemulsion droplets to increase in size as the monolayer becomes less curved with decreasing HLB values. For sufficiently low HLB values, a transition from Winsor I to Winsor III is expected, and further decreasing HLB would result in the formation of Winsor II microemulsion. In addition to the volume of the hydrophilic and hydrophobic parts, respectively, the presence of electric charge on the surfactant head group also contributes to increase the HLB value and promotes a positive (oil-in-water) curvature of the droplets.
It is possible to tune the HLB value of ionic surfactants by means of adding salt to an aqueous phase and HLB decreases with increasing electrolyte concentration. Similarly, the HLB value of non-ionic ethylene oxide-based surfactants is found to be considerably temperature-sensitive and, as a result, it is possible to observe the sequence of transitions oil-in-water ! bicontinuous ! water-in-oil microemulsion by means of increasing the temperature [4]. The amount of oil dissolved in the microemulsion increases during this sequence of transitions and, at a certain temperature denoted the phase inversion temperature (PIT), the microemulsion phase contains equal amounts of oil and water. It is also possible to tune the curvature of a surfactant monolayer by means of adding a cosurfactant that is mixed into the layer [5].
In contrast to conventional (macro)emulsions, microemulsions are thermodynamically stable systems. This means that the size of the droplets may fall within a wide range, depending on the surfactant HLB value, from small micelle-like droplets of about 1 nm to about 100 nm, above the size of which the droplets usually transform into a macroscopic bicontinuous structure. In contrast, kinetically stabilized emulsion droplets are usually larger falling in the range 10 nm to 100 μm.
In this chapter, we present a comprehensive, yet simple, theory that rationalizes the structural behaviour of spherical microemulsion droplets. The theory is based on conventional solution thermodynamics combined with bending elasticity theory and it predicts several experimentally available quantities, such as droplet size and polydispersity, interfacial tension and solubilisation capacity. We only consider the case of rather rigid microemulsion droplets that are spherically shaped. More flexible droplets, consisting of interfacial monolayers with bending rigidities approaching zero (see further below), may assume a more spheroidal or ellipsoidal shape [6] or undergo undulatory fluctuations of the droplet interfaces [7].

Thermodynamics of self-assembly
Microemulsions are thermodynamically stable equilibrium structures. Hence, the theoretical treatment of microemulsion droplets necessarily needs to take into account thermodynamics of self-assembling surfactant molecules to form an interfacial monolayer that encapsulates the discrete phase (cf. Figure 1). The latter process may be considered in terms of a set of multiple equilibrium reactions where N is the free surfactant molecules dissolved in the continuous phase (Bancroft's rule states that free surfactants are always present in the continuous phase) that form a microemulsion droplet M N . The change in entropy for the process of self-assembling free surfactant molecules to the interface of a single droplet that encloses a certain amount of oil (o/w microemulsions) or water (w/o microemulsions) in accordance with Eq. (1) may thus be written as where k is the Boltzmann's constant. For the case of o/w microemulsions, the free surfactant molecules are present in the aqueous continuous phase with a volume fraction φ free , whereas oil molecules are brought from the adjacent macroscopic excess oil phase to form droplets encapsulated by surfactant monolayers with volume fraction φ N . The expression in Eq. (2) is also valid for the reverse case of w/o microemulsions in which the surfactant is dissolved in the continuous oil phase and water is brought from an excess aqueous phase. ΔS agg is always a negative quantity and, as a result, the process of self-assembly must be unfavourable in the absence of other contributions favouring the self-assembly process. Equation (2) can be understood as the difference between the entropy of mixing self-assembled surfactant monolayers and solvent, and the entropy of mixing the corresponding amount of free surfactant molecules with solvent.
In order to allow for the spontaneous self-assembly of surfactant molecules, additional contributions to the overall free energy change must be added. These contributions are related to the formation of the surfactant monolayer and can be collected into a single quantity denoted by Δμ. The most important contribution to Δμ comes from the hydrophobic effect, i.e. the largely entropic favourable process of reducing interfacial contact between surfactant hydrophobic tails and water molecules as free surfactants are self-assembled in the microemulsion-droplet interfaces. From a molecular point of view, the hydrophobic effect is believed to be the result of organization of water molecules adjacent to the hydrophobic moieties of the surfactant molecules. We may now write the total change in Gibbs free energy of the self-assembly process as Introducing the following free energy parameter [8].
Ε N NΔμ À NkTlnφ f ree (4) and combining Eqs. (2)-(4) gives the following set of equilibrium conditions one for each surfactant aggregation number N of the interfacial monolayers. By setting ΔG N = 0, we account for a situation where surfactant molecules are reversibly exchanged between droplet interface and as free molecules in the continuous phase. Equation (5) may now be rewritten as Summing up the different volume fractions in Eq. (6) gives the total volume fraction φ s of surfactant self-assembled in the droplets, i.e.
which may be evaluated from the following integral approximation We have been able to approximately set the lower limit in the integral to zero since the largely curved droplets with low aggregation numbers are too energetically unfavourable to contribute to the integral (see further below). Equation (8) gives the full size distribution of surfactants self-assembled in microemulsion-droplet interfaces; in so far an mathematical expression for the function E(N) is available. In the following section, we will demonstrate that such an expression may be obtained by means of considering the properties of bending a surfactant monolayer.

Bending elasticity
The free energy of an arbitrarily shaped surfactant monolayer may be calculated taking into account bending elasticity properties. The curvature at a single point on a surfactant monolayer (most conveniently defined at the hydrocarbon/water interface [cf. Figure 1]) may be defined by considering two perpendicular curves on the interface with radii of curvature, R 1 and R 2 , respectively. Hence, each point at the aggregate interface may be distinguished by two principal curvatures defined as c 1 = 1/R 1 and c 2 = 1/R 2 . By means of changing variables from c 1 and c 2 to the mean curvature H = ½(c 1 + c 2 ) and Gaussian curvature K = c 1 Ác 2 , respectively, the following fruitful expression for the free energy per unit area may be introduced [9] γðH, The Helfrich expression in Eq. (9) is a second-order expansion with respect to H and K, respectively, and introduces three important parameters, spontaneous curvature (H 0 ), bending rigidity (k c ) and saddle-splay constant (k c Þ. γ 0 is the free energy per unit area at H = H 0 and K = 0. The total free energy of a surfactant monolayer may be obtained by means of integrating Eq. (9) over the entire interfacial area A of the monolayer giving For a spherically curved monolayer, an expression for E is obtained by simply introducing the proper values H = 1/R and K = 1/R 2 in Eq. (9) and multiplying γ with the interfacial area 4πR 2 .
As a result, the following expression for E as a function of droplet radius R is obtained for a surfactant monolayer at the interface of a spherical microemulsion droplet where γ p is the interfacial tension of a strictly planar interface, i.e. γ p = γ(H = K = 0). The free energy in Eq. (11) may be considered as a sum of the contribution of forming a planar surfactant monolayer out of free surfactant (=4πR 2 γ p ) and the contribution of bending the monolayer into spherical shape.
The three bending elasticity constants H 0 , k c and k c influence various properties of self-assembled surfactant monolayers and are expected to assume values that reflect the chemical properties of the constituent surfactant or surfactant mixture as well as environmental properties such as temperature and electrolyte concentration. They may be interpreted as thermodynamic parameters that may be determined from experiments [10][11][12][13] as well as calculated from a suitable molecular model [14][15][16] by means of minimizing the free energy per molecule of a surfactant interfacial layer at given values of H and K.

Spontaneous curvature
The spontaneous curvature H 0 represents the sign and magnitude of the preferential curvature of a single surfactant monolayer. Since smaller microemulsion droplets must be more curved than larger ones, H 0 is expected to have a large impact on the droplet size. We have defined the mean and spontaneous curvature of o/w microemulsion droplets as positive, whereas w/o droplets have a negative spontaneous curvature (cf. Figure 1). H 0 depends on the chemical structure of a surfactant molecule in so far it is expected to increase with an increasing HLB value [14,16]. This means that H 0 rapidly decreases as the tail length is increased for a surfactant with given head group. Likewise, ionic surfactants usually have larger spontaneous curvature than non-ionic surfactants with a head group of similar size, and adding electrolyte is expected to reduce H 0 of an ionic surfactant monolayer. As a result, in the presence of both water and oil, surfactants with high HLB values and positive H 0 tend to form o/w microemulsion droplets whereas surfactants with low HLB values and negative H 0 tend to form w/o microemulsions.

Bending rigidity
The bending rigidity k c is a measure of the ability of a surfactant monolayer to resist deviations from a uniform mean curvature H = H 0 . k c must always be a positive quantity in order to realize a minimum of γ as a function of mean curvature in the Helfrich expression in Eq. (9). Large bending rigidities are expected to favour surfactant monolayers with small deviations from a homogenous curvature, i.e. rigid and monodisperse droplets with a uniform shape. This means that microemulsion droplets are expected to be strictly spherical for sufficiently large values of k c . Likewise, a bicontinuous microemulsion is expected to form an ordered cubic phase as the bending rigidity is high whereas a more disordered sponge phase forms at low k c values [17].
From a molecular point of view, it has been demonstrated that k c assumes a maximum value for some optimal value of the HLB of a surfactant [16]. The flexibility of the surfactant tails appears to have a crucial impact on k c which is found to be significantly larger for surfactants with a flexible rather than rigid tail [14]. As a matter of fact, it has been demonstrated that k c must always be equal to zero, thus preventing self-assembly of non-ionic surfactants made up of hydrophobic and hydrophilic parts that both are rigid [16]. For surfactants with a flexible tail, the bending rigidity is expected to increase with increasing tail length [10,14]. Moreover, k c has been found to become substantially reduced in surfactant mixtures [18][19][20] and the magnitude of the reduction has been demonstrated to increase with increasing asymmetry between two surfactants with respect to head group charge number, volume of surfactant tail, etc [15].

Saddle-splay constant
The saddle-splay constant k c is related to the Gaussian curvature K. As implied by its name, high positive values of k c influence the curvature of an interface so as to favour a saddle-like structure, with negative Gaussian curvature, that is curved in opposite directions perpendicular to one another. Examples of surfactant assemblies with saddle-like structure are torus-like micelles and bicontinuous microemulsions. Hence, the latter is expected to be favoured by large (positive) values of k c whereas a planar lamellar Winsor III microemulsion is favoured by low and negative values of k c [17].
Moreover, according to the Gauss-Bonnet theorem, the last integral in Eq. (10) equals where g represents the number of handles or holes present in a surfactant monolayer. A microemulsion droplet consists of one geometrically closed monolayer and, as a result, g = 0.
As a consequence, the last term in Eq. (10) equals 4πk c , the quantity of which does not depend on the size of the microemulsion droplet. Since k c contributes with a size-independent term to E, the value of k c indirectly determines the size of the droplets. Positive values of k c indicate that the total free energy in the system increases as one single droplet is split up to form two or many smaller droplets. Likewise, negative values of k c imply a decreasing total free energy as droplets are split up to increase their number while maintaining the surfactant concentration fixed. In either case, increasing values of k c favour larger microemulsion droplets.
From a molecular point of view, k c is expected to be negative for a pure symmetric surfactant that tends to curve equivalently in all lateral directions [14]. However, introducing some kind of asymmetry of the surfactant may result in positive values of k c . For instance, the nonuniform distribution of charges in monolayers formed by ionic gemini surfactants may result in a positive saddle-splay constant [13,21]. Likewise, mixtures of two or more asymmetric surfactants, in particular mixtures of oppositely charged surfactants, may result in positive values of k c [22].

Microemulsion droplets
The size distribution of surfactant monolayers in Eq. (8) may be rewritten so as to give where φ s (R) is the volume fraction distribution of surfactant monolayers and may be evaluated from the geometrical relation Nv ¼ 4πR 2 ξ for a surfactant monolayer with thickness ξ and surfactant molecular volume v. However, the self-assembled monolayer is excluded from the interior volume occupied by the discrete phase of oil molecules (in the case of o/w droplets). As a result, it is more appropriate to consider the size distribution based on the volume fraction of entire droplets φ d , including both the surfactant monolayers at the droplet interfaces and solubilised oil molecules in the interior core of the droplets.
The geometrical relations φ d ðRÞ ∝ 4πR 2 dR and φ s ðRÞ ∝ 8πRξdR imply that the volume fraction distribution of droplets equals Now we may substitute Eq. (11) into Eq. (13) and combine with Eqs. (14) and (15) to arrive at the full size distribution of spherical microemulsion droplets from which the total volume fraction of droplets (including both surfactant and the discrete phase) may be evaluated so as to give The size distribution according to Eq. (16) is plotted in Figure 2 for some different values of volume fraction of droplets. In accordance with Eq. (16), the formation of finite-sized microemulsion droplets may be considered as the result of bending elasticity properties as well as self-association entropy effects. The latter contribution always favours small droplets.
Since the driving force towards smaller droplets due to the entropy of self-assembly increases in magnitude with decreasing droplet volume fraction, the average size of the droplets depends on φ d and the droplets increase in size with increasing surfactant (or droplet) concentration [cf. further below].
Equation (17) may be solved so as to give where we have introduced the dimensionless parameter [6,23] as well as the quantity Except for values of α close to zero, the following two approximations of Eqs. (18) and (20), respectively, are accurate  The entropy parameter k s may be added to the sum of twice the bending rigidity and the saddle-splay constant to make up the effective bending constant k ef f ¼ 2k c þ k c þ k s . Hence, k s has a similar impact on microemulsion droplet properties, such as average size and polydispersity (see further below), as k c and k c , although it is not (primarily) related to bending elasticity properties. From a physical point of view, it is a measure of the entropy of self-assembly and, since this process is unfavourable from an entropic point of view, k s is always found to be a negative property (cf. Figure 3). Moreover, the entropy of selfassembly increases in magnitude as the concentration of surfactant decreases. As a result, k s becomes a strong function of volume fraction of surfactant and its magnitude increases (that is k s becomes more negative) as a solution of microemulsion droplets is diluted.

Average droplet size
The following expression for the (volume weighted) average droplet radius may be derived from the size distribution in Eq. (16) where the quantity only depends on α. For sufficiently large values of α, Eq. (24) may be simplified so as to give Figure 3. The self-assembly entropy parameter k s plotted against the volume fraction φ s of self-assembled surfactant, in accordance with Eq. (20). The three bending elasticity constants k c , k c and H 0 assume identical values as in Figure 2.
Combining Eqs. (23) and (25) gives the following rather simple expression for the average droplet radius According to Eqs. (23) and (26), the average droplet size depends on k eff as well as on k c H 0 and in Figure 4 we have plotted 〈R〉 against k eff for a certain value of k c H 0 corresponding to a gemini surfactant that forms elongated micelles in water [13]. It is seen that the more simple equation (26) is an excellent approximation in virtually the entire regime of k eff values. As expected, 〈R〉 is found to increase considerably with increasing values of k c as well as decreasing values of H 0 . The dependence on the more shape-related bending rigidity is more complex.
For the case where k c ≫ k c + k s , Eq. (26) simplifies to 〈R〉 = 1/H 0 and droplet size only depends on spontaneous curvature. For small k c values, however, the size of the droplets is explicitly determined by the quantity k c times H 0 .
In addition to the three bending elasticity constants, the droplet size also depends on the selfassembly entropy parameter k s . The always negative parameter k s reduces the size of the droplets and its effect becomes stronger as the concentration of droplets decreases. As a result, microemulsion droplets grow in size with increasing surfactant concentration (cf. Figure 2). We may note that the expression 〈R〉 ¼ ð2k c þ k c Þ=2k c H 0 , according to which 〈R〉 does not  depend on surfactant concentration, was derived by Safran [24] for the case where self-assembly entropy effects were neglected. This is the optimal value of 〈R〉, obtained for the case where only bending elasticity properties are considered. As a matter of fact, the Safran equation is recovered as a special case of Eq. (26) in the limit 2k c þ k c ! ∞, as k s becomes negligible. Since k s is always negative, the size of the droplets must be smaller than the value predicted by the Safran equation (cf. Figure 2).

Polydispersity
Like average droplet size, microemulsion droplet polydispersity may be predicted from the present model. From the size distribution function in Eq. (16), we may derive the following expression for the average of the squared radius Hence, the polydispersity in terms of the relative standard deviation equals in accordance with Eqs. (23) and (27) and where the definition of x is given in Eq. (24). Equation (29) may be simplified so as to give the following excellent approximation for the polydispersity as a function of α = 4πk eff /kT (cf. Figure 5).
The polydispersity according to Eqs. (29) and (30), respectively, is plotted in Figure 5. Notably, the droplet polydispersity does not at all depend on spontaneous curvature, but is seen to be a sole function of k ef f ¼ 2k c þ k c þ k s . This means that k eff may be directly determined from experimentally available polydispersities. It is seen that σ R /〈R〉 decreases with increasing values of the effective bending constant and it approaches σ R =〈R〉 ¼ ffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi kT=8πk ef f p in the limit of large k eff . The latter expression (with k s = 0) was first derived by Milner and Safran for the case where self-association entropy effects were neglected [7]. Since the entropy parameter k s depends on surfactant concentration [cf. Eq. (22)], σ R /〈R〉 according to Eq. (30) is predicted to decrease with increasing surfactant concentration.

Interfacial tension
The (planar) interfacial tension between the water and oil phases in a microemulsion in equilibrium with excess solvent phase is an experimentally available quantity. Hence, it has been observed that unusually small interfacial tensions are, in general, generated in microemulsion systems. The values are found to decrease in magnitude as the droplet size increases and may approach as low values as 10 À3 -10 À2 mNm À1 in non-ionic surfactant microemulsion systems at the phase inversion temperature (PIT) [10,25]. From our present model we obtain the planar interfacial tension γ p from the definition of k eff in Eq. (19), i.e.
Similar to droplet size, the interfacial tension depends on the three bending elasticity constants as well as the entropy parameter k s . In Figure 6 we have plotted γ p against the effective bending constant k eff for some different values of the spontaneous curvature H 0 . It is seen that γ p decreases significantly as k eff is increased or H 0 is reduced. For spontaneous curvatures approaching the limit H 0 = 1/4ξ, where bilayer structures start to form in aqueous solvents [26], γ p may drop well below 1 mNm À1 . Notably, γ p decreases with increasing surfactant concentration and droplet volume fraction since it, according to Eq. (31), depends on the entropy parameter k s . By means of combining Eqs. (26) and (31), we may eliminate k c H 0 to deduce the following expression that relates interfacial tension and droplet size According to Eq. (32), the interfacial tension must decrease in magnitude as the size of the droplets increases (cf. Figure 7). This behaviour agrees very well with experimental observations that surfactants forming larger microemulsion droplets are, in general, found to have smaller interfacial tensions. γ p is seen to drop below 1 mNm À1 as the radius exceeds a few nanometres for some typical value of 2k c þ k c . Moreover, the interfacial tension is seen to fall below about 0.01 mNm À1 as the droplet size approaches orders of magnitudes equal to 1 μm. These values agree very well with observations of minimum interfacial tensions obtained for microemulsions formed by temperature sensitive non-ionic surfactants at PIT, where the microemulsion contains equal amounts of oil and water [10].
Since k eff is directly related to the droplet polydispersity, it may, in principle, be eliminated by means of combining Eqs. (30) and (32) to yield a relation between the three experimentally available quantities 〈R〉, σ R /〈R〉 and γ p . The following approximate expression (valid for polydispersities lower than about σ R /〈R〉 = 0.1) is obtained by means of inserting in Eq. (32) Figure 6. The planar interfacial tension γ p of microemulsions plotted against effective bending constant k eff in accordance with Eq. (31) for some different values of the spontaneous curvature H 0 .
The relation between interfacial tension, droplet radius and polydispersity may also be expressed in terms of the dimensionless parameter [27] w 8π Combining Eqs. (26), (30) and (31), the following approximate expression for w may be derived In the limit k eff ! ∞ (equivalent to σ R /〈R〉! 0), Eq. (35) may be further simplified so as to give w = 1, the expression of which has been derived by Gradzielski et al. [27] for the case where selfassociation entropy effects were neglected. It is necessary that γ p assumes low values for a system of small droplets, with a large overall interfacial area, to form. The ultra-low values of interfacial tensions experimentally observed in microemulsion systems appear as a consequence of the comparatively large interfacial area between oil and water that is formed in colloidal systems consisting of small droplets. Microemulsions are thermodynamically equilibrated systems and the  disadvantageous (positive) free energy of increasing the interfacial area is exactly balanced by the favourable (negative) free energy of mixing oil droplets with solvent water. The entropy of mixing increases with the number of dissolved droplets, equivalent to a reduction in droplet size, and the corresponding contribution to the free energy becomes increasingly more negative with decreasing droplet size. As a result, the interfacial tension must decrease in magnitude with increasing droplet size in accordance with Eq. (33).

Solubilisation
A notable property of microemulsion systems is the ability to dissolve hydrophobic components which otherwise are insoluble in an aqueous solvent. A quantitative measure of the solubilisation capacity may be considered as the ratio of the volume of oil (V oil ) that is solubilized and the volume of surfactant (V surf ) present in the solution. Hence, we may define the solubilisation capacity as follows Since the surfactant is exclusively located at the droplet interface whereas the oil molecules are confined to the interior of the droplets, σ becomes strongly dependent on droplet size. The volumes of solubilized oil and aggregated surfactant, respectively, are obtained from the following geometrical relations valid for spherical droplets with an average radius 〈R〉 and an interfacial surfactant monolayer with thickness ξ (cf. Figure 1) where N drop is the total number of droplets. Combining Eqs. (36)-(38) gives the following expression for the solubilisation capacity as a function of the single quantity δ defined as The second equality in Eq. (40) is obtained by means of employing the expression in Eq. (26) for the average droplet radius.
In Figure 8, we have plotted the average droplet radius 〈R〉 against the volume fraction of surfactant self-assembled in the droplet interfaces. It is seen that the droplet size increases with an increasing surfactant concentration as a result of the self-association entropy effects discussed above. This means that we obtain an enhanced solubilisation effect as the amount of added surfactant is increased. The number of droplets, as well as the average size of the droplets, increases as the surfactant concentration is raised and both effects contributes to an increased solubilisation capacity.
The solubilisation capacity σ, as defined in Eq. (39), is plotted against the surfactant concentration in Figure 9. It is seen that σ increases considerably with increasing surfactant concentration. σ strongly depends on the different bending elasticity constants. The bending rigidity k c and saddle-splay constant k c have been found to be of the order of magnitude kT for many common surfactants and the possibility to tune the values of them is usually limited [11]. The spontaneous curvature H 0 , on the other hand, strongly depends on the surfactant HLB value and may assume values from close to zero, for amphiphilic molecules forming bilayers, to much higher values for surfactants forming small micelles. In Figure 9 we have chosen to compare σ for three different values of H 0 . The smallest value ξH 0 = 1.55 À1 = 0.154 > 0.25 corresponds to a (gemini) surfactant that form elongated micelles in an aqueous solution and the solubilisation capacity is found to be rather small, with σ considerably lower than unity. However, σ is found to increase substantially as the spontaneous curvature is decreased to ξH 0 = 0.25, which corresponds to the point of transition from micelles to bilayers in an aqueous solvent in the absence of oil [26]. Further decreasing H 0 to ξH 0 = 6.5 À1 = 0.645, i.e. to values where rather large bilayer structures are expected to form, the solubilisation capacity σ exceeds unity, which means that the amount of oil that is dissolved exceeds the amount of surfactant present in the solution.
As a matter of fact, these results agree very well with observations that the solubilisation capacity of decane (oil) in mixed tetradecyldimethylamine oxide (TDMAO)-hexanol microemulsion droplets was found to be significantly larger at high hexanol compositions, where vesicles are formed in the absence of oil, as compared to compositions where mixed micelles form in the absence of oil [5].

Summary
A comprehensive theoretical model for the formation of microemulsion droplets has been presented. The theory is based on thermodynamics of self-assembling surfactant molecules aggregated in the droplet interfaces, combined with bending elasticity properties of the surfactant monolayer as taken into account by the three parameters spontaneous curvature, bending rigidity and saddle-splay constant. It relates properties depending on surfactant chemical structure with experimentally available properties of the droplets, such as average size, polydispersity, interfacial tension and solubilisation capacity. It has recently been demonstrated that all three constants H 0 , k c and k c may be determined from measurements of growth behaviour of surfactant micelles in aqueous solutions (in the absence of oil). [12,13,22] Hence, from the present model it is possible to predict, for instance, the solubilisation capacity of oil of a certain surfactant by studying the growth behaviour of micelles that is formed by the surfactant in the absence of oil.
Likewise, it is possible to determine the various bending elasticity constants by means of measuring the average droplet radius, polydispersity and/or interfacial tension of a microemulsion system. We have included the important contribution of entropy of self-assembly and, as a result, we obtain quantitative expressions that differ from previously derived ones [7,24,27] where selfassociation entropy effects were omitted. Nevertheless, the previously derived expressions for average droplet radius, polydispersity and interfacial tension are all recovered as special cases from our model in the limit 2k c +k c ! ∞ . Taking into account self-association entropy effects, quantified in terms of the parameter k s , have the impact on the model to predict smaller values of the droplet size as compared to the case where the entropy of self-assembly is neglected. Moreover, k s explicitly depends on surfactant concentration (k s increases with increasing concentration), implying that the droplet size is predicted to increase, whereas the polydispersity decreases, with increasing surfactant concentration. Moreover, the model predicts ultra-low values of the interfacial tension that agrees excellently with experimental measurements [10,25]. Likewise, the capacity to solubilize a hydrophobic component (oil) is predicted to increase dramatically by means of reducing the spontaneous curvature (equivalent to decreasing the HLB value) of the surfactant, in agreement with experiments. In addition, self-association entropy effects are shown to generate an enhanced solubilisation capacity as the size of the microemulsion droplets is found to increase with increasing surfactant concentration.