the description of strong electron correlation effects is still nowadays a challenging problem for both density - functional theory ( dft ) and wavefunction theory ( wft ) communities . even though kohn
 sham dft ( ks - dft ) is in principle exact , standard approximate exchange
correlation functionals usually do not enable an adequate description of multi - configurational systems .
on the other hand , the `` gold standard '' single - reference _ coupled cluster _ ( cc ) method should also be able to model near - degeneracies but high - order excitations would then be necessary . in practice multi - reference ( mr )
perturbation theories such as the _ second order complete active space perturbation theory _ ( caspt2 )
@xcite or the _ n - electron valence state perturbation theory _ ( nevpt2 )  @xcite are usually employed . despite their success
, these methods may suffer from the perturbative description of the short - range dynamical correlation which is sometimes not accurate enough . for these reasons
multi - configurational extensions of both dft and cc have been investigated for many years and despite significant efforts and encouraging results , it is still unclear how the incorporation of a _ complete active space self - consistent field _ ( casscf ) into dft and cc should be performed . a recent review article by bartlett and coworkers  @xcite gives an extensive picture of the developments in mrcc .
we shall focus here on dft and restrict ourselves to the description of the ground state .
+ the rigorous combination of casscf with dft is already difficult at the formal level due to the infamous double counting problem .
specific functionals must indeed be developed for complementing the casscf energy that already contains some correlation effects  @xcite .
the latter are usually referred to as static correlation effects . only the complementary correlation , known as dynamical correlation , should therefore be assigned to the density functional .
turning such a scheme into a practical computational method is not an easy task since approximations used for standard functionals can not be applied straightfowardly in this context .
one way to overcome the double counting problem consists in separating correlation effects in the coordinate space , as initially proposed by savin  @xcite .
the resulting range - separated dft methods combine rigorously long - range wft with short - range dft .
even though range separation allows for a multi - configurational description of the electron density , it can not completely isolate static correlation from dynamical correlation simply because the former is usually not a purely long - range effect , even in a dissociated molecule .
this point will be discussed further in the following .
+ let us mention that savin earlier combined configuration interaction ( ci ) with dft by separating correlation effects in the natural orbital space  @xcite
. orbitals with occupation numbers larger than a given threshold @xmath0 were correlated explicitly , at the ci level , while the remaining correlation energy was described with a complementary @xmath0-dependent density functional .
more recently , gutl and savin  @xcite proposed alternative ci - dft schemes where the correlation energy is still split in the orbital space but through the introduction of gap shift or cutoff parameters . in all these hybrid ci - dft schemes , a ci energy is complemented by a functional of the electron density . in this work
we propose to revisit at the formal level the separation of correlation effects in the orbital space . as a major difference with the approaches discussed previously ,
orbitals occupation will be used as basic variable rather than the electron density .
the idea originates from condensed matter physics where strongly correlated electrons are usually described with model hamiltonians ( such as the hubbard hamiltonian ) rather than the true physical one  @xcite . as discussed
further in the following , adapting such an approach to quantum chemistry is appealing since static correlation is usually defined in the orbital space .
molecular orbitals can then be viewed as sites with an occupation to be determined .
as shown in this work , a wft description of static correlation can be combined rigorously with an orbital - occupation functional modeling of dynamical correlation .
a potential drawback of such a scheme lies in the fact that , in contrast to range - separated dft , the complementary functional loses its universality since it depends on the molecular orbitals .
the theory is derived for a particular choice of orbitals , namely those obtained by diagonalizing the non - interacting hamiltonian ( kinetic and nuclear potential energy operators only are considered ) .
its generalization is left for future work .
+ the paper is organized as follows : after a short introduction to the double counting problem in multi - configuration dft ( sec .
[ sec:2blecounting ] ) , the exact multi - determinant extension of standard dft based on range separation is briefly reviewed ( sec .
[ subsec : srdft ] ) and its application to the dissociated h@xmath1 molecule is discussed in sec .
[ subsec : h2_min_basis ] .
we then leave the framework of dft by using orbitals occupation as basic variable rather than the electron density .
this approach is first introduced for model hamiltonians such as the hubbard hamiltonian ( sec .
[ sec : soft ] ) . by isolating an impurity site and separating its on - site repulsion from the interactions on the remaining sites ( referred to as the bath ) , an exact _ site - occupation embedding theory _ ( soet ) is derived and compared with the _ density matrix embedding theory _ ( dmet ) of knizia and chan  @xcite .
illustrative results are then given for the simple two - site hubbard model in sec .
[ sec : two_site_hubb ] . finally , we show in sec .  [
sec : soft_in_qc ] that an exact combination of ci with orbital occupation functional theory is obtained when applying soet to a quantum chemical hamiltonian .
conclusions and perspectives are given in sec .
[ sec : conclusions ] .
in wft the exact ground - state energy of an electronic system can be obtained variationally as follows , @xmath2 where @xmath3 is the kinetic energy operator , @xmath4 denotes the two - electron repulsion operator and @xmath5 is the nuclear potential operator . in a regular casscf calculation , the minimization in eq .
( [ eq : exactvp_wft ] ) is restricted to linear combinations of slater determinants that belong to a given _ complete active space _ ( cas ) .
the orbitals are also optimized variationally .
the cas is obtained by distributing a given number of active electrons in selected active orbitals that become consequently partially occupied .
the selection of active electrons and orbitals is usually based on chemical intuition .
doubly occupied orbitals are referred to as inactive orbitals .
note that the latter are not frozen in a casscf calculation .
the remaining orbitals , that are not occupied , are the virtuals .
for convenience we will denote @xmath6 the space of trial casscf wavefunctions .
the superscript @xmath7 specifies all the restrictions in the minimization ( number of active electrons and active orbitals for example ) .
the casscf energy can therefore be written as @xmath8 note that , in practice , a trial casscf wavefunction @xmath9 will be parametrized in second quantization as follows  @xcite , @xmath10 where @xmath11 is the basis of determinants for the cas .
the singlet excitation operator @xmath12 that allows for orbital rotations is defined as @xmath13 where @xmath14 and @xmath15 can be inactive , active or virtual orbitals .
consequently , the minimization in eq .
( [ casscf_ener_exp ] ) is performed over both @xmath16 and @xmath17 parameters .
+ the energy difference @xmath18 is usually referred to as dynamical correlation energy . the latter is routinely described within multi - reference perturbation theories such as caspt2 or nevpt2 .
a combined casscf - dft approach should ideally be exact and variational .
dynamical correlation would be described with a density functional @xmath19 $ ] that would complement the casscf energy .
this can be formulated rigorously when considering the universal levy
lieb functional  @xcite , @xmath20 =   \underset{\psi\rightarrow n}{\rm min}\langle \psi\vert \hat{t}+{\hat{w}_{\rm ee}}\vert\psi\rangle,\end{aligned}\ ] ] where the minimization is performed over all wavefunctions with density @xmath21 , and the following partitioning  @xcite , @xmath22 = f^{m}[n]+\overline{e}^{{{\rm dyn}}}_{\rm c}[n ] .
\end{aligned}\ ] ] the ll functional associated with the casscf space @xmath6 is defined as @xmath23 =   \underset{\psi\in\mathcal{s}^{m}\rightarrow n}{\rm min}\langle \psi\vert \hat{t}+{\hat{w}_{\rm ee}}\vert\psi\rangle,\end{aligned}\ ] ] where the minimization is restricted to casscf wavefunctions in @xmath6 with density @xmath21 . according to the hohenberg
kohn ( hk ) theorem  @xcite , the exact ground - state energy can be obtained variationally as follows , @xmath24 + ( v_{\rm ne}\vert
n ) \big\},\end{aligned}\ ] ] where the minimization is performed over electron densities @xmath25 that integrate to a fixed number @xmath26 of electrons . the notation @xmath27 has been used .
combining eq .
( [ ll_cas_fun ] ) with eq .
( [ eq : hk2_theo ] ) leads to @xmath28+\overline{e}^{{{\rm dyn}}}_{\rm c}[n]+ ( v_{\rm ne}\vert n ) \big\}.\end{aligned}\ ] ] consequently , for any trial casscf wavefunction @xmath9 in @xmath6 , @xmath29}}&\geq & f^{{{m}}}[n_{\psi}]+{\overline{e}^{{{\rm dyn}}}_{\rm c}[n_{\psi}{]}}+(v_{\rm ne}\vert n_\psi ) \nonumber\\ & \geq & e,\end{aligned}\ ] ] where @xmath30 is the electron density obtained from the trial casscf wavefunction @xmath9 .
thus we conclude that , provided that the exact ground - state density can be reproduced by a casscf wavefunction , the exact ground - state energy can be written as @xmath31 } } \big\}.\end{aligned}\ ] ] in order to turn eq .
( [ eq : var_casdft_ener ] ) into a practical computational method , approximate complementary correlation functionals should be developed .
obviously standard correlation functionals can not be used otherwise correlation effects would be double counted .
this double counting problem arises also when the size of the cas varies . in the exact theory
, the functional should vary with the cas so that correlation effects are transferred from the casscf to the dft part of the energy and the total energy remains constant and equal to the exact ground - state energy .
deriving cas - dependent density functionals is a complicated task simply because the electron density is defined in the coordinate space while a cas is defined in the orbital space . in order to rigorously overcome such difficulties ,
savin  @xcite proposed to separate correlation effects in the coordinate space , thus leading to the so - called range - separated dft methods .
this approach is briefly reviewed in sec .
[ sec : mc - srdft ] .
we then propose an alternative approach where correlation effects are separated in the orbital space , like in a regular casscf calculation .
in this context , the orbitals occupation will be used as basic variable rather than the electron density .
in the standard ks - dft scheme  @xcite , the universal ll functional in eq .
( [ eq : universal_ll_fun ] ) is decomposed into the non - interacting kinetic energy functional @xmath32= \underset{\psi\rightarrow n}{\rm min}\langle \psi\vert \hat{t}\vert\psi\rangle $ ] and the hartree - exchange - correlation ( hxc ) energy functional , @xmath33=t_{\rm s}[n]+e_{\rm hxc}[n].\end{aligned}\ ] ] consequently , the hk variational principle can be reformulated as follows , @xmath34 \big\},\end{aligned}\ ] ] where the minimizing ks determinant @xmath35 fulfils the self - consistent equation @xmath36\,\hat{n}({\bf r})\bigg ) \vert\phi^{\rm
ks}\rangle=\mathcal{e}^{\rm ks}\vert\phi^{\rm ks}\rangle.\end{aligned}\ ] ] therefore , within ks - dft , the two - electron interaction is fully described by a density functional . as originally shown by savin  @xcite
, it is in fact possible to describe only a part of the two - electron repulsion within dft and leave the remaining part to wft .
this is achieved by separating the interaction in the coordinate space into two complementary contributions .
the range separation based on the error function , @xmath37 has for example been used extensively in the last decade ( see ref .
@xcite and the references therein ) , but any separation like the simpler linear one  @xcite can be considered .
range separation , that is controlled by the @xmath38 parameter , is appealing as it enables to isolate the coulomb hole and assign it to a density functional while long - range correlation can be described in wft .
this hybrid range - separated wft - dft approach can be derived rigorously from the alternative partitioning of the ll functional , @xmath39= f^{\rm lr,\mu}[n]+\overline{e}^{\rm sr,\mu}_{\rm hxc}[n],\end{aligned}\ ] ] where the long - range ll functional equals @xmath40= \underset{\psi\rightarrow n}{\rm min}\langle \psi\vert \hat{t}+{\hat{w}^{\rm lr,\mu}_{\rm ee}}\vert\psi\rangle , \end{aligned}\ ] ] and @xmath41 $ ] is the complementary @xmath38-dependent short - range hxc density functional . let us stress that , for any @xmath38 value , this functional is universal since it depends on the electron density only .
combining eq .
( [ eq : ll_range_separation ] ) with eq .
( [ eq : hk2_theo ] ) leads to the exact ground - state energy expression @xmath42 \right\},\end{aligned}\ ] ] where the minimizing wavefunction @xmath43 fulfils the self - consistent equation @xmath44 \hat{n}(\mathbf{r}).\end{aligned}\ ] ] note that ks - dft and pure wft are recovered when @xmath45 and @xmath46 , respectively .
+ since the long - range interaction is now described explicitly , the auxiliary wavefunction @xmath43 , that has exactly the same density as the physical system , is multi - determinantal .
casscf can therefore be applied in this context in conjunction with local or semi - local short - range functionals  @xcite . even though encouraging results were obtained with such a range - separated casdft scheme , better functionals
are still needed for the method to be reliable  @xcite . as discussed further in sec .
[ subsec : h2_min_basis ] , one major problem is that static correlation is not a purely long - range correlation effect .
+ returning to the exact theory , an adiabatic connection ( ac )  @xcite expression for the short - range functional can be obtained from the auxiliary equations , @xmath47 and the density constraint @xmath48 for @xmath49 .
indeed , according to eqs .
( [ eq : ll_range_separation ] ) and ( [ eq : ll_longrange ] ) , @xmath50&=&f[n]-f^{\rm lr,\mu}[n ] \nonumber\\ & = & \int^{+\infty}_\mu \dfrac{\ddroit\mathcal{e}^\nu}{\ddroit \nu}\ddroit\nu\ + ( v^\mu - v^{+\infty}\vert n),\end{aligned}\ ] ] thus leading , according to the hellmann  feynman theorem , to @xmath51 & = & \int^{+\infty}_\mu \langle { \psi}^{\nu}\vert\partial \hat{w}^{\rm lr,\nu}_{\rm ee}/\partial \nu \vert{\psi}^{\nu } \rangle .
\end{aligned}\ ] ] we obtain for @xmath45 a range - separated ac expression for the conventional hxc functional , @xmath52 & = & \int^{+\infty}_0 \langle { \psi}^{\nu}\vert\partial \hat{w}^{\rm
lr,\nu}_{\rm ee}/\partial \nu \vert{\psi}^{\nu } \rangle .
\end{aligned}\ ] ] it is readily seen from eqs .
( [ eq : ac_srfun ] ) and ( [ eq : ac_rs_convhxc ] ) that conventional functionals can not be used straightfowardly in multi - determinant range - separated dft otherwise there would be double counting of long - range correlation effects through the wft treatment . by connecting the non - interacting ks system ( @xmath53 ) to the long - range interacting one ( @xmath54 ) as follows , @xmath55-t_{\rm s}[n ] \nonumber\\ & = & { e}^{\rm lr,\mu}_{\rm hxc}[n],\end{aligned}\ ] ] where @xmath56 $ ] is the purely long - range hxc functional , we obtain from eqs .
( [ eq : ac_srfun ] ) and ( [ eq : ac_rs_convhxc ] ) the following expression , @xmath50={e}_{\rm hxc}[n]-{e}^{\rm lr,\mu}_{\rm hxc}[n],\end{aligned}\ ] ] that has been used by toulouse  @xcite for developing approximate local and semi - local short - range functionals .
we should stress that all the formalism briefly reviewed in this section is in fact general and can be applied in a different context , for example in dft for model hamiltonians , as proposed in sec .
[ subsec : soet ] .
we consider in this section the h@xmath1 molecule in a slater minimal basis consisting of the @xmath57 and @xmath58 atomic orbitals localized on the left and right hydrogen atoms , respectively  @xcite .
the basis functions are identical with @xmath59 . for large bond distances the bonding and anti - bonding molecular orbitals are equal to @xmath60 and @xmath61 , respectively .
the ground state will then be written in the basis the two slater determinants @xmath62 and @xmath63 . since the latter differ by a double excitation , they are not coupled by one - electron operators .
therefore , when approaching the dissociation limit , the matrix representation of the auxiliary long - range hamiltonian in eq .
( [ timeindepsrdfteq ] ) reduces to @xmath64=\left [ \begin{array}{c c } { e}^\mu & k^{{\mu } } \\ { k^{{\mu } } } & { e}^\mu \end{array } \right],\end{aligned}\ ] ] where diagonal elements are identical since atomic orbitals do not overlap , and the coupling term equals @xmath65 while the second long - range two - electron integral on the right - hand side of eq .
( [ eq : kmu_exp ] ) reduces to the regular one @xmath66 and becomes zero in the dissociation limit , the first term is an `` on - site '' integral computed with the long - range interaction .
thus we conclude that the coupling between @xmath62 and @xmath63 determinants is determined by the contribution at short range of the long - range interaction . in this case
static correlation , that is also referred to as left - right correlation , can obviously not be interpreted as a purely long - range correlation effect . nevertheless ,
if the error function is used , we see from its taylor expansion for small @xmath67 , @xmath68 that the coupling term can be expanded as follows @xmath69 consequently , as already pointed out by gori - giorgi and savin  @xcite , even an infinitesimal @xmath38 value ensures that @xmath70 is not strictly equal to zero , thus providing the correct multi - configurational description of the dissociated h@xmath1 molecule in the ground state : @xmath71 this simple example illustrates how difficult it is to describe static correlation in the coordinate space , in contrast to short - range dynamical correlation that is connected with the coulomb hole .
we propose in the rest of this paper an alternative approach where correlations are separated in the orbital space . in order to overcome the double counting problem ,
the orbitals occupation will be used rather than the electron density . for clarity
, this approach will be introduced first for model hamiltonians .
as mentioned previously , it is convenient to work in the orbital space rather than the coordinate space when it comes to separate static and dynamical correlation effects .
a change of paradigm is then necessary in order to avoid the double counting problem .
the basic variable in dft is the electron density @xmath72 that is defined in the coordinate space .
as we want to reformulate dft in the orbital space , orbitals occupation @xmath73 seems to be the variable of choice .
this is known in condensed matter physics as _ site - occupation functional theory _ ( soft )  @xcite .
the latter is nothing but the formulation of dft for model hamiltonians such as the hubbard hamiltonian . after a short introduction to the ks - soft scheme in sec .
[ subsec : ks - soft ] , we will show in sec .
[ subsec : soet ] how formal analogies with range - separated dft can lead in this context to an exact embedding theory for model hamiltonians .
let us consider the hubbard hamiltonian with an external potential @xmath74 : @xmath75 where @xmath76 is the hopping integral , @xmath77 denotes the on - site two - electron repulsion , @xmath78 and @xmath79 .
the site - occupation operator equals @xmath80 and , for a given wavefunction @xmath9 , the occupation of site @xmath81 is defined as @xmath82 . for simplicity
we will consider a fixed number of electrons @xmath26 and only discuss ground - state properties .
as shown by gunnarsson and schnhammer  @xcite , the hk theorem can be adapted to the hamiltonian in eq .
( [ eq : hubb_hamil_ext_pot ] ) .
there is indeed a one - to - one correspondence between the external potential @xmath83 and the ground - state sites occupancy @xmath84 .
consequently , for fixed @xmath76 and @xmath77 parameters , the exact ground - state energy can be obtained from the following variational principle , @xmath85 where @xmath86 . the analog of the universal hk functional is a function of the sites occupation that can be written within the ll constrained - search formalism as @xmath87    since the hk theorem still holds when @xmath88 , a ks formulation of soft ( ks - soft ) is obtained from the following partitioning , @xmath89 where the non - interacting kinetic energy functional is defined in analogy with ks - dft as @xmath90 note that , in this context , the non - interacting ks system is _ not _ described by sites that are strictly singly , doubly or not occupied , @xmath91 in other words @xmath92 is not a single slater determinant in the basis of the site orbitals .
indeed , @xmath92 is the ground state of the non - interacting hamiltonian @xmath93 where @xmath94 is _ non - local _ in the orbital space .
let us stress that there is no hk theorem without a non - local term like @xmath94 included into the hamiltonian .
for example , for two electrons , the two hamiltonians that differ by more than a constant , @xmath95 and @xmath96 where @xmath97 and @xmath98 , yield the same ground - state occupancies @xmath99 and @xmath100 .
( 1,0 ) ; ( -0.2,0 ) node [ left]@xmath101 ; ( 0.4,-0.5 )  ( 0.4,0.5 ) ; ( 0.6,-0.5 )  ( 0.6,0.5 ) ; ( 0,1.4 )
( 1,1.4 ) ; ( -0.2,1.4 ) node [ left]@xmath102 ;    ( 0,0 )  ( 1,0 ) ; ( -0.2,0 ) node [ left]@xmath101 ; ( 0.4,-0.5 )  ( 0.4,0.5 ) ; ( 0.6,-0.5 )  ( 0.6,0.5 ) ; ( 0,2.0 )
( 1,2.0 ) ; ( -0.2,2.0 ) node [ left]@xmath103 ;    this is an important difference with ks - dft where the quantity to reproduce from singly or doubly occupied ks orbitals is the exact electron density @xmath72 . in ks - soft , the quantity to reproduce
is the exact sites occupation @xmath21 that may be fractional due to the non - zero hopping integral @xmath104 .
+ by analogy with ks - dft , an exact ac expression for the hxc site - occupation functional in eq .
( [ eq : hxcsofun_def ] ) can be obtained from the following auxiliary equations , @xmath105 with the site - occupation constraint @xmath106 fulfilled for @xmath107 , thus leading to @xmath108 according to the hellmann
feynman theorem .
note that within the mean - field approximation and in the particular case of a singlet ground state , the hxc functional is simplified as follows , @xmath109 the latter expression is usually referred to as the hartree energy in condensed matter physics  @xcite .
the exchange energy is then considered to be zero .      in the spirit of _ dynamical mean - field theory _
@xcite and dmet  @xcite , we propose to isolate one site , referred to as impurity and for which correlation effects will be described explicitly , while the remaining sites ( the bath ) are treated within soft .
such an approach , that will be referred to as _ site - occupation embedding theory _
( soet ) in the following , is formally very similar to range - separated dft .
the main difference is that the interactions are separated here in the orbital space rather than in the coordinate space . for convenience the impurity and bath sites
will be labelled as @xmath110 and @xmath111 , respectively . by analogy with eqs .
( [ eq : ll_range_separation ] ) and ( [ eq : ll_longrange ] ) , soet can be derived from the following partitioning of the ll functional , @xmath112 where the impurity ll functional equals @xmath113 and @xmath114 denotes the complementary hxc bath functional that describes the bath as well as the coupling between the bath and the impurity .
note that , in contrast to range - separated dft , the separation of correlation effects is , in soet , not controlled by a single continuous parameter .
it rather relies on the selection of an impurity site . as discussed in the following more impurity sites
can actually be considered .
+ we now return to the single - impurity case . since , according to eqs .
( [ eq : vp_soft ] ) , ( [ eq : fimp_plus_ebath ] ) and ( [ eq : llimpfun ] ) , any normalized trial wavefunction @xmath9 fulfils @xmath115 the exact ground - state energy can be expressed as @xmath116 where the minimizing wavefunction @xmath117 fulfils the self - consistent equation @xmath118\hat{n}_i\bigg)\vert{\psi}^{\rm
imp}\rangle \nonumber \\ & & = { \mathcal{e}}^{\rm imp } \vert{\psi}^{\rm imp}\rangle .
\end{aligned}\ ] ] thus we obtain an exact embedding scheme where the bath is described by a site - occupation functional .
its connection with dmet will be discussed in the following .
let us first focus on the complementary bath functional . from the ks decomposition @xmath119 where @xmath120 describes the repulsion on the impurity that is embedded into a non - interacting bath , we obtain with eqs .
( [ eq : hxcsofun_def ] ) and ( [ eq : fimp_plus_ebath ] ) , @xmath121 this expression is convenient for developing approximate functionals . for example , within the mean - field approximation , the hxc functional for the impurity is simplified as follows @xmath122 thus leading , when combined with eq .
( [ mf_approx_hxc ] ) , to @xmath123 better approximate functionals could be developed from the local density approximation for a luttinger liquid  @xcite , in complete analogy with the development of short - range functionals from the uniform electron gas  @xcite .
work is currently in progress in this direction .
+ returning to the exact theory , a more explicit expression for the complementary bath functional can be obtained within the ac formalism when considering , in analogy with eq .
( [ eq : ac_full_scale ] ) , the following auxiliary equations @xmath124 where @xmath125 for @xmath126 , @xmath127 , and @xmath128 . in such an ac
the on - site repulsion is switched on continuously site after site .
a graphical representation is given in fig .
[ fig : ac_graph_rep ] .
c    the embedded impurity system is recovered when @xmath129 and @xmath130 or , equivalently , when @xmath131 and @xmath132 . from eqs .
( [ eq : llimpfun ] ) and ( [ eq : ll_imp_ksdecomp ] ) , the site - occupation constraint @xmath133 and the hellmann  feynman theorem , we obtain @xmath134 by introducing the multiple - impurity ll functional @xmath135 the ( fully - interacting ) hxc functional becomes @xmath136 since @xmath137 , thus leading to the alternative expression ( see eq .
( [ eq : hxc_fun_acexp ] ) ) , @xmath138 combining eqs .
( [ eq : bathfun_as_diff ] ) , ( [ eq : impfun_ac_formula ] ) and ( [ eq : hxc_fun_ac_alternative ] ) finally gives the exact expression @xmath139 if the exact complementary bath functional was known , its derivative could be computed and , by solving eq .
( [ eq : sc_eq_imp ] ) self - consistently , we would obtain the exact embedded impurity wavefunction @xmath117 that has the same sites occupation @xmath21 as the physical system described by the hubbard hamiltonian in eq .
( [ eq : hubb_hamil_ext_pot ] ) : @xmath140 if @xmath21 is known , the exact embedding potential @xmath141 can in fact be reached without knowing the complementary bath functional . indeed , according to the variational principle , for any local potential @xmath142 , the embedded impurity wavefunction fulfils @xmath143 where @xmath144 is the ground - state energy of @xmath145 .
thus we obtain from eqs .
( [ eq : llimpfun ] ) , ( [ eq : sc_eq_imp ] ) and ( [ eq : nimp_equals_n ] ) , the legendre  fenchel transform expression  @xcite @xmath146 where @xmath21 is the input sites occupation and @xmath147 the embedding potential to be optimized . in practice approximations to @xmath21 could be used for example by performing a mean - field description of the physical system . in this case ,
the ( approximate ) embedding potential would enable to reproduce the mean - field occupation in both the impurity and the bath .
+ thus a connection with dmet can be established ( see eq .
( 5 ) in ref .
a substantial difference though is that , in dmet , the potential to be optimized is the mean - field interaction operator @xmath147 that defines the mean - field approximation to the hubbard hamiltonian , @xmath148 and the sites occupation to be reproduced by the ground state of @xmath149 is the one obtained from the ground state of the dmet impurity hamiltonian , @xmath150 where @xmath151 .
the terms @xmath152 and @xmath153 are defined in ref .
@xcite and depend indirectly on @xmath147 . note that for simplicity dmet equations are written here for a uniform external potential .
+ at a given iteration @xmath154 of a dmet calculation , the mean - field potential is set to @xmath155 and the impurity hamiltonian generates the ground - state sites occupation @xmath156 .
if it exists a local potential @xmath157 such that the ground state @xmath158 of @xmath159 has the same sites occupation then , according to the variational principle , any trial potential @xmath147 fulfils @xmath160 where @xmath161 is the ground - state energy of @xmath149 , thus leading to the legendre
fenchel transform expression : @xmath162 if a maximum exists , it then corresponds to @xmath163 .
the updated occupation @xmath164 is then obtained from @xmath165 , thus generating a new potential until convergence is reached .
in this respect , dmet can be considered as a ks _ optimized effective potential _ ( oep ) scheme since the mean - field hamiltonian is in this context nothing but a non - interacting hamiltonian .
the oep is used to model the interaction on the impurity .
a formal analogy can actually be made with the range - separated ks - oep approach proposed in ref .
@xcite , where the long - range interaction is described with an oep while the short - range interaction is modeled with a density functional .
the oep is then obtained from a density constraint , exactly like in dmet .
+ let us finally mention that an alternative and exact formulation of dmet is obtained when introducing the following partitioning of the ll functional @xmath166    where @xmath167    in connect to connect dmet with soet one may fix @xmath147 and let @xmath153 vary independently as a potential to be optimized .
if it exists a potential @xmath168 such that the ground state @xmath169 of @xmath170 has its site occupation equal equal to the one generated by @xmath149 , that we denote @xmath171 , then we have according to the variational principle , @xmath172    at the first iteration of a dmet calculation , the latter will be fixed and for any new trial potential @xmath147 , the variational principle leads to @xmath173    let us finally mention that , in practical soet calculations , the embedding potential would be obtained self - consistently by using an approximate ( local ) hxc functional for the bath and solving eq .
( [ eq : sc_eq_imp ] ) .
the computational implementation of soet is left for future work .
an illustration is given in sec .
[ sec : two_site_hubb ] with the simple two - site hubbard hamiltonian .
+    returning to soet , let us finally mention that the theory can be extended to more impurity sites simply by using the following partitioning of the ll functional , @xmath174 where @xmath175 . in this case , the ac formula for the complementary bath functional becomes @xmath176
soet is applied in this section to the simple two - site hubbard model with a uniform external potential .
the analytical construction of the ac is presented in sec .
[ subsec : ac_for_two_sites ] .
the resulting integrand expressions are then analyzed in sec .
[ subsec : results ] .
let us consider a two - electron system described by the two - site hubbard hamiltonian @xmath177 since we are interested here in the singlet ground state only , the matrix representation of the hamiltonian can be reduced to the basis of the `` doubly - occupied site '' states @xmath178 with @xmath110 or @xmath179 , and @xmath180 that corresponds to singly - occupied sites , thus leading to @xmath181 & = & \left [ \begin{array}{c c c } u_0+v_0-v_1 & 0 & -\sqrt{2}t\\ 0 & u_1+v_1-v_0 & -\sqrt{2}t\\   -\sqrt{2}t & -\sqrt{2}t & 0   \end{array } \right ] + ( v_0+v_1 )   .\end{aligned}\ ] ] for simplicity we choose for the physical hamiltonian @xmath182 and @xmath183 .
consequently the site occupation is uniform : @xmath184 we obtain by diagonalization the well - known expressions for the ground - state energy @xmath185 and the corresponding wavefunction @xmath186 where @xmath187 since @xmath188 , thus leading to the more explicit expression @xmath189 note the hellmann
feynman theorem , @xmath190 that will be used in the following .
the double occupancy is the same for both sites and equal to @xmath191 along the ac described in eq .
( [ eq : ac_impurities ] ) and fig .
[ fig : ac_graph_rep ] , symmetry in the on - site repulsions is broken in the auxiliary hamiltonian operator .
this is clear for the embbeded impurity where @xmath192 and @xmath193 .
nevertheless , symmetry can be restored in the hamiltonian matrix simply by adjusting the local potential components @xmath194 and @xmath195 such that @xmath196 thus leading to @xmath197 the latter condition defines uniquely ( up to a constant ) the local embedding potential that gives a uniform site occupation . in this case
the effective on - site repulsion is the same on each site and is simply expressed as @xmath198 therefore the ac can be constructed analytically as follows for the impurity ( @xmath129 ) , @xmath199 and for the bath ( @xmath131 ) , @xmath200 consequently the wavefunctions along the ac for the bath and the impurity are @xmath201 and @xmath202 respectively .
the corresponding double occupancies , that are nothing but hxc integrands per unit of @xmath77 for the impurity and the bath ( see eqs .
( [ eq : impfun_ac_formula ] ) and ( [ eq : ac_comp_bath_exp ] ) ) , respectively , can then be expressed , according to eqs .
( [ eq : explicit_exp_2bleocc ] ) and ( [ eq : same_2bleocc ] ) , as @xmath203 in order to derive the expression for the bath integrand in eq .
( [ eq : ac_comp_bath_exp ] ) , we express the embedded wavefunction along the ac as @xmath204 where the coefficient @xmath205 is obtained by diagonalizing the hamiltonian matrix in eq .
( [ hmatrix ] ) and using eqs .
( [ eq : hubb_gs_ener ] ) and ( [ eq : ac_p1 ] ) , thus leading to @xmath206 since @xmath207 we obtain the alternative expression @xmath208 and @xmath209 note that the embedded impurity wavefunction is obtained when @xmath130 along the ac for the impurity , and @xmath132 along the ac for the bath : @xmath210 since , according to eqs .
( [ hmatrix ] ) and ( [ eq : auxhamil_ac_bath ] ) , @xmath211 it comes from the site occupation constraint in eq .
( [ eq : so_constraint_two_sites ] ) that @xmath212 moreover , according to eqs .
( [ hell  feyn_th_2sites ] ) and ( [ eq : wf_ac_bath ] ) , the double occupancy can be rewritten along the ac for the bath as @xmath213 thus giving with eq .
( [ eq : ac_comp_bath_exp ] ) , @xmath214 or , more explicitly , @xmath215 as expected , the physical energy @xmath216 is recovered when adding contributions from the impurity and the bath : @xmath217    let us finally mention that , in the conventional ac ( see eq .
( [ eq : ac_full_scale ] ) ) , the on - site repulsion is scaled on both the impurity and the bath sites . in this case
the interaction on the impurity is not separated from the repulsions in the bath .
the corresponding ac will therefore be defined as @xmath218 and @xmath219 as a result the ac integrand can be expressed as @xmath220 thus leading to the explicit expression @xmath221 note that , according to eq .
( [ hell  feyn_th_2sites ] ) , the hxc integrand can also be rewritten as @xmath222 finally by rewriting the double occupancy along the ac for the impurity ( see eq .
( [ eq:2ble_occ_imp_explicit_exp ] ) ) as follows , @xmath223 we see from eqs .
( [ eq : hxc_fun_acexp ] ) , ( [ eq : hxc_fun_ac_alternative ] ) , ( [ eq:2ble_occ_ac_bath ] ) and ( [ eq : convac_hellth_int ] ) that , as expected , both the conventional hxc integrand and the sum of bath and impurity integrands lead to the same hxc energy after integration over [ 0,1 ] : @xmath224 or , more explicitly , @xmath225    * * * * * * * * * * * * * * * * * * * * * * * * + the embedded impurity wavefunction is actually obtained for @xmath132 , @xmath226 so that since sites are both singly occupied , it comes @xmath227    when defining the hartree
exchange energy as the mean field one , that is @xmath228 for a singly occupied site , the bath correlation energy can be written as @xmath229 where the bath correlation integrand equals @xmath230    ac for the impurity : @xmath231      c    double occupancies along the ac for the impurity and the bath ( see eqs .
( [ eq:2ble_occ_imp_explicit_exp ] ) and ( [ eq:2ble_occ_bath_explicit_exp ] ) ) , that can be identified as hxc integrands per unit of @xmath77 , have been plotted for various correlation regimes .
results are shown in figs .
[ fig:2site_hubb_imp_ac ] and [ fig:2site_hubb_bath_ac ] . for analysis purposes
, comparison is made between the sum of these two integrands and the conventional hxc integrand in eq .
( [ eq : conv_ac_int_explicit_exp ] ) .
results are shown in fig .
[ fig:2site_hubb_bath_comparacs ] .
+ let us first focus on the impurity .
as expected from accurate quantum chemical calculations of the ac for the h@xmath1 molecule at equilibrium  @xcite , the integrand varies linearly with the interaction strength @xmath232 in the weakly correlated regime ( @xmath233 ) .
the integrand gains curvature as @xmath234 increases . in the strongly correlated regime ( @xmath235 ) , the double occupancy of the impurity becomes zero for large interaction strengths ( see the bottom panel in fig .
[ fig:2site_hubb_imp_ac ] ) .
it only varies in the vicinity of the non - interacting case ( @xmath132 ) with a large negative slope .
this pattern was found by teale  for the stretched h@xmath1 molecule when computing the ac at the _ full configuration interaction _ ( fci ) level in a large basis set  ( see fig .
5 ( a ) in ref .
one should mention that , along the ac for the impurity , the on - site repulsion in the bath is set to zero .
comparison with ref .
@xcite is then more relevant when considering the conventional ac integrand for the two sites . as shown in fig .
[ fig:2site_hubb_bath_comparacs ] , all correlation regimes observed numerically for h@xmath1 along the ac and the bond breaking coordinate are qualitatively well reproduced by applying soft to the two - site hubbard hamiltonian .
+    c    turning to the ac for the bath , we first notice in the top panel of fig .
[ fig:2site_hubb_bath_ac ] that the double occupancy of the bath remains essentially linear in @xmath232 for all correlation regimes , in contrast to the ac for the impurity .
this is due to the fact that the on - site repulsion is already switched on for the impurity .
scaling is only used for the repulsion on the bath site so that the embedded impurity ( @xmath132 ) can be connected with the physical system ( @xmath130 ) .
as expected , the two systems exhibit different double occupancies even though they have the same sites occupation ( see the bottom panel in fig .  [ fig:2site_hubb_bath_ac ] ) . variations with the correlation strength are consistent with those obtained for the 1d - hubbard model ( see case @xmath236 in fig .  2 of ref .
+    c    let us finally discuss the separation of interactions between the impurity and the bath . as mentioned previously and illustrated in fig .
[ fig:2site_hubb_bath_comparacs ] , both the conventional and bath+impurity ac integrands give the same total hxc energy after integration over @xmath237 $ ] but they of course differ along the ac . in the weakly correlated regime ,
both the bath and the impurity contribute significantly to the hxc energy .
note that , in practical calculations , we want to describe the embedded impurity only , meaning that the contributions of the bath to the hxc energy will be described by a site - occupation functional .
it becomes clear from this simple example that the latter can contribute significantly and that approximations beyong the mean field ( which is exact only when @xmath238 ) will be needed .
returning to the ac , in the strongly correlated regime , the contribution of the bath becomes negligible for all interaction strengths .
note that , for @xmath132 , the impurity and conventional hxc integrands are equal to @xmath239 and @xmath240 , respectively .
since both are non - zero only for small interaction strengths and behave similarly ( with a large negative slope ) , they both integrate to the same hxc energy simply because the effective on - site repulsion along the ac is smaller for the impurity alone ( @xmath241 ) than for the two sites ( @xmath242 ) .
we propose in this section to apply soft to quantum chemical hamiltonians .
the formulation of a ks scheme in this context is discussed in sec .
[ subsec : ks - soft_qc ] .
we then show in sec .
[ subsec : cassoft ] that soet can lead to an alternative casdft method .      in order to adapt soft to quantum chemistry ,
let us first consider the orthonormal basis of molecular orbitals @xmath243 obtained for non - interacting electrons , so that @xmath244 with @xmath245 this is obviously not a good starting point for conventional quantum chemical calculations but it is convenient for deriving a ks - soft scheme in this context . as explained further in the following , the theory can in principle be adapted to different choices of orbitals . using eq .
( [ eq : orbitals_soft_qc ] ) leads to the following second - quantized expression for the molecular hamiltonian @xmath246 , @xmath247 where @xmath248 denotes the fluctuation potential .
the hartree  fock ( hf ) potential operator can be decomposed as follows , @xmath249 where @xmath250 denotes a regular two - electron integral . the first term on the right - hand side of eq .
( [ eq : uhfdecomp ] ) is local in the orbital space while the second one is non - local and thus enables hopping between occupied and unoccupied orbitals , hence the notation @xmath251 in analogy with the kinetic energy operator in the hubbard hamiltonian .
therefore the molecular hamiltonian becomes @xmath252 with @xmath253 .
this expression enables a direct comparison with the hubbard hamiltonian in eq .
( [ eq : hubb_hamil_ext_pot ] ) : the sites are now molecular orbitals and their energies @xmath254 play the role of the external potential . + following sec .
[ subsec : ks - soft ] , we introduce the ll functional @xmath255 so that the exact ground - state energy , that is the fci energy in this context , can be written as @xmath256 let us consider the ks partitioning , @xmath257 where @xmath258 and , in contrast to eq .
( [ eq : hxcsofun_def ] ) , no hartree and exchange energy contributions have been introduced .
indeed , the fluctuation potential induces correlation effects only since , when written with normal ordering  @xcite , it generates double excitations only .
consequently , within ks - soft , the fci energy is obtained variationally as follows , @xmath259 where the minimizing wavefunction @xmath260 fulfils the self - consistent equation @xmath261\hat{n}_p\bigg)\vert{\psi^{\rm
ks}}\rangle = { \mathcal{e}^{\rm ks}}\vert{\psi^{\rm ks } } \rangle.\end{aligned}\ ] ] eqs .
( [ eq : fciener_var_soft ] ) and ( [ eq : ks - soft_eq_qc ] ) should in principle enable to recover the fci energy without treating electron correlation explicitly .
orbital rotations would therefore be sufficient .
one should of course investigative potential @xmath262-representability problems for example by computing legendre  fenchel transforms .
approximate functionals may also be developed from model hamiltonians such as the hubbard hamiltonian .
let us finally note that different orbitals could be used in eq .
( [ eq : ks - soft_eq_qc ] ) .
one would need to adapt the correlation functional in order to recover the correct fci energy .
work is in progress in these directions .
following knizia and chan  @xcite , we propose in this section to apply soet to the molecular hamiltonian in eq .
( [ eq : mol_hamil_like_hubb ] ) .
the impurity sites will be the active orbitals ( @xmath263 ) whose selection is usually based on chemical intuition , while the inactive and virtual orbitals correspond to the bath .
a graphical representation is given in fig .
[ fig : active_space_graph_rep ] . as already mentioned for soet , the separation of correlation effects is not controlled here by a single continuous parameter like in range - separated dft . in the spirit of a regular casscf calculation
, it rather relies on the selection of active orbitals .
this ensures that only correlation effects within the active space are described explicitly in wft while the remaining correlation , including core correlation , is modeled by an orbital - occupation functional .
@xmath264 separating correlation effects in the active orbital space @xmath265 from the remaining ones leads to the alternative partitioning of the ll functional , @xmath266 where @xmath267 and @xmath268 is the reduction of the fluctuation potential operator to the active orbital space : @xmath269 note that the hf potential in the second term on the right - hand side of eq .
( [ eq : fluctu_pot_active ] ) is here calculated with the active orbitals that are doubly occupied , by analogy with eq .
( 31 ) in ref .
@xcite . combining eqs .
( [ eq : fciener_vp_n ] ) , ( [ eq : part_f_imp_dyn ] ) and ( [ eq : fimp_ll_qc ] ) leads to the exact variational expression @xmath270 where the minimizing wavefunction @xmath271 fulfils the self - consistent equation @xmath272\hat{n}_p\bigg)\vert{\psi_{\mathcal a}^{\rm imp}}\rangle = { \mathcal{e } _ { \mathcal a}^{\rm imp}}\vert{\psi_{\mathcal a}^{\rm imp } } \rangle.\end{aligned}\ ] ] the method will be referred to as cassoft .
interestingly the impurity hamiltonian is very similar to the one proposed by zgid and chan in ref .
a major difference is that the embedding potential @xmath273 should in principle enable to recover the exact orbitals occupancy in both the bath and the impurity and , through the functional , the fci energy .
the impurity hamiltonian can also be viewed as an embedded version of dyall s hamiltonian  @xcite . @xmath274k .
g. dyall , j. chem .
102 , 4909 ( 1995 ) .
+ @xmath275c .
angeli , r. cimiraglia , and j .-
malrieu , j. chem .
117 , 9138 ( 2002 ) .
a numerical validation of eq .
( [ eq : sc_cassoft ] ) for small molecules can be achieved when computing from the input fci orbitals occupation @xmath21 the legendre
fenchel transform @xmath276 where @xmath277 is the ground - state energy of @xmath278 .
work is in progress in this direction .
+ let us stress that , in contrast to the short - range density functional used in range - separated dft , the complementary bath orbital - occupation functional @xmath279 is not a universal functional .
it will _ a priori _ depend on the molecular orbital basis in which cassoft is formulated .
it also depends on the active orbital space .
even though range - separated dft is much simpler in that respect , the choice of an optimal range - separation parameter is , in terms of accuracy , not universal partly because static correlation is not a purely long - range effect .
local and semi - local approximations to the short - range density functional are usually not accurate enough for modeling strongly multi - configurational systems ( see ref .
@xcite and the references therein ) . working in the orbital space will at least allow for a proper separation of static and dynamical correlation effects .
one may consider the development of approximate functionals to be a cumbersome task in this context . a similar challenge is to some extent encountered in _ natural orbital functional theory _
( noft )  @xcite .
recent advances in the field might be useful for such developments .
this is left for future work .
+ returning to exact cassoft , an ac formula can actually be derived for the complementary bath correlation functional @xmath279 by introducing orbital - dependent active spaces @xmath280 where the indices @xmath281 and @xmath282 refer to the lowest and highest active orbitals in energy , respectively .
let @xmath283 denote the active space of interest @xmath265 .
we then consider the auxiliary equations @xmath284 where the orbitals occupation constraint @xmath285 is fulfilled for @xmath107 , @xmath286 and @xmath287 .
the virtual and inactive increment operators are defined as @xmath288 respectively .
the explicit expression for @xmath289 is deduced from eq .
( [ eq : fluctu_pot_active ] ) .
the superscripts `` @xmath290 '' and `` @xmath291 '' refer to the incorporation into the active space @xmath292 of the virtual @xmath293 and inactive @xmath294 orbitals , respectively .
a graphical representation is given in fig .
[ fig : active_space_graph_rep ] .
c    the embedded active electrons are recovered along the ac when @xmath295 , @xmath296 and @xmath132 .
according to eqs .
( [ eq : part_f_imp_dyn ] ) and ( [ eq : fimp_ll_qc ] ) , the complementary bath correlation functional can be expressed as @xmath297 thus leading , by analogy with eq .
( [ eq : hxc_fun_ac_alternative ] ) , to @xmath298 it is readily seen from eq .
( [ eq : ac_bath_cassoft ] ) that the complementary bath correlation functional describes correlation effects that ( i ) couple inactive or active orbitals with the virtual orbitals ( first term on the right - hand side ) and that ( ii ) couple inactive orbitals with active orbitals ( second term ) .
the correlation effects within the active space are treated explicitly .
note that cassoft is in principle exact ( in a given basis set ) and it is free from double counting problems .
the reason is that it relies on an exact separation of the ll functional in the orbital space .
the exact formulation of multi - configuration density - functional theory has been discussed .
the infamous double counting problem can be avoided when separating correlation effects either in the coordinate space or in the orbital space . in the latter case
, orbitals occupation should be used as basic variable rather than the electron density .
this approach has been applied to the hubbard hamiltonian , thus leading to an exact _ site occupation embedding theory _ ( soet ) .
the connection with _ density matrix embedding theory _ ( dmet ) has been discussed and an adiabatic connection ( ac ) formula has been derived for the complementary bath functional .
the ac has been constructed analytically for the simple two - site hubbard model .
the computational implementation of soet as well as the development of approximate local bath functionals is left for future work .
+ we then proposed to apply soet to a quantum chemical hamiltonian , thus showing that multi - configuration methods can be merged rigorously with orbital - occupation functionals .
the method is referred to as cas _
site occupation functional theory _ ( cassoft ) . in this
context , impurity sites correspond to the active orbitals while inactive and virtual orbitals are the bath . an ac formalism has also been derived .
the latter should be useful for developing approximate functionals .
a connection with multi - reference perturbation theory may be achieved from a perturbation expansion of the ac integrand .
note also that , in regular dft , the uniform electron gas played a crucial role in the development of electron density functionals . by analogy ,
model hamiltonians such as the hubbard hamiltonian could be used for developing local and semi - local orbital - occupation functionals . finally , following knizia and chan  @xcite , it would be interesting to explore connections between cassoft and _ density matrix renormalization group _ ( dmrg ) methods  @xcite .
work is in progress in these directions .
the author would like to thank the editorial board for the kind invitation to submit a `` new views '' paper to mol .
ef is grateful to markus reiher and stefan knecht for their kind invitation to give a seminar at eth in august 2014 and their helpful comments on this work .
ef thanks masahisa tsuchiizu , vincent robert , laurent mazouin , andreas savin , lucia reining and bernard amadon for stimulating discussions on strongly correlated electrons .
the author finally acknowledges financial support from the phc program sakura 2969uk , the labex `` chemistry of complex systems '' and the anr ( mcfunex project ) .


in recent years , a considerable amount of research has focused on the study of graph structures arising from technological , biological and sociological systems .
graphs are the tool of choice in modeling such systems since the latter are typically described as a set of pairwise interactions .
important examples of such datasets are the internet graph ( vertices are routers , edges correspond to physical links ) , the web graph ( vertices are web pages , edges correspond to hyperlinks ) , social networks ( vertices are humans , edges correspond to friendships ) , information networks like facebook and linkedin ( vertices are accounts , edges correspond to online friendships ) , biological networks ( vertices are proteins , edges correspond to protein interactions ) , math collaboration network ( vertices are mathematicians , edges correspond to collaborations ) and many more .    towards the end of the 90s , a series of papers observed that the classic models of random graphs introduced by erds and rnyi @xcite and gilbert @xcite did not explain the empirical properties of real - world networks @xcite .
typical properties of such networks include skewed degree distributions @xcite , large clustering coefficients @xcite and small average distances @xcite , a phenomenon typically referred to as `` small worlds '' .
skewed degree distributions have widely been modeled as power laws : the proportion of vertices of a given degree follows an approximate inverse power law , i.e. , the proportion of vertices of degree k is approximately @xmath14 .
understanding the properties of real world networks has attracted considerable research interest in the recent years @xcite .
a large body of work has focused on finding models which generate graphs which mimick real world networks .
most of the existing work on network formation models falls into the stochastic and the strategic categories .
kronecker graphs @xcite , the cooper - frieze model for the web graph @xcite , the aiello - chung - lu model @xcite , protean graphs @xcite and numerous other models belong to the former category .
the strategic approach has its origins in the work of boorman @xcite , aumann @xcite and aumann and myerson @xcite and a variety of models fall into this category , e.g. , the fabrikant - koutsoupias - papadimitriou model @xcite . despite the large amount of work on models which generate power law graphs ,
a considerably smaller amount of work has focused on generative models for planar graphs .
planar graphs have been studied mainly in the context of transportation networks and occur in numerous real - world problems : modelling city streets @xcite , crowd simulations @xcite , river networks , railway and road maps printed circuits , chemical molecules , see also @xcite and references therein .
in this work , we prove fundamental properties of random apollonian networks ( rans ) @xcite , a popular model of planar graphs with power law properties @xcite .
we use the symbol @xmath15 to denote the random graph at time @xmath2 .
the details of the model are exposed in section  [ sec : ran ] .
specifically our main results are the following theorems :    let @xmath16 be the @xmath0 highest degrees of the random apollonian network @xmath15 at time @xmath2 where @xmath0 is a fixed positive integer .
also , let @xmath17 be a function such that @xmath4 as @xmath5 .
then _ whp _
@xmath18    and for @xmath19    @xmath20 [ thrm : thrm1 ]    the growing function @xmath17 can not be removed , see @xcite .
using theorem  [ thrm : thrm1 ] we relate the highest @xmath0 degrees and eigevalues .
let @xmath21 be the largest @xmath0 eigenvalues of the adjacency matrix of @xmath15 . then _ whp _
@xmath22 [ thrm : thrm2 ]    also , we prove the following theorem for the diameter .    the diameter @xmath10 of @xmath15 satisfies in probability @xmath23 where @xmath12 is the unique solution greater than 1 of the equation @xmath13 .
[ thrm : ternary ]    the outline of the paper is as follows : in section  [ sec : ran ] we describe the model and discuss existing work on rans . in sections
[ sec : degree ] and  [ sec : eigen ] we prove theorems  [ thrm : thrm1 ] and  [ thrm : thrm2 ] respectively . in section  [ sec :
diameter ] we give a simple proof for the asymptotic growth of the diameter and prove theorem  [ thrm : ternary ] . finally in section  [ sec : waiting ] we investigate another property of the model .
finally , in section  [ sec : concl ] we conclude by proposing several problems for future work .
apollonius of perga was a greek geometer and astronomer noted for his writings on conic sections .
apollonius introduced the problem of space filling packing of spheres which in its classical solution , the apollonian packing see @xcite for remarkable properties of apollonian circles- , exhibits a power law behavior .
specifically , the circle size distribution follows a power law with exponent of about 1.3 @xcite .
apollonian networks ( ans ) were introduced in @xcite and independently in @xcite .
zhou et al . @xcite
introduced the random apollonian networks ( rans ) and gave formulae for their order , size , degree distribution ( using heuristic arguments ) , clustering coefficient and diameter .
high dimensional rans were introduced in @xcite .
the degree distribution of rans was shown to follow a power law in @xcite .
other properties have also been analyzed , e.g. , average distance of two vertices @xcite , properties of the connectivity profile @xcite .
zhou et al .
@xcite proposed a simple rule that generates a random two - dimensional apollonian networks with very large clustering coefficient .
rans are planar 3-trees , a special case of random @xmath0-trees @xcite .
the general result of the degree distribution of random @xmath0-trees was proved by cooper & uehara and gao @xcite . in rans in
contrast to the general model of random @xmath0 trees the random @xmath0 clique chosen at each step has never previously been selected .
for example in the two dimensional case any chosen triangular face is being subdivided into three new triangular faces by connecting the incoming vertex to the vertices of the boundary .
random @xmath0-trees due to their power law properties have been proposed as a model for complex networks , see , e.g. , @xcite and references therein .
recently , a variant of @xmath0-trees , namely ordered increasing @xmath0-trees has been proposed and analyzed in @xcite .
an example of a two dimensional ran is shown in figure  [ fig : fig1 ] .
the ran generator takes as input a parameter @xmath24 equal to the number of iterations the algorithm will perform and runs as follows :    * let @xmath25 be a triangle @xmath26 embedded on the plane , e.g. , as an equilateral triangle . * _ for _ @xmath27 to @xmath24 : * * sample a face @xmath28 of the planar graph @xmath29 uniformly at random . * * insert the vertex @xmath30 inside this face ( e.g. , in the barycenter of the corresponding triangle ) and draw the three edges @xmath31 , @xmath32,@xmath33 , e.g. , as straight lines .
it s worth pointing out why we expect this model to yield a power law degree distribution : consider any vertex @xmath34 at some time @xmath2 .
vertex @xmath34 sees a wheel graph around it ( except for vertices 1,2,3 who see a wheel modulo one edge ) as figure  [ fig : fig2 ] shows . since we pick a face uniformly at random the probability that @xmath34 increases its degree by 1 is ( roughly ) proportional to its degree .
therefore , we expect this process to result in a graph , which is obviously planar , and which _ whp
_   occurs with high probability _ whp _  if @xmath35}}=1 $ ] .
] has properties similar to those generated by a preferential attachment process @xcite .
therefore , it should not come as a surprise that existing techniques developed in the context of preferential attachment models can be adapted to solve problems concerning the structure of random apollonian networks .
finally , we shall make use of the following formulae for the number of vertices ( @xmath36 ) , edges ( @xmath37 ) and faces ( @xmath38 ) at time @xmath2 in a ran @xmath15 : @xmath39 note that a ran is a maximal planar graph since for any planar graph @xmath40 .
we decompose the proof of theorem  [ thrm : thrm1 ] into several lemmas which we prove in the following .
the proof of theorem  [ thrm : thrm1 ] follows from lemmas  [ lem : lemma2 ] ,  [ lem : lemma3 ] ,  [ lem : lemma4 ] ,  [ lem : lemma5 ] ,  [ lem : lemma6 ] .
we partition the vertices into three sets : those added before @xmath41 , before @xmath42 and after @xmath42 where @xmath43 and @xmath44 .
we define a supernode to be a collection of vertices and the degree of the supernode the sum of the degrees of its vertices .
let @xmath45 denote the degree of vertex @xmath46 at time @xmath2 . and let @xmath47denote the rising factorial function .
then , for any positive integer
@xmath0    @xmath48 } \leq \frac { ( k+2)!}{2 } \big ( \frac{2t}{s } \big)^{\frac{k}{2}}. \label{eq : rising}\ ] ]    [ lem : lemma1 ]    we distinguish two types of vertices .
this analysis will be omitted in the subsequent lemmas . ] since as it can also be seen in figure  [ fig : fig2] the three initial vertices @xmath49 have one face less than their degree whereas all other vertices have degree equal to the number of faces .
+ note that @xmath50 and that @xmath51 . by conditioning
successively we obtain    @xmath52 }                 & = { { \mathbb e}\left [ { { { \mathbb e}\left [ { d_t(s)^{(k ) } | d_{t-1}(s)}\right ] } } \right]}\\                                     & = { { \mathbb e}\left [ { ( d_{t-1}(s))^{(k ) } \big ( 1-\frac{d_{t-1}(s)}{2t-1 } \big ) +   ( d_{t-1}(s)+1)^{(k ) } \frac{d_{t-1}(s)}{2t-1}}\right ] } \\                                      & =   { { \mathbb e}\left [ {   ( d_{t-1}(s))^{(k ) } \big ( 1-\frac{d_{t-1}(s)}{2t-1 } \big ) + ( d_{t-1}(s))^{(k ) } \frac{d_{t-1}(s)+k}{d_{t-1}(s ) } \frac{d_{t-1}(s)}{2t-1}}\right]}\\                                     & =   { { \mathbb e}\left [ { ( d_{t-1}(s))^{(k)}}\right]}\big ( 1+\frac{k}{2t-1 } \big ) = ... = 3^{(k ) } \prod_{t'=s+1}^t ( 1+\frac{k}{2t'-1 } ) \\                                     & \leq   3^{(k ) }   \exp{\big ( \sum_{t'=s+1}^t \frac{k}{2t'-1 } \big ) } \leq 3^{(k ) } \exp { \big ( k \int_{s}^{t } \ !
\frac{\mathrm{d}x}{2x-1 }    \ ,    \big ) }   \\                                      & \leq \frac{(k+2)!}{2 } \exp { \big ( k \log { \frac{2t-1}{2s-1 } } \big ) } \leq \frac { ( k+2)!}{2 } \big ( \frac{2t}{s } \big)^{\frac{k}{2}}. \end{aligned}\ ] ]     + note that initially the degree of any such vertex is 2 . for any @xmath53    @xmath52 }                 & = { { \mathbb e}\left [ { { { \mathbb e}\left [ { d_t(s)^{(k ) } | d_{t-1}(s)}\right ] } } \right]}\\                                     & = { { \mathbb e}\left [ { ( d_{t-1}(s))^{(k ) } \big ( 1-\frac{d_{t-1}(s)-1}{2t-1 } \big ) +   ( d_{t-1}(s)+1)^{(k ) } \frac{d_{t-1}(s)-1}{2t-1 } } \right]}\\                                      & =   { { \mathbb e}\left [ { ( d_{t-1}(s))^{(k ) } \big ( 1+\frac{k}{2t-1 } \big ) - ( d_{t-1}(s))^{(k ) } \frac{k}{(2t-1)d_{t-1}(s ) }   } \right]}\\                                     & \leq { { \mathbb e}\left [ { ( d_{t-1}(s))^{(k)}}\right ] } \big ( 1+\frac{k}{2t-1 } \big ) \leq \ldots \leq \frac { ( k+2)!}{2 } \big ( \frac{2t}{s } \big)^{\frac{k}{2}}.\end{aligned}\ ] ]    the degree @xmath54 of the supernode @xmath55 of vertices added before time @xmath41 is at least @xmath56 _ whp _  .
[ lem : lemma2 ]    we consider a modified process @xmath57 coupled with the ran process , see also figure  [ fig : fig3 ] .
specifically , let @xmath58 be the modified degree of the supernode in the modified process @xmath57 which is defined as follows : for any type of insertion in the original ran process note there exist three types of insertions with respect to how the degree @xmath54 of the supernode ( black circle ) gets affected , see also figure  [ fig : fig3] @xmath58 increases by 1 .
we also define @xmath59 .
note that @xmath60 for all @xmath61 .
let @xmath62 and @xmath63}}$ ] .
the following technical claim is proved in the appendix  [ sec : appendix ] .
@xmath64 [ claim : lemma2 ]
let @xmath65 denote the event that the supernode consisting of the first @xmath41 vertices has degree @xmath58 in the modified process @xmath57 less than @xmath56 .
note that since @xmath66 it suffices to prove that @xmath67 } } =   o(1)$ ] .
using claim  [ claim : lemma2 ] we obtain    @xmath68 } } & \leq \sum_{r=0}^ { t_0^{1/4}\sqrt{t } -(6t_0 + 6 ) } { r+6t_0 + 5 \choose 6t_0 + 5 }   \big(\frac{2t_0 + 3}{2t+1}\big)^{3t_0 + 3 } e^{-\frac{3}{2}-2t_0+\frac{2t_0^{1/4}}{3 } } \\
& \leq t_0^{1/4}t^{1/2 } \frac{\big ( t_0^{1/4}t^{1/2 } \big)^{6t_0 + 5}}{(6t_0 + 5 ) ! } \big(\frac{2t_0 + 3}{2t+1}\big)^{3t_0 + 3 } e^{-\frac{3}{2}-2t_0+\frac{2t_0^{1/4}}{3 } } \\   & \leq \big ( \frac{t}{2t+1 } \big)^{3t_0 + 3 }   \frac{t_0^{3t_0/2 + 3/2 } ( 2t_0 + 3)^{3t_0 + 3}}{(6t_0 + 5)^{6t_0 + 5 } } e^{4t_0 + 7/2 + 2/3 t_0^{1/4}}\\ & \leq 2^{-(3t_0 + 3 ) } \frac{e^{4t_0 + 7/2 + 2/3 t_0^{1/4}}}{(6t_0 + 5)^{\tfrac{3}{2}t_0+\tfrac{1}{2 } } }   = o(1).\end{aligned}\ ] ]    no vertex added after @xmath42 has degree exceeding @xmath69 _ whp _  . [
lem : lemma3 ]    let @xmath70 denote the event that some vertex added after @xmath42 has degree exceeding @xmath69 .
we use a union bound , a third moment argument and lemma  [ lem : lemma1 ] to prove that @xmath71}}=o(1)$ ] .
specifically    @xmath72 } } & \leq   \sum_{s = t_1}^t { { { \bf{pr}}\left[{d_t(s ) \geq t_0^{-2}t^{1/2}}\right ] } } = \sum_{s = t_1}^t { { { \bf{pr}}\left[{d_t(s)^{(3 ) } \geq ( t_0^{-2}t^{1/2})^{(3)}}\right ] } } \\                        & \leq   t_0^{6}t^{-3/2 } \sum_{s = t_1}^t { { \mathbb e}\left[{d_t(s)^{(3)}}\right ] } \leq 5!\sqrt{2 } t_0 ^ 6 \sum_{s = t_1}^t s^{-3/2 }             \leq 5!2\sqrt{2 } t_0 ^ 6 t_1^{-1/2 } = o(1 ) .
\end{aligned}\ ] ]    no vertex added before @xmath42 has degree exceeding @xmath73 _ whp _  .
[ lem : lemma4 ]    let @xmath74 denote the event that some vertex added before @xmath42 has degree exceeding @xmath73 .
we use again a third moment argument and lemma  [ lem : lemma1 ] to prove that @xmath75}}=o(1)$ ] .
@xmath76 } } & \leq \sum_{s=1}^{t_1 } { { { \bf{pr}}\left [ { d_t(s ) \geq t_0^{1/6}t^{1/2}}\right ] } } = \sum_{s=1}^{t_1 } { { { \bf{pr}}\left [ { d_t(s)^{(3 ) } \geq ( t_0^{1/6}t^{1/2})^{(3 ) } } \right ] } } \\              & \leq t_0^{-1/2}t^{-3/2 } \sum_{s=1}^{t_1 } { { \mathbb e}\left [ { d_t(s)^{(3 ) } } \right ] }   \leq t_0^{-1/2}t^{-3/2 } \sum_{s=1}^{t_1 } 5!\sqrt{2 } \frac{t^{3/2}}{s^{3/2 } } \\              & \leq 5!\sqrt{2}\zeta(3/2)t_0^{-1/2 } = o(1)\end{aligned}\ ] ]    where @xmath77 .
the @xmath0 highest degrees are added before @xmath42 and have degree @xmath78 bounded by @xmath79 _ whp _  .
[ lem : lemma5 ]    for the upper bound it suffices to show that @xmath80 .
this follows immediately by lemmas  [ lem : lemma3 ] and  [ lem : lemma4 ] .
the lower bound follows directly from lemmas  [ lem : lemma2 ] ,  [ lem : lemma3 ] and  [ lem : lemma4 ] .
assume that at most @xmath81 vertices added before @xmath42 have degree exceeding the lower bound @xmath82 .
then the total degree of the supernode formed by the first @xmath41 vertices is @xmath83 .
this contradicts lemma  [ lem : lemma2 ] .
finally , since each vertex @xmath84 has degree at most @xmath85 the @xmath0 highest degree vertices are added before @xmath42 _ whp _  .
the @xmath0 highest degrees satisfy @xmath86 _ whp _  .
[ lem : lemma6 ]    let @xmath87 denote the event that there are two vertices among the first @xmath42 with degree @xmath82 and within @xmath88 of each other .
by the definition of conditional probability and lemma  [ lem : lemma3 ]    @xmath89 } }   & = { { { \bf{pr}}\left[{\mathcal{a}_4|\bar{\mathcal{a}_3}}\right ] } } { { { \bf{pr}}\left[{\bar{\mathcal{a}_3}}\right ] } } + { { { \bf{pr}}\left[{\mathcal{a}_4|\mathcal{a}_3}\right]}}{{{\bf{pr}}\left[{\mathcal{a}_3}\right ] } } \leq { { { \bf{pr}}\left[{\mathcal{a}_4|\bar{\mathcal{a}_3}}\right ] } } + o(1 ) \end{aligned}\ ] ]    it suffices to show that @xmath90}}=o(1)$ ] . note that by a simple union bound    @xmath89 } } & \leq \sum_{1 \leq s_1 < s_2 \leq t_1 } \sum_{l=-\frac{\sqrt{t}}{f(t)}}^{\frac{\sqrt{t}}{f(t ) } } p_{l , s_1,s_2 } = o\big(t_1 ^ 2 \frac{\sqrt{t}}{f(t ) } \max { p_{l , s_1,s_2 } } \big)\end{aligned}\ ] ]    where @xmath91}}$ ] .
we consider two cases and we show that in both cases @xmath92 .
+ note that at time @xmath42 there exist @xmath93 edges in @xmath94 .
@xmath95 } } \tag{2}\\                 & \leq t_0^{1/6}t^{1/2 }   \sum_{d_1,d_2=3}^{4t_1 } { 2t_0^{1/6}t^{1/2 } \choose d_1 - 1 } { 2t_0^{1/6}t^{1/2 } \choose d_2 - 1 } \big(\frac{2t_0 + 3}{2t+1}\big)^{(d_1+d_2)/2 } e^{\frac{3}{2}+t_1+\frac{2t_0^{1/6}}{3 } } \tag{3}\\                 & \leq t_0^{1/6}t^{1/2 } \sum_{d_1,d_2=3}^{4t_1 }   ( 2t_0^{1/6}t^{1/2})^{d_1+d_2 - 2 }   \big(\frac{2t_0 + 3}{2t+1}\big)^{(d_1+d_2)/2 } e^{2t_1 } \\                 & \leq   t_0^{1/6}t^{1/2 } e^{2t_1 } t_1 ^ 2 ( 2t_0^{1/6}t^{1/2})^{8t_1 - 2 } \big(\frac{2t_0 + 3}{2t+1}\big)^{4t_1 } \\                & = t_0^{4t_1/3 + 1/6}t^{-1/2 } e^{2t_1 } t_1 ^ 2 2^{8t_1 } ( 2t_0 + 3)^{4t_1 } \big ( \frac{t}{2t+1 } \big)^{4t_1 } \\                 & = o\big ( \frac{f(t)}{t_1 ^ 2\sqrt{t } } \big)\end{aligned}\ ] ]    note that we omitted the tedious calculation justifying the transition from ( 2 ) to ( 3 ) since calculating the upper bound of the joint probability distribution is very similar to the calculation of lemma  [ lem : lemma2 ] .
+   + notice that in any case @xmath96 share at most two faces ( which may change over time ) .
note that the two connected vertices @xmath97 share a common face only if @xmath98 .
the other case is treated in the same manner . ] .
consider the following modified process @xmath99
: whenever an incoming vertex `` picks '' one of the two common faces we do nt insert it .
we choose two other faces which are not common to @xmath97 and add one vertex in each of those . notice that the number of faces increases by 1 for both @xmath97 as in the original process and the difference of the degrees remains the same .
an algebraic manipulation similar to case 1 gives the desired result .
in complex network analysis the spectrum of the adjacency matrix is an important aspect of the network or more generally of the network model with several applications , e.g. , triangle counting @xcite . to the best of our knowledge
the spectrum of random apollonian networks has been studied only experimentally @xcite . here
, we prove theorem  [ thrm : thrm2 ] . we decompose the proof of the main theorem in lemmas  [ lem : lemma7 ] ,  [ lem : lemma8 ] ,  [ lem : lemma9 ] ,  [ lem : lemma10 ] .
having computed the highest degrees of a ran in section  [ sec : degree ] , eigenvalues are computed by adapting existing techniques @xcite , see also @xcite for a closely related analysis , and taking into account the special properties of the model .
we decompose the proof of theorem  [ thrm : thrm2 ] in four lemmas .
specifically , in lemmas  [ lem : lemma7 ] ,  [ lem : lemma8 ] we bound the degrees and co - degrees respectively . having these bounds
, we decompose the graph into a star forest and show in lemmas  [ lem : lemma9 ] and  [ lem : lemma10 ] that its largest eigenvalues , which are @xmath100 , dominate the eigenvalues of the remaining graph .
we partition the vertices into three set @xmath101 .
specifically , let @xmath102 be the set of vertices added after time @xmath103 and at or before time @xmath104 where    @xmath105    in the following we use the recursive variational characterization of eigenvalues @xcite .
specifically , let @xmath106 denote the adjacency matrix of a simple , undirected graph @xmath107 and let @xmath108 denote the @xmath109-th largest eigenvalue of @xmath106 . then    @xmath110    where @xmath111 ranges over all @xmath112 dimensional subspaces of @xmath113 .
we shall use the following lemma in our proof , specifically in the proof of lemma  [ lem : lemma9 ] .    for any @xmath114 and
any @xmath17 with @xmath4 as @xmath5 the following holds _ whp _  : for all @xmath46 with @xmath115 , for all vertices @xmath116 , then @xmath117 .
[ lem : lemma7 ]    set @xmath118 .
we use lemma  [ lem : lemma1 ] , a simple union bound and markov s inequality to obtain :    @xmath119 } } & \leq \sum_{s = f(t)}^t \sum_{r=1}^s   { { { \bf{pr}}\left [ { d_s(r)^{(q ) } \geq ( s^{1/2+\epsilon } r^{-1/2})^{(q ) }   } \right]}}\\ & \leq \sum_{s = f(t)}^t \sum_{r=1}^s   { { { \bf{pr}}\left [ { d_s(r)^{(q ) } \geq ( s^{-(q/2+q\epsilon ) } r^{q/2 } ) } \right ] } } \\   & \leq \sum_{s = f(t)}^t \sum_{r=1}^s   \frac{(q+2)!}{2 } \big ( \frac{2s}{r } \big)^{q/2 } s^{-q/2 } s^{-q\epsilon } r^{q/2 }   \\ & = \frac{(q+2)!}{2 } 2^{q/2 } \sum_{s = f(t)}^t s^{1-q\epsilon } \\ & \leq \frac{(q+2)!}{2 } 2^{q/2 }    \int_{f(t)-1}^t \ !
x^{1-q\epsilon }   \ , \mathrm{d}x     \\ & \leq \frac{(q+2)!}{2(q\epsilon-2 ) } 2^{q/2 } ( f(t)-1)^{2-q\epsilon } = o(1).\end{aligned}\ ] ]    let @xmath120 be the set of vertices in @xmath121 which are adjacent to more than one vertex of @xmath122
. then @xmath123 _ whp _  .
[ lem : lemma8 ]    first , observe that when vertex @xmath46 is inserted it becomes adjacent to more than one vertex of @xmath122 if the face chosen by @xmath46 has at least two vertices in @xmath122 .
we call the latter property @xmath124 and we write @xmath125 when @xmath46 satisfies it . at time
@xmath42 there exist @xmath126 faces total , which consist of faces whose three vertices are all from @xmath122 . at time
@xmath127 there can be at most @xmath128 faces with at least two vertices in @xmath122 since each of the original @xmath126 faces can give rise to at most 3 new faces with at least two vertices in @xmath129 .
consider a vertex @xmath130 , i.e. , @xmath127 . by the above argument , @xmath131 } } \leq \frac{6t_1 + 3}{2t+1}$ ] .
writing @xmath132 as a sum of indicator variables , i.e. , @xmath133 and taking the expectation we obtain    @xmath134 } & \leq \sum_{s = t_2}^t \frac{6t_1 + 3}{2t+1 } \leq ( 6t_1 + 3 ) \int_{t_2}^t \ !
( 2x+1)^{-1 } \ , \mathrm{d}x   \\            & \leq ( 3t^{\tfrac{1}{8}}+\tfrac{3}{2 } )    \ln{\frac{2t+1}{2t_2 + 1}}= o(t^{1/7})\end{aligned}\ ] ]    by markov s
inequality :    @xmath135 } } \leq \frac { { { \mathbb e}\left[{|s_3'|}\right ] } } { t^{1/6 } }   = o(1).\ ] ]    therefore , we conclude that @xmath136 _ whp _  .
let @xmath137 be the star forest consisting of edges between @xmath122 and @xmath138 .
let @xmath16 denote the @xmath0 highest degrees of @xmath107
. then @xmath139 _ whp _  .
[ lem : lemma9 ]    it suffices to show that @xmath140 for @xmath141 .
note that since the @xmath0 highest vertices are inserted before @xmath42 _ whp _  , the edges they lose are the edges between @xmath122 and the ones incident to @xmath120 and @xmath142 and we know how to bound the cardinalities of all these sets . specifically by lemma  [ lem : lemma8 ] @xmath136 _ whp _  and by theorem  [ thrm : thrm1 ] the maximum degree in @xmath143 is less than @xmath144 , @xmath145 for @xmath146 respectively _ whp _  . also by theorem  [ thrm : thrm1 ] @xmath147 .
hence , we obtain    @xmath148    to complete the proof of theorem  [ thrm : thrm2 ] it suffices to prove that @xmath149 is @xmath150 where @xmath151 .
we prove this in the following lemma . the proof is based on bounding maximum degree of appropriately defined subgraphs using lemma  [ lem : lemma7 ] and then using standard inequalities from spectral graph theory @xcite .
@xmath152 _ whp _  .
[ lem : lemma10 ]    from gershgorin s theorem @xcite the maximum eigenvalue of any graph is bounded by the maximum degree .
we bound the eigenvalues of @xmath153 by bounding the maximum eigenvalues of six different induced subgraphs .
specifically , let @xmath154 $ ] , @xmath155 where @xmath156 $ ] is the subgraph induced by the vertex set @xmath111 and @xmath157 is the subgraph containing only edges with one vertex is @xmath111 and other in @xmath158 .
we use lemma  [ lem : lemma9 ] to bound @xmath159 and lemma  [ lem : lemma8 ] for the other eigenvalues .
we set @xmath160 .
@xmath161 @xmath162 @xmath163 @xmath164 @xmath165 @xmath166    therefore _ whp _
we obtain    @xmath167
as shown in appendix b of @xcite using physicist s methodology the diameter of a ran is asymptotically upper bounded by a logarithmic factor @xmath168 .
similarly , zhang et al .
@xcite show that the diameter scales logarithmically . in the following we give a simple proof that the diameter @xmath10 of the graph @xmath15 created @xmath2 steps is @xmath169 _ whp _  and in theorem  [ thrm : ternary ] we give a refined upper bound for the diameter .    .
each face is labelled with its depth.,scaledwidth=30.0% ]    we begin with a necessary definition for the proof of claim  [ claim : diameter ] .
we define the _ depth of a face _ recursively .
initially we have one face @xmath170 its depth is @xmath171 . for each new face
@xmath172 created by picking a face @xmath173 , we have @xmath174 . an example is shown in figure  [ fig : depth ] , where each face is labelled with its corresponding depth .
the diameter @xmath10 of the graph @xmath15 created @xmath2 steps is @xmath169 _ whp _  .
[ claim : diameter ]    note that if @xmath175 is the maximum depth of a face then @xmath176 .
hence , we need to upper bound the depth of a given face after @xmath2 rounds .
let @xmath177 be the number of faces of depth @xmath0 at time @xmath2 , then :    @xmath178 } & = \sum_{1 \leq t_1<t_2<\ldots
< t_k\leq t } \prod_{j=1}^k \frac{1}{2t_j+1 } \leq \frac{1}{k ! } ( \sum_{j=1}^t \frac{1}{2j+1 } ) ^k \leq \frac{1}{k ! } ( \frac{1}{2 } \log{t})^k \leq    ( \frac{e\log{t}}{2k})^{k+1}\end{aligned}\ ] ]    by the first moment method we obtain @xmath179 _ whp _  and hence @xmath180 _ whp _  .        the depth of a face can be formalized via a bijection between random ternary trees and rans . using this bijection we prove theorem  [ thrm : ternary ] which gives a refined upper bound on the asymptotic growth of the diameter .
consider the random process which starts with a single vertex tree and at every step picks a random leaf and adds three children to it .
let @xmath158 be the resulting tree after @xmath2 steps .
there exists a natural bijection between the ran process and this process , see @xcite and also figure  [ fig : ternary ] .
the depth of @xmath158 in probability is @xmath181 where @xmath12 is the unique solution greater than 1 of the equation @xmath13 , see broutin and devroye @xcite , pp .
284 - 285 .
note that the diameter @xmath10 is at most twice the height of the tree and hence the result follows .
the above observation , i.e. , the bijection between rans and random ternary trees can not be used to lower bound the diameter . a counterexample is shown in figure  [ fig : counter ] where the height of the random ternary tree can be made arbitrarily large but the diameter is 2 .
albenque and marckert proved in @xcite that if @xmath182 are two i.i.d .
uniformly random internal vertices , i.e. , @xmath183 , then the distance @xmath184 tends to @xmath185 with probability 1 as the number of vertices @xmath186 of the ran grows to infinity .
however , an exact expression of the asymptotic growth of the diameter to the best of our knowledge remains an open problem .
finally , it is worth mentioning that the diameter of the ran grows faster asymptotically than the diameter of the classic preferential attachment model @xcite which _ whp _   grows as @xmath187 , see bollobs and riordan @xcite .
consider the three faces created after the insertion of the first point .
let s call the face which gets the first vertex @xmath188 and the other two faces @xmath189 , see also figure  [ fig : fig1](b ) . also for simplicity of the notation let the configuration of figure  [ fig : fig1](b ) correspond to time @xmath190 .
let @xmath191 equal the number of steps until a new vertex picks face @xmath192 or @xmath193 .
clearly , @xmath194 . what is the expectation @xmath195}$ ] ? for any @xmath196    @xmath197 } } & = \prod_{j=1}^t \frac{3 + 2(j-1)}{5 + 2(j-1 ) } = \frac{3}{2t+3}.\end{aligned}\ ] ]    using now the fact that @xmath195 } = \sum_{t=1}^{+\infty } { { { \bf{pr}}\left[{x \geq t}\right]}}= 1 + \sum_{t=1}^{+\infty } { { { \bf{pr}}\left[{x > t}\right]}}$ ] and substituting we obtain that @xmath195 } = + \infty$ ] .
in this work we studied several aspects of random apollonian networks , namely the highest degree vertices , the spectrum , the diameter and the `` waiting times '' .
there are various research directions which we plan to address in future work .
indicatively , we report the following : ( a ) in contrast to the preferential attachment model , ran have small separators .
specifically , by the lipton - tarjan theorem @xcite we know the existence of a good separator of size @xmath198 .
this is in accordance with the empirical observation that real world networks have good separators @xcite .
we plan to investigative the minimum size separator in future work .
( b ) consider the natural process where an adversary with probability @xmath170 at every step deletes a vertex from the network , e.g. , uniformly at random or in a maliscious way .
under which conditions do we obtain a connected ran ?
( c ) what is the exact asymptotic expression for the diameter of @xmath15 ?
we would like to thank luc devroye and alexis darrasse for pointing out references @xcite and @xcite respectively .
99    aiello , w. , chung , f. , lu , l. : _ a random graph model for power law graphs _ experimental mathemathics , 10(1 ) , 53 - 66 , 2001    albenque , m. , marckert , j.f . :
_ some families of increasing planar maps _ electronic journal of probability , 13 , pp .
1624 - 1671    albert , r. , barabsi , a. , jeong , h. : _ diameter of the world wide web .
_ nature 401 , pp .
103131 , 1999    andrade , j.s .
, herrmann , h.j . , andrade , r.f.s .
, da silva , l.r . : _
apollonian networks : simultaneously scale - free , small world euclidean , space filling , and with matching graphs _ phys .
94 ( 2005 ) 018702 .
andrade , r.f.s . ,
miranda , j.g.v . : _
spectral properties of the apollonian network _
physica a , 356 , ( 2005 )    aumann , r. : _ acceptable points in general cooperative @xmath186-person games _ contributions to the theory of games iv , pp .
287 - 324 , princeton university press , 1959    aumann , r. , myerson , r. : _ endogenous information of links between players and coalition : an application to the shapley value _ shapley value , pp .
175 - 191 , cambridge university press , 1988    bala , v. , goyal , s. : _ a noncooperative model of network formation _ econometrica , econometric society , vol .
68(5 ) , pp . 1181 - 1230 , 2000    barabsi , a. , albert , r. : _ emergence of scaling in random networks _ science 286 , pp . 509512 , 1999    barthlemy , m. , flammini , a. : _ modeling urban street patterns _ physical review letters , 100(13 ) , 138702(4 ) , 2008    blandford , d. , blelloch , g. , kash , i. : _ compact representations of separable graphs _
proceedings of the fourteenth annual acm - siam symposium on discrete algorithms ( soda 03 ) , pp .
679 - 688 , 2003    bodini , o. , darrasse , a. , soria , m. : _ distances in random apollonian network structures _
available at arxiv http://arxiv.org/abs/0712.2129    bollobs , b. , riordan , o. : _ the diameter of a scale - free random graph _
combinatorica , 2002    bonato , a. , hadi , n , horn , p. , pralat , p. : _ models of online social networks _
internet mathematics , to appear ( 2011 )    boorman , s. : _ a combinatorial optimization model of transmission of job information through contact networks _ bell journal of economics , 6:216 - 249 , 1975    boyd , d.w . : _ the sequence of radii of the apollonian packing _ mathematics of computation , vol .
249 - 254 , 1982    broder , a. , kumar , r. , maghoul , f. , raghavan , p. , rajagopalan , s. , stata , r. , tomkins , a. , wiener , j. : _ graph structure in the web _ in proc . of the 9th intl .
world wide web conference , pp .
amsterdam : north - holland publishing company , 2000    chung graham , f. : _ spectral graph theory _ american mathematical society 1997    chung graham , f. , lu , l. : _ complex graphs and networks _ american mathematical society , no . 107 ( 2006 )
chung , f. , lu , l , vu . v.h . : _ spectra of random graphs with given expected degrees
_ proceedings of the national academy of sciences of the united states of america , 100 , 63136318    cooper , c. , frieze , a. : _ a general model of web graphs _
random structures & algorithms , volume 22 issue 3,pp . 311335 , 2003    cooper , c. , uehara , r. : _ scale free properties of random @xmath0-trees _ mathematics in computer science , 3(4 ) , pp .
489496 , 2010    darrasse , a. , soria , m. : _ degree distribution of random apollonian network structures and boltzmann sampling _ 2007 conference on analysis of algorithms , aofa 07 , dmtcs proceedings .
darrasse , a. , hwang , h .- k . ,
bodini , o. , soria , m. : _ the connectivity - profile of random increasing k - trees _ available at arxiv : http://arxiv.org/abs/0910.3639    broutin , n. , devroye , l. : _ large deviations for the weighted height of an extended class of trees _
algorithmica , 46 , pp .
271 - 297 , 2006    doye , j.p.k , massen , c.p . :
_ self - similar disk packings as model spatial scale - free networks _
e 71 ( 2005 ) 016128 .
eisenstat , d. : _ random road networks : the quadtree model _
available at arxiv    erds , p. , rnyi , a. : _ on random graphs i _
publicationes mathematicae debrecen 6 , pp .
290 - 297 , 1959    erds , p. , rnyi , a. : _ on the evolution of random graphs _ magyar tud .
akad . mat .
kutat int .
kzl 5 , pp .
17 - 61 , 1960    fabrikant , a. , koutsoupias , e. , papadimitriou , c. : _ heuristically optimized trade - offs _ proceedings of the 29th international colloquium on automata , languages , and programming ( icalp ) , pp . 110122 , 2002    faloutsos , m. , faloutsos , p. , faloutsos , c. : _ on power - law relationships of the internet topology _ proceedings of the conference on applications , technologies , architectures , and protocols for computer communication ( sigcomm ) , 1999    flaxman , a. , frieze , a. , fenner , t. : _ high degree vertices and eigenvalues in the preferential attachment graph _
internet mathematics , 2(1 ) , 2005    gao , y. : _ the degree distribution of random k - trees _ theoretical computer science , 410(8 - 10 ) , 2009 .
gao , y. , hobson , c. : _ .
random k - tree as a model for complex networks _ in the fourth workshop on algorithms and models for the web - graph ( waw2006 ) , 2006    gilbert , e.n .
: _ random graphs _
annals of mathematical statistics , 30 , pp .
1141 - 1144 , 1959    graham , r.l . , lagarias , j.c .
, mallows , c.l . , wilks , a.r .
, yan , c.h . : _
apollonian circle packings : number theory _
j. number theory 100 , no . 1 , 1 - 45 , mr1971245    guy , s. , chhugani , j. , curtis , s. , dubey , p. , lin , m. , manocha , d. : _ pledestrians : a least - effort approach to crowd simulation _ eurographics / acm siggraph symposium on computer animation , 2010 .    hayashi , y. , matsukubo , j. : _ a review of recent studies of geographical scale - free networks _
available at arxiv : cond / physics/0512011    kloks , t. : _ treewidth : computations and approximations _
springer - verlag , 1994    kolountzakis , m.n , miller , g.l . , peng , r. , tsourakakis , c.e . : _
efficient triangle counting in large graphs via degree - based vertex partitioning _
7th workshop on algorithms and models for the web graph ( waw 10 ) , 2010    lipton , r.j . ,
tarjan , r. : _ a separator theorem for planar graphs _ siam journal on applied mathematics 36 : 177189 , 1979    leskovec , j. , faloutsos , c. : _ scalable modeling of real graphs using kronecker multiplication _ machine learning proceedings of the twenty - fourth international conference ( icml 2007 ) , corvalis , oregon , usa , june 20 - 24 , 2007    masucci , a.p . , smith , d. , crooks , a. , batty , m. : _ random planar graphs and the london street network _
european physics journal b ( 2009 )    mihail , m. , papadimitriou , c. : _ on the eigenvalue power law _ random 02 , pp .
254 - 262 ( 2002 )    panholzer , a. , seitz , g. : _ ordered increasing @xmath0-trees : introduction and analysis of a preferential attachment network model _ dmtcs proc . ,
aofa10 , pp .
549 - 564 , 2010    pralat , p. , wormald , n. : _ growing protean graphs _
internet mathematics , 4(1 ) , pp . 1 - 16 , 2007    strang , g. : _ linear algebra and its applications _ brooks cole , 2005    tsourakakis , c.e . : _
fast counting of triangles in large real networks without counting : algorithms and laws _ eighth ieee international conference on data mining ( icdm 08 ) , 2008    watts , d. , strogatz , s. : _ collective dynamics of small - world networks _ nature , 1998 , * 393 * , 440 - 442    wu , z .- x .
, xu , x .- j . ,
wang , y .- h .
: _ comment on `` maximal planar networks with large clustering coefficient and power - law degree distribution '' _ physical review , e 73 , 058101 ( 2006 )    zhang , z.z . ,
comellas , f. , fertin , g. , rong , l.l . : _
high dimensional apollonian networks _ available at arxiv http://arxiv.org/abs/cond-mat/0503316    zhou , t. , yan , g. , wang , b.h . : _
maximal planar networks with large clustering coefficient and power - law degree distribution _ phys .
e 71 ( 2005 ) 046141 .
let @xmath199 be a vector denoting that @xmath58 increases by 1 at @xmath200 for @xmath201 .
we upper bound the probability @xmath202 of this event in the following.note that we consider the case where the vertices have same degree as the number of faces around them , figure  [ fig : fig2](b ) .
the other case is analyzed in exactly the same way , modulo a negligible error term .
@xmath203 \bigg [ \prod_{k=0}^r \prod_{j=\tau_k+1}^{\tau_{k+1}-1 } \big(1-\frac{d_0+k}{2j+1}\big ) \bigg ] \\            & = d_0 ( d_0 + 1 ) \ldots ( d_0+r-1)\bigg [ \prod_{k=1}^r \frac{1}{2\tau_k+1 } \bigg ]           \exp{\bigg ( \sum_{k=0}^r \sum_{j=\tau_k+1}^{\tau_{k+1}-1 } \log{\big(1-\frac{d_0+k}{2j+1}\big ) }   \bigg ) } \\           & = \frac{(d_0+r-1)!}{(d_0 - 1)!}\bigg [ \prod_{k=1}^r \frac{1}{2\tau_k+1 } \bigg ]           \exp{\bigg ( \sum_{k=0}^r \sum_{j=\tau_k+1}^{\tau_{k+1}-1 } \log{\big(1-\frac{d_0+k}{2j+1}\big ) }   \bigg ) } \\\end{aligned}\ ] ]                                          @xmath225 \\    & \leq c(r , d_0,t_0,t ) \sum_{t_0-\lceil \tfrac{d_0}{2 } \rceil + 1 \leq \tau_1 ' \leq .. \leq \tau_r ' \leq t-\lceil \tfrac{d_0+r}{2 } \rceil }   \prod_{k=1}^r   \big [ ( 2\tau_k'+3)^{-1/2 } \big ( 1+\frac{1}{\tau_k'+\frac{d_0+k}{2 } + 1/2}\big)\big ] \\
& \leq \frac{c(r , d_0,t_0,t)}{r ! } \bigg ( \sum_{t_0-\lceil \tfrac{d_0}{2 } \rceil}^{t-\lceil \tfrac{d_0+r}{2 } \rceil } ( 2\tau_k'+3)^{-1/2}+\tfrac{1}{\sqrt{2 } } ( \tau_k'+3/2)^{-3/2 } \bigg)^r \\    & \leq \frac{c(r , d_0,t_0,t)}{r ! }   \bigg ( \int_{0}^{t-\frac{d+r}{2 } } \ ! \big [   ( 2x+3)^{-1/2 } + \tfrac{1}{\sqrt{2 } } ( x+3/2)^{-3/2 } \big ] \ , \mathrm{d}x\bigg)^r \\     & \leq \frac{c(r , d_0,t_0,t)}{r ! } \big(\sqrt{2t+3-(d_0+r)}+2/3\big)^{r }   \\     & \leq \frac{c(r , d_0,t_0,t)}{r ! } ( 2t)^{r/2 } e^{-\tfrac{r}{2 } \tfrac{d_0+r-3}{2 t } } e^{\frac{2r}{3\sqrt{2t-(d_0+r)+3 } } } \\    & \leq { d_0+r-1 \choose d_0 - 1 } \big(\frac{2t_0 + 3}{2t+1}\big)^{d_0/2 } \big [ \big ( 1-\frac{d_0}{2t_0 + 3}\big)^{-(1-\frac{d_0}{2t_0 + 3 } ) } \big]^{t_0 + 3/2 } \times \\     & \big(\frac{2t}{2t+1}\big)^{r/2 } \exp{\big ( -\frac{d_0}{2 } + \frac{(d_0+r)^2}{8t+4}+ \frac{(d_0+r)^3}{4(2t+1)^2 } - \frac{r(d_0+r-3)}{4t}+\frac{2r}{3\sqrt{2t+3-(d_0+r ) } }   \big ) } \\ \end{aligned}\ ] ]