Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: one of the main industrial applications of bitumen is as a binder in asphalt pavement  [ ] . bitumen links together mineral aggregates and filler particles to form a cohesive asphalt on the road surface . over the years , a chemical reaction takes place in bitumen increasing the number of heavy aromatic molecules  [ ] . this process is called chemical aging . the aromatic molecules in bitumen , especially the asphaltene molecules , tend to align to form nanoaggregates . the aging reaction leads also to an increase in the nanoaggregate size  [ ] correlated with an unwanted increase in bitumen viscosity and brittleness  [ ] . the change in bitumen mechanical properties , which go from `` liquid - like '' to more `` solid - like '' during chemical aging finally results in cracks in the pavement and road deterioration . to prevent or reverse the effects of chemical aging , a first step is to gather more knowledge about the nanoaggregate structure and stability in bitumen . much progress has been made over the last 50 years in the experimental and numerical literature to determine the structure of the nanoaggregates and the conditions under which they are formed . one can cite in particular the design of the yen - mullins model  [ , ] describing the nanoaggregate structure , the determination of the critical nanoaggregate concentration in different solvents  [ , ] , the evaluation of the average size and polydispersity of a nanoaggregate  [ ] on the experimental side . on the numerical side , different stable conformations of nanoaggregates were identified depending on the asphaltene structure using molecular mechanical calculations  [ ] and molecular dynamics ( md ) simulations  [ ] . md simulations were also used to determine the molecular orientation inside the nanoaggregates  [ ] , and the effects of solvent and presence of other molecules on the nanoaggregate structure  [ ] . however , analytical models for the thermodynamics stability and dynamics of the nanoaggregates in relation to their structure are still quite rare , to the notable exception of ref .  [ ] . the purpose of this paper is precisely to suggest simple and generic models , which can reproduce md results on the nanoaggregate stability . we present molecular dynamics results concerning the nanoaggregate size in the stationary regime . in bitumen , nanoaggregates are composed of asphaltene molecules , the most heavy and aromatic fraction in bitumen , but also of lighter aromatic molecules like resin and resinous oil  [ ] . we give results on the probability of having a nanoaggregate containing a given number of aromatic molecules and the probability of having a nanoaggregate containing a given number of asphaltene molecules . the two probabilities are shown to differ qualitatively . we first model the simpler case , where only asphaltene molecules are counted . the focus in this case is to establish a mechanism for the nanoaggregation of asphaltene molecules . from this aggregation mechanism , a master equation and a statistical mechanics model are derived and provide a simple theoretical framework for interpreting the simulation results on nanoaggregate stability . the results of the master equation approach are also compared to the aggregation dynamics observed in the md simulations . then , the more complicated case where all aromatic molecules are counted is discussed in terms of statistical mechanics arguments . the md simulations carried out in this work are based on the four - component united - atom - unit model developed in ref . [ ] in the framework of the cooee project  [ ] . the model is shown to reproduce a generic bitumen reasonably well  [ ] . the simulations were performed on graphic - processor - units ( gpu ) using the roskilde university molecular dynamics ( rumd.org ) package  [ ] . the paper is organized as follows . in sec . [ sec : results ] , we provide simulation details , give a definition of a linear nanoaggregate , and present the md results about the probability of having a nanoaggregate of a given size . section  [ sec : mastereq ] is devoted to model the results on the stability of linear aggregates where only asphaltene molecules are counted with a master equation approach . in sec . [ sec : statmech ] , a related statistical mechanics model is described and shown to reproduce the md results on asphaltene nanoaggregate stability . the more complicated case where all aromatic molecules are counted in a nanoaggregate is also discussed in sec . [ sec : statmech ] in terms of statistical mechanics arguments . [ sec : discussion ] includes a comparison of our results to existing experimental and numerical results and discuss the limit of our model . finally , sec . [ sec : conclu ] contains a summary and a conclusion . before presenting the molecular dynamics ( md ) results on nanoaggregate stability , we mention a few details about the simulations and define precisely a nanoaggregate .      as mentioned in the introduction , the simulation method and molecular potentials are described in details in ref . only information necessary to understand the study on nanoaggregate stability carried out in the present paper is given here . the simulated system contains four types of molecule , chosen to resemble the sara classification  [ ] : a saturated hydrocarbon , a resinous oil molecule , which is denoted aromatic in the sara scheme , a resin molecule and an asphaltene molecule . the molecular structures chosen are shown in fig . [ fig : molecule ] . the main system studied in this paper contains @xmath0 saturated hydrocarbons , @xmath1 resinous oil molecules , @xmath1 resin molecules and @xmath1 asphaltene molecules , which corresponds to @xmath2 united atom units . the methyl ( ch@xmath3 ) , methylene ( ch@xmath4 ) , and methine ( ch ) groups are represented by the same united atom unit of molar mass 13.3 g@xmath5mol@xmath6 and the sulfur atoms are represented by a united atom unit with a different molar mass 32 g@xmath5mol@xmath6 . the potential between the united atom units contains four terms : an intermolecular potential , corresponding to a lennard - jones potential with parameters @xmath7 @xmath8 and @xmath9 k , where @xmath10 is the boltzmann constant and three terms for the intramolecular potential . these three terms describe the bond length between two connected particles , the angle between three consecutive particles , and the dihedral angle between four consecutive particles . the parametrization of the intramolecular potential is described in details in refs  [ ] . the simulations are performed in the canonical ensemble ( nvt ) at a constant temperature @xmath11 k and a constant density . the density @xmath12 kg@xmath5l@xmath6 is chosen to obtain an average pressure around the atmospheric pressure . a nos - hoover thermostat is used . the time step is @xmath13 fs and the duration of the simulations is @xmath14 ns . eight independent simulations are performed at the same state point . the molecular dynamics package rumd  [ ] is used to perform the calculation . +      it was shown in a previous work  [ ] that the asphaltene , resin , and resinous oil molecules tend to align with respect to each other . they align at a distance of around @xmath15  , close to the minimum of the lennard - jones potential between the molecules . the lennard - jones potential used in the md simulations mimics the @xmath16-stacking interaction observed between aromatic molecules experimentally . the alignment of aromatic molecules in bitumen is the basis of the nanoaggregate formation . the definition of a nanoaggregate is described in detail in ref . it is based on the following rule : two aromatic molecules are nearest neighbors in the same nanoaggregate if they are `` well - aligned '' and `` close enough '' . a nanoaggregate is composed of all molecules connected by this rule . moreover , the asphaltene molecule chosen in this model has two parts , a flat head and a flat body which can rotate with respect to each other . for this reason , aromatic molecules can align in the direction of an asphaltene body or in the direction of an asphaltene head , thus creating branched nanoaggregates . these branches can link together purely linear nanoaggregates . we believe this is one mechanism explaining the formation of clusters of nanoaggregates , also observed experimentally  [ ] . in this paper , we focus on purely linear nanoaggregates . it means that any aromatic molecules linked to an asphaltene head will not be considered as part of this asphaltene nanoaggregate . we define a linear aggregate as composed of asphaltene bodies , resin and resinous oil molecules , and not asphaltene heads because heads are the smallest parts and probably lead to the smallest interaction energy . figure  [ fig : linearaggregate ] shows conformations of two molecules corresponding to limiting cases of the nanoaggregate definition used in this paper . figures  [ fig : linearaggregate](a ) and ( b ) show two conformations of two asphaltene molecules where these molecules are considered as nearest neighbors . conversely , figs  [ fig : linearaggregate](c ) to ( e ) show three conformations where the two molecules are not considered as nearest neighbors . a picture of a linear nanoaggregate obtained from the md simulations is shown in fig . [ fig : linearaggregate](f ) . note that this nanoaggregate is not typical as a resin and a resinous oil molecule are aligned on the same side of an asphaltene body . ( color online ) ( a ) to ( e ) scheme of limiting cases illustrating whether a molecule is or is not the nearest neighbor of the first molecule . the first molecule is an asphaltene molecule and is black . the second molecule is green if it is the neighbor of the first one and red otherwise . ( a ) head - to - head conformation . ( b ) head - to - tail conformation . ( c ) non - aligned molecules . ( d ) molecule aligned to the head of the asphaltene molecule . ( e ) molecules far from each other and aligned . modified with permission from c. a. lemarchand , t. b. schroder , j. c. dyre , and j. s. hansen , j. chem . 139 , 124506 ( 2013 ) . copyright 2013 american institute of physics . ( f ) snapshot of a linear nanoaggregate , obtained in molecular dynamics . asphaltene molecules are in blue , resin molecules in red and resinous oil molecules in green . , title="fig : " ]   ( color online ) ( a ) to ( e ) scheme of limiting cases illustrating whether a molecule is or is not the nearest neighbor of the first molecule . the first molecule is an asphaltene molecule and is black . the second molecule is green if it is the neighbor of the first one and red otherwise . ( a ) head - to - head conformation . ( b ) head - to - tail conformation . ( c ) non - aligned molecules . ( d ) molecule aligned to the head of the asphaltene molecule . ( e ) molecules far from each other and aligned . modified with permission from c. a. lemarchand , t. b. schroder , j. c. dyre , and j. s. hansen , j. chem . 139 , 124506 ( 2013 ) . copyright 2013 american institute of physics . ( f ) snapshot of a linear nanoaggregate , obtained in molecular dynamics . asphaltene molecules are in blue , resin molecules in red and resinous oil molecules in green . , title="fig : " ]      sections  [ sec : simulationdetail ] and  [ sec : defnanoaggregate ] gave necessary information about the simulations and the definition of linear nanoaggregates . the results of the molecular dynamics simulations concerning the nanoaggregate size and stability can now be presented . there are different possibilities to quantify the size of a linear nanoaggregate . for example , it can be quantified as the total number of aromatic molecules or as the number of asphaltene molecules which it contains . the first definition makes use of the fact that nanoaggregates are not only composed of asphaltene molecules , but also of smaller aromatic molecules . the second definition accounts for the fact that asphaltene molecules are the largest aromatic molecules and probably the most important in the nanoaggregate stability . this is why experimentalists often define the nanoaggregate size as the latter  [ ] . the molecular dynamics simulations enable us to study the consequences of both definitions on the nanoaggregate stability . we quantify the stability as the probability of having a nanoaggregate of a given size in the stationary regime . the probability @xmath17 of having a nanoaggregate containing @xmath18 aromatic molecules is defined as : @xmath19 where @xmath20 is the number of linear aggregates containing @xmath18 aromatic molecules and @xmath21 is the total number of linear aggregates . in a similar way , the probability @xmath22 of having a nanoaggregate with @xmath18 molecules of type @xmath23 is defined as @xmath24 where @xmath25 is the number of linear aggregates containing @xmath18 molecules of type @xmath23 and @xmath26 is the total number of linear aggregates containing at least one molecule of type @xmath23 . the type @xmath23 can be @xmath27 for asphaltene , @xmath28 for resin , @xmath29 for resinous oil or @xmath30 for resin and resinous oil . the probabilities @xmath31 and @xmath32 are plotted in fig . [ fig : pdfall ] ( b ) and ( c ) , respectively . the probabilities @xmath33 and @xmath34 are plotted in fig . [ fig : pdfall ] ( d ) . two main points can be noticed in fig . [ fig : pdfall ] . the probability @xmath35 of having a nanoaggregate containing @xmath18 aromatic molecules seems to have two slopes in a log - lin scale , _ i. e. _ can be described by a biexponential . on the contrary , when only asphaltene molecules are counted , the probability @xmath36 has only one slope , _ i. e. _ can be characterized as a simple exponential . for the probability @xmath37 of having a nanoaggregate containing @xmath18 resin or resinous oil molecules , the biexponential shape seems to be prevailing . it is the same for the probability @xmath38 of having a nanoaggregate containing @xmath18 resin molecules . for the probability @xmath39 of having a nanoaggregate containing @xmath18 resinous oil molecules , the monoexponential shape seems to be recovered . the transition between the two slopes is not sharp for the probability @xmath38 . to highlight the existence of the two slopes in this case in fig . [ fig : pdfall ] ( d ) , a dashed line was drawn in continuation of the line corresponding to the first slope . it departs further and further away from the second slope . the rest of the paper is devoted to understand the molecular dynamics results presented in fig . [ fig : pdfall ] . more specifically , the paper is aimed at unraveling some of the dynamical and thermodynamical origins lying behind the monoexponential and biexponential behaviors . a dynamical approach will first be proposed for the monoexponential behavior . in fig .  [ fig : pdfall ] ( b ) , the probability of having a nanoaggregate containing @xmath18 asphaltene molecules was shown to be monoexponential . the aim of this section is to establish a simple aggregation mechanism reproducing this behavior . the mechanism will be used to derive a master equation which provides a dynamical framework for understanding the monoexponential behavior . a simple mechanism which can be written in the asphaltene aggregation problem is @xmath40 where @xmath41 denotes an aggregate containing @xmath18 asphaltene molecules . the assumptions lying behind this mechanism are the following : only one asphaltene molecule at a time can attach to or detach from an existing aggregate , the medium is homogeneous and the aggregation dynamics is not limited by diffusion . the same mechanism was already suggested in ref . [ ] to model asphaltene aggregation from a thermodynamical point of view , but the possibility to obtain a master equation from this mechanism and consequently to model the dynamics of asphaltene aggregation was left aside .    to obtain a master equation linked to this aggregation reaction , we make supplementary assumptions : the reaction rate constants do not depend on the size @xmath18 of the aggregate and the attachment rate does not depend on the probability of having free asphaltene molecules in the medium . under these assumptions , the master equation linked to the aggregation reaction  eq . can be written as : @xmath42 where @xmath43 is the attachment rate constant , @xmath44 the detachment rate constant and @xmath45 the probability of having an aggregate containing @xmath18 asphaltene molecules at time @xmath46 , given that its initial size was @xmath47 . the probability denoted @xmath31 in sec . [ sec : mdresults ] corresponds to the stationary state to eqs .   and   and depends neither on time @xmath46 nor on the initial condition @xmath47 . the master equation eq . was solved analytically in the 1950 s  [ , ] for different boundary conditions . this master equation and its generalized version in which the rate constants depend on the aggregate size @xmath18 have also been extensively used to model different stochastic processes such as alcohol clusters  [ , ] , biological adhesion clusters  [ ] or aggregation in freeway traffic  [ ] , to cite only a few . the expression of this master equation deserves some further clarification in connection to the specific aggregation reaction used . following the law of chemical kinetics , the attachment rate linked to the aggregation mechanism   should be @xmath48 , because the probability of having a collision between two small spherical molecules is usually proportional to the product of their concentrations . we can give empirical reasons explaining why it is not proportional to the probability @xmath49 of having a free asphaltene molecule here . we assume that the attachment rate is not due to the collision between two small spherical molecules , but limited by the steric hindrance around the aggregate and the orientation of the free asphaltene molecule @xmath50 . the attachment reaction is successful if the free asphaltene molecule is well oriented and placed at a reactive end of the aggregate . the number of free asphaltene molecules is higher than the number of any aggregate of a given size at any time in our md simulations ( not shown ) . as the number of free asphaltene molecules is high enough , there is always a free asphaltene molecule nearby an aggregate but it is not always well - oriented and prevent other molecules from approaching due to the density of the system . thus , the attachment rate is independent of the free asphaltene concentration and the attachment rate constant @xmath43 is averaged over possible orientations . the boundary condition  eq .   needs also some explanation . it corresponds to the fact that the system is closed . in other words , the total number of asphaltene molecules is kept constant . as the free asphaltene molecules are involved in each aggregation reaction  , the derivative with respect to time of the probability @xmath51 to have a free asphaltene molecule depends on infinitely many terms . it leads to @xmath52 these series are then simplified using the fact that at every time @xmath46 : @xmath53 , to give the boundary condition  eq .  .    to test the validity of this master equation for the present problem , two different quantities are evaluated : the stationary probability of having a nanoaggregate of a given size , as shown in sec . [ sec : mdresults ] and the aggregation dynamics , quantified as the time evolution of the fraction of aggregated asphaltene molecules .      the stationary distribution of the master equation  eq .   with the boundary condition  eq . is easy to derive . in the stationary state , the probability of having a nanoaggregate containing @xmath18 asphaltene molecules depends neither on time nor on the initial size of a nanoaggregate . for this reason , the stationary probability will be denoted @xmath31 . it satisfies the following equation : @xmath54 the solution of this equation is a geometrical law  [ , ] . by induction and using the fact that @xmath55 , one can show that the solution of this equation is : @xmath56 in a log - lin scale , the stationary probability predicted by the master equation approach is a straight line of the form : @xmath57 it is a monoexponential distribution , as observed in fig .  [ fig : pdfall ] ( b ) , obtained in md . the exponential distribution can be used to fit the molecular dynamics data , as shown in fig .  [ fig : aonlybb ] . the value of the dynamical parameter @xmath58 is , in this case , @xmath59 . the master equation framework provides a dynamical interpretation of the parameter @xmath58 as the ratio between the attachment and detachment rate constants . figure  [ fig : aonlybb ] also shows the probability of having a nanoaggregate containing @xmath18 asphaltene molecules obtained in md simulations for two different system sizes and the corresponding exponential fits . the results for a system containing @xmath1 asphaltene molecules and for a system containing @xmath60 times less molecules agree surprisingly well , giving the values @xmath59 and @xmath61 respectively . it shows that , while finite size effects are present , they are not very important when it comes to the number of asphaltene molecules in a nanoaggregate . more specifically , the relative probability of small nanoaggregates is accurately described even in the small system . the probability of having larger nanoaggregates is not so well described in the small system , as expected , but does not affect much the value of the parameter @xmath58 since these nanoaggregates are rare even in the large system . it is important to note , however , that physical properties due to large nanoaggregates or to nanoaggregates filling up the box in one direction , such as residual stresses , could be affected by finite size effects . as shown above , the stationary state predicted by the master equation agrees with the md results on the probability of having a nanoaggregate containing @xmath18 asphaltene molecules . if the master equation  eq .   correctly describes the aggregation process , it should also reproduce the aggregation dynamics . checking this fact is the purpose of this section .    to compare the prediction of the master equation approach and the molecular dynamics results , we quantified the aggregation dynamics in the following way . we used the fraction @xmath62 of aggregated asphaltene molecules versus time . the fraction @xmath62 of aggregated asphaltene molecules is defined as the ratio of the number of asphaltene molecules in all aggregates containing at least two asphaltene molecules and the total number of asphaltene molecules in the system . the fraction @xmath62 verifies : @xmath63 where @xmath64 is the fraction of free asphaltene molecules . the fraction of free asphaltene molecules is defined as @xmath65 where @xmath66 is the number of asphaltene molecules in aggregate containing one asphaltene molecule and @xmath67 is the total number of asphaltene molecules in the system . the fraction @xmath64 can be expressed in terms of the probability @xmath68 of having an asphaltene aggregate of size @xmath69 at time @xmath46 and of the average size of an asphaltene aggregate @xmath70 , given that all aggregates had size 1 initially . indeed , we know that : @xmath71 where @xmath72 is the number of aggregates containing one asphaltene molecule and @xmath73 is the total number of aggregates containing at least one asphaltene molecule , as in eq . [ eq : padef ] . we also know that : @xmath74 where @xmath75 is again the total number of asphaltene molecules and @xmath73 is the total number of aggregates containing at least one asphaltene molecule . it is easy to show that @xmath76 and finally that @xmath77    to obtain the dynamics predicted by the master equation  eq .   with the specific boundary condition  eq . , a numerical implementation of the scheme was carried out . in the numerical implementation the time step @xmath78 is a hundred times smaller than the inverse of the detachment rate constant @xmath79 . an aggregate can attach to a single molecule with the probability @xmath80 and release a single molecule with the probability @xmath81 . an aggregate of size @xmath69 can attach to a molecule but can not release one . the total number of molecules is kept constant in the numerical implementation of the master equation . the total number of molecules is chosen to be @xmath82 to reduce the statistical noise . initially , in the numerical implementation of the master equation , all the aggregates are of size @xmath69 . the fraction @xmath64 of free asphaltene molecules can be obtained through this numerical implementation . it depends _ a priori _ on the two rate constants @xmath43 and @xmath44 . however , the ratio @xmath83 is known from the stationary state result . it leaves us with one dynamical parameter , say @xmath43 , to fit . the initial state in the molecular dynamics simulations is not as well defined as in the master equation approach . to reach the desired density , a first md simulation where the system is compressed is performed . during the compression period , the asphaltene molecules begin to aggregate , but the data can not be recorded . the data are recorded just after the compression period , in a state where small asphaltene aggregates are already formed . the fraction of free asphaltene molecules in that state is @xmath84 . it is a second fitting parameter . the curve predicted by the numerical implementation of the master equation was shifted in time so that time @xmath85 corresponds to @xmath84 as in the md simulations . figure  [ fig : dynamics ] shows the fraction of aggregated asphaltene molecules versus time in the molecular dynamics simulations and in the master equation approach . both results agree well , indicating that the master equation approach correctly describe the aggregation process of asphaltene molecules . the value of the dynamical parameter @xmath43 is found to be @xmath86 s@xmath6 . it is worth mentioning that the inverse rate constant @xmath87 s is much larger than the upper limit @xmath88 s , needed for an asphaltene molecule to diffuse of a distance equals to the average distance @xmath89 nm between the centers of mass of @xmath90 aromatic molecules in a homogeneous system of volume @xmath91 nm@xmath92 minus the intermolecular distance @xmath93 @xmath8 in an aggregate . this characteristic time is evaluated using the diffusion coefficient of a single asphaltene molecule in the docosane solvent : @xmath94 m@xmath95.s@xmath6 and the formula @xmath96 for three - dimensional diffusion . this time is an upper limit because only the distance between centers of mass is considered , whereas molecules are extended in space and can be close to each other even if the distance between their centers of mass is larger than @xmath97 . we can conclude from that fact that the nanoaggregation process is not limited by diffusion .    in fig . [ fig : dynamics ] , there seems to be a discrepancy between the master equation approach and the md results at long times . this could indicate the existence of another aggregation process , taking place at a longer time scale . checking carefully the existence of this second process requires further investigation .    the master equation approach related to the aggregation reaction  eq .   is able to reproduce both the stationary behavior and the dynamics of the asphaltene aggregation process at intermediate time scales . it indicates that the assumptions made to derive the master equation  eq . are relevant . the master equation approach gives a dynamical framework to interpret the monoexponential distribution observed in md . even if the same master equation has already been used in the general field of clustering processes  [ ] , it is the first time , to our knowledge , that it is applied to asphaltene aggregation in molecular dynamics simulations . section  [ sec : mastereq ] provided a dynamical interpretation of the monoexponential behavior of the probability of having a nanoaggregate containing @xmath18 asphaltene molecules based on a master equation description on a simple aggregation reaction . the aim of the present section is to describe a statistical mechanics model based on the same aggregation scheme and consistent with the monoexponential behavior . in this section , we will also look into possible thermodynamical interpretations of the biexponential behavior , obtained when all aromatic molecules are counted .      to formulate a simple statistical mechanics model linked to the aggregation reaction  eq .  , some assumptions are needed . we make the following common assumptions :    1 . the nanoaggregates are linear . 2 .   the system is homogeneous and isotropic on average and the nanoaggregates are rigid . 3 .   the aggregation process occurs at a characteristic time much larger than the equilibration of the pressure , temperature and solvent molecules . the system is dilute enough to consider no interaction between the nanoaggregates except through the aggregation reaction . in other words , the solution of nanoaggregates of different sizes in a solvent is ideal . the free energy of a nanoaggregate depends linearly on its size . the idea of considering a mixture of free molecules and dimers as an ideal mixture with an aggregation reaction occurring at a larger time scale dates back to the beginning of the 20@xmath98 century  [ ] as is nicely explained in the recent review  [ ] . the set of assumptions was then completed to include the case of linear aggregates of any size and widely used to describe rodlike micelles  [ ] , linear polymer chains  [ ] , and discotic liquid crystals  [ ] . we propose in this section our own version of the derivation applied to the case of linear asphaltene aggregates and obtain an analytical formula for the probability of having an aggregate of a given size . the analogy between linear assembly in discotic liquid crystals and asphaltene aggregates was also made in ref  [ ] , which derives very similar equations to the ones presented in this section , but do not present a successful comparison with their coarse - grained simulations of asphaltene - resin nanoaggregates at high density . the assumption of having an ideal mixture can be relaxed to take into account more complicated interactions between the aggregates such as excluded volume interactions  [ ] , the gain in entropy when the chain breaks  [ ] and interactions due to the flexibility of the aggregates [ ] , to cite only a few . for the sake of simplicity , these interactions are neglected in the case of linear asphaltene aggregates . the validity of the assumptions made in this section will be discussed in sec . [ sec : assumption ] . concurrently to these approaches expressing the probability of having an aggregate of a given size , another theory of aggregation based on statistical thermodynamics was developed by wertheim  [ ] . this theory considers molecules of a reference liquid , typically a lennard - jones fluid , with a finite number of binding sites . the potential modeling the interaction between the binding sites is designed to take into account steric hindrance . the theory then counts the number of molecules with no bond on a given binding site instead of the number of aggregates of a given size . it is a clever way to count physically meaningful graphs . the outcome of the theory is an expression for the equilibrium pressure and the concentration of free molecules versus the composition of the system . these predictions were successfully checked numerically for strongly associating fluid  [ ] . however , this theory does not provide an expression for the probability of having an aggregate of a given size . for this reason , it is not considered in further details in this section but is mentioned for the sake of completeness .    in the framework of the common assumptions stated above , the free energy of the system can be written as  [ ] : @xmath99 where @xmath100 is the number of nanoaggregates containing @xmath18 asphaltene molecules , @xmath101 the volume of the system , @xmath102 the effective free energy of an aggregate , @xmath10 the boltzmann constant and @xmath103 the temperature . the first term in eq .   corresponds to the free energy of an ideal mixture of ideal gas and the second term corresponds to the energy of the aggregates . the effective free energy @xmath102 is to be understood as the energy of an asphaltene aggregate when the degrees of freedom due to solvent molecules have been integrated out . the expression eq .   for the free energy of the system holds for each value of the number @xmath100 of aggregates of a given size , because a state of local equilibrium acting on entropy , pressure and solvent molecules is assumed for each step along the aggregation reactions . the total equilibrium of the system depends now only on the equilibrium of the aggregation reactions  . the condition for chemical equilibrium of each aggregation reaction @xmath104 is @xmath105 where @xmath106 is the chemical potential of an asphaltene aggregate of size @xmath18 . using the definition of the chemical potential @xmath107 , we can obtain after some calculations the well - known mass action law : @xmath108 where @xmath109 is the equilibrium constant of the reaction @xmath104 . the equilibrium constant is given by : @xmath110 we now make use of assumption 5 , stating that the free energy of a nanoaggregate depends linearly on its size . this assumption amounts to write the effective free energy @xmath102 as : @xmath111 where @xmath112 is the free energy of a single asphaltene molecule and @xmath113 is the effective free energy between two asphaltene molecules . using this form in the expression of the equilibrium constant  eq . gives : @xmath114 thus , the equilibrium constant does not depend on the size @xmath18 of the considered nanoaggregate under the set of assumptions considered . by induction , it is now easy to show from  eq . that @xmath115 to compare the model to the simulations results it is more useful to obtain the probability @xmath31 of having a nanoaggregate of size @xmath18 , which is defined in eq.[eq : padef ] @xmath116 where @xmath117 is the total number of asphaltene nanoaggregates . to express @xmath31 , we make use of the conservation of the total number @xmath75 of asphaltene molecules : @xmath118 having this in mind , one can show ( see appendix  [ ap : paramp ] ) that @xmath119 where @xmath120 and @xmath121 we have now recovered the exponential distribution observed in the md simulations and fig .  [ fig : pdfall ] ( b ) . the md simulations provide a value of the parameter @xmath58 characterizing the exponential distribution : @xmath122 . according to the thermodynamical interpretation  eq . , it leads to the value @xmath123 the simple statistical mechanics model provides an interpretation for the parameter @xmath58 in terms of an effective free energy @xmath113 between two asphaltene molecules . within the assumptions of this simple model , the effective free energy is the interaction energy between two asphaltene molecules when the degrees of freedom related to solvent molecules are integrated out . the physical meaning of free energy @xmath113 depends on its definition in the framework of the assumptions made here , but also on the validity of these assumptions . this will be discussed in the next section . the validity of each assumption and its consequences on the physical meaning of the free energy @xmath113 will now be listed .    1 . _ linear dependence of @xmath102 on the size @xmath18 of the nanoaggregate . _ + it is very easy to check that the potential energy , and not the free energy , of a nanoaggregate in vacuum depends linearly on its size . the potential energy of a linear nanoaggregate of size @xmath18 is plotted versus @xmath18 for nanoaggregates in vacuum in fig . [ fig : potenergya ] . the origin of the energy is set arbitrarily to zero in this figure . the nanoaggregate used to plot this figure is a linear nanoaggregate containing five asphaltene molecules found in one simulation . to find the potential energy of an aggregate of size @xmath124 , only the first @xmath18 asphaltene molecules in the aggregate were considered in vacuum . this figure clearly shows that the potential energy of a nanoaggregate in vacuum is linear and the slope @xmath125 corresponding to the interaction energy between two asphaltene molecules is equal to @xmath126 in practical terms , it means that when an asphaltene molecule is added to a nanoaggregate of size @xmath18 , this asphaltene molecule interacts with the energy @xmath127 @xmath128 with the molecule at the end of the nanoaggregate but do not interact with the other ones , which are further away . + the fact that the potential energy of a nanoaggregate in vacuum depends linearly on its size is a good indication that it might also be the case for the effective free energy @xmath102 of a nanoaggregate in a solvent . in that case , the effective free energy @xmath113 corresponding to the slope of @xmath102 versus the size @xmath18 of the nanoaggregate is well defined . + we remind the reader that the potential energy of a nanoaggregate in vacuum and the effective free energy @xmath102 are not identical and the slopes @xmath113 and @xmath125 are indeed quite different . within the framework of the statistical model presented in sec . [ sec : thermomodel ] , the two causes of this difference are : the interaction between the asphaltene molecules and the solvent and the entropic effects due to the effective integration of the degrees of freedom of the solvent . to quantify the effective free energy @xmath113 properly , a possibility is to set up umbrella sampling simulations controlling the distance between the center of mass of two asphaltene molecules in a bath of docosane molecules at the same temperature and pressure as the one used in the md simulations . the free energy of this system can be derived in terms of the distance between the two asphaltene molecules . the effective free energy @xmath113 would then be the difference between the free energy of such a system when the two asphaltene molecules are far away and the free energy of the system when the two asphaltene molecules are aligned and close . implementing umbrella sampling simulations is beyond the scope of this paper . moreover , the free energy @xmath113 found with this method might be slightly different from the value @xmath129 @xmath128 , found using the statistical mechanics model of sec . [ sec : thermomodel ] , because of the other assumptions made to derive this model . _ asphaltene nanoaggregates . _ + to derive the statistical mechanics model in sec . [ sec : thermomodel ] , we only considered the aggregation of asphaltene molecules , whereas resin and resinous oil molecules are also part of the nanoaggregates . the degrees of freedom related to the position of resin and resinous oil molecules should be integrated in an effective way to obtain the free energy @xmath113 , just as it was done for the degrees of freedom related to the solvent . the fact that the monoexponential behavior is valid for the probability of having a nanoaggregate containing @xmath18 asphaltene molecules means that this effective integration can be done . linear nanoaggregates . _ + the nanoaggregates considered in the simulations are only some linear portions of bigger branched nanoaggregates . the existence of the branches could also produce some degrees of freedom to be integrated to obtain the effective free energy @xmath113 . 4 .   _ dilute limit . _ + the assumption stating that the solution of nanoaggregates is ideal neglects the interaction between the nanoaggregates . some of these interactions , for example between the nanoaggregates @xmath41 , @xmath130 and @xmath50 , are later taken into account through the aggregation reactions   corresponding to the asphaltene aggregation . but in the md simulations , the asphaltene molecules interact in a more complicated way . they can not for example be too close to each other , due to short range repulsion . this could be taken into account through excluded volume interactions . it is known  [ , ] , that taking into account excluded volume interactions preserves the monoexponential behavior . consequently , excluded volume interactions might play a role in the value of the effective free energy @xmath113 . _ time scale of the aggregation reactions . _ + to derive the statistical mechanics model , we assumed that the aggregation reactions occurred on a time scale much larger than the time scale associated with the equilibrium of pressure , entropy , and solvent molecules for a given number of each nanoaggregate . it is very difficult to predict the effect of the relaxation of this fundamental assumption . one way to check the assumption , however , is to compare once again the characteristic time of the aggregation reaction @xmath87 s and the upper limit @xmath88 s for the diffusion of an asphaltene molecule . a factor @xmath131 exists between the two characteristic times which should be enough to ensure the validity of the assumption . _ rigidity of the nanoaggregates and isotropy of the system . _ + [ sec : relaxrigidity ] the assumptions of rigid nanoaggregates and isotropic system can be considered together , because if they are both relaxed , they lead to a new term in the free energy of the system  [ ] . this term depends on a persistence length @xmath132 and reads as  [ , ] @xmath133^{1/2 } \nabla[\rho_n(\mathbf{u})]^{-1/2},\ ] ] where @xmath134 is the number density of asphaltene molecules part of a nanoaggregate of size @xmath18 oriented in the direction @xmath135 with respect to some reference direction and @xmath136 is the corresponding solid angle . this term is helpful to describe the nematic phase of liquid crystals , where long range order is seen . the addition of this term in the free energy expression leads to a biexponential behavior  [ ] , which is not observed for the probability of having a linear nanoaggregates with @xmath18 asphaltene molecules . consequently , the assumption of rigid nanoaggregates and isotropic system probably holds , at least in an effective way , for the probability of having a linear nanoaggregate with @xmath18 asphaltene molecules and does not participate in the value of the effective free energy @xmath113 . the addition of this flexible term to the free energy is potentially of interest to explain the biexponential behavior of the probability of having a nanoaggregate containing @xmath18 aromatic molecules .    to summarize , one can say that the statistical mechanics model developed in sec . [ sec : thermomodel ] provides a thermodynamical interpretation of the monoexponential behavior based on the effective free energy between two asphaltene molecules . this effective free energy should be understood as the interaction energy between two asphaltene molecules when all the degrees of freedom related to solvent molecules , resin and resinous oil molecules and branched nanoaggregates have been integrated out . the picture emerging from the study of the probability of having a linear nanoaggregate containing @xmath18 asphaltene molecules is that asphaltene bodies interact with an effective energy through the aggregation reaction  . this leads to a monoexponential behavior . we know that the probability of having a linear nanoaggregate with @xmath18 aromatic molecules has a different behavior . it is biexponential as shown in fig . [ fig : pdfall ] ( a ) . this section is devoted to identifying possible statistical mechanics explanations for this different behavior . the first explanation that comes to mind is that when all molecules are counted , different interaction energies are involved . each interaction energy taken separately would lead to a specific monoexponential behavior and the combinations of several interaction energies could lead to a bi- or multiexponential behavior .    to test this idea , we set up one - dimensional lattice monte - carlo simulations . a linear lattice of @xmath137 sites is created . each site can contain one asphaltene molecule , one resin molecule , one resinous oil molecule or nothing . one site can not contain two molecules . the total numbers of asphaltene molecules , resin molecules , and resinous oil molecules are constant . a monte - carlo move consists in exchanging the content of two sites providing that the content is different . this condition makes the molecules indiscernible . the system is initialised with the largest possible aggregate where all asphaltene molecules are next to each other , then comes resin molecules and then resinous oil molecules . ten million ( @xmath138 ) monte - carlo moves are realised using the metropolis algorithm . the algorithm converges quite quickly despite its elementary implementation . the potential energy of a nanoaggregate is calculated in the following way : when two asphaltene molecules are next to each other the interaction @xmath139 between two asphaltene molecules is added , when a resinous oil molecule is next to an asphaltene molecule or another resinous oil molecule the interaction energy @xmath140 is added , finally when a resin molecule is next to any other molecule the interaction energy @xmath141 is added . as two molecules can not be on the same site , effective excluded volume interactions are created . the assumptions underlying the establishment of the one - dimensional lattice monte - carlo simulations are similar to the ones made in sec . [ sec : thermomodel ] except for excluded volume interactions :    1 . the nanoaggregates can only be linear because the system is one - dimensional . 2 .   the nanoaggregates are rigid . the system is one dimensional , so isotropy is not a criterion . 3 .   the only energies involved are those related to the aggregation process . all potential energies related to interaction with and within solvent molecules are averaged out . no kinetic energy is involved . 4 .   the system is not dilute and excluded volume interactions are taken into account . 5 .   the energy of a pure asphaltene nanoaggregate depends linearly on its size . one of the main advantages of lattice monte - carlo simulations compared to the analytical approach is to take into account resin and resinous oil molecules and not only asphaltene molecules . we checked that in the dilute limit , when only asphaltene molecules are present in the monte - carlo simulations , the same monoexponential behavior with the same value for the parameter @xmath58 as the one predicted by the analytical approach is recovered .    in the lattice monte - carlo simulations , it is possible to obtain a biexponential behavior for the probability of having @xmath18 molecules in a nanoaggregate as can be seen in fig . [ fig:1dlatmc]a . the biexponential behavior is characterized by the presence of two straight lines with different slopes in a log - lin scale . the biexponential behavior occurs in the lattice monte - carlo simulations when the interaction energies @xmath139 between asphaltene molecules on the one hand and @xmath140 and @xmath141 with resin and resinous oil molecules on the other hand are substantially different . for example , the choices @xmath142 @xmath128 , @xmath143 @xmath128 and @xmath144 @xmath145 at the same temperature and same volume as the md simulations give a biexponential behavior very close to the one obtained in md in fig . [ fig : pdfall ] ( a ) . it is shown in fig . [ fig:1dlatmc ] ( a ) . however , the biexponential behavior does not have the same causes as in the md simulations . in the monte - carlo simulations , the biexponential behavior is due to the fact that pure resin or resinous oil nanoaggregates and pure asphaltene nanoaggregates are formed . there are very few mixed nanoaggregates . thus , the first slope in the biexponential behavior is due to resin and resinous oil nanoaggregates and the second slope is due to asphaltene nanoaggregates . one consequence of this fact is that the probability of having a nanoaggregate with @xmath18 resin or resinous oil molecules has a monoexponential behavior with a slope very close to the first slope of the biexponential behavior . in the same way , the probability of having a nanoaggregate with @xmath18 asphaltene molecules has a monoexponential behavior with a slope very close to the second slope of the biexponential behavior . this can be seen in fig . [ fig:1dlatmc](a ) . on the contrary , in the md simulations , the two slopes of the biexponential behavior do not correspond to two different slopes in two different monoexponential behaviors . this can be seen in fig . [ fig:1dlatmc](b ) . the existence of mixed nanoaggregates can be checked directly in the md simulations . [ fig:1dlatmc](c ) displays the ratios @xmath146 , @xmath147 , and @xmath148 of asphaltene , resin , and resinous oil molecules respectively versus the size of the aggregate . for a given molecule type @xmath149 , the ratio @xmath150 is defined as : @xmath151 where @xmath152 is the number of molecules of type @xmath149 in the aggregate and @xmath18 the total number of molecules in the aggregate . [ fig:1dlatmc](c ) shows that there is indeed a change in the nanoaggregate composition with their size . the ratio of asphaltene molecules increases versus the size of the aggregates until it reaches an approximately constant value for aggregates of size @xmath153 . at the same time , the ratio of resin molecules decreases versus the size of the aggregates and reaches an approximately constant value for aggregates of size @xmath153 , while the ratio of resinous oil molecules is roughly constant . however , for any size the nanoaggregates contain all molecule types . we can conclude that , the biexponential behavior in md is probably not only due to the difference in effective interaction energies involved . another possible explanation of the biexponential behavior obtained for the probability of having a nanoaggregate containing @xmath18 aromatic molecules , is that it derives from the nanoaggregate flexibility .    as mentioned in sec . [ sec : assumption ] , if the aggregates are flexible and the system is not isotropic , an extra term in the free energy of the system should be considered ( see  eq .  ) . this leads to a biexponential behavior . the existence of the biexponential behavior can be explained qualitatively in this way : bent nanoaggregates , less stable than straight ones , tend to be smaller and are responsible for the first slope of the biexponential behavior ; on the contrary straight nanoaggregates are larger and give rise to the second slope  [ ] . this additional term is of course well - suited to explain the biexponential behavior observed for the probability of having a nanoaggregate with @xmath18 aromatic molecules . however , we are not convinced that it is the main reason explaining the biexponential behavior , but there are indications that it might play a role . first , some nanoaggregates are bent . a picture of a bent nanoaggregate can be seen in fig . [ fig : flexchain](a ) . to quantify the variation of the rigidity of an aggregate with the distance inside this aggregate , we computed the orientation correlation function @xmath154 , where @xmath155 is the unit vector normal to the molecule @xmath156 in a given linear aggregate , @xmath157 stands for the @xmath158 neighbor of molecule @xmath156 in the same linear aggregate , and @xmath159 corresponds to the average over different nanoaggregates and over time . [ fig : oriencorrel ] shows the variation of the orientation correlation function @xmath154 with @xmath160 . the error bars corresponding to the standard deviation over eight independent simulations are very large for this plot and are not shown for the sake of visibility . considering the large errors , only a qualitative discussion on the average trend is possible . we can see that there is , on average , an initial decrease of the orientation correlation with the number of neighbors . this trend shows that nanoaggregates are not perfectly rigid . for larger distances and consequently larger nanoaggregates , the orientation correlation seems to plateau around the value @xmath161 . this last trend matches the qualitative picture of small aggregates being more bent than larger aggregates . the change in the nanoaggregate composition observed in fig . [ fig:1dlatmc](c ) could explain the change in the nanoaggregate rigidity with the nanoaggregate size : aggregates containing many asphaltene molecules tend to be longer and more rigid . second , the system is not strictly isotropic . the isotropy was quantified using the nematic order parameter @xmath162  [ ] . to define the nematic order parameter , the following order tensor needs to be defined for each nanoaggregate : @xmath163 where @xmath18 is the number of molecules in the considered linear nanoaggregate , @xmath164 and @xmath165 cartesian coordinates , @xmath156 the index of a molecule inside the aggregate , @xmath155 the unit vector normal to molecule @xmath156 , and @xmath166 the kronecker delta . the nematic order parameter @xmath162 is the largest eigenvalue of the order tensor @xmath167 averaged over different nanoaggregates and time . the nematic order tensor is equal to @xmath69 in a system where all aggregates are perfectly aligned and to @xmath168 in a perfectly isotropic system . in the md simulations , the nematic order parameter is equal to @xmath169 where the error is the standard deviation corresponding to eight independent simulations . the value of the nematic order parameter indicates a system closed to being isotropic but not quite . it can be due to the fact that some linear nanoaggregates are branches of bigger aggregates . they can be connected through asphaltene heads and the angle between an asphaltene head and an asphaltene body is fixed by the dihedral potential and does not vary this much from one asphaltene molecule to another  [ ] . it can also be due to steric hindrance : long nanoaggregates can not interpenetrate each other and consequently tend to align . the facts that nanoaggregates are flexible and that the system is not perfectly isotropic are consequently a plausible explanation for the biexponential behavior of the probability of having a nanoaggregate containing a given number of aromatic molecules . ( color online ) picture of a bent nanoaggregate , obtained in molecular dynamics . the color code is the same as in fig . [ fig : linearaggregate ] . ]    the fact that the nanoaggregates are not strictly rigid and that the system is not strictly isotropic was also valid when we looked at the probability of having a nanoaggregate containing @xmath18 asphaltene molecules . however , the fact that this probability has a monoexponential behavior means that the additional term  eq .   can be neglected in this case . it could be because the effective persistence length is longer when only asphaltene molecules are considered . our bitumen model is quite simplified compared to a real bitumen . in this section , we will compare our results to other md results and to experimental results . we will also discuss the new idea that such a simplified model can bring to the field of asphaltene nanoaggregation and the specific aspects of asphaltene nanoaggregation which are neglected in this simplified model . our model relies on the fact that the @xmath16-stacking interaction is the origin of the nanoaggregate formation . in our md simulations , the @xmath16-stacking is modeled through the lennard - jones potential . recent molecular dynamics simulations of asphaltene molecules in toluene  [ ] reported that the @xmath16-stacking interaction is indeed the most important one to explain the nanoaggregate formation , which justifies our model .    in our md simulations , the average number of asphaltene molecules in a linear nanoaggregate is given by the exponential distribution . we can compare the average and standard deviation predicted by the exponential distribution to experimental results . to do that , we considered a renormalized distribution of the nanoaggregate size , where the aggregate of size @xmath69 are left out . in other words , free asphaltene molecules are not considered here as a nanoaggregate . the corresponding probability reads : @xmath170 in this case , the average number of asphaltene molecules in linear nanoaggregates is @xmath171 the average number of aromatic molecules in linear aggregates can be computed directly from the md results and is @xmath172 . these results are in agreement with the general consensus that nanoaggregates contain less than @xmath173 molecules  [ ] . our simulation also offers a precise quantification of the polydispersity of the nanoaggregate size . the standard deviation of the distribution given in eq .   is : @xmath174 again , the standard deviation on the number of aromatic molecules in linear nanoaggregates can be computed from the md results and is @xmath175 . these results are compatible with previous md results reporting for example that : for asphaltene nanoaggregates in toluene , the number of asphaltene molecules in the largest aggregate varies from 2 to 18 depending on the asphaltene structure  [ ] ; for asphaltene nanoaggregates in vacuum , the number of asphaltene molecules in any aggregate varies from 1 to 5 and the precise distribution depends both on temperature and the asphaltene structure  [ ] . on the experimental side , recent laser - based mass spectrometry experiments  [ ] were able to obtain not only the average size of a nanoaggregate but also the polydispersity . the aggregation numbers are found to range roughly from @xmath176 to @xmath177 or from @xmath177 to @xmath178 depending on the bitumen chemical composition . it is also compatible with our results . moreover , our simulations bring the idea , already suggested in ref .  [ ] without a successful match to simulation results at high density , that the simple monoexponential distribution is a typical distribution of the number of asphaltene molecules in linear aggregates . this distribution is consequently a good basis to model more complicated cases . one of these more complicated cases is considered here : when all aromatic molecules are counted in the nanoaggregates the probability of having a nanoaggregate of a given size becomes biexponential . the biexponential distribution can be seen as a modification of the exponential distribution when the flexibity of the aggregates and the anisotropy of the system are taken into account . many more specific aspects relevant for bitumen science can be considered and are not treated here . for example , the effect of having branched nanoaggregates on the distribution of the nanoaggregate size could be addressed . furthermore , it seems that the presence of long aliphatic chains in the asphaltene molecules modifies the typical structure of a nanoaggregate  [ ] enhancing t - shaped geometry ( @xmath16-@xmath179 interaction ) and offset @xmath16-stacked geometry ( @xmath179-@xmath179 interaction ) compared to @xmath16-@xmath16 geometry  [ ] . it would be very interesting to consider the effect of adding asphaltene and resin molecules with long aliphatic chains in our simulations on the shape of the nanoaggregate size distribution compared to the simple exponential distribution . in conclusion , we have shown that the aggregation reactions   give a good description of the stability of linear asphaltene nanoaggregates as observed in molecular dynamics ( md ) simulations . from these aggregation reactions , we derived a master equation able to reproduce the monoexponential behavior of the stationary probability of having a nanoaggregate containing a given number of asphaltene molecules in md . the parameter of the monoexponential behavior is interpreted as the ratio between the attachment and detachment rate constants of a single asphaltene molecule to a nanoaggregate . the master equation approach is also able to reproduce the aggregation dynamics . from the same aggregation reactions , we also derived a statistical mechanics model giving a thermodynamics interpretation of the monoexponential behavior . the main parameter is then the effective free energy between two asphaltene molecules , when the degrees of freedom corresponding to solvent molecules , resin , and resinous oil molecules and branched nanoaggregates are integrated out . finally , a possible thermodynamic explanation for the biexponential behavior , observed for the stationary probability of having a nanoaggregate of @xmath18 aromatic molecules in md , is the flexibility of these nanoaggregates .    to continue this work on bitumen nanoaggregate two directions are possible and equally interesting . a first direction is to consider a simpler system without resin and resinous oil and even without any possibility of branching . then , the integration of the degrees of freedom related to the solvent molecules could be done and the effective energy could be evaluated in an independent way . this direction would lead to a quantitative understanding of the effective free energy between two asphaltene molecules . the second and opposite direction is to add more molecule types to resemble a real bitumen . for example asphaltene molecules without a head and with long alkyl chains could be added and the consequences of this addition on the probability of having a nanoaggregate of a given size investigated . interesting md simulations have been performed recently using many molecules types  [ ] and could serve as an inspiration . another interesting route is to quantify the evolution of the nanoaggregate size distribution with the composition of the bitumen mixture . this work is sponsored by the grants 1337 - 00073b and 1335 - 00762b of the danish council for independent research @xmath180 technology and production science . it is in continuation of the cooee project ( co@xmath4 emission reduction by exploitation of rolling resistance modeling of pavements ) , sponsored by the danish council for strategic research . the centre for viscous liquids dynamics `` glass and time '' is supported by the danish national research foundation s grant dnrf61 . the probability @xmath31 of having a nanoaggregate with @xmath18 asphaltene molecules is defined as : @xmath181 where @xmath182 is the number of aggregates with @xmath18 asphaltene molecules and @xmath117 the total number of asphaltene nanoaggregates . according to eq . , @xmath183 so that one can express the total number of asphaltene nanoaggregates : @xmath184 to obtain the last expression eq . , we assume that the number of asphaltene molecules is very large , so that the sum goes to infinity and that the ratio @xmath185 is smaller than 1 , _ i. e. _ the volume @xmath101 is big enough for the aggregates to develop given the total number of asphaltene molecules and the equilibrium constant . inserting the expression eq . of the total number of asphaltene aggregates back into eq .   gives : @xmath186 it is a geometrical law of the form @xmath187 and the parameter @xmath58 can be identified as : @xmath188 an expression for @xmath72 , the number of asphaltene aggregates of size @xmath69 , is given by the conservation law : @xmath189 , where @xmath75 is the total number of asphaltene molecules . it leads to : @xmath190 the last expression eq . is a quadratic equation , whose solutions are : @xmath191 where @xmath192 the solution @xmath193 is the physical one , because in the limit where there is no interaction between the asphaltene molecules and even a large repulsion , all molecules should be in aggregates of size @xmath69 . in mathematical terms , it gives : @xmath194 only @xmath193 satisfies this last equation . it leads finally to : @xmath195 which is the same result as eq .  . another useful expression is the one giving the equilibrium constant @xmath196 in terms of the parameter @xmath58 : @xmath197    finally , it is noteworthy that the minimization of the free energy of the system as given in eq . subject to the conservation condition eq .   leads to the same result . the minimization can be done using a lagrange multiplier to guarantee the conservation of the total number of asphaltene molecules . in this calculation , the free energy @xmath102 of an aggregate should be expressed right away as @xmath198 , where the arbitrary origin of the energy @xmath199 should be equal to @xmath200 . in this case , @xmath113 represents the energetic penalty of having a free end . j. f. branthaver , j. c. petersen , r. e. robertson , j. j. duvall , s. s. kim , p. m. harnsberger , t. mill , e. k. ensley , f. a. barbour , and j. f. schabron , tech . shrp - a-368 , strategic highway research program ( 1993 ) . o. c. mullins , h.  sabbah , j.  eyssautier , a.  e. pomerantz , l.  barr , a.  b. andrews , y.  ruiz - morales , f.  mostowfi , r.  mcfarlane , l.  goual , r.  lepkowicz , t.  cooper , j.  orbulescu , r.  m. leblanc , j.  edwards , and r.  n. zare , energy & fuels * 26 * , 3986 ( 2012 ) .
Abstract:

Answer: asphaltene and smaller aromatic molecules tend to form linear nanoaggregates in bitumen . over the years bitumen undergoes chemical aging and during this process , the size of the nanoaggregate increases . this increase is associated with an increase in viscosity and brittleness of the bitumen , eventually leading to road deterioration . this paper focuses on understanding the mechanisms behind nanoaggregate size and stability . we used molecular dynamics simulations to quantify the probability of having a nanoaggregate of a given size in the stationary regime . to model this complicated behavior , we chose first to consider the simple case where only asphaltene molecules are counted in a nanoaggregate . we used a master equation approach and a related statistical mechanics model . the linear asphaltene nanoaggregates behave as a rigid linear chain . the most complicated case where all aromatic molecules are counted in a nanoaggregate is then discussed . the linear aggregates where all aromatic molecules are counted seem to behave as a flexible linear chain .    corresponding author : c. a. lemarchand , e - mail : clairel@ruc.dk