typical proton - conducting polymer electrolyte membranes ( pem ) for fuel cell applications consist of a perfluorinated polymeric backbone and side chains that are terminated by strongly acidic sulfonic head groups .
they provide the transport medium for the protons generated during the anodic oxidation reaction of the fuel .
pem fuel cells offer numerous advantages as power generators for emission - free vehicles , portable applications , and individual homes @xcite .
when such ionomers are exposed to water or humid air , the acidic groups terminating the side chains dissociate into so@xmath1 groups on the chains , and protons in the aqueous sub - phase .
the sub - phase extends through the membrane , creating a bi - continuous nano - phase separated network of aqueous pores and polymer .
the percolation structure and its ability to support proton diffusion have to meet strict requirements for a pem to be considered as an ` ideal ' membrane for fuel cells : a ) it must form a chemically , mechanically , thermally and morphologically stable structure over a wide range of water content , b ) show steady performance and tolerance of elevated temperatures up to 150@xmath2c at high proton conductance , and c ) be impermeable to gases and contaminating ions .
development of new pem membranes that could meet all these demands is a primary task for many research groups involved in fuel - cell studies .
one extensively used pem material is nafion@xmath0 , a perfluorinated polymer from dupont .
the integrity and structural stability of nafion membranes is provided by the polytetrafluoroethylene ( ptfe ) backbones , but it is the hydrophilic clusters of side chain material that facilitate the transport of ions in the membrane .
the microstructure of these clusters consists of an inter - facial region of solvated perfluoroether side chains separating the polymer matrix from the more bulk - like water found in the pores . in the presence of water ,
one finds a nano - separation of hydrophobic and hydrophilic domains . according to the early work of gierke _ et al .
_ @xcite these pores are inverted micellar spheres with diameters around 4 nm , and are interconnected by channels with a diameter and length of about 1 nm .    recent experimental investigations suggest that the spherical shape and uniform spacing of these pores are serious oversimplifications @xcite .
new small - angle x - ray scattering ( saxs ) and small - angle neutron scattering ( sans ) experiments by gebel _
@xcite for nafion membranes in different swelling states have given a different picture for the morphology .
their results suggest that elongated bundles of hydrophobic backbone material with diameters of about 4 nm and lengths greater than 100 nm are surrounded by sulfonate groups and water molecules . in other words ,
hydrophobic material is surrounded by hydrophilic domains .
the fact that controversy continues to surround the nature of the hydrated morphology and the shape of the domains suggests that there is still room for theoretical investigation and computer simulations to address this issue . _
ab initio _ electronic structure calculations of polymeric fragments in contact with water @xcite and quantum molecular dynamics studies on model pem systems @xcite have provided a basis for understanding the energetically favorable conformations of nafion chains and the molecular ingredients of the conduction process @xcite . for systems large enough to contain several examples of hydrophilic inclusions and their interconnecting networks ,
however , treatment of the polymer in an _ ab initio _ manner , i.e. a full electronic treatment with molecular orbital theory , is not computationally feasible .
the objective of this paper is to develop a deeper understanding of the driving mechanism behind the nanophase aggregation of sulfonate head groups , and to test the extent to which coarse - graining methods can be applied .
we attempt this by using a simple dipolar head group model for ionomer side chains .
we show that the dipole - dipole correlations strongly depend on the molar fraction of acidic groups and dielectric permittivity of the membrane .
globular micelles formed in the case of zero dipole moment evolve to percolated cluster structures for a non - zero dipole moment of head groups when an artificial attraction between the end groups of side chains is introduced .
accounting for the dipole moment of sulfonate head groups thus changes the conformational structure of clusters from a gierke model of spherical aggregations to a gebel model of elongated aggregations . the determination of the parameters necessary for success in modeling nanophase aggregation will help us understand the essential role played by the electrostatic interactions of the head group in determining microstructure formation in pem materials .    the rest of this paper is organized as follows . after a brief review of segregation models in section
ii , we present our dipolar model for side chains in section iii . section iv is devoted to our simulation results for dry pem ionomers .
we conclude with a discussion of results and planned future work in section v.
significant experimental and theoretical efforts have been put into the characterization of the microstructure of dry and swollen nafion membranes . although the fact of microphase segregation and its influence on the thermodynamic and transport properties of nafion membranes are well established , some uncertainty regarding the morphology of the subphases in swollen nafion remains .
it has been generally assumed that the hydrophilic phase is continuous , in accord with the experimental evidence of high conductivity and water permeability .
the conventional model of microphase segregation in nafion membranes , put forward by gierke _
@xcite , where the protons and water diffuse between globular clusters through cylindrical narrow channels , was further modified by mauritz and rogers @xcite and eikerling _ et al .
also , ioselevich _ et al . _
@xcite developed a lattice based micelle - channel model to describe a connecting structure of a water bridges between the cages of the micelles .
however , there is as yet no direct experimental evidence for the existence of the channels connecting the clusters in gierke s model .
other models explaining the properties of hydrated nafion exploit bilayer , lamellar , and sandwich - like structures for nanophase separation in pem membranes .
gebel @xcite describes the nafion membrane as an aggregation of polymeric chains forming elongated objects ( simplified as cylinders ) , embedded in a continuous ionic medium . at larger scales ,
those aggregates form bundles with definite orientational order .
this new picture of multi - scale structure can explain membrane swelling as a continuous process from the dry state to a colloidal suspension .
scanning calorimetry experiments on cast nafion films @xcite led to the conclusion that a continuous water phase exists only when the relative humidity of the membrane exceeds 90% .
this conclusion contradicts the simulation results of khalatur _ et al . _
@xcite , where formation of continuous channels at very low hydration , and even in almost dry membranes , is reported .
this is certainly consistent with the conductivity measurements in nafion films with very little water present @xcite .
in contrast to all the above - mentioned network models , vishnyakov _ et al . _
@xcite showed that the water clusters do not form a continuous hydrophilic subphase .
instead , the cluster - size distribution evolves in time as a consequence of the formation and rupture of temporary bridges between clusters .
molecular dynamics ( md ) studies , using both classical force fields and an _ ab initio _ approach , appear to support the idea of irregularly shaped ionic clusters of sulfonic groups connected by smaller channels of similar composition embedded in the hydrophobic covalent matrix .
however , the characteristic length scales of the morphological conformations suggested by experimental data are currently inaccessible to atomistic modeling . to probe the morphology at the nanometer scale ,
a mesoscale model is required .
in this work we are concentrating on the effects of electrostatic interactions among the elements of the side chains , and so in order to avoid the excessively slow processes resulting from chain entanglement we consider a segmented system in which the side chains are detached from the backbone .
the later introduction of very strong interactions between the ends of the side chains and the backbone , which have the effect of knitting together the segments into a single macromolecule , do not seem to affect our results significantly .    a generic acidic side chain of nafion is commonly denoted as cf@xmath3ocf@xmath4cf(cf@xmath3)o(cf@xmath4)@xmath4so@xmath3h . in order to expedite these simulations ,
we adopt a simplified model for the side chain by using a united - atom presentation of the cf@xmath4 and cf@xmath3 groups @xcite and the so@xmath1 sulfonate head groups @xcite .
these groups are modeled as lennard - jones ( lj ) monomers .
we employ the same united - atom approach for backbone chain segments .    in the framework of our dipole model
all the partial charges on the ether oxygen , carbon and fluorine atoms of the side chain are set to zero .
the fluorocarbon groups of the backbone skeleton are also assigned zero partial charges .
the electrostatic charge is located entirely on the sulfonate head group , so@xmath1 , and on the hydrogen ion h@xmath5 permanently bound to it .
the dipole moment produced when the proton and the sulfonate group , modeled as lj monomers , carry the full formal charges of @xmath6 and @xmath7 respectively , was reduced by a factor @xmath8 in some of our simulations . though our model is not applicable for the measurement of proton diffusion in pem membranes , it is a good starting point for a step - by - step exploration of the nanophase morphology in dry pem materials .
the overall inter - monomer interaction in our model system is a combination of electrostatic coulomb interactions between each pair of charged particles , 12 - 6 lennard - jones ( lj ) potentials between neighboring monomers , and intrachain stretch and bend energies .
the lennard - jones interaction parameters for the side chain are chosen to agree in most instances with the nafion model of paddison @xcite .
we make a further simplification for the lj parameters of the side chain by assigning a universal value of @xmath9 , chosen to be 3.5    for all monomers . the total potential energy of the system is then given as @xmath10 where @xmath11 and @xmath12 are the lennard - jones , electrostatic , bond - stretching ( two - body term ) , and angle - bending ( three - body term ) components of the total energy , and are given by @xmath13 @xmath14 @xmath15
@xmath16 we omit consideration of the dihedral potential term ( four - body term ) , as this would be expected to make only a small contribution to the energetics of the single chain conformations in these dense polymeric systems .
( this assumption was tested in sample simulations in which this term was restored , and found to have little effect on our results . )
the following parameters for the polymer chains were used : an equilibrium bending angle @xmath17 , an equilibrium bond length @xmath18 ( 0.44@xmath9 ) , a bending force constant @xmath19 kcal / mol deg@xmath20 , and a stretching force constant @xmath21 kcal / mol @xmath22 .
the last two constants correspond to the thermal energy @xmath23 at room temperature for a bond - length displacement of 0.1  and bond - angle displacement of @xmath24 from the equilibrium position .
the bond length between the so@xmath1 head - group monomer and the hydrogen ion was chosen to be @xmath9 , and hence 0.35 nm .
the lj interaction coefficient @xmath25 was chosen to be 3@xmath23 for hydrophobic - hydrophobic ( hh ) interactions , and 1@xmath23 for both hydrophobic - hydrophilic ( hp ) and hydrophilic - hydrophilic ( pp ) interactions . in the latter case
the lj potential was modified to be purely repulsive by truncating at its minimum , where @xmath26 , and raising it by the addition of an amount @xmath25 .
the reason for this choice is that the interaction between protons is not well represented by van der waals forces , which arise from virtual dipole - dipole interactions , and also that such attractive forces are small in comparison to the coulomb interaction .
the normalized charge @xmath27 ( @xmath28 ) in eq.(3 ) is + 1 for protons and -1 for sulfonate groups .
we used this coarse - grained model to perform molecular dynamics simulations of ensembles of ionomer segments .
because chain entanglement greatly slows the rate of equilibration in polymers , we adopted the strategy of temporarily cutting the system into more mobile segments , which were later reassembled . in order to increase the diffusional rate of our model ionomer we detached the side chains from the backbone skeleton
, an approach previously adopted by vishnyakov @xcite , and also cut the backbone skeleton into segments @xcite . in total the system then consists of @xmath29 sidechain segments and @xmath30 backbone segments .
taking into account that the sulfonic acid groups are hydrophilic while the ethers and fluorocarbon groups are hydrophobic @xcite , we adopted the following general representation for a coarse - grained architecture of side chains and backbone segments .
each side chain contains @xmath31 monomers of which @xmath32 are hydrophilic and @xmath33 are hydrophobic .
its architecture is written as : @xmath33h+@xmath32p .
the backbone segments are fully hydrophobic with an architecture @xmath34h . in the majority of systems studied in this paper ,
the side chain architecture was 8h+2p and the backbone segment architecture was 10h .
the simulations were then performed in several stages .
the initial configuration was first prepared by using monte - carlo techniques to grow a mixture of polymer chain segments with fixed bending angle @xmath35 and bond length @xmath36 . after an initial equilibration
the constraints on @xmath35 and @xmath36 were removed .
subsequently , md runs were performed for time periods between 50ps and 500ps at constant volume @xmath37 and constant temperature @xmath38 until the system of polymer segments was fully equilibrated .
the system temperature was controlled by coupling the system to a langevin thermostat with a friction coefficient @xmath39 and gaussian white - noise force of strength @xmath40 .
the equations of motion were integrated using the velocity verlet algorithm with a time step of 0.5 fs .
we imposed standard periodic boundary conditions to our system , filling space with translational replications of a fundamental cell .
long - range electrostatic interactions were calculated using the standard lekner summation algorithm @xcite .    to check
whether equilibrium had been reached we examined the proton pair distribution function @xmath41 in every picosecond time interval , and continued the run until no significant further change was observed . to avoid
mistaking the trapping of the simulated system into a metastable glassy state for true equilibrium , we repeated each run with several different initial configurations . for each computer experiment
we typically gathered system statistics over time intervals of a few nanoseconds .
we believe this to be sufficient to allow for the self - organization of the side chains and backbones into clustered structures .    in the next stage of the simulations
the segments were reassembled into the branched chain characterizing the original nafion .
this was achieved by the simultaneous introduction of a specific mutually attractive force between the ends of each backbone segment , which united them into a single chain , and a similar force acting between the tail monomer of the detached side chains and the fifth monomer of every backbone segment .
this was applied in the form of a lj potential of strength @xmath42 . to avoid the formation of star - like branched polymers
, only a single occupancy of the backbone attachment sites was permitted .
the simulation was then resumed , and run until a new equilibrium was attained .
it was found that the bonds on the ends of all the backbone segments were occupied , and that a large majority ( typically around 95% ) of the side chains were reattached to a backbone segment .
this fraction was a function of the ratio @xmath43 of the number @xmath29 of side chains ( and hence of sulfonate groups ) to the number @xmath44 of backbone monomers .
the quantity @xmath43 thus measures the concentration of sulfonate groups relative to the concentration of backbone monomers .
the few side - chains that did not bond to the special sites on the backbone attached themselves more loosely to the backbone through the weaker regular hydrophobic - hydrophobic lj interaction potential .
we also introduce the absolute molar concentration of sulfonate groups , @xmath45 , with @xmath46 being avogadro s number and @xmath47 the volume @xmath48 of the simulation cell , as an alternative parameter with which to describe the sulfonate concentration , independently of the amount of backbone polymer . while it is true that in real dry nafion , no significant free volume exists , and so the total number of monomers per unit volume is roughly fixed at the time of synthesis , in a simulation , however , it is possible to change the density of the material .
it is of interest in understanding the mobility of ions in nafion to know what structures would form in the presence of imposed free volume , were such a reduction in density of the material possible .    in order to observe the effects of the dipolar nature of the sulfonate - proton group on the structural characteristics of side chain aggregates
, we calculated the normalized pair distribution function @xmath41 of head group protons ( which are the end group monomers of the side chains ) .
the function @xmath41 indicates the relative probability of finding two protons at a separation distance @xmath49 averaged over the equilibrium trajectory of the simulated system , so that @xmath50 here @xmath51 is the number of protons located in a shell of thickness @xmath52 at the distance @xmath49 from a fixed proton , and @xmath53 is the total number of protons in the system of volume @xmath54 .
in addition to this spherically averaged quantity , it is also useful to note some indications of the local anisotropy in order to distinguish between the gierke and gebel models . to this end
we also evaluated the partial structure factors @xmath55^{2}+\left[\sum_{i=1}^{n_{+}}\sin q_{j}r_{i}\right]^{2}\right\rangle .\ ] ] here @xmath56 denotes the @xmath57 , @xmath58 and @xmath59 components of the vector @xmath60 , with @xmath61 , and the sums can proceed over the entire sample or over smaller volumes .
.list of key variables . [ cols="<,^ " , ]
the goal of these extensive coarse - grained molecular - dynamics simulations was to investigate the relative importance of the various model parameters on the aggregation of the sulfonate ( acidic ) head groups . to achieve this we performed a series of simulation runs in which the key parameters were varied as summarized in table ii .
these parameters included the absolute concentration and the relative concentrations of side chains and backbone polymer , as well as the strength of the electrostatic interactions .
this last quantity is characterized by a coulomb coupling parameter defined as @xmath62 .
here @xmath63 , which is known as the bjerrum length , is about 56  nm for a dielectric constant @xmath64 and @xmath65 k. the coupling parameter @xmath66 was varied between 2 and 160 by varying the system dielectric permeability @xmath67 between 80 and 1 , with higher @xmath66 representing stronger electrostatic correlations in the system . in a typical simulation , a few thousand side chains
are confined in a cubic box whose side is of length @xmath68 , with @xmath69 nm being the monomer diameter .
each simulation run was started from a randomly generated assembly of side chains and backbone segments within the system volume . however , in order to avoid packing problems , and to make this process effective for high monomer densities , a preliminary simulation of the same number of monomers was first performed in a larger box of side @xmath70 , and this was then followed by a gradual reduction of the system volume until the desired molar concentration was achieved .      in this first simulation
the constraints imposed by the attachment of the side chains to the backbone were removed in order to concentrate on the effects of the electrostatic interactions .
we then analyzed the proton - proton pair distribution function @xmath41 calculated for an assembly of unattached side chains for three different dielectric constants @xmath67 of the ionomer medium , and for two different molar concentrations .
the higher concentration , for which @xmath71
mol / l , is characteristic of typical dry nafion samples .
figure [ fig4 ] shows the proton - proton pair distribution function for the parameters described in table ii as run  1 and run  7 .
the first observation is that , as expected , small @xmath67 , or , equivalently , high @xmath66 , facilitates clustering of the head groups , as shown in the red dot - dashed lines in figure [ fig4 ] .
the graph of @xmath41 shows well - defined maxima at @xmath72,2.0 , and 3.4 .
a detailed investigation shows that this is a consequence of the formation of three simple multiplet structures : compact and linear multiplets , each created by two head groups , and a branched multiplet composed of three head groups , as shown in fig .
[ clusters ] .    for @xmath73
there is no correlation between head groups on different side chains except for steric hindrances , as seen in the continuous blue lines in figure [ fig4 ] .
a more surprising result is the effective repulsion between acidic head groups at intermediate dielectric permittivities , as seen in the green dashed lines in fig .
it appears that , in the absence of strong correlations , each head - group dipole prefers to stay in the bulk of the sample in order to increase its polarization energy .
this is similar to the well known salt - depletion effect in colloidal physics , when salt pairs are effectively repelled from the surface of a neutral colloid @xcite .
the other unexpected result is the suppression of the correlations between the end groups at high molar sulfonate concentrations @xmath74 , as seen in the dot - dashed line and symbols in figure [ fig4](b ) .
this can be traced back either to the reduction of the debye screening length , which effectively shields the dipoles from each other , or to the enhancement of the hydrophobic - hydrophobic attraction in the system .
the latter emerges from the hydrophobic tails of the side chain segments .    to decide which of these two possible mechanisms was responsible , we added into the system described in fig .
[ fig4](b ) a small number of backbone segments in order to boost the effective hydrophobicity ( or association energy ) in our model system .
one backbone segment of 10 monomers was added for each side chain unit , making @xmath75 .
after equilibration , the side chains were then loosely attached to the backbone segments and the backbone segments were mostly united in the manner previously described , whereupon the system was further equilibrated .
the results , shown as run  2 in figure [ fig6 ] , clearly indicate that the backbone segments have a minor impact on the head group distribution function .
the system has almost the same microphase structure as observed for the original side chain system .
a fivefold increase in the backbone concentration in figure [ fig6](b ) has only a small effect on @xmath41 .
therefore , it is the dipole - dipole interaction that tends to destroy the correlations when the molar fraction of sulfonates is high .
this effect takes place when the screening length becomes less than the mean separation distance between the dipoles .
this has some relevance to the search for ionomer membranes with superior transport properties : whereas for better proton conductance one needs a higher sulfonate molar fraction , the latter hinders sulfonate aggregation , which is vital for water - sulfur channel formation in wet ionomers .    typical snapshots of the simulation box , shown in figure [ fig5 ] , make visible the process of microphase separation in ionomer systems .
the tips of the side chains , which correspond to end group protons , are shown as balls , whereas the rest of the side chains are plotted as lines .
the backbone segments , drawn also as lines , mostly cluster together to create large aggregates .
the side chains either anchor onto the surface of the backbone aggregate , or create their own small clusters .
figure [ fig3d - density ] presents a 3-dimensional density distribution of monomers for run  5 .
the areas of low ( high ) density colored in blue ( red ) show the channel structure of the aggregates .
the polymer matrix is composed of hydrophobic tails of side chains and backbone segments .
the hydrophilic head groups form the walls of these channels .
though this 3d density picture differs from the long - bundle structural model of nafion @xcite , it supports the model of a dynamic system of non - spherical clusters with temporary inter - cluster bridges @xcite .      in the previous section , in which we varied the strength of the coulomb interactions
, we observed that in the systems with the higher sulfonate concentration , the coulomb correlations were suppressed .
this suggested that a study of the morphology of aggregates over a wider range of concentrations would be helpful to understand the physics behind sulfonate cluster formation .
simulation results for @xmath41 for systems with sulfonate concentrations @xmath74 varying between 1.6 and 2.5
mol / l are shown in fig .
[ fig9 ] for runs 8 and 10 - 12 .
results for the much wider range of 0.16 to 4.8
mol / l are shown in fig .
[ fig9a ] , but for these it was necessary to delete the backbone segments entirely in order to accommodate the very large concentration of side chains at the upper part of the range constituting run  13 .
the dipole - dipole correlations are seen to depend strongly on the molar concentration of acidic groups . for a typical dilute case ,
run  10 in fig .
[ fig9 ] , the first peak of @xmath41 reaches a value @xmath76 , indicating a strong condensation of nearly all the side chains .
a possible explanation for this effect can be drawn from the clustering features of pure dipolar liquids . in a dilute system , where the average interchain distance is larger than the gyration radius @xmath77 of the side chain , the lj interaction between the neighboring chains is negligible .
the main contribution to the energy in this case comes from the electrostatic dipole - dipole interaction of head groups .
although this interaction decays with a @xmath78 distance dependence , it is not weak over distances comparable with the simulation box length @xmath79 , and is capable of forcing the chains to aggregate .
however , for a dense system , where the lj contribution to the energy becomes comparable to the dipolar term , the side chains start to lose their rotational entropy , and thus can not reorient themselves in order to minimize the electrostatic energy of the aggregate . while the heights of the maxima of @xmath41 gradually decrease as
the system becomes more tightly packed , their positions do not change and remain fixed at @xmath80 and@xmath81 .
the line with circles in fig .
[ fig9a ] corresponds to the case where the concentration of backbone segments was double the case for run  8 , and corresponds to run  9 .
it is seen that the proton pair distribution function is hardly affected by the existence of the extra hydrophobic segments .
however this conclusion is not necessarily true for a system with the more attractive lj parameters for hh and hp interactions that will be considered in section d. both the snapshot picture in fig .
[ fig15 ] and the 3d monomer density in fig .
[ fig8 ] show microphase separation .
another way to examine the head group clustering is to calculate the partial structure factors @xmath82 ( @xmath83 , @xmath58 , @xmath59 ) of the protons .
because these did not exhibit any marked anisotropy , an averaged @xmath84 is plotted in fig .
[ figsk ] .
the position of the peak at @xmath85 is associated with a correlation distance , @xmath86 in the sample structure via the bragg relationship @xmath87 . a sharp growth at small @xmath88 of the values of @xmath84 in fig .  [ figsk ] results from microscopic aggregation of hydrophilic and hydrophobic monomers in the system .
to summarize the tendencies of dipole - dipole correlations presented in this and the previous section , we conclude that the hydrophilic head group aggregation takes place preferentially at low dielectric permittivity @xmath67 and low densities of sulfonate groups .
although the bulk of this work describes the conformational structures of side chains of length @xmath90=10 , it is worth while investigating systems with shorter side chains , since recent experiments @xcite show a high proton transfer rate in ionomer membranes with shorter pendant segments .
this section contains simulation results for different side chain lengths , @xmath91 .
the length of the hydrophilic part ( the number of hydrophilic monomers ) is fixed at @xmath92 and the length of the hydrophobic tail ( the number of hydrophobic monomers ) @xmath33 is varied between 0 and 8 ( which corresponds to the variation of total side chain length @xmath90 between 2 and 10 ) in order to reveal how the side chains hydrophobic tail affects the pair distribution function of protons .
no backbone material was included in these systems .
simulation results for run  14 are shown in fig .
[ fig12 ] .
note that the case @xmath93 corresponds to a system of dipolar dumbbells .
even a short hydrophobic tail with @xmath94 strongly suppresses the aggregation of head groups .
thus the hydrophobicity of side chains , or the length of hydrophobic tail , is a strong obstacle to head group clustering . in other words ,
the dipole - dipole attraction between sulfonates of longer side chains is not sufficient to produce a nanophase separation .
this conclusion is in contrast to the suggestion of eisenberg @xcite that each ionic group in the vicinity of a cluster of sulfonate groups is drawn in by dipolar attraction to form a depleted zone free of charge .
one intriguing question in fuel cell technology is how the morphology of nanophase separation in pem membranes depends on the polymeric structure of nafion . in other words ,
what is the optimal comb - like structure for the ionomer to be a good material for a fuel cell ?
using our segmental model for nafion we can mimic different comb - like nafion molecular conformations by changing the parameter @xmath43 , which denotes the ratio of the number of sulfonate groups to the number of backbone monomers .
the impact of @xmath43 on the distribution function @xmath41 was shown to be negligible for @xmath95 for hh and @xmath96 for both hp and pp interactions in fig .
[ fig6 ] for runs 1 - 6 . in order to explore more generally the effect of changing @xmath43
, we exaggerate the hydrophilic and hydrophobic interactions by tripling the interaction constants @xmath25 . the results for runs 1 - 3 presented in fig .
[ fig16 ] for @xmath97 for hh , and @xmath95 both for hp and pp interactions , now show that the addition of backbone chain segments brings more conformational disorder into the head group correlations .
this is also evident from the system snapshots for run  3 in fig .
[ fig18 ] for the two different strengths of lennard - jones interaction parameters . in other words , a strong interaction between hydrophobic and hydrophilic monomers
can disrupt the side chain aggregation , causing all the side chains to be evenly distributed on the surface of backbone aggregations , as in fig .
[ fig18](b ) .
this distribution will have a globular , lamellar , or cylindrical - rod form .      to complete this study of the effects of interactions on morphology
, we look at the effect of both decreasing and increasing the interaction between head groups . we have already seen in section iiia some of the consequences of reducing the coulomb interaction by increasing the dielectric constant @xmath67 . in the absence of interactions other than the dipole - dipole interaction ,
this is equivalent to a reduction in the strength of the dipole moments . in fig .
[ fig1 ] we show the effect on the structure factor @xmath84 of the protons of a reduction in the strength of the dipole moments by various factors @xmath8 ( run  8) .
when the parameter @xmath8 is small , there is no phase separation in the ionomers , and the function @xmath84 has a liquid - like structure with a single maximum .
the behavior of @xmath84 is changed for an unreduced dipole strength @xmath8=1 in a manner similar to the behaviour of @xmath84 in the less dense system shown in fig .
[ figsk ] : there is a strong peak at @xmath98=4.2 , and a sharp growth at small @xmath98 .
the latter , again , is an indication of a strong microscopic separation of hydrophilic and hydrophobic monomers in the system .
it is clear that the effective attraction between the acidic groups of side chains can not be accounted for by the sulfonate - hydrogen dipole moment alone .
the sulfonate end group has its own internal dipole moment due to the charge distribution between the sulfur atom and the three oxygen atoms bound to it .
partial charges on the other side chain monomers can also contribute to the effective attraction between head groups . to take all these contributions into account at a phenomenological level we introduce an artificial attractive interaction between the end groups of the side chains , in an approach similar to that employed by shirvanyanz @xcite .
for this purpose we use a yukawa - type attractive potential , @xmath99 here @xmath100 is an effective yukawa charge for the head group sulfonates . for the sake of simplicity ,
we restrict ourselves to the case of constant debye screening length @xmath101 .
in fact , the parameters @xmath100 and @xmath102 determine the stability , size and shape of the aggregates forming in the system .
the simple attractive yukawa term can serve as a convenient basis for interpreting all the important features that can be observed in real systems or predicted by statistical mechanical theories .    in the absence of electrostatic interactions in the system ,
i.e. when the head group dipole moment vanishes and @xmath103 , the ionomer system becomes spatially inhomogeneous
. the pair distribution functions @xmath41 for runs 7 - 9 , plotted in fig .  [ fig19 ] for @xmath104 ( case ( a ) ) and @xmath105 ( case ( b ) ) show strong structuring .
the aggregation is substantially stronger for a system missing backbone polymer ( run  7 ) than for a system with a full complement of backbone monomers ( run  9 ) , as seen from a comparison of the respective lines in fig .
[ fig19 ] .
there is a clear tendency to micellar arrangement , which is characterized by some short - range order .
snapshot pictures , given in fig .
[ fig23 ] , show formation of a space - filling web - like network of cluster - like aggregates surrounded by neutral polymer sections .
the size of these clusters grows as the density of side chains increases .
a comparison of figs .
[ fig23](a ) and [ fig23](b ) reveals that the effective yukawa charge @xmath100 , or in other words , the strength of the the screened coulomb attraction , has a strong impact on the shape of the clusters .
the structure of the system for a moderate yukawa attraction of @xmath106 between the end groups is mobile : the aggregates can diffuse , decay and reappear again .
they can also create temporary bridge - like connections .
therefore the nanophase separation of the head groups , shown in fig .
[ fig23](a ) can be viewed as an example for a gierke - like sulfonate aggregation . however , a strong attraction between the end groups of the side chains , the case @xmath105 in fig .
[ fig23](b ) , results in a frozen network of globular clusters of diameter roughly 1.5 nm in the polymer matrix .
this cluster picture resembles the ion segregation phenomenon proposed by eisenberg @xcite .
when the dipole moment of the sulfonate head groups is restored , in this case with a value of @xmath107 , the globular micelles form super - clusters of diameter 3 - 4 nm , as seen is fig .
[ fig24](b ) for @xmath105 . on the other hand
, this nonzero dipole moment facilitates the formation of percolated cluster structures for @xmath104 , as in fig .
[ fig24](a ) .
the inclusion of the dipole moment of the sulfonate head groups thus changes the conformational structure of clusters from a gierke model of spherical aggregations to a gebel model of elongated aggregations . a snapshot picture of the same system but with no artificial attraction between the head groups is given in the fig .
[ fig20 ] .
it is clear that the introduction of a screened coulomb attraction results in the formation of compact and interconnected clusters of sulfonate head  groups .
in this study we have explored the issue of nanophase separation in dry nafion@xmath0 films by means of simulations within a simple dipolar - head - group model for the side chains .
our procedure involved detaching the side chains from the backbone and cutting the backbone into segments to speed equilibration , and then reassembling the macromolecule by means of a strong imposed attractive force between the cut ends of the backbone , and between the non - ionic ends of the side chains and the midpoints of the backbone segments .
as we vary the interaction parameters and dielectric permittivity of the medium we find different structures to be formed by the head groups .
small permittivities facilitate the clustering of head groups with a well defined maximum in the proton - proton pair distribution function @xmath41 that indicates three simple multiplet structures : compact and linear multiplets each created by two head groups , and a branched cluster composed of three head groups . at higher permittivity no correlations between head groups are obtained . an effective repulsion between acidic head groups at intermediate dielectric permittivities
was revealed from our data .
we assume that , in the absence of strong correlations , each head - group dipole prefers to stay in the bulk of the sample in order to increase its polarization energy .    for the ionomer interaction parameters considered in this work the microphase structure of the head groups
was found to have little dependence on the concentration of backbone segments . since the molar fraction of sulfonates in a fuel cell membrane is known to be a crucial parameter for its performance , this factor was also studied .
we find that for a dilute case a strong condensation of nearly all side chains takes place .
the main contribution to the hamiltonian in this case comes from a long - range dipole - dipole interaction of head groups which forces the chains to aggregate .
however , for a dense system , where the hydrophobic lj contribution to the hamiltonian becomes comparable with the dipolar term , the side chains can not reorient themselves in order to minimize the electrostatic energy of the aggregate .
while on the one hand , for high proton conductance one needs higher sulfonate molar fractions @xmath74 , on the other hand we find that a high value of this parameter @xmath74 hinders the sulfonate aggregation that is vital for water - sulfonate channel formation in ionomers .
overall , our simulations show that the dipole - dipole interaction is the main driving force for the aggregation of sulfonates .
when the dipole moment of sulfonates is small , so that @xmath108 , there is no phase separation in the ionomers .
in contrast , for @xmath109 a microscopic aggregation of hydrophilic and hydrophobic monomers takes place in the system .
the most significant feature learned from the simulations is the fact that the dipolar - head - group approximation is not sufficient to drive the compact side chain cluster formation suggested by experimental measurements .
an additional artificial attraction between the sulfonates , introduced as a screened yukawa interaction potential , boosts formation of nanometer - sized clusters dispersed throughout the system volume .
the final size of these clusters , and the geometry of the cluster network they create , depend on the dipole moment @xmath8 of the sulfonate head groups for fixed yukawa parameters . as the strength of @xmath8 increases from zero to one , the conformational structure of clusters changes from a gierke model of spherical aggregations to a gebel model of elongated aggregations .
more work is required to develop a deeper understanding of the mechanisms behind the aggregation picture observed experimentally .
because the shape and size distribution of the aqueous pores in wet membranes is not completely known , it will be useful to extend our preliminary results for the dry nafion considered in this paper to systems with water present and with a more detailed chain structure . in future work
, we plan to extend the model presented here to more complex models .
in particular we will employ partial charges and lennard - jones parameters @xcite for side chain and backbone monomers of nafion - like oligomers .
the issues of how the partial charges of ionomer molecules , in fixed and free proton models for side chains , affect the nanophase structure of the pem materials will be addressed .
a complete understanding of the role of head group charges on the cluster conformations of side chains can provide helpful guidelines for understanding existing membranes and designing new membrane materials for fuel cells with improved characteristics .
this work was supported by the us department of energy under grant de - fg02 - 05er46244 , and was made possible by use of facilities at the case its high performance computing cluster and the ohio supercomputing center .
98 m. i. perry , t. f. fuller , j. electrochemical society * 149 * , s59 ( 2002 ) . t. d. gierke , g. e. munn , f. c. wilson , j. polym .
phys . ed .
* 19 * , 1687 - 1704 ( 1981 ) .
w. y. hsu , t. d. gierke , macromolecules * 15 * , 101 - 105 ( 1982 ) ; _ ibid _ j. memb .
* 13 * , 307 - 326 ( 1982 ) .                                                       for ( a ) run  1 ( low density ) , and ( b ) run  7 ( high density ) for @xmath67=1 ( dot - dashed line ) , 15 ( dashed line ) and 80 ( full line ) .
simulation results for @xmath67=1 are scaled down by a factor of four to fit the figure dimensions .
the @xmath110 ( symbols ) in ( b ) , given for comparison , corresponds to the @xmath64 line in ( a ) .
[ fig4],title="fig:",scaledwidth=60.0% ]   for ( a ) run  1 ( low density ) , and ( b ) run  7 ( high density ) for @xmath67=1 ( dot - dashed line ) , 15 ( dashed line ) and 80 ( full line ) .
simulation results for @xmath67=1 are scaled down by a factor of four to fit the figure dimensions .
the @xmath110 ( symbols ) in ( b ) , given for comparison , corresponds to the @xmath64 line in ( a ) .
[ fig4],title="fig:",scaledwidth=60.0% ]       for run  2 and different dielectric permittivities .
@xmath67=1 ( dot - dashed line ) , 15 ( dashed line ) and 80 ( full line ) .
results for run  1 , not shown here , are very close to the results for @xmath67=1 .
simulation results for @xmath67=1 are scaled down by a factor of four to fit the figure dimensions .
( b ) proton - proton pair distribution function @xmath41 for different amounts of backbone polymer ( runs 2 and 6 ) .
results for run 3 - 5 ( not shown here ) that correspond to intermediate amount of backbone polymer , lay between these two curves .
[ fig6 ] , title="fig:",scaledwidth=60.0% ]   for run  2 and different dielectric permittivities .
@xmath67=1 ( dot - dashed line ) , 15 ( dashed line ) and 80 ( full line ) .
results for run  1 , not shown here , are very close to the results for @xmath67=1 .
simulation results for @xmath67=1 are scaled down by a factor of four to fit the figure dimensions .
( b ) proton - proton pair distribution function @xmath41 for different amounts of backbone polymer ( runs 2 and 6 ) .
results for run 3 - 5 ( not shown here ) that correspond to intermediate amount of backbone polymer , lay between these two curves .
[ fig6 ] , title="fig:",scaledwidth=60.0% ]         for different box sizes @xmath79 , and hence sulfonate concentrations @xmath74 , for runs 8 and 10 - 12 .
the first peak of @xmath41 for run  10 ( full line , blue in colored figure ) has a value @xmath76 , out of the scale of the plot .
[ fig9],scaledwidth=90.0% ]     for different sulfonate molar concentrations @xmath74 , run  13 .
( a ) : 0.16
mol / l @xmath111 1.6 mol / l .
( b ) : 1.6
mol / l @xmath111 4.8 mol / l .
the case @xmath112 ( line with circles in ( b ) ) is for run  9 .
the first peak of @xmath41 for @xmath113 ( @xmath114 ) mol / l ( full blue line ( dashed green line ) in ( a ) ) has a value @xmath115 ( @xmath116 ) , which is out of the scale of the plot .
[ fig9a],title="fig:",scaledwidth=65.0% ] +   for different sulfonate molar concentrations @xmath74 , run  13 .
( a ) : 0.16 mol / l @xmath111 1.6
( b ) : 1.6
mol / l @xmath111 4.8 mol / l . the case @xmath112 ( line with circles in ( b ) ) is for run  9 .
the first peak of @xmath41 for @xmath113 ( @xmath114 ) mol / l ( full blue line ( dashed green line ) in ( a ) ) has a value @xmath115 ( @xmath116 ) , which is out of the scale of the plot .
[ fig9a],title="fig:",scaledwidth=65.0% ]    mol / l ( case ( a ) ) and @xmath117 mol / l ( case ( b ) ) head groups .
[ fig15],title="fig:",scaledwidth=65.0% ] + mol / l ( case ( a ) ) and @xmath117 mol / l ( case ( b ) ) head groups .
[ fig15],title="fig:",scaledwidth=65.0% ]           for run  1 ( solid line ) , run  2 ( dashed line ) and run  3 ( dot - dashed line ) .
thick lines : @xmath97 for hydrophobic - hydrophobic ( hh ) interaction , @xmath95 for hydrophobic - hydrophilic ( hp ) and hydrophilic - hydrophilic ( pp ) interactions .
thin lines : @xmath95 for hh , @xmath96 for hp and pp interactions .
[ fig16],scaledwidth=80.0% ]    2 . case ( a ) : @xmath95 for hh , @xmath96 for hp and pp interactions .
case ( b ) : @xmath97 for hh and @xmath95 for hp and pp interactions . in the case
( b ) the sulfonate groups are mostly on the surface of clusters.,title="fig:",scaledwidth=60.0% ] + 2 .
case ( a ) : @xmath95 for hh , @xmath96 for hp and pp interactions . case ( b ) : @xmath97 for hh and @xmath95 for hp and pp interactions . in the case
( b ) the sulfonate groups are mostly on the surface of clusters.,title="fig:",scaledwidth=60.0% ]       for runs 7 - 9 and @xmath104 ( a ) and @xmath105 ( b ) .
@xmath100 is a fictitious charge on the head groups which causes an artificial yukawa like attraction between sulfonates .
thick ( thin ) lines are for @xmath103 ( @xmath107 ) .
note that the scales on @xmath58-axis for short - range and long - range behavior of @xmath41 ( shown as an inset ) are different .
for other details see the text .
[ fig19],title="fig:",scaledwidth=70.0% ] +   for runs 7 - 9 and @xmath104 ( a ) and @xmath105 ( b ) .
@xmath100 is a fictitious charge on the head groups which causes an artificial yukawa like attraction between sulfonates .
thick ( thin ) lines are for @xmath103 ( @xmath107 ) .
note that the scales on @xmath58-axis for short - range and long - range behavior of @xmath41 ( shown as an inset ) are different .
for other details see the text .
[ fig19],title="fig:",scaledwidth=70.0% ]    2 , @xmath8=0 , @xmath100=0.1 ( case ( a ) ) and @xmath100=0.5 ( case ( b ) ) .
[ fig23],title="fig:",scaledwidth=70.0% ] + 2 , @xmath8=0 , @xmath100=0.1 ( case ( a ) ) and @xmath100=0.5 ( case ( b ) ) .
[ fig23],title="fig:",scaledwidth=70.0% ]    2 , @xmath8=0.5 , @xmath100=0.1 ( case ( a ) ) and @xmath100=0.5 ( case ( b ) ) .
[ fig24],title="fig:",scaledwidth=70.0% ] + 2 , @xmath8=0.5 , @xmath100=0.1 ( case ( a ) ) and @xmath100=0.5 ( case ( b ) ) .
[ fig24],title="fig:",scaledwidth=70.0% ]


the proliferation of mobile devices and the surge of a myriad of multimedia applications has resulted in an exponential growth of the mobile data traffic . in this context
, wireless caching has emerged as a powerful technique to overcome the backhaul bottleneck , by reducing the backhaul rate and the delay in retrieving content from the network .
the key idea is to store popular content closer to the end users . in @xcite , a novel system architecture named _
femtocaching _ was proposed .
it consists of deploying a number of small base stations ( bss ) with large storage capacity , in which content is stored during periods of offpeak traffic .
the mobile users can download content from the small bss , which results in a higher throughput per user . in @xcite , it was proposed to store content directly in the mobile devices .
users can then retrieve content from neighboring devices using device - to - device ( d2d ) communication or , alternatively , from the serving bs .    in both scenarios
, content may be stored using an erasure correcting code , which brings gains with respect to uncoded caching @xcite .
the use of erasure correcting codes establishes an interesting link between distributed caching for content delivery and distributed storage ( ds ) for reliable data storage .
the key difference is that in the wireless network scenario , data can be downloaded from the storage nodes ( the small bss or the mobile devices ) but also from a serving macro bs , which has always the content available .
therefore , the reliability requirements in ds for reliable data storage can be relaxed . in @xcite ,
the placement of content encoded using a maximum distance separable ( mds ) code to small bss was investigated and it was shown that a careful placement allows to significantly reduce the backhaul rate . in @xcite , for the scenario where content is stored directly in the mobile devices , the repairing of the lost data when a device storing data leaves the network was considered . assuming instantaneous repair , the communication cost of data download and repair was investigated . in @xcite , a repair scheduling where repair is performed periodically was introduced and analytical expressions for the overall communication cost of content download and data repair as a function of the repair interval were derived . using these expressions , the overall communication cost entailed by storing content using mds codes , regenerating codes @xcite , and locally repairable codes @xcite was evaluated in @xcite and it was shown that storing content using erasure correcting code can reduce the overall communication cost with respect to the scenario where content is downloaded solely from the bs .    in this paper
, we consider a similar cellular network scenario as the one in @xcite , where content is stored in a number of mobile devices using an erasure correcting code .
mobile devices roam in an out of a cell according to a poisson random process .
however , as opposed to @xcite , where the download of a single file is considered , here we consider that users may request files , of different popularity , from a library of files .
our focus is on the delay of retrieving content from the network , which was not considered in @xcite .
we derive analytical expressions for the download delay if content is stored in the mobile devices using mds codes and show that mds - coded distributed caching can significantly reduce the download delay with respect to the case where content is solely downloaded from the bs and the case where uncoded caching is used .
the download delay of a single file was analyzed in  @xcite .
the remainder of the paper is organized as follows .
the system model is introduced in section  [ s : system_model ] .
the average download delay incurred when mds - coded distributed caching is used is analyzed in sections  [ s : download ] and section  [ s : d2d ] . section  [ sec : numericalresults ] presents and discusses numerical results and finally some conclusions are drawn in section  [ sec : conclusions ] .
_ notation .
_ the probability density function ( pdf ) of a random variable @xmath0 is denoted by @xmath1 and the expectation with respect to @xmath0 is denoted by @xmath2 .
probability is denoted by @xmath3 and @xmath4 represents the all - ones vector of length @xmath5 .
we denote by @xmath6 the stationary distribution of an @xmath7 queueing system described by a poisson birth - death process with arrival rate @xmath8 and departure rate per node @xmath9 , which is given by @xmath10 where @xmath11 .
we consider a single cell in a cellular network where @xmath12 mobile devices , referred to as nodes , request files , each of size @xmath13 bits , from a library of @xmath14 files .
the files have different popularities and accordingly have a given probability to be requested .
depending on the placement strategy , some files are encoded and stored into @xmath15 mobile devices , referred to as storage nodes .
for ease of language , the set of storage nodes is referred to as the ds network and nodes not storing any content are referred to as _
regular _ nodes .
a copy of each encoded file is also available at the bs serving the cell .
a node requesting a file attempts to retrieve it from the storage nodes using d2d communication , and , if the file can not be completely retrieved from the ds network , the bs assists in providing the missing data . in order to increase the system efficiency ,
we allow multiple d2d communications to coexist if they are sufficiently far apart in space .
therefore , we divide the cell in @xmath16 virtual clusters and assume that the size of the cluster and the transmit power are properly chosen such that only one d2d communication can be established between any two nodes in the cluster and the interference across different clusters can be neglected .
a similar model is considered in  @xcite .
_ data allocation and coding strategy .
_ we adopt a deterministic allocation strategy , where the @xmath17 most popular files are cached in a distributed fashion in @xmath18 storage nodes in the cell , according to the storage capacity of the devices .
these files are partitioned into @xmath19 packets , called symbols , of @xmath20 bits each and are encoded into @xmath18 coded symbols using an @xmath21 mds erasure correcting code of rate @xmath22 .
we use the same code for every file in order to simplify the analysis .
we assume that each storage node stores a single symbol for each of the @xmath23 most popular files .
overall , @xmath24 symbols are stored in @xmath18 storage nodes and no two storage nodes store the same symbol .
we model the popularity of the files in the library using the time - invariant zipf distribution  @xcite .
accordingly , the probability that the @xmath5th file is requested is @xmath25 where parameter @xmath26 regulates the relative popularity of the files . in the following , the set of @xmath23 files stored in the cache of the mobile devices will be referred to as the ds library .
we assume that the mobile devices are free to move inside the cell .
we consider a uniform spatial distribution of the nodes in the cell , and hence there are @xmath27 devices per cluster on average and among them @xmath28 storage nodes .
we focus on a single cluster in isolation , and assume that the devices roam in and out of it .
the arrival , departure and request model of the nodes are borrowed from @xcite .
the considered scenario is shown in fig .
[ f : cluster ] .    _ arrival - departure model .
_ we assume that nodes arrive to the cluster according to a poisson random process with exponential independent , identically distributed ( i.i.d . ) random inter - arrival times @xmath29 with pdf @xmath30 where @xmath31 is the expected arrival rate and @xmath32 is time , measured in time units ( t.u . ) .
the nodes stay in the cluster for an i.i.d . exponential random lifetime @xmath33 with pdf @xmath34 where @xmath35 is the expected departure rate per node .
we assume that @xmath36 , which implies that the expected number of nodes in the cluster is @xmath37 .
this model corresponds to an @xmath7 queuing model and the probability that there are @xmath5 nodes in the cluster is @xmath38 .
the arrival of storage nodes to the cluster can also be described as a poisson random process . in particular , the inter - arrival times @xmath39 of the set of storage nodes has pdf @xmath40 the related lifetime is described by  ( [ e : departure ] ) and the probability that there are @xmath5 storage nodes in the cluster is @xmath41 .. more precisely , the only guarantee is that there are on average @xmath42 storage nodes per cluster , but there are no constraints on their instantaneous number , which can even exceed @xmath18 . on the other hand , the probability of having a high number of storage devices in one cluster is generally very low . for @xmath43 , we have @xmath44 , @xmath45 and @xmath46 .
the same consideration holds for the total number of mobile devices . ]     mobile devices , and @xmath42 storage nodes among them ( red circles ) , caching one different coded symbol for each of the most popular files .
a device requesting a file ( pink circle ) , must collect @xmath19 symbols .
it attempts to recover them by using the ds network if the requesting file is stored in the devices .
it uses the bs to collect the symbols that it is not able to download from the devices .
the download of a symbol from a storage node takes @xmath47 t.u . , and from the bs @xmath48 t.u .. ,width=332 ]    _ ds network update .
_ we assume that the nodes storing content that arrive to the cluster from neighboring clusters are not immediately available for download , but the bs serving the cell keeps track of them and periodically updates and broadcasts to all mobile devices the list of storage nodes in the cell every @xmath49 t.u .. in the sequel , parameter @xmath49 is referred to as the update interval and the set of storage nodes in the list broadcasted by the bs as the ds list .    _ data delivery . _
nodes request the file at random times with i.i.d .
random inter - request time @xmath50 with pdf @xmath51 where @xmath52 is the expected request rate per node .
we focus on the download process .
the node that requests a file attempts to retrieve it from the ds network using d2d communication .
thanks to the mds property , an encoded file can be reconstructed by accessing any @xmath19 encoded symbols .
if the file can not be completely retrieved from the ds network , the bs assists in providing the missing coded symbols .
the download of a coded symbol from a storage node incurs @xmath47 t.u . and from the bs @xmath48 t.u .. we assume that @xmath53 due to the congestion of the bs - to - node link and the fact that d2d communication occurs over a better channel due to the reduced distance between the involved nodes .
we further assume that only one d2d link at a time can be established , and that the d2d communication does not interfere with the communication between the bs and the nodes .
we say that the d2d network is _ idle _ if there is no active d2d communication in the cluster . if the d2d network is not idle when one node requests the file , the whole file is downloaded from the bs .
moreover , to simplify the analysis , we assume that multiple bs - to - node links can coexist .
we investigate the average time that is required to retrieve one file from the wireless network , referred to as the download delay . if a requested file is stored in the ds library , the requesting node attempts to retrieve it from the ds network using d2d communication , otherwise the file is entirely downloaded from the bs .
therefore , we introduce the binary random variable ( rv ) @xmath54 which describes the hitting of the ds cache , i.e. , @xmath55 if a file of the ds library is requested and @xmath56 otherwise . moreover , the d2d network can be used only if it is idle , i.e. , if there are no active d2d communications .
accordingly , we introduce the binary rv @xmath57 that describes the status of the d2d network . @xmath58 if the network is idle and @xmath59 otherwise .
if the d2d network is idle , the requesting node tries to collect the necessary coded symbols from the nodes of the ds list provided by the bs using d2d communication .
if the requesting node is a storage node of the ds list , it needs to download @xmath60 symbols , otherwise @xmath19 symbols must be downloaded .
we thus introduce the binary rv @xmath61 , which represents the type of request , i.e. , @xmath62 for requests originating from a node that belongs to the ds list and @xmath63 for the other requests .
the download from the storage nodes can be either fully successful or partially accomplished , in which case the requesting node turns to the bs to recover the missing symbols . on the other hand ,
if the d2d network is not idle and the requested file is stored in the ds library , the node downloads @xmath19 or @xmath60 symbols from the bs , depending on the type of node . from the discussion above , the average file download delay , @xmath64 , may be formalized as    the average file download for the cellular network described in section  [ s : system_model ] where the @xmath23 most popular files are stored in the mobile devices using an @xmath21 mds code is @xmath65 where @xmath66 is the average number of coded symbols downloaded per request using d2d communication and @xmath67 , referred to as the average d2d download delay , is the corresponding delay .
the computation of @xmath66 , @xmath67 and @xmath68 is addressed in section  [ s : d2d ] .
the probability of hitting the cache can be expressed as @xmath69 where the probabilities @xmath70 are given in  ( [ e : zipf ] ) .
it follows that @xmath71 if @xmath72 .
the next step is the computation of the probability that the d2d network is idle .
let @xmath73 be the status of the network at the time of the @xmath74th request .
it follows @xmath75 in order to compute @xmath76 , we introduce the rv @xmath77 that denotes the time instant of the @xmath78th request .
also , let @xmath79 be the time during which the d2d network is occupied by the @xmath78th request .
the d2d network is idle at the time of the @xmath74th request if none of the previous requests is still using d2d communication .
therefore , @xmath80 and @xmath81 assuming that if the d2d network is not idle at time @xmath82 is because of the @xmath83th request , the product in ( [ e : p_i exact ] ) reduces to the term involving the @xmath83th request only , i.e. , @xmath84 since the requests are i.i.d . with inter - request time distributed as in ( [ e : request time ] ) and on average there are @xmath37 nodes in the cluster , we can compute @xmath85 and ( [ e : p_i approx1 ] ) can be written as @xmath86 if @xmath87 , @xmath88 and @xmath89 in  ( [ e : idle_ell ] ) , we used the fact that the probability of hitting the cache and the average d2d download delay are independent of the request index ( if @xmath74 is sufficiently large ) , as it is proven in lemma  [ l : tdd_ell ] in section  [ s : d2d ] . substituting ( [ e : idle_ell ] ) in ( [ e : idle ] ) and after some simple calculations , we obtain @xmath90
in this section , we consider the computation of the average d2d download delay @xmath67 and the average number of coded symbols @xmath66 downloaded per request using d2d communication .
we assume that a node can not download in parallel from multiple nodes , but it serially tries to download the coded file symbols from the nodes in the ds list .
when a node requests the file , if the d2d network is idle and the requested file belongs to the ds library , it randomly chooses one of the storage nodes from the list supplied by the bs .
after each downloaded symbol , the requesting node randomly chooses another storage node from the ds list and still alive .
we assume that a requesting node that has collected fewer than the @xmath19 symbols necessary to reconstruct the file turns to the bs when all the reference storage nodes left or when the download of a symbol fails , even if other storage nodes are available . to simplify the analysis
, we assume that both cases ( the failed symbol download and the absence of storage nodes ) incur @xmath47 t.u .
, even if the node could contact the bs earlier .
we also assume that the download from the d2d network fails if the requesting node itself leaves the cluster before collecting @xmath19 symbols . in this case
, the download is also completed from the bs .    to derive the average d2d download delay , we introduce three rvs describing the number of nodes of different type that are present in the cluster at the instant of a request : the number of storage nodes of the ds list , the total number of storage nodes ( belonging or not to the list , the latter corresponding to the storage nodes that arrive to the cluster after the ds list update and that have not left the cluster at the time of the request ) , and the number of regular nodes . in particular , we denote by @xmath91 the rv that describes the number of storage nodes of the ds list when a request arrives .
we describe by the rvs @xmath92 and @xmath93 the total number of storage nodes and the number of regular nodes at the instant of a request , respectively .
moreover , we denote by @xmath94 the rv that represents the total number of storage nodes ( belonging or not to the ds list ) at the beginning of the update interval of length @xmath49 . in the following three lemmas , we give a probabilistic description of the above rvs .
[ l : px1 ] the probability that there are @xmath95 storage nodes of the ds list at the time of a request is @xmath96 where @xmath97 is the probability that @xmath98 is equal to @xmath99 , given that @xmath100 storage nodes are in the cluster at the beginning of the update interval of length @xmath49 , and is @xmath101 where @xmath102 and @xmath103 .    the proof is given in appendix  [ app : prooflem1 ] .
[ l : pq ] the probability that there are @xmath104 storage nodes in the cluster at the time of a request is given by @xmath105    the proof follows the same lines as the proof of lemma  [ l : px1 ] .    [ l : pv ] the probability that there are @xmath106 regular nodes in the cluster at the time of a request is given by @xmath107    the proof follows the same lines as the proof of lemma  [ l : px1 ] .    based on the above lemmas
, we can compute the probability that the request originates from a storage node of the ds list and the probability of having a given number of storage nodes in the ds list at the time of the request conditioned to the type of request .    using bayes rule , the probability that there are @xmath108 storage node of the ds list alive at the time of a request , conditioned to the type of request , is given by @xmath109 the probability @xmath110
is given in lemma  [ l : px1 ] .
we now compute @xmath111 and @xmath112 .
we start with the probability of having one request from the ds list conditioned to the number of storage nodes in the ds list at the time of the request . for @xmath113
, it can be written as @xmath114 clearly , the condition @xmath115 implies that the request can not originate from a storage node of the ds list , therefore @xmath116 .
we approximate the probability that there are @xmath117 storage nodes at the instant of the request , given the number of nodes of the ds list alive , by using the steady state probability of a poisson birth - death process with arrival rate @xmath118 and departure rate @xmath119 .
in particular , we compute @xmath120 where the expectations @xmath121 and @xmath122 are obtained starting from the probabilities  ( [ e : pq ] ) and ( [ e : px1 ] ) , respectively .
the number of regular nodes @xmath123 is independent of the number of storage nodes at the instant of the request , therefore we finally have @xmath124 where @xmath125 is given in lemma  [ l : pv ] .
note that in the expression above , with some abuse of notation , we used equal sign to avoid carrying all the way the approximation sign due to the approximation introduced in ( [ e : pq_approx ] ) .
starting from this result , we compute the probability that the request originates from the ds list as @xmath126 where @xmath127 is given in lemma  [ l : px1 ] .
the probability @xmath128 is easily computed as @xmath129 .
similarly , we have @xmath130 . following the same approach for the proof of lemma  [ l : px1 ] , it is easy to show that @xmath131 and @xmath132 are independent of the specific request ( when @xmath74 grows large ) , where @xmath133 is the binary rv describing the type of the @xmath74th request .
the case @xmath134 represents the case of instantaneous update , where the nodes contact directly the bs when they request a file and receive the list of the storage nodes through a dedicated link .
for instantaneous update , the number of storage nodes at the instant of the request and the type of request is described by the following probabilities @xmath135 the probability that there are @xmath95 storage nodes at the time of the request given the type of request can be computed by replacing the above probabilities in  ( [ e : px1cond ] ) .    in order to describe the d2d download , let @xmath136 be the binary rv that describes the success of the download at the first attempt .
more precisely , @xmath137 represents the successful download of the coded symbol from the first contacted storage node . if the download is not successful from the first contacted storage node , @xmath138 .
similarly , we denote by @xmath139 the binary rv describing the download at the @xmath78th attempt and we denote by @xmath140}$ ] , @xmath141 the random vector ( @xmath142 ) . in the following , in lemmas  [ l : s1 ] ,
[ l : sk ] , and [ l : sj ] , we derive the probability that no symbols can be downloaded from the d2d network , @xmath143 , the probability that the content is fully recovered from the ds network , @xmath144}={\mathbf{1}}_{k - i } | r\ } $ ] , and the probability that it is only partially recovered , @xmath145}={\mathbf{1}}_{j } , s_{j+1}=0 | r\}$ ] , respectively .
[ l : s1 ] the probability that no symbols are downloaded through d2d communication , conditioned to the type of request , is given by @xmath146 where @xmath147 is given in  ( [ e : px1cond ] ) and @xmath148 with @xmath102 .
the proof is given in appendix  [ app : prooflem4 ] .
[ l : sk ] the probability that the file can be completely retrieved from the ds network , i.e. , the probability that k symbols are downloaded through d2d communication when @xmath63 , or @xmath60 when @xmath62 , conditioned to the type of request , is given by @xmath149 } = { \mathbf{1}}_{k - i } |r = i\ } =   e^{-(k - i)\mu { t_\mathrm{d } } }    \sum_{g = 1}^\infty \sum_{d=0 } ^{g } \frac{g - d}{g } \gamma_{k - i}(g , d , i )   , i=0,1,\ ] ] where @xmath150 is defined by the recursion @xmath151 for @xmath152 and @xmath153 , with initial condition @xmath154 and where @xmath155 the function @xmath156 is given in  ( [ e : theta ] ) , and @xmath157 is given in  ( [ e : px1cond ] ) .    the proof is given in appendix  [ app : prooflem5 ] .
[ l : sj ] the probability of consecutively download @xmath158 symbols and to fail the download of the @xmath159th one is @xmath160}={\mathbf{1}}_{j } , s_{j+1}=0 |r = i \ } =   \gamma_{j+1}(0,0,i)a_{j+1 } +   \sum_{g=1}^\infty \sum_{d=0}^d \big (   \frac{d}{g } \gamma_{j+1}(g , d , i)a_{j+1 } + \frac{g - d}{g } \gamma_{j}(g , d , i)b_{j+1 }   \big )    \nonumber\ , , \ ] ] where @xmath161 , @xmath162 , and @xmath163 is given in  ( [ e : gamma ] ) .
the proof follows the same lines as the proof of lemma  [ l : sk ] .    finally , the average d2d download delay and the average number of downloaded symbols from ds network are given in the following theorem .
[ t : t_ni ] consider the network described in section  [ s : system_model ] , where an @xmath21 mds erasure correcting code is employed and where there are @xmath42 storage nodes in the cluster on average .
let @xmath47 be the time to download a symbol through d2d communication .
the average d2d download delay and the corresponding average number of downloaded symbols are given by @xmath164 where @xmath165 is the probability that the request comes from a storage node of the ds list , and @xmath166}={\mathbf{1}}_{k - i } | r = i\ }   + \sum_{j=1}^{k-1-i } j    \pr \{{\mathbf{s}}_{[j]}={\mathbf{1}}_{j } , s_{j+1}=0| r = i \ } \ , , \nonumber\\ { \overline{t}_i   } = & { t_\mathrm{d}}\big ( { \eta_i } +   c_{k , i}\pr \{\ s_{1}=0 | r = i\ } + \sum_{j=0}^{k-1-i }    \pr \{{\mathbf{s}}_{[j]}={\mathbf{1}}_{j } , s_{j+1}=0 | r = i\ } \big)\ , , \qquad \qquad i=0,1 \nonumber   \end{aligned}\ ] ] where @xmath167 for @xmath168 and @xmath169 otherwise .
the average d2d download delay is obtained as the sum of the average d2d delays in the case of requests originated from the ds list and of requests originating from the other nodes , weighted by the probabilities @xmath170 and @xmath171 , respectively .
the same approach is used for the corresponding average number of downloaded symbols . according to our model
, the requesting node completes the download of @xmath172 symbols from the ds network in @xmath173 t.u . with probability @xmath144}={\mathbf{1}}_k
| r = i \ } $ ] , while the partial download of @xmath174 symbols happens with probability and incurs @xmath175 t.u .. for @xmath168 , in the computation of the average d2d download delay , we also consider the case where download from the ds network completely fails .
the corresponding probability is @xmath176 and the delay is @xmath47 .
when the request originates from the ds list and @xmath177 , i.e. @xmath178 , no symbols need to be downloaded , therefore @xmath179 and @xmath180 are equal to zero .    [
l : tdd_ell ] the average d2d download time for the @xmath74th request , @xmath181 , is independent of the specific request if the index @xmath74 is sufficiently large .    similarly to the average d2d download delay , @xmath181 is @xmath182 where @xmath183}={\mathbf{1}}_{k - i } | r^{(\ell)}=i\ } +   c_{k , i}\pr \{\
s^{(\ell)}_{1}=0 | r^{(\ell)}=i\ } \nonumber\\                & + \sum_{j=0}^{k-1-i } ( j+1 )   \pr \{{\mathbf{s}}^{(\ell)}_{[j]}={\mathbf{1}}_{j } , s^{(\ell)}_{j+1}=0 | r^{(\ell)}=i\ } \big)\ , , \qquad \qquad i=0,1 \nonumber \end{aligned}\ ] ] the lemma follows from the fact that the probabilities in the expressions above are independent of @xmath74 , when @xmath74 grows large .
in this section , we evaluate the average download delay when content is cached using mds codes for a cluster with @xmath184 nodes on average , departure rate @xmath185 , and request rate @xmath186 . we compare the average file download delay @xmath64 of the considered network with mds - coded ds with the delay of the traditional scenario where the content is solely downloaded from the bs , denoted by @xmath187 and with uncoded caching . in the following , with no loss of generality
, we set @xmath188 t.u ..    we recall that each file that is cached in the ds library is divided into @xmath19 symbols , and encoded using an ( @xmath189 ) mds code , where @xmath18 is the number of storage nodes in the cell .
the code parameters are chosen such that @xmath190 . in this way ,
the _ average storage overhead _ in a cluster , @xmath191 , is positive , which increases the probability that the content is downloaded through d2d communication only .
alternatively , we can use the same @xmath192 mds code for each cluster , but in this case the bs must continuously restore the initial state of reliability of the ds network when storage nodes leave the clusters  @xcite .
we first consider the special case where @xmath72 , therefore the probability of hitting the cache @xmath193 is @xmath194 . in figs .
[ f:0.1][f:0.001 ] , we show the gain that can be achieved using mds - coded caching , by reporting the ratio between @xmath195 and @xmath64 as a function of the update interval @xmath49 .
the infinite series involved in the computation of @xmath64 are truncated to a given value @xmath32 , chosen according to @xmath196 when involving the number of storage nodes and to @xmath197 when involving the number of nodes in general .
we fix the ratio @xmath198 to be @xmath199 , and consider several mds codes .
moreover , we also consider an uncoded scenario where one storage node on average in the cluster stores the uncoded files . in figs .
[ f:0.1 ] , [ f:0.01 ] , and [ f:0.001 ] , @xmath47 is @xmath200 , @xmath201 , and @xmath202 times , respectively , smaller than @xmath48 . in the figures ,
the solid lines correspond to the analytical closed - form expressions derived in the previous sections and markers correspond to simulation results .
it is observed that the analytical expressions predict very well the actual performance , which shows the goodness of the approximations introduced in  ( [ e : approxi ] ) and ( [ e : pq_approx ] ) .
the results clearly show that mds - coded ds can greatly improve the performance in terms of content download delay with respect to the case where content is downloaded from the bs , provided that the update interval , @xmath49 , is sufficiently small .
for example , for @xmath203 and @xmath204 , a speed - up factor of around @xmath205 in the download is achieved with respect to the case of downloading from the bs using a @xmath206 mds code .
interestingly , the results also show that the performance improves when @xmath19 increases .
in particular , simple replication ( repetition coding ) is very inefficient and much better performance are achieved using larger mds codes ( of the same rate ) .
we now consider the more general case where only part of the library of files is cached in the devices .
we assume that the library has a size of @xmath207 files , and each storage node stores one symbol for each of the @xmath23 most popular files .
we assume that the each file is of size @xmath201 mb , which corresponds to a @xmath200-minutes video .
we further assume that each storage node makes @xmath208 gb available for caching . in this case , the number of cached files @xmath23 , that corresponds to the number of symbols cached by each storage node , increases by increasing @xmath19 , since the size of one encoded symbols is @xmath201mb@xmath209 and @xmath23 is clearly obtained by dividing the storage capacity of the node by the symbol size . in fig .
[ f : sigma ] , we show the download speedup factor @xmath210 as a function of the parameter @xmath26 , which regulates the relative popularity of the files .
for example , a large value of @xmath26 represents the case where few popular files are responsible for the majority of the download traffic .
the figure refers to the case where @xmath211 and @xmath212 .
the results confirm the gain that can be achieved by mds - coded distributed caching , but highlight another important aspect .
when the whole library is not cached in the network , i.e. , @xmath213 , the gain reduction is not negligible , especially for low values of @xmath26 .
this fact suggests that the adopted deterministic allocation , where the @xmath23 most popular files are equally stored in the ds network , is suboptimal .    .
solid lines show analytical results and markers simulation results . ]    .
solid lines show analytical results and markers simulation results . ]    .
solid lines show analytical results and markers simulation results . ]    . @xmath214 and @xmath212 . ]
in this paper , we considered the cache of popular content in the mobile devices of a cellular network using maximum distance separable erasure correcting codes to speed - up content delivery .
we derived analytical expressions for the average download delay and showed that mds - coded distributed caching may dramatically reduce the download delay with respect to the traditional case where content is always downloaded from the base station .
we denote by @xmath215 the number of storage nodes available for download at the time of the @xmath74th download request , i.e. , the number of storage nodes of the ds list that have not left the cluster at the time of the request .
we compute @xmath127 by averaging over an infinite number of requests , @xmath216 similarly , let @xmath217 be the number of storage nodes at the beginning of the update interval wherein the @xmath74th request arrives , denoted by @xmath218 .
we have @xmath219 in  @xcite , it was shown that the probability @xmath220 does not depend on @xmath74 ( when @xmath74 grows large ) , and is given by  ( [ e : px_j ] ) .
its derivation is based on the observation that the number of storage nodes available for download in the update interval is described by a poisson death process .
the probability @xmath221 can be written as @xmath222 where in the second equality we used bayes rule . in  ( [ e : y ] ) , @xmath223 is the number of storage nodes at the beginning of the update interval , which is described by a birth - death poisson process , thus @xmath224 .
the probability @xmath225 is the probability that there is at least one request in @xmath218 .
it depends on the inter - request time , which in turn depends on the number of nodes in the cluster .
therefore , we compute @xmath226 where @xmath227 is the probability that the inter - request time is shorter than @xmath49 when @xmath228 nodes are present in the cluster . similarly , we compute @xmath229 as @xmath230 where @xmath231 is the probability that there are @xmath228 nodes in the cluster , given that there are @xmath232 storage nodes .
since these probabilities are independent of the specific request , we conclude that @xmath221 is also independent of @xmath74 . substituting ( [ e : y ] ) into ( [ e : x1l ] ) , we observe that @xmath233 is also independent of @xmath74 and using  ( [ e : x1 ] ) we prove the lemma .
we compute the conditional probability that no symbols are downloaded by averaging over an infinite number of requests , @xmath234 where @xmath235 is the rv describing the download of the first symbol for the @xmath74th request .
we now consider the computation of @xmath236 .
the recovery of the first symbol fails with probability @xmath194 if the requesting node leaves the cell before completing the download .
it also fails if the requesting node stays in the cluster but no storage nodes are available or if it chooses to download from a storage node which departs before @xmath47 t.u . from the start of the download
let @xmath237 be the departure time of the node which places the @xmath74th request and let @xmath238 be the event that the node which places the @xmath74th request stays in the network for more than @xmath239 t.u . from the start of the download .
the corresponding probability does not depend on @xmath74 and is easily computed as @xmath240 .
similarly , the probability that the requesting node departs before @xmath47 t.u . from the start of the download is @xmath241 .
therefore , the conditional probability that the @xmath74th request fails the first symbol download is @xmath242 let @xmath243 be the number of storage nodes useful for download for the @xmath74th request , i.e. , the number of storage nodes in the ds list that have not left the cluster at the time of the @xmath74th request , excluding the requesting node itself if it belongs to the ds list .
let @xmath244 the number of departures in @xmath47 t.u . among the @xmath245 storage nodes .
the probability @xmath246 can be written as @xmath247 equation  ( [ e : s1cond ] ) is obtained by observing that @xmath248 the probability that the download of the first symbol fails conditioned to @xmath249 and @xmath250 is independent of the type of request , and that @xmath251 the number of useful storage nodes and the number of departures in @xmath47 t.u . among them is independent of the departure of the requesting node .
the probability @xmath252 is equal to 1 if there are no useful storage nodes , i.e. , @xmath253 .
otherwise , it equals the probability to choose one of the @xmath254 storage nodes that leave the cell in @xmath47 t.u .
, i.e. , @xmath255 , with @xmath256 .
we observe that the number of departures of useful storage nodes conditioned to their number is independent of the type of request and that the number of storage nodes useful for download @xmath249 is related to @xmath257 by @xmath258 since when the request originates from a storage node of the ds list , the requesting node itself is not counted among the useful storage nodes . therefore , the probability @xmath259 can be written as @xmath260 we denote by @xmath156 the probability @xmath261 , given in  ( [ e : theta ] ) .
its derivation is similar to that of @xmath262  @xcite .
the probability @xmath263 is independent of the specific request and is given by  ( [ e : px1cond ] ) .
after simple manipulations , we finally obtain @xmath264 since the probabilities involved in  ( [ e : ps1 ] ) are all independent of @xmath74 , the lemma is proved .
to evaluate the probability of complete download from the ds network , we start with the following limit , @xmath265 } = { \mathbf{1}}_{k - i}| r = i \}= \lim_{l \to \infty }    \frac{1}{l } \sum_{\ell=1}^l \pr\ {   { \mathbf{s}}_{[k - i]}^{(\ell)}={\mathbf{1}}_{k - i}|r^{(\ell)}=i \}\ , , \ ] ] where @xmath266}^{(\ell)}= (   s_1^{(\ell ) } , \dots , s_j^{(\ell)})$ ] and @xmath267 describes the successful symbol download at the @xmath5th attempt of the @xmath74th request .
we consider the @xmath74th request and , similarly to the proof of lemma  [ l : s1 ] , we will find that this probability is independent of @xmath74 .
we denote by the rv @xmath268 the number of storage nodes useful for download at the time of the @xmath78th attempt of the @xmath74th request , i.e. , the storage nodes of the ds list not yet contacted , excluding the requesting node if it belongs to the ds list .
we denote by the rv @xmath269 the number of departures in @xmath47 t.u . among the @xmath270 nodes .
we also denote by @xmath271 the event that the node which places the @xmath74th request stays in the network for more than @xmath272 t.u . from the start of the download .
the corresponding probability does not depend on @xmath74 and is given by @xmath273 .
the probability of complete download is zero if the requesting node departs before @xmath274 t.u . from the start of the download , therefore we can write @xmath275}^{(\ell)}={\mathbf{1}}_{k - i}|r^{(\ell)}=i \}=e^{-(k - i)\mu { t_\mathrm{d } } } \pr\ {   { \mathbf{s}}_{[k - i]}^{(\ell)}={\mathbf{1}}_{k - i}|\mathcal{a}_{k - i}^{(\ell ) } , r^{(\ell)}=i   \}\,.\ ] ] the probability @xmath276}^{(\ell)}={\mathbf{1}}_{k - i}|\mathcal{a}_{k - i}^{(\ell ) } , r^{(\ell)}=i   \}$ ] can be written as @xmath277}^{(\ell)}={\mathbf{1}}_{k - i}|\mathcal{a}_{k - i}^{(\ell ) } , r^{(\ell)}=i   \}=\nonumber\\ & = \sum_{g } \sum_{d } \pr\{s^{(\ell)}_{k - i}=1|   \mathcal{a}_{k - i}^{(\ell ) } , g_{k - i}^{(\ell)}=g , d_{k - i}^{(\ell)}=d   \ } \pr\ {   g_{k - i}^{(\ell)}=g , d_{k - i}^{(\ell)}=d,{\mathbf{s}}_{[k - i-1]}^{(\ell)}={\mathbf{1}}_{k - i-1 } |   \mathcal{a}_{k - i}^{(\ell ) } , r^{(\ell)}=i \}\nonumber \\ & = \sum_{g=1}^\infty \sum_{d=0}^{g } \frac{g - d}{g } \gamma^{(\ell)}_{k - i}(g , d , i)\ , .
\label{e : gamma1 } \end{aligned}\ ] ] the last equality is obtained by observing that the probability @xmath278 for @xmath279 and @xmath280 equals the probability to choose one of the storage nodes of the ds list that remains in the cluster , i.e. , @xmath281 .
moreover , we have defined @xmath282}^{(\ell)}={\mathbf{1}}_{j-1 } |   \mathcal{a}_{j}^{(\ell ) } , r^{(\ell)}=i \}$ ] , that can be computed as @xmath283}^{(\ell)}={\mathbf{1}}_{j-2 } |   \mathcal{a}_{j-1}^{(\ell ) } , r^{(\ell)}=i \}\ , , \nonumber\end{aligned}\ ] ] for @xmath284 .
we also define @xmath285 , which is equal to one if @xmath286 and @xmath287 , and zero otherwise .
the condition @xmath286 follows from the fact that the number of useful storage nodes after a successful symbol download is equal to the number of useful storage nodes still alive , , minus the storage node just used .
the condition @xmath287 comes from the fact that the @xmath288th symbol download is assumed to be successful , i.e. , @xmath289 .
it is easy to prove by induction that the probability in  ( [ e : gamma1 ] ) does not depend on @xmath74 . by defining @xmath290
, we obtain the recursion  ( [ e : gamma ] ) , with initial condition  ( [ e : gammaic ] ) .
the probabilities @xmath291 , @xmath158 , are equal to zero for @xmath292 .
k.  shanmugam , n.  golrezaei , a.  g. dimakis , a.  f. molisch , and g.  caire , `` femtocaching : wireless content delivery through distributed caching helpers , '' _ ieee trans .
theory _ , vol .
59 , no .  12 , pp . 84028413 , dec .
n.  golrezaei , p.  mansourifard , a.  f. molisch , and a.  g. dimakis , `` base - station assisted device - to - device communications for high - throughput wireless video networks , '' _ ieee trans .
wireless commun .
_ , vol .  13 , no .  7 , pp
36653676 , july 2014 .
j.  pedersen , a.  graell i amat , i.  andriyanova , and f.  brnnstrm , `` repair scheduling in wireless distributed storage with d2d communication , '' in _ proc .
. theory work .
( itw ) _ , jeju island , korea , 2015 .
j.  pedersen , a.  graell i amat , i.  andriyanova , and f.  brnnstrm , `` distributed storage in mobile wireless networks with device - to - device communication , '' _ ieee trans . commun .
_ , vol .
64 , no .  11 , pp . 48624878 , nov .
2016 .
a.  piemontese and a.  graell i amat , `` mds - coded distributed storage for low delay wireless content delivery , '' in _ proc .
turbo codes & iterative inform .
( istc ) _ , brest , france , sep .