when dissolved in aqueous solvent , lipids self - assemble into an amazing variety of different structures as a function of concentration and temperature @xcite .
most prominent is the lamellar phase , which consists of a stack of lipid bilayers separated by layers of water .
in fact , the lipid bilayer constitutes the basic building block of biological membranes ; in eucaryotic cells , it not only envelopes the cell itself , but also separates its nucleus and its organelles from the cytosol , and builds up the endoplasmic reticulum , the golgi apparatus and the uni - lamellar vesicles which serve as reaction chambers and transport vehicles .
the predominance of the lamellar phase at ambient temperatures derives from the fact that in contrast to surfactants , which usually have large headgroups and form micelles , lipids have rather bulky hydrocarbon chains .
a so - called _ bilayer lipid _ has a chain region cross section which is of similar size as the area per head group ; thus the spontaneous curvature is small and the lipid forms a lamellar phase .
if temperature is increased , the size of the head groups hydration shell decreases , spontaneous curvature increases and the lipid monolayers tend to curve towards the water regions . since the tendency of lipid monolayers to curve
is frustrated in the lamellar phase , it is often found to transform into a inverse hexagonal phase which consists of a two - dimensional array of inverse cylindrical micelles .
spontaneous curvature can also be increased by changing molecular architecture , by adding lipids with bulkier chains or by increasing salt concentration ( since salt screens electrostatic repulsion between charged headgroups ) @xcite .
lipids with spontaneous curvature are often called _ non - bilayer lipids_. for large spontaneous curvature , a cubic phase of inverse spherical micelles is sometimes observed @xcite .
the inverse phases have been denoted _ type ii ( water - in - oil ) phases _ in contrast to _ type i ( oil - in - water ) phases _ which dominate surfactant phase behavior @xcite . although there are some lipids which also form type i phases , we focus here on type ii phases since these are the ones formed by the majority of biologically relevant lipids .
one of the intriguing aspects of the polymorphism of lipid - water mixtures is the existence of yet another structural type , which often occurs between the lamellar and the inverse hexagonal phase . like the lamellar phase , these _ inverse bicontinuous cubic phases _ ( ibcps ) consist of lipid bilayers .
however , now a _
lipid bilayer extends throughout the whole of space , and divides it into two disconnected but interpenetrating labyrinths of cubic symmetry , which both are filled with water .
thus the structure can also be considered to be built from water channels ( that is inverse cylindrical micelles ) as is the inverse hexagonal phase , only that now the channels meet in vertices with a coordination number related to the topology of the cubic phase .
although the ibcps unite the structural elements of the lamellar and the inverse hexagonal phases , in contrast to them they are optically isotropic and highly viscous . like the inverse hexagonal phase ,
the ibcps are stabilized by spontaneous curvature .
therefore they are formed by non - bilayer lipids and their stability can be controlled by changing the lipids spontaneous curvature .
the generic occurrence of ibcps in lipid - water mixtures and their property of dividing space into two interwoven aqueous compartments lends itself readily to speculations on their biological relevance , e.g.  on their possible occurrence in the endoplasmic reticulum or the golgi apparatus @xcite .
a recent survey of a vast number of tem pictures suggests that ibcps are indeed ubiquitous in biological systems @xcite .
the abundance of non - bilayer lipids in biological systems is a longstanding puzzle , and several studies have explored the possible relationship between lipid polymorphism and biological membrane function ( for reviews see refs .
@xcite ) . in particular
, it was shown that certain cells like _ e. coli _ regulate the relative amount of bilayer and non - bilayer lipid in order to control membrane morphology and function of membrane proteins .
_ e. coli _ mutants lacking the non - bilayer lipid phosphatidylethanolamine were found to be severely impaired in regard to protein transport across the plasma membrane @xcite .
luzzati and coworkers have pointed out that if any ibcp had a biological function , it should be stable at physiological temperatures and in excess water @xcite .
ibcps are also relevant for biotechnological applications , e.g.  they have been used recently as artificial matrix which enables membrane proteins to crystallize in a three - dimensional array @xcite .    since ibcps consist of a lipid bilayer separating two labyrinths of water , is it reasonable to assume that the two monolayers are arranged in a way which is _ locally _ symmetrical about their mid - surface . in differential geometry ,
a surface which curves to both sides in the same way is known as a _ minimal surface _ ,
i.e.  a surface for which the mean curvature vanishes at each of its points
. therefore the mid - surfaces of the ibcp s lipid bilayers are often modeled as _ triply periodic minimal surfaces _
( tpms , sometimes also denoted ipms for _ infinite periodic minimal surfaces _ ) of cubic symmetry @xcite .
indeed it has been confirmed by a thorough analysis of electron density maps derived from x - ray data that the mid - surfaces of these structures can be very close to minimal surfaces @xcite .
note that the stronger requirement of the monolayers _ global _ symmetry about the mid - surface leads to a certain subclass of cubic tpms which is called _ balanced _ since then the two water labyrinths have to be congruent to each other .
since the seminal work of schoen @xcite , it is known that there exists a large variety of different cubic tpms . before schoen s
work , only p , d and c(p ) were known from the 19th century work of schwarz and his students ; then schoen described g , f - rd , i - wp , o , c - to and c(d ) .
more cubic tpms have been found later by karcher and polthier @xcite , fischer and koch @xcite and others , but most of them seem to be too complicated as to be of physical relevance . in figs .
[ figureparallelsurfaces1 ] and [ figureparallelsurfaces2 ] we depict surfaces corresponding to the seven different tpms treated in this work . until now
, only g , d and p have been established to be realized as ibcp in lipid - water mixtures @xcite .
one system for which all three of these phases is stable is 2:1 lauric acid ( la)/dilauroyl phosphatidylcholine ( dlpc ) and water ( with p coexisting with excess water ) @xcite .
the only other system for which g , d and p are known to be stable is didodecylphosphatidylethanolamine and water @xcite ; however , here g exists with hardly no hydration at all , and the regions for d and p are very difficult to access experimentally , thus this system should not be compared with the theoretical results for fully hydrated phases presented here . in all other systems with ibcps ,
only one or two of them are stable ( for example , for monoolein and water , g and d are stable , with d coexisting with excess water @xcite ) . in this work
we present a model which can explain why only g , d and p are observed in lipid - water mixtures and why their actual occurrence seems to depend on the specific system under investigation .
the phase diagram predicted by our model agrees well with the experimental phase diagram for 2:1 la / dlpc and water .
exact ( _ weierstrass _ ) representations are only known for g , d , p and i - wp @xcite .
however , it has been shown in recent studies @xcite that these and other cubic tpms can be generated as isosurfaces of a density field , which minimizes the free - energy functional of a simple ginzburg - landau model for amphiphilic systems @xcite .
using the fourier representations of ref .
@xcite , it is now possible to investigate also the physical properties of the ibcps , for whose minimal mid - surfaces no exact representations are known . in particular , the question of whether these structures might correspond to stable ibcps in lipid - water mixtures can now be addressed quantitatively . in this work
we study g , d , p , i - wp , s , c(p ) and f - rd , of which g , d , p , s and c(p ) are balanced .    the most important contribution to the free energy of amphiphilic interfaces is the curvature energy , which dominates all other contributions if the radii of curvature are large compared with molecular length scales @xcite . for the inverse phases in lipid - water mixtures , it is generally assumed that another important contribution to the free energy is the stretching energy of the hydrocarbon chains @xcite . in order to attain optimal free energy ,
the neutral surfaces of the two monolayers should realize at the same time constant mean curvature @xmath0 and constant distance @xmath1 to the minimal mid - surface . since this is not possible geometrically , the free energy will always be frustrated @xcite .
the relative importance of interface ( bending ) and bulk ( stretching ) contributions in self - assembled interfaces is a long - debated issue .
for amphiphilic systems , for which the free energy contributions come mainly from the interfaces , the main problem is of a mathematical nature : the overall solution can be investigated only in a model which allows the interfacial positions to adjust themselves freely . in contrast , in diblock copolymer systems bulk contributions are much more important and models based on random polymer coils and interacting monomer densities are more adequate .
self - consistent mean - field theory for diblock - copolymer systems suggests that in this case , interfacial tension and stretching are equally important @xcite .
recently , this framework has been adapted to treat lipid - water mixtures @xcite .
this approach is complementary to the membrane approach ; it is most useful for the investigation of systems , in which the internal structure of the bilayer is not uniform .
however , self - consistent field theory for large molecules is not expected to describe amphiphilic systems well at large lattice constants , when the membranes internal structure is essentially preserved and the main contributions to the free energy depends on their shape . in particular , only in the membrane approach does it becomes possible to address the role of thermally activated membrane undulations .    in order to study interface - dominated systems
, it has been suggested in ref .
@xcite to consider the two classes of surfaces , which can completely relax one of the two relevant contributions .
_ parallel surfaces _ to the minimal midsurface realize the optimal chain length @xmath1 and are frustrated only in their bending energy , while _ constant - mean - curvature surfaces _ ( cmc - surfaces ) realize the given spontaneous curvature , but are frustrated in their stretching energy .
the mathematical properties of parallel surfaces are well known in terms of the properties of the underlying minimal surface ; the parallel - surface model has therefore been the main tool for the investigation of ibcps @xcite . up to now
, the main difficulty of this approach has been to take into account the variation of gaussian curvature over the minimal mid - surface . in earlier work
, the curvature energy was therefore expanded in moments of the distribution of the gaussian curvature . in ref .
@xcite this was done up to eighth order for d , and in ref .
@xcite up to second order for p , d and g. cmc - surfaces @xcite have the advantage that they offer an unified picture for all relevant phases  lamellar , inverse bicontinuous cubic , inverse hexagonal and inverse micellar cubic phases are modelled by planes , cubic cmc - surfaces , cylinders and spheres , respectively .
however , since no exact representations are known , the stretching frustration for ibcps has so far been evaluated numerically for d only @xcite .
we will show below that in the case of short chains ( or large lattice constants ) , stretching energy is prohibitively large , and the two monolayers of a ibcp in lipid - water systems can be modeled by parallel surfaces . in the case of long chains ( or small lattice constants ) , we expect stretching and bending to be equally important ; this agrees with the results for diblock copolymers , which are in the long chain limit .
our reasoning implies that cmc - surfaces are not a good representation of ibcps in lipid - water systems for physical reasons .
however , for small distances to the minimal mid - surface , their mathematical properties are very similar to the ones of parallel surfaces @xcite , thus many results from the parallel surface model will carry over to the cmc - model .
in general , however , the parallel surface model is more appropriate both for physical and mathematical reasons . moreover , there is strong experimental evidence for it , both from a detailed reconstruction of electron densities @xcite and from analysis of swelling data @xcite .
in contrast to earlier work , we evaluate the curvature energy in the parallel surface model without any expansion and to high numerical accuracy by first calculating the distribution of gaussian curvature over the minimal mid - surface . for g ,
d , p and i - wp this can be done from their weierstrass - representations .
we also consider the structures s , c(p ) and f - rd by calculating the distribution from the representations which we obtained from ginzburg - landau theory .
our main result is that the free - energy densities of the structures g , d and p are considerably lower than those of the other investigated structures due to their narrow distribution of gaussian curvature over the minimal midsurface .
we show that this result persists when thermal membrane undulations and higher order terms in the bending energy are considered , and argue that potential additional contributions to the free energy are unlikely to change it .
this explains why only g , d and p have been observed in lipid - water mixtures .
in fact , the calculated phase diagram agrees nicely with the experimental one for 2:1 la / dlpc and water , for which these three phases coexist .
we show that due to the existence of a bonnet transformation between g , d and p , these phases coexist along a triple line in our model .
simultaneously , p coexists with an excess excess water due to a mechanism called _
emulsification failure_. the bonnet transformation also implies that the free - energy densities of g , d and p scale as a function of concentration with a universal geometrical quantity , which we term _ topology index_. since the topology index decreases from g to d to p , the gyroid g is most prominent , followed by smaller regions of stability for d and p at higher water concentrations .
any additional contribution to the free energy , which introduces a new length scale , is expected to change the delicate balance between these three phases ; this includes stretching contributions , van der waals and electrostatic interactions .
since such additional effects are specific to a given experimental system , this result explains qualitatively why usually only one or two of these phases are observed in lipid - water mixtures .
the paper proceeds as follows . in sec .
[ sec : model ] , we specify the free - energy expression which has to be evaluated for each type of minimal mid - surface , and discuss some immediate consequences of the model .
the calculation of the distribution of gaussian curvature is explained in sec .
[ sec : weierstrass ] for p , d , g and i - wp from their weierstrass representation .
results for the gaussian - curvature distributions for s , c(p ) and f - rd are obtained in sec .
[ sec : ginzburglandau ] from the representations derived recently from a simple ginzburg - landau model .
these results are combined in sec .
[ sec : phasediagrams ] to predict phase behavior and to numerically calculate phase diagrams . in sec .
[ sec : symmetry_breaking ] , we discuss which physical mechanisms will break the bonnet symmetry between p , d and g. finally , we compare our results with experiments and conclude in sec .
[ sec : discussion ] .
the bending energy of one lipid monolayer is described by the _ canham - helfrich hamiltonian _
@xcite @xmath2 where the integration extends over the monolayer s neutral surface , which is at a distance @xmath1 from the minimal midsurface of the bilayer . in the framework of the curvature model ,
bending and stretching contributions decouple at the neutral surface @xcite ; this allows to disregard area - stretching contributions .
the neutral surface has been shown experimentally to be located close to the polar - apolar interface @xcite . in the curvature model ( [ helfrichhamiltonian ] ) , the neutral surface is characterized by mean curvature @xmath3 , gaussian curvature @xmath4 and differential area element @xmath5 ; the superscript @xmath1 indicates that the neutral surface is a parallel surface .
the model parameters are the spontaneous curvature @xmath0 , the bending rigidity @xmath6 and the saddle - splay modulus @xmath7 of the monolayers . for lipids ,
the bending rigidity @xmath6 for a monolayer is of the order @xmath8 .
usually the saddle - splay modulus @xmath7 is assumed to have a small negative value ; in fact the curvature model without any additional constraints is well - defined only for @xmath9 @xcite .
the spontaneous curvature @xmath0 can be considered to result from the mismatch between optimal head group area , volume incompressibility of the chains and optimal chain length . for
non - ionic lipids , it depends mainly on temperature , which changes the head group hydration and therefore optimal head group area .
since lipids typically form lamellar and inverse phases , where the monolayers prefer to bend towards their polar sides , we define positive curvature to be towards the water regions .
we assume that spontaneous curvature scales with temperature as @xmath10 , since such a linear relationship seems to hold very well for ternary surfactant systems @xcite .    the stretching energy can be assumed to be harmonic about the average chain length @xmath1 , so that @xmath11 where @xmath12 is the local chain length .
the stretching modulus @xmath13 is itself a function of the average chain length @xmath1 .
the relevant scaling law has to be obtained from a microscopic model ; here , we model the lipid monolayer as a gaussian polymer brush grafted to the polar - apolar interface , but more realistic models give similar results @xcite .
the stretching energy per chain reads @xmath14 , where @xmath15 is the ( microscopic ) kuhn length and @xmath16 the number of kuhn segments . enforcing the volume constraint @xmath17 ( where @xmath18 is the grafting density ) , and transforming from energy per chain to energy per area , results in a potential of the form with @xmath19 .
thus the stretching modulus @xmath13 scales inversely with average chain length @xmath1  stretching becomes very difficult for short chains @xcite .
although the bending energy is well known not to depend on microscopic details , for our purpose it is important to note that the bending rigidity @xmath6 also depends on the average chain length @xmath1 .
starting from the stretching energy of and taking into account how the chain volume varies with curvature , results in the bending energy of with @xmath20 @xcite .
the result that the bending rigidity @xmath6 scales with the third power of film thickness @xmath1 is well known also from the elasticity of thin solid sheets .
we now can estimate the relative importance of bending and stretching contributions .
it is known that for both cmc - surfaces and parallel surfaces , the normalized standard deviations of @xmath12 and @xmath3 , respectively , are almost independent of hydrocarbon volume fraction @xcite , so that @xmath21 and @xmath22 , where @xmath23 is the lattice constant
. therefore @xmath24 .
thus , for small chain length @xmath1 ( large lattice constant @xmath23 ) , the stretching term is much larger and the chains can be assumed to be of constant length .
this corresponds to the well - known result that the curvature energy is the relevant contribution to the free energy as long as the curvature radii ( that is the lattice constants ) are much larger than any molecular length @xcite . only for large chain length @xmath1 ( small lattice constant @xmath23 ) ,
the bending energy becomes comparable and both terms have to be considered .
experimentally , this corresponds to the fact that for small water content ( i.e. chain length @xmath1 of comparable size as lattice constant @xmath23 ) , the bicontinuous cubic phases become unstable with respect to the hexagonal phase , whose geometrical properties require larger chain stretching . note that the relative importance of stretching and bending is asymmetric : since the chain length @xmath1 can not be larger than the lattice constant @xmath23 , there is no regime in which bending can be neglected .
thus cmc - surface will never be a good approximation for monolayers in lipid - water systems .
however , the stretching frustration can be relaxed by swelling the bilayer with oil ; therefore cmc - surfaces are good approximations for amphiphilic monolayers in ternary systems with oil , water and amphiphile @xcite .
similar considerations of the relative importance of stretching and bending can be found also in the context of large membrane inclusions ( like transmembrane proteins ) where chains have to stretch in order to decrease the hydrophobic mismatch at the inclusion boundary @xcite .
for the rest of this paper , we proceed within the framework of the parallel surface model .
we assume that the interfaces between the tail and head group regions of the two lipid monolayers in a icbp are located at distances @xmath25 away from a given cubic tpms ; typically @xmath1 has values between one and two nanometers .
note that for a tpms the mean curvature satisfies @xmath26 everywhere ( where @xmath27 and @xmath28 are the two principal curvatures ) , while the gaussian curvature @xmath29 varies as function of the position and is only restricted to satisfy @xmath30 .
the differential area element @xmath31 , the mean curvature @xmath32 and the gaussian curvature @xmath33 on the parallel surfaces then follow as functions of @xmath1 and the quantities @xmath5 and @xmath4 on the minimal mid - surface @xcite , with @xmath34 note that these formulae are special cases of steiner s theorem of integral geometry ; they are therefore no approximations for small @xmath1 , but exact .
since minimal surfaces have @xmath30 , positive mean curvature is defined here to correspond to positive @xmath1 .
the effective bending energy of the lipid bilayer is only a function of the gaussian curvature @xmath4 of the minimal midsurface since its mean curvature @xmath3 vanishes . in an expansion for small chain length @xmath1 ,
it can be expressed as an infinite series in powers of @xmath4 .
using and yields up to order @xmath35 @xmath36 therefore the effective saddle - splay modulus for the lipid bilayer , @xmath37 , is not just the sum of the monolayer s moduli , but is corrected to higher positive values due to the presence of the spontaneous curvature @xmath0 @xcite .
we conclude that as long as @xmath38 , the ibcps are favored over the lamellar phase , since the preferred curvature of the monolayers translates into a topological advantage of saddle - type bilayer structures .
the third term in favors saddle - type structures in any case , and also removes the degeneracy which minimal surfaces might experience under deformations which preserve the minimal surface property .
this has been noted in a seminal work by bruinsma @xcite ; in his notation , we have @xmath39 for the elastic modulus of the @xmath40-term .
the volume fraction occupied by the hydrocarbon can be calculated as @xmath41 where we have used the gauss - bonnet theorem @xmath42 . @xmath43 and @xmath44 denote surface area and euler characteristic of the given cubic tpms in the conventional unit cell with unit lattice constant . for the following ,
it is useful to introduce the _ topology index _ @xmath45 .
this quantity is independent of scaling and choice of unit cell and characterizes a given tpms s topology in an universal way ; the higher its value , the smaller is the porosity and the larger the specific surface area of a tpms .
it is the only quantity which characterizes a two - dimensional minimal surface in three - dimensional space independently of lattice constant and choice of unit cell , and its relevance for the description of bicontinuous cubic phases in amphiphilic systems has been discussed before @xcite .
the topology index of the various tpms is of order @xmath46 and is given in tab .
[ tableindex ] .
the gyroid g is expected to have the highest value since it divides space into labyrinths with threefold coordinated vertices .
s has a similar value for @xmath47 , since it features a mixture of three- and fourfold coordinated vertices . then comes the double diamond structure with fourfold coordinated vertices , while all other structures have higher values . in tab .
[ tableindex ] we collect the values of @xmath43 , @xmath44 and @xmath47 for the cubic tpms investigated in this work as obtained from our numerical representations in the framework of a simple ginzburg - landau theory ( see below ) .
exact values are known for g , d , i - wp , p and c(p ) and given in tab .
[ tablenumbers ] . comparing numerical and exact values shows that for most structures , our numerical results are quite accurate . due to numerical limitations for the more complicated structures , the values for s and f - rd are less precise .
in particular , we expect s to have a smaller value for @xmath47 than g.    equation ( [ hydrocarbonvolume ] ) can be inverted numerically to give the dimensionless function @xmath48 , the lattice constant @xmath23 in units of the chain length @xmath1 , as a function of hydrocarbon volume @xmath49 . for small @xmath49
we find @xmath50 in fig .
[ figurelatticeconstant ] we plot the first - order approximation @xmath51 of for all structures considered ; it can be hardly distinguished from the full curve over the full range of @xmath49 . by visual inspection of our numerical representations
( see below ) , we find that for all structures considered , parallel surfaces do not self - intersect as long as @xmath52 ; ibcps can accommodate large amounts of hydrocarbon since their geometry is so close to that of the lamellar phase .
a typical value for @xmath1 is @xmath53 , so that for @xmath54 we have @xmath55 , @xmath56 , @xmath57 , @xmath58 , @xmath59 , @xmath60 and @xmath61 .
note that at a given @xmath49 , the ratio of the lattice constants of two coexisting ibcps is simply the ratio between their scaled surface areas @xmath43 .
in figs .  [ figureparallelsurfaces1 ] and [ figureparallelsurfaces2 ] we depict for each of the structures considered one of the two monolayers for the hydrocarbon volume fraction @xmath54 .
then the volume fraction of each of the two labyrinths is @xmath62 .
we are now in a situation to rewrite the curvature energy from for our purpose .
we consider a typical experimental situation , in which temperature is controlled at constant volume . as explained above , changing temperature @xmath63 amounts to changing spontaneous curvature @xmath0 . under the assumption of incompressibility of both components
, there is only one degree of freedom for composition , which we take to be the hydrocarbon volume @xmath49 ( which is equivalent to the amount of interfacial area ) . dividing the free energy by the constant volume
, we arrive at the free energy per _ unit volume _ as a function of @xmath49 .
however , @xmath49 is not a controlled quantity , since at two phase coexistence , it is not clear how it will distribute into the two phases . the controlled quantity can be considered to be the chemical potential coupled to @xmath49 ( that is the derivative of free energy  or free energy volume density  for hydrocarbon volume fraction  or interfacial area ) , and therefore phase transitions have to be calculated by using the maxwell construction for the free energy density as a function of @xmath49 ( which amounts to requiring the same chemical potential for both phases ) .
we want to remark parenthetically that the free energy per _ unit area _ can not be used to calculate phase coexistences , since an exchange of area , controlled by the chemical potential of the lipid , has to be allowed between coexisting phases .
we write the free energy volume density dimensionless by using a factor @xmath64 . for the lamellar phase the hydrocarbon volume fraction @xmath65 ,
thus its dimensionless free - energy density follows from to be @xmath66 the larger @xmath49 , the smaller @xmath23 and the more frustrated curvature energy per volume accumulates .    for ibcps
, we again use and in order to write the dimensionless free - energy density as a function of hydrocarbon volume fraction @xmath49 and the properties of the bilayers minimal mid - surface : @xmath67 where again we have used the gauss - bonnet theorem and defined @xmath68 . as mentioned above
, the curvature model without any additional constraints is well - defined for @xmath69 and bicontinuous phases are favored for @xmath70 .
@xmath43 , @xmath71 and @xmath44 refer to a conventional unit cell . using to lowest order , @xmath51 , and defining @xmath72 , we rewrite as @xmath73 in eq .
( [ freeenergydensitybicont2 ] ) and below , spontaneous curvature @xmath0 is measured in units of @xmath74 .
since for the minimal midsurface both @xmath71 and @xmath44 are negative , @xmath75 is a positive quantity which varies over the surface .
note that @xmath76 due to the gauss - bonnet theorem . the free - energy density @xmath77 of the lamellar phase , eq .
( [ freeenergydensitylamellar ] ) , is a special case of @xmath78 , with @xmath79 . in this work
, we will treat the free - energy density without any further approximations .
previously , similar expressions have been expanded in small @xmath49 ( which is equivalent to small @xmath1 ) @xcite , similarly as discussed in .
in fact , each term @xmath80 in an expansion for small @xmath1 corresponds to a term @xmath81 in an expansion for small hydrocarbon volume @xmath49 . since @xmath82 ,
higher - order terms in would result in a free - energy density which has a similar structure as . as we will see below
, this means that our results will remain valid in the case that higher - order terms are included in the curvature energy of the monolayers .
it is well known that thermal fluctuations can contribute significantly to the free - energy of amphiphilic systems . since we consider the regime where the stretching energy dominates , undulations are the favored fluctuation modes . for the lamellar phase , they give rise to steric repulsion @xcite , which in our case can be written as @xmath83 where @xmath84 from monte carlo simulations @xcite and field - theoretical calculations @xcite .
for an ibcp , steric repulsion is hardly relevant since the lateral correlation length @xmath85 is of the order of the lattice constant @xmath23 .
therefore , the perpendicular correlation length @xmath86 , and few membrane collisions should occur .
however , fluctuations of the lipid bilayer lead to a renormalization of @xmath6 and @xmath7 @xcite for all phases .
since @xmath6 multiplies the average mean curvature squared , which vanishes for the mid - plane , only the renormalization of @xmath87 at length scale @xmath88 has to be taken into account .
we identify the typical length scale of a cubic structure with @xmath89 .
this implies @xmath90 , so that @xmath91 in eq .  ( [ freeenergydensitybicont2 ] )
gets renormalized , with @xmath92\ .\ ] ] since the renormalization increases @xmath7 , this effect increases the topological advantage of the saddle - shaped structures as does the spontaneous curvature .
note that for the lamellar phase , these effects might favor the formation of wormholes @xcite . here , we neglect this aspect , as well as the acoustic modes of the bicontinuous structures at large wavelengths , which have been discussed by bruinsma @xcite .
note also that the thermal fluctuations also give a contribution to the free energy , which is linear in the membrane area .
we have omitted this term here , since it can be absorbed into the chemical potential of the lipid .
having defined the model , we now turn to its phase behavior .
we first note that the dimensionless mean curvature averaged over the parallel surface follows from as @xmath93 where again @xmath94 and the gauss - bonnet theorem has been used .
it is plotted in fig .
[ figuregeometry1]a .
note that all curves fall on a universal curve when the hydrocarbon volume fraction @xmath49 is scaled with the topology index @xmath47 .
it reaches the value @xmath0 at @xmath95 thus , @xmath0 is realized in a sequence given by increasing topology index ; in particular , the gyroid structure g reaches @xmath0 at the highest value of @xmath49 .
we conclude that if all considered structures were stable , they would appear in the sequence g - s - d- i - wp - p etc with increasing water concentration .    the stability of the different phases is determined by the free - energy density given in .
the second ( topological ) term is easy to understand .
the topological properties of the parallel and mid - surfaces are the same : for negative saddle - splay modulus @xmath7 ( @xmath96 ) , the different ibcps are favored according to their values of the topology index @xmath47 . here
the g - structure performs best since it has the lowest porosity .
the first term in is more complicated , since the mean - curvature properties of the parallel surfaces translate into the gaussian - curvature properties of the minimal mid - surface in a quite complicated way .
as discussed above , the spontaneous curvature of the monolayers implies a topological advantage of the bicontinuous phases . for a more detailed analysis , which includes all terms of the expansion
, we note that the first term in measures the standard deviation of the parallel surface s mean curvature from the spontaneous curvature .
since it follows from that the distribution of mean curvature over the parallel surface , @xmath97 , is proportional to the distribution of gaussian curvature over the minimal surface , @xmath98 , we expect that those structures will be more favorable which have small standard deviations for their gaussian curvature distributions @xcite . in order to quantify this concept
, we now turn to the distributions of gaussian curvature which later will allow us to evaluate without any further approximations .
weierstrass representations are known for p , d and g @xcite and i - wp @xcite . for each of these tpms
, a fundamental domain can be identified , so that the rest of the surface follows by replicating it with the appropriate space group symmetries ( @xmath99 , @xmath100 , @xmath101 and @xmath99 , respectively ) .
the weierstrass representation is a conformal mapping of certain complicated regions within the complex plane onto the fundamental domain : @xmath102 where @xmath103 are the internal ( and conformal ) coordinates of the minimal surface .
the geometrical properties of such a surface follow as @xmath104 with @xmath105 .
obviously the ( isolated ) poles of @xmath106 correspond to the flat points ( @xmath107 ) of the minimal surface .
only few choices of @xmath106 yield embedded minimal surfaces .
the ones for d and p have been known since the 19th century from the work of schwarz : for d it is @xmath108 .
p corresponds to the same region and follows simply by the bonnet transformation @xmath109 with @xmath110 .
equation ( [ weierstrassdarstellung2 ] ) implies that p and d have the same metric and the same distribution of gaussian curvature .
however , since they map differently into embedding space , they have different space groups and lattice constants .
the gyroid g was discovered in 1970 by schoen @xcite as another bonnet transformation of d , with @xmath111 .
the weierstrass representation for i - wp was found only recently @xcite .
if one of its poles is chosen to be at infinity , one has @xmath112 . in fig .
[ figureweierstrass ] we show the fundamental domains in the complex plane of the different weierstrass representations . in tab .
[ tablenumbers ] we collect the values for scaled surface area @xmath43 and lattice constant @xmath23 as they follow from the weierstrass representations .
we also give the number @xmath16 of replications of the fundamental domain needed to build up the surface in one conventional unit cell and the values for the euler characteristic @xmath44 which follows from the topology generated in this process .
note that for d one can choose another conventional unit cell which is contained @xmath113 times in the one chosen here .
this has been done in refs .
@xcite ; then @xmath44 , @xmath114 and @xmath4 scale as @xmath115 , @xmath116 and @xmath117 , respectively .
in particular , one then has @xmath118 and @xmath119 .
our choice is motivated by the fact that it is the relevant one also for ternary systems .
we define the distribution function of gaussian curvature , @xmath120 over some surface parametrized by internal coordinates @xmath103 , where in our case @xmath121 and @xmath122 are determined by . for the generating functions @xmath123 detailed above ,
it is not possible to calculate @xmath124 analytically .
therefore , we evaluate it numerically by considering a histogram @xmath125 for a set of discrete gaussian curvature values @xmath126 spaced equidistantly with @xmath127 .
we cover the @xmath103-region corresponding to the fundamental domain with a square grid of @xmath128 points .
for each point , we calculate @xmath5 and @xmath4 from and add @xmath5 to the @xmath129 which corresponds to the @xmath130 with @xmath131 $ ] . in order to obtain results for one conventional unit cell , the values @xmath132
are then multiplied by the number @xmath16 of replications and the values @xmath126 are scaled with @xmath133 ( both @xmath16 and @xmath23 are given in tab .
[ tablenumbers ] ) . in practice
we use @xmath134 and @xmath135 ; the latter value amounts to about @xmath136 bins for the histogram . by this procedure , we have @xmath137 the moments @xmath138 of the distribution @xmath139 follow as @xmath140 in particular , @xmath141 and @xmath142 .    the procedure described here
was carried out to obtain @xmath143 for d and i - wp . since
p and g are related to d by a bonnet transformation , their distributions can be obtained simply by appropriately rescaling .
for example , for the transformation @xmath144 , the @xmath130 have to be rescaled with @xmath145 and the @xmath129 with @xmath146 , as given in tab .
[ tablenumbers ] . in fig .
[ figuredistributions1]a , we plot the four different distributions obtained .
p , d and g feature one peak around a ( relatively large ) k - value , that means they are quite uniformly curved , although they have flat regions , too ( which correspond to the poles of @xmath106 ) .
i - wp is bimodal , that is it has prominently both flat and strongly curved regions .
the strong curvature of d is somehow artificial due to our choice of unit cell .
if one corrects for this by a factor @xmath147 for the k - values , it becomes clear that i - wp has a much broader distribution than p , d and g. this can also be inferred from tab .
[ tablemoments ] , where we list the first moments @xmath138 of the distributions for d and i - wp .
the ones for d are identical to the ones given in ref .
@xcite after appropriate rescaling @xmath138 with @xmath148 due to the different choice in unit cell ( the numerical deviation is less than 0.02 % even for @xmath149 ) .
the moments for p and g can be calculated from the ones for d due to the bonnet - transformation between them ; e.g.  @xmath150 with the values for @xmath16 and @xmath23 given in tab .
[ tablenumbers ] .
in particular , @xmath151 , @xmath152 and @xmath153  as can be verified from tab .
[ tablenumbers ] .
it follows immediately that the quantity @xmath72 defined for is invariant under a bonnet transformation .    in order to quantify the variance of the different distributions ,
we calculate standard deviation , @xmath154 ( here @xmath155 means area average ) . this quantity is independent of scaling , choice of unit cell and bonnet - transformation , and given in the last row of tab .
[ tablemoments ] . as expected , its value for i - wp is higher than for p / d / g .
in order to investigate the curvature properties of s , c(p ) and f - rd , we use the representations as isosurfaces to fields @xmath156 which we obtained recently as local minima of a simple ginzburg - landau theory for ternary amphiphilic systems @xcite .
the same model can also be employed for binary systems if the scalar order parameter field @xmath156 is interpreted to be the local concentration difference between water in the first labyrinth and water in the second labyrinth .
then the mid - surface of the lipid bilayer is identified with the surface @xmath157 .
the free - energy functional is given by @xmath158 = \int\ d\textbf{r}\ \left\ { ( \delta \phi)^{2 } + g(\phi ) ( \nabla \phi)^{2 } + f(\phi ) \right\}\ .\ ] ] a reasonable choice for @xmath124 and @xmath159 has been found to be @xmath160 for given model parameters @xmath161 , the free - energy functional is minimized in fourier space by implementing the correct space group symmetry and minimizing for fourier amplitudes and lattice constant with conjugate gradients . for balanced tpms , an additional black and white symmetry has to be implemented .
the resulting fields @xmath156 correspond to the various local minima .
we showed in ref .
@xcite that by tuning @xmath161 appropriately , the @xmath162-isosurfaces can be made to become very close to minimal surfaces , and generated representations for nine different tpms : g , d , p , s , c(p ) , c(d ) , c(y ) , i - wp and f - rd .
k - distribution can be obtained from these representations as follows . by triangulating the isosurfaces with the marching cube algorithm
, we numerically obtain a metric on the surface . from the fourier representations of the fields
@xmath156 , it is possible to calculate exactly the k - values on the vertices of the triangulation . using a histogram technique like in the weierstrass case ,
it is thus possible to obtain @xmath143 even for structures like s , c(p ) and f - rd , which might be of physical relevance but for which no exact representations are known .
the distributions obtained in this way for @xmath163 , @xmath164 and @xmath165 are plotted in fig .
[ figuredistributions1]b .
all three of them are roughly trimodal and feature much stronger curved parts than the ones discussed above .
this reflects their complicated shape , which makes them more difficult to access numerically than the ones which have been generated by weierstrass representations ; it also explains the difficulty to obtain smooth data for the k - distributions . in tab .
[ tablemoments ] , we give their first moments @xmath138 and the measure @xmath166 of their variance .
all of them have broader distributions than i - wp , in the sequence s - f - rd - c(p ) .
therefore we can conclude already from the distributions that the ibcps should become less favorable in the sequence p / d / g - i - wp - s - f - rd - c(p ) .
given the distribution of gaussian curvature as histogram @xmath143 , it is straightforward to evaluate the first term of the curvature energy , , of the bicontinuous phases . since it is an area average over a complicated function of @xmath71 , we have to replace @xmath167 and @xmath71 by @xmath168 and @xmath130 , respectively . in fig .
[ figuregeometry1]b , we plot the quantity @xmath169 as a function of hydrocarbon volume fraction @xmath49 for @xmath170 .
in contrast to @xmath171 , compare , no simple formula exists , and it has to be calculated numerically as @xmath172 since @xmath173 is invariant under a bonnet transformation , @xmath174 is the same for p , d and g , except for a rescaling of @xmath49 with @xmath47 . in particular , the curves reach the same minimal values and the minima are shifted to the right with increasing topology index @xmath47 .
thus the structures p , d and g perform equally well , but at different values of @xmath49 . from fig .
[ figuregeometry1]a , we see that all other structures can not reach this low value of frustration since they have much broader distributions than p , d and g ( compare tab .  [ tablemoments ] ) .
note that the scaling with @xmath175 is valid in for all structures and in only for those related by a bonnet transformation .
except for a factor of @xmath176 , the free - energy density ( [ freeenergydensitybicont2 ] ) for @xmath177 is identical with @xmath178 ( see ) . in fig .
[ figurefreeenergies ] , we show the free - energy density of all considered phases for @xmath170 , @xmath177 and @xmath179 . > from the graph , we infer the phase sequence l - g - d - p with increasing water content ; the other phases are not stable since they can not reach that small an amount of bending frustration as can p , d and g. as discussed above , implies that the sequence g - d - p is determined by the topology index .
since the free - energy density for p rises again for small @xmath49 , it would be more favorable for it to coexist with excess water .
this mechanism is known as _
emulsification failure _ in surfactant systems @xcite and allows to extend the sequence to include a p+w coexistence .
the emulsification failure can be interpreted as a maxwell construction of p with the pure water phase with @xmath180 at @xmath181 .
it follows from and the scaling of @xmath178 with @xmath175 that the free - energy density @xmath78 for p , d and g can be written in the scaling form @xmath182 with a universal function @xmath183 . for small @xmath91 , this function has @xmath184 , a minimum at @xmath185 , and diverges at @xmath186 .
a maxwell construction shows that these structures always lie on a triple line , irrespective of the values for @xmath0 , @xmath91 and @xmath6 ( compare inset of fig .
[ figurefreeenergies ] ) .
this means that if we calculated the whole phase diagram with the maxwell construction , d and p would be stable only along a line . in the following we identify phase transitions between g , d and p with the intersections of the free - energy curves , since under experimental conditions , small additional contributions to the free energy
are certain to destroy the delicate free - energy balance resulting from the bonnet transformation , and then will lead to extended one - phase regions for d and p ( or remove one or both of them from the overall phase diagram ) .    in
we use this reasoning to construct phase diagrams as a function of @xmath49 and @xmath0 for different values of @xmath91 and @xmath187 . with increasing water concentration @xmath188
we always obtain the sequence l - g - d - p - w+p .
the most stable ibcp is the gyroid structure g close to the lamellar phase at high values of @xmath49 since it has both a narrow distribution of gaussian curvature and the highest topology index . at lower concentrations of @xmath49 ,
one finds narrow regions of stability for d and p since they have narrow distributions , too , but increasingly lower values of @xmath47 .
increasing @xmath0 shifts the loci of least frustration to higher @xmath49 as evident from fig .
[ figuregeometry1]a . decreasing @xmath7 ( i.e. increasing @xmath91 ) disfavors the ibcps and increases the stability region for the lamellar phase , as can be seen by comparing a and [ figurephasediagram]b , where @xmath91 is increased from @xmath189 to @xmath190 , respectively .
however , for @xmath191 there always will be a region of stability for the ibcps , even for @xmath192 .
this stands in marked contrast to similar modelling for ternary systems where bicontinuous cubic phases are expected to disappear for @xmath193 @xcite .
the explanation lies in the fact that changing the spontaneous curvature @xmath0 of the monolayers amounts to changing the effective saddle - splay modulus of the bilayer , compare eq .
( [ eq : effectivebilayerbending ] ) .
thus , increasing @xmath0 can counteract the effect of making @xmath194 more negative .
the bicontinuous phases are also favored by lowering @xmath6 without changing @xmath91 , since this increases the fluctuation effects which favor ibcps ( positive renormalization of @xmath7 ) and disfavor the lamellar phase ( increase of steric repulsion ) .
this can be seen by comparing b and [ figurephasediagram]c , where @xmath195 is decreased from @xmath196 to @xmath197 , respectively , while @xmath91 is kept constant at @xmath190 .
it is clear from our discussion in sec .  [
sec : model ] that all contributions to the free energy , which can be written in the scaling form @xmath198 , will preserve the triple - line coexistence of g , d and p. obviously , an interaction which introduces an additional length scale , in addition to @xmath1 and @xmath89 , is needed to break the bonnet symmetry and to lift the triple - line coexistence . here
we want to discuss van der waals interactions and chain stretching .
the non - retarded van der waals interaction at distance @xmath199 has the general form @xmath200 where @xmath201 is the hamacker constant , and @xmath202 and @xmath203 are the number densities of atoms in the two interacting bodies .
the total interaction energy between two colloidal particles of volumes @xmath204 and @xmath205 is then given by @xmath206 unfortunately , the energy density of the van der waals interaction can not be calculated easily in a cubic bicontinuous phase . however , it is clear from eq .
( [ eq : vdw_energy ] ) that it is the embedding of the tpms in three - dimensional space , which is important for the van der waals energy .
distances between different parts of the bilayer are not determined uniquely by the @xmath4-distribution , which describes the internal geometry of a surface .
this is evidenced by the self - intersections of the parallel surfaces for @xmath207 , which do _ not _ affect any of the considerations of the curvature energy described in the previous sections , but would give a divergent contribution to the van der waals energy . thus , we conclude that the van der waals interaction must break the bonnet symmetry .
we illustrate this point in appendix  [ app : fcc ] , where we approximate the tpms by a fcc array of spherical shells .    in order to estimate the frustration energy due to chain stretching ,
it has been suggested in ref .
@xcite to calculate the variance of the distance between constant - mean - curvature companions of a given tpms .
for the d - surface , the variance @xmath208 has been shown @xcite to be almost independent of @xmath49 , with @xmath209 and @xmath210 , where eq .
( [ latticeconstant ] ) has been employed in the second equality .
similar results can be expected for g and p , with prefactors @xmath211 and @xmath212 , respectively .
we introduce a function @xmath213 , which is defined such that @xmath214 , @xmath215 and @xmath216 for the three values of @xmath43 for the p , d and g phases .
then , the stretching - energy density has the scaling form @xmath217 where @xmath218 is the stretching energy from and has been used for @xmath219 . note that this energy density has been made dimensionless with the same factor as the curvature energy in sec .
[ sec : model ] , i.e. it is normalized to the free - energy density of the lamellar phase [ compare eq .
( [ freeenergydensitylamellar ] ) ] .
these results make the scaling argument of sec .
[ sec : model ] more precise .
equation  ( [ deltal ] ) shows very clearly that the lattice constant @xmath23 appears as new length scale , and that the bonnet symmetry is broken in general by the stretching contributions .
it is interesting to consider the possible functional dependencies of @xmath213 in more detail .
first , if @xmath220 is a monotonically increasing ( decreasing ) function of @xmath43 , then p is stabilized ( destabilized ) by chain stretching due to its small value of @xmath43 , while d is disfavored ( favored ) . in the case of increasing @xmath220
, d would not appear in the phase diagram , while in the case of decreasing @xmath220 , p can not be present .
second , if @xmath221 const , the bonnet symmetry is restored .
finally , @xmath220 can of course be a non - monotonic function of @xmath43 , which stabilizes the phase with the smallest value of @xmath220 .    in any case
, the difference between the three phases can be expected to be small , since the prefactor @xmath222 ( and presumably also @xmath211 and @xmath212 ) is very small ; nevertheless , the stretching energy might account for small differences , like the ones needed to produce the phase diagram of 2:1 la / dlpc and water , which shows no triple - line symmetry , but something close to it .
a more precise statement requires the calculation of @xmath219 for g and p in the framework of the cmc - model .
other estimates of the stretching energy have been discussed in ref .
@xcite ; they also indicate a breaking of the bonnet symmetry .
in this work , we have investigated inverse bicontinuous cubic phases ( ibcps ) in lipid - water mixtures .
we systematically included all ibcps which might be expected to be of physical relevance , and found that only the ones with a narrow distribution of gaussian curvature over the minimal mid - surface are stable , i.e.  p , d and g. while this result has been anticipated in ref .
@xcite , we were able to prove it quantitatively for the first time by calculating the distributions of gaussian curvature over the minimal mid - surfaces of all relevant structures ; from these we calculated the variance @xmath166 of the @xmath4-distributions given in which is the same for p , d and g  due to bonnet transformations  and higher for all other structures .
the bonnet transformation also implies that the free - energy densities for p , d and g have the scaling form @xmath223 with hydrocarbon volume fraction @xmath49 and topology index @xmath47 .
the topology index is a dimensionless inverse euler characteristic which is independent of scaling and choice of unit cell ; the higher its value , the less porous the structure . since it is highest for the gyroid structure g , it has the largest region of stability .
d and p have progressively lower values for @xmath47 , and have therefore smaller regions of stability at lower values for @xmath49 .
we also used the concept of the emulsification failure ( which is well established in surfactant science ) to explain the stability of ibcps in excess water .
a maxwell construction with the excess phase shows that p , d , g and excess water lie on a triple line for all values of @xmath0 and @xmath91 as a consequence of the bonnet transformation .
this bonnet symmetry should not be interpreted as a degeneracy between g , d and p , since their different topologies ensures that their regions of stability are arranged according to their respective values of the topology index .
it should be noted that the free energy per _ unit area _ ( rather than per unit volume ) has the scaling form @xmath224 ; this might be interpreted as some kind of degeneracy ( compare refs .
@xcite ) , since now all three structures attain the same minimal values .
our results agree nicely with experimental observations . at fixed temperature , it is always the sequence g - d - p which is observed with increasing water content , with g having the largest region of stability and only d or p coexisting with excess water @xcite .
for example , for monoolein and water at @xmath225 , with increasing water concentration the sequence l - g - d has been found , where g has the largest region of stability and d coexist with excess water @xcite .
our phase diagram for @xmath226 and @xmath179 ( a ) agrees particularly well with the one for 2:1 la / dlpc and water @xcite ( compare ) .
we do not treat the inverse hexagonal phase , which experimentally occurs at large @xmath49 and large @xmath0 , but we showed above that ibcps will occur at large @xmath0 and _ small _ @xmath49 , and argued that at large @xmath49 stretching energy should become comparable to bending energy .
this , in fact , favors the inverse hexagonal phase , whose geometry requires chain stretching .    finally , we can discuss the role of temperature for 2:1 la / dlpc and water by using the ansatz @xmath227 , with a balanced temperature @xmath228 and a phenomenological constant @xmath229 . using @xmath230 , we find @xmath231 and @xmath232 . since the chain melting temperature is known to be @xmath233 , we conclude that the balanced temperature @xmath228 is well below the main transition .
this seems to be the generic case for lipid systems ; it might be related , e.g. , to the anomalous swelling of pc - bilayers when the main transition is approached from above @xcite .
it is interesting to note that the constant @xmath229 obtained here for 2:1 la / dlpc and water is very close to the corresponding value @xmath234 for monolayer of @xmath235 in ternary microemulsions @xcite .    in this paper , we only considered inverse bicontinuous phases with cubic symmetry
. a similar analysis for non - cubic phases should be straight - forward , since here weierstrass representations are known for the main structures , especially for the rhombohedral and tetragonal variants of g , d and p @xcite .
however , we expect that these phases will not be more favorable than g , d and p for the same reason as for the other cubic phases . in order to make further progress in understanding their phase behavior
, it seems more important to investigate the different mechanisms which might break the bonnet symmetry between g , d and p. in our view , two calculations are needed now : ( i ) of the van der waals energy of g , d and p and ( ii ) of the variance of the distance distribution of a constant - mean - curvature model of g and p.
for a rough estimate of the energy density due to the van der waals interaction in cubic bicontinuous phase , we use a fcc array of spherical shells of radius @xmath236 , thickness @xmath237 and lattice constant @xmath23 . here
, the lattice constant @xmath23 corresponds to the conventional unit cell and is taken to be the same as for the corresponding bicontinuous phase , and @xmath236 is determined by the requirement of equal area @xmath114 per unit cell , so that @xmath238 .
thus we have @xmath51 again .
the distance of closest approach between two neighboring spheres is @xmath239 ; in order to avoid contact between them , @xmath240 is necessary , which implies @xmath241 ; this inequality is satisfied for all tpms studied here , compare tab .
[ tableindex ] . for the crystal ,
we consider the interaction between nearest - neighbor spheres only .
the van der waals energy between two hollow spheres is described by the girifalco potential , which is often used to model the van der waals interaction between two buckyballs @xmath242 @xcite . for our purpose
it is sufficient to use the derjaguin approximation for the girifalco potential , which describes the leading contributions for close approach correctly , and shows that the van der waals attraction for two hollow spheres scales as @xmath243 .
we thus find for the van der waals energy density @xmath244 this energy density has again been made dimensionless with the same factor as the curvature energy in sec .
[ sec : model ] .
its scaling shows that _ two _ ratios of length scales , @xmath245 and @xmath246 , determine the van der waals energy .
it is intersting to note that the curvature - energy density in this model of spheres scales as @xmath247 ^ 2\ .\ ] ] since @xmath248 ( for fixed @xmath44 ) , is consistent with the scaling function ( [ eq : scaling ] ) .
this indicates that our identification of the sphere radius and lattice constant gives a consistent approximation for both van der waals and bending energies .                                  ,  e. ;  rummel ,  g. ;  rosenbusch ,  j.  p. ;  ,  e.  m. _ science _ * 1997 , * _ 277 , _ 16761681 .
rummel ,  g. ;  hardmeyer ,  a. ;  widmer ,  c. ;  chiu ,  m.  l. ;  nollert ,  p. ;   locher ,  k.  p. ;  pedruzzi ,  i. ;  landau ,  e.  m. ;  rosenbusch ,  j.  p. _ j. struct .
biol . _ * 1998 , * _ 121 , _ 8291 .
templer ,  r.  h. ;  seddon ,  j.  m. ;  warrender ,  n.  a. ;  syrykh ,  a. ;   huang ,  z. ;  winter ,  r. ;  erbes ,  j. _ j. phys . chem .
b _ * 1998 , * _ 102 , _ 72517261 .
templer ,  r.  h. ;  seddon ,  j.  m. ;  duesing ,  p.  m. ;  winter ,  r. ;   erbes ,  j. _ j. phys .
b _ * 1998 , * _ 102 , _ 72627271 .                                                                              .euler characteristic @xmath44 , scaled surface area @xmath43 ( both in the conventional unit cell ) and topology index @xmath249 for all tpms considered , as obtained from local minima of the free - energy functional of a simple ginzburg - landau theory .
exact values are known for g , d , i - wp , p and c(p ) and given in .
the structures are ordered according to decreasing @xmath47 ; the value for s is probably too large , but difficult to improve numerically . [ cols="<,<,<,<,<,<,<,<",options="header " , ]