Instruction: Write an abstract for the following acadamic paper.

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

Paper: band structures are of fundamental interest to the solid - state physicist as they reveal important properties of crystals . @xcite the determination of band structures on _ ab initio _ level is therefore an important issue in theoretical solid state physics . nowadays predominantly density functional theory  ( dft )  @xcite with the local density approximation  ( lda ) or its variants are employed to treat this problem . @xcite dft focuses on the ground - state electron density and yields ground - state properties . in several cases , the interpretation of the kohn - sham orbital energies as quasiparticle energies also turns out to be successful and particular attention has been paid to the  lda , due to its numerical simplicity . however , in insulating materials , the lda tends to underestimate the band gap , e.g. , ref .  . as is well known , it is difficult to improve band structure calculations based on density functional theory . this is in particular the case when electron correlations are strong such that in quite some cases  @xcite the one - electron picture becomes inappropriate and corrections are needed . nevertheless , a number of ingenious methods have been devised to improve band structure results based on density functional theory . for example , the calculations have been supplemented by a gw  treatment , e.g. , refs  , or the so - called lda+@xmath1 method . @xcite other improvements concern the use of optimized effective potentials  @xcite such as the exact exchange potential  @xcite instead of its local , simplified form . time - dependent dft is another route to treat excited states more rigorously . furthermore , lda calculations have been coupled with various extensions of the coherent potential approximation  ( cpa ) like the dynamical cpa  @xcite or dynamical mean field theory  ( dmft ) . @xcite also faddeev s method  @xcite of treating the three - particle @xmath2-matrix has been applied in conjunction with lda band structure calculations . @xcite the same holds true for different forms of the projection operator technique . @xcite yet one must face the fact that the various approximations remain uncontrolled and therefore have to be reconsidered from case to case . when one is aiming at controlled approximations , two different routes offer themselves . one is bound to quantum monte carlo  ( qmc ) calculations , e.g. , ref .  . this is a simple and straightforward way of dealing with the many - body problem . yet it suffers from one serious shortcoming , the so - called sign problem when treating fermions . therefore , we are pursuing the second route . it starts from hartree - fock band structures and improves them by including correlations by means of many - particle techniques similar to the ones employed for molecules . @xcite thereby one can take advantage of program packages like wannier which has been written by shukla  _ et al . _ @xcite or crystal which originates from torino . @xcite both provide hartree - fock bands as well as wannier orbitals for crystals . these orbitals are especially suited for an _ ab initio _ treatment of electron correlations because of their local character .  @xcite    a number of correlation calculations have been performed . @xcite a local hamiltonian approach  @xcite was shown to improve hartree - fock energy bands substantially . suhai used toyozawa s electronic polaron model  @xcite to repartition mp2 pair energies to estimate quasiparticle band structures . @xcite by inserting the orbital energies into the energy - dependent self - energy , this model was shown to be a special case of the outer valence green s functions  ( ovgfs )  @xcite that were derived in terms of crystal orbitals by liegener . @xcite by applying the approximation to the self - energy of igarashi  _ et al . _ ,  @xcite quasiparticle band structures were obtained by albrecht  _ et al._.  @xcite    in the present paper , we derive the algebraic diagrammatic construction  ( adc ) scheme  @xcite for crystals . the adc scheme has proven to be superior to the ovgf method in molecular studies ,  @xcite and numerous works have been carried out over the last two decades , e.g. , refs .  , including studies of oligomers and clusters chosen to model infinite chains or crystals . @xcite the adc scheme is a method to approximate the feynman - dyson perturbation series for the self - energy and contains sums of certain proper and improper diagrams to infinite order . @xcite the basic properties of the adc scheme and their derivations , as applied to molecules , are found in refs .  . among those size - extensivity is an important one since it is crucial when solids are considered . furthermore , the adc method is known to be robust and facilitates also to study strong electron correlations due to the efficient and stable evaluation of the one - particle green s function in terms of a hermitian eigenvalue problem . in molecules strong correlations for instance occur when inner valence electrons are treated . @xcite with the crystal orbital adc ( co - adc ) devised here , the spectral representation of the one - particle green s function is obtained . we suggest a formulation of the expressions in wannier representation instead of crystal momentum representation in order to be able to exploit the local character of the correlation hole around an electron especially in the case of crystals with a large unit cell .    as a first application of the co - adc method to a three - dimensional crystalline solid , we have chosen lithium fluoride which occurs in nature as the mineral griceite . lif not only has a wide range of technological applications like in x - ray monochromators or in filters for ultraviolet radiation , e.g. , ref . ( and references . therein ) , but is also interesting for a number of fundamental physical reasons . it is considered to be the `` most ionic substance '' and a prototypical insulator which manifests in its very large fundamental band gap of@xmath3 . @xcite piacentini  @xcite estimates the band gap of lif to be  @xmath4 . this value is refined in ref . to  @xmath5 . shirley  _ et al . _ @xcite communicate  @xmath6 from the dissertation of himpsel . ] @xmath6 that is the largest one found in nature apart from exotic systems . some authors even consider  lif to be comparable to a he  ne rare - gas solid . @xcite its optical spectra show strong excitonic effects , which complicates the experimental determination of the band gap  @xcite by optical spectroscopy . similarly , many - particle effects have to be accounted for in the measurement of the widths of the f@xmath0 and f@xmath7  valence bands . lif has been thoroughly studied both experimentally and theoretically . @xcite poole  _ et al . _ survey early experimental and theoretical data . @xcite recent studies of lif comprise density functional theory calculations at the lda level supplemented by an inclusion of correlation effects in terms of the gw  approximation . @xcite _ ab initio _ investigations comprise hartree - fock studies of ground - state properties  @xcite and of the band structure . @xcite an accurate treatment of electron correlations for the ground - state properties of  lif ( and other alkali halides ) has been carried out by doll and stoll . @xcite quasiparticle band structures have been obtained by kunz  @xcite on the basis of the electronic polaron model  @xcite and by albrecht  @xcite with the green s function approach of igarashi  _ et al._.  @xcite    the article is structured as follows : section  [ sec : band ] puts forward an _ ab initio _ foundation for band structures based on the one - particle green s function . in secs . [ sec : dynamic ] and [ sec : static ] , we derive the co - adc approximation to the self - energy in terms of wannier orbitals . section  [ sec : poledyson ] combines the results of the previous sections to formulate the band structure problem as a hermitian matrix eigenvalue problem . criteria for the truncation of lattice sums are discussed in sec . [ sec : onion ] . in sec . [ sec : lif_crystal ] , the application of co - adc theory to lithium fluoride is described . conclusions are drawn in sec . [ sec : conclusion ] . the nonrelativistic hamiltonian of a crystal in atomic units  @xcite reads @xmath8 \\      & & { } + \frac{1}{2 } \ , { \sum\limits}_{\scriptstyle m , n=1 \atop \scriptstyle m \neq n}^n \frac{1}{|\vec r_m -        \vec r_n| } \\      & & { } \stackrel{\textstyle + \frac{1}{2 } \ , { \sum\limits}_{i , j=1}^{n_0 } { \sum\limits}_{a , b=1}^m        \frac{z_{\vec r_i a } z_{\vec r_j b}}{|\vec r_{\vec r_i a } -        \vec r_{\vec r_j b}| } } { \qquad\scriptstyle i \neq j \lor a \neq b \hfill } \ ; ,    \end{array}\ ] ] where fixed nuclei are assumed . here @xmath9  is the number of electrons in the crystal , @xmath10  denotes the number of unit cells , and @xmath11  is the number of nuclei per unit cell . @xmath12  stands for the charge of nucleus  @xmath13 in unit cell  @xmath14 , and @xmath15  is the distance between the @xmath16-th electron and the @xmath13-th nucleus in unit cell  @xmath17 . finally , @xmath18  represents the distance between the @xmath19-th and @xmath16-th electron and @xmath20  denotes the distance between the nuclei  @xmath13 , @xmath21 of charge  @xmath12 , @xmath22 in unit cells  @xmath14 , @xmath17 . the schrdinger equation with the hamiltonian  ( [ eq : operator_hamilton ] ) can be solved for the ground state in a restricted , closed - shell hartree - fock approximation ,  @xcite yielding a set of spin bloch orbitals  @xmath23 where @xmath24  is the crystal momentum and @xmath25  is a band index . @xcite we use these orbitals to represent the hamiltonian  ( [ eq : operator_hamilton ] ) in second quantization , yielding the so - called bloch or crystal momentum representation . @xcite it is decomposed according to mller and plesset  @xcite into the hartree - fock part  @xmath26 and the residual interaction  @xmath27 , a suitable form to apply perturbation theory :    @xmath28    the energy bands ( or bloch orbital energies ) are denoted by  @xmath29 , and the operators  @xmath30 ( @xmath31 ) create ( annihilate ) electrons in  @xmath23 . the negative of the hartree - fock potential is given by  @xmath32 } \ > n_{\vec k \ , i}$ ] , the occupation numbers  @xmath33 being unity for occupied and zero for unoccupied ( or virtual ) bloch orbitals ( and @xmath34 ) . generally , we use the band indices  @xmath35 , @xmath36 , @xmath19 , @xmath16 ,  to denote occupied hartree - fock bands , @xmath37 , @xmath38 , @xmath39 , @xmath40 ,  for unoccupied bands , and @xmath25 , @xmath41 , @xmath42 , @xmath43 ,  for bands which are occupied or unoccupied . the two - electron integrals in eq . ( [ eq : bloch_hamilton ] ) are defined with respect to bloch orbitals by  @xcite @xmath44 the antisymmetrized two - electron integrals are @xmath45 } = v_{\vec k_1 \ , p \ ; \vec k_2 \ , q \ ; \vec k_3 \ , r \ ; \vec        k_4 \ , s } - v_{\vec k_1 \ , p \ ; \vec k_2 \ , q \ ; \vec k_4 \ , s \ ; \vec        k_3 \ , r}\ ] ] and consist of a coulomb term  @xmath46 and an exchange term  @xmath47 . all crystal momentum vectors in eqs . ( [ eq : bloch_hamilton ] ) , ( [ eq : twoint ] ) , and ( [ eq : antitwoint ] ) are restricted to the first brillouin zone . furthermore , we assume that the first brillouin zone is discretized by the born von krmn boundary conditions , or equivalently , that numerical integrations in reciprocal space are carried out in terms of a monkhorst - pack net . @xcite in both cases , integrations over the brillouin zone are approximated by finite sums . the one - particle green s function ( or particle propagator ) in terms of bloch orbitals  @xcite reads @xmath48 { \left|\right.\!\psi_0^n\!\left.\right > } \ ; , \ ] ] with wick s time - ordering operator  @xmath49 . @xmath50  is the exact ground state of the @xmath9  particle system . due to translational symmetry of the hamiltonian  ( [ eq : operator_hamilton ] ) , the one - particle green s function depends only on one crystal momentum ,  @xcite i.e. , @xmath51 . similarly , the time independence of the hamiltonian  ( [ eq : operator_hamilton ] ) causes the particle propagator  ( [ eq : k - prop ] ) to depend only on the time difference  @xmath52 . @xcite fourier transforming  @xmath53 with respect to  @xmath52 yields the one - particle green s function in energy space  @xmath54 which can be recast in terms of the spectral or lehmann representation as  @xcite    [ eq : p - spectrum ] g_pq(k , ) & = & _ n \{n+1 } + & & + _ n \{n-1 } + & = & g_pq^+(k , ) + g_pq^-(k , ) . here @xmath16  numerates the excited states of the @xmath55  particle system . the negative of the pole positions in eq . ( [ eq : p - spectral ] ) is given by either the electron affinities  @xmath56 or the ionization potentials @xmath57 where @xmath58  are the energies of the different excited states of the @xmath59  particle system . the summands  @xmath60 are required to ensure the convergence of the fourier transformation . whenever such a factor occurs , we implicitly take the limit  @xmath61 . @xcite    the pole strengths in eq . ( [ eq : p - spectral ] ) are given by the transition amplitudes  @xcite @xmath62 where @xmath63  denote excited states of the @xmath59  particle system with crystal momenta  @xmath64 . the pole strengths can be interpreted in terms of the spectral intensities observed in photoelectron spectroscopy experiments , similarly to the molecular case discussed in ref .  . to obtain electron affinities and ionization potentials from the band structure of three - dimensional crystals , one has to specify the energy of the added ( removed ) electron at the fermi level of the neutral crystal . @xcite this is usually done by introducing chemical potentials  @xmath65 which are added to the pole positions of the one - particle green s function , i.e. , @xmath66 and @xmath67 .    the ( fundamental ) band gap is the smallest difference between the energies for removing an electron from and attaching an electron to an @xmath9  particle system : @xmath68 here @xmath13  stands for the largest electron affinity and @xmath69  designates its smallest ionization potential . correspondingly , @xmath70  denote the energies of the electron attachment ( removal ) state with the lowest ( highest ) energy and @xmath71  represents the ground - state energy . the self - energy  @xmath72 with respect to the residual interaction  @xmath27 is defined by the dyson equation  @xcite    [ eq : dyson ] g(k , ) = g^0(k , ) & + & g^0(k , ) ( k , ) g(k , ) + = g^0(k , ) & + & g^0(k , ) ( k , ) g^0(k , )   + & + &  ,    which can be solved formally by @xmath73^{-1 } \ ; , \ ] ] where @xmath74\ ] ] is the free green s function with respect to the mller - plesset partition  ( [ eq : bloch_hamilton ] ) . in fact , for crystals with a band gap , the occupation numbers in eq . ( [ eq : freeprop ] ) are independent of  @xmath24 , i.e. , @xmath75 and @xmath76 . the self - energy can be decomposed into an @xmath77  independent part , the static self - energy  @xmath78 , and an @xmath77  dependent part , the dynamic self - energy  @xmath79 ,  @xcite @xmath80 where @xmath81 . the dynamic self - energy is considered first in the ensuing sec . [ sec : dynamic ] as the static self - energy is determined using the former one in sec . [ sec : static ] . the partition  ( [ eq : bloch_hamilton ] ) of the hamiltonian has the same form and meaning as the mller - plesset partition for canonical molecular orbitals  @xcite used frequently in molecular physics . therefore , all equations derived on the basis of the molecular mller - plesset partition are in complete analogy to the equations in the crystalline case . the only difference between molecules and crystals is the occurrence of composite indices in the hamiltonian  ( [ eq : bloch_hamilton ] ) , consisting of a crystal momentum vector  @xmath24 and a band index  @xmath25 which vary independently . given an equation in terms of molecular orbitals , it can be written immediately in terms of bloch orbitals by replacing all molecular orbital indices by composite indices . afterwards one can exploit ( translational ) symmetry to simplify the equation . the dynamic self - energy  @xmath82 is represented in terms of the @xmath83/@xmath84-propagator ( two - particle - one - hole / two - hole - one - particle - propagator ) and possesses the spectral representation  @xcite @xmath85 where @xmath86  denote the pole positions of the @xmath83/@xmath84-propagator and the @xmath87  are termed dyson amplitudes . the retarded dynamic self - energy  @xmath88 and the advanced dynamic self - energy  @xmath89 are associated with excitations of the @xmath59  particle system , respectively . yet the pole positions  @xmath86 do _ not _ directly correspond to the energies of physical states as the ( dynamic ) self - energy is only defined in conjunction with the dyson equation  ( [ eq : dyson ] ) . the algebraic diagrammatic construction scheme is a stable and efficient method to determine the spectral representation of the dynamic self - energy  ( [ eq:2h1p - spectral ] ) . the construction starts by making the ansatz @xmath90 which is termed general algebraic form or adc  form . the vectors  @xmath91 are called modified coupling amplitudes for the crystal orbital  @xmath92 . both @xmath93 and @xmath94 are hermitian matrices , @xmath94 being referred to as modified interaction matrices  @xcite while @xmath93 are assumed to be diagonal . explicit expressions for these matrices are given below in eq . ( [ eq : k - co - adc2 ] ) . the adc  form  ( [ eq : adc - form ] ) reproduces the analytic structure of the feynman - dyson perturbation series for the dynamic self - energy as can be seen from its expansion into a geometric series : @xmath95^n \nonumber \\    \label{eq : adcanalytic }    & & { } \times ( \omega \ , { { \hbox{\boldmath{$\mathbbm{1}$}\unboldmath}}}- { \hbox{\boldmath{$k$}\unboldmath}}^{\pm}(\vec k))^{-1 }      \ , \vec u_q^{\pm}(\vec k ) \ ; .\end{aligned}\ ] ] the @xmath16-th order approximation to the adc  form  ( [ eq : adc - form ] ) is constructed by inserting the perturbation expansions @xmath96 in terms of the residual interaction  @xmath97 into eq . ( [ eq : adcsuminfinity ] ) , and comparing the resulting expressions with the diagrammatic expansion of the dynamic self - energy up to @xmath16-th order . this constitutes a scheme which is denoted by  adc(@xmath16 ) . note that the analytic structure imposed by eq . ( [ eq : adcanalytic ] ) occasionally requires one to associate linear combinations of the analytic expression of several diagrams with particular terms in the above expansion .    in adc(2 ) and adc(3 ) , each entry of the  @xmath91 and the @xmath98 in eq . ( [ eq : adc - form ] ) is characterized by one or two , respectively , arrangements of two particles and one hole  ( @xmath83 ) for  @xmath99 and two holes and one particle  ( @xmath84 ) for  @xmath89 . in brief notation these arrangements are referred to as @xmath83- and @xmath84-configurations . all such configurations that can be formed with the one - particle basis set underlying a hartree - fock calculation constitute the configuration space in adc(2 ) and adc(3 ) . the configuration space is enlarged in adc(4 ) by @xmath100- and @xmath101-configurations . @xcite in general , every second order in adc , the configuration space enlarges by the next higher excitation class . restricting our discussion to @xmath83- and @xmath84-configurations , we can rewrite eq . ( [ eq : adc - form ] ) as    @xmath102    the prime on the summation symbol indicates that the summation only runs over indices with @xmath103 and @xmath104 to avoid double counting of contributions . the corresponding ansatz for  @xmath105 is formally identical to eq . ( [ eq : adc - form - expanded ] ) apart from the changed occupation numbers of the band indices  @xmath106 , @xmath107 , @xmath108 , @xmath109 , @xmath110 , and @xmath111 . the summation variables  @xmath112 , @xmath113 , @xmath114 , and @xmath115 are crystal momentum vectors like  @xmath24 . the quantities @xmath116 , @xmath117 , @xmath93 , and @xmath94 are invariant under translations by an arbitrary lattice vector , a property which is taken into account for the internal summations in the adc  form  ( [ eq : adc - form - expanded ] ) . explicit molecular adc equations up to fourth order are given in ref .  . here we show only the adc equations for the dynamic self - energy in terms of crystal ( bloch ) orbitals up to second order  [ co - adc(2 ) ] :    [ eq : k - co - adc2 ] u^+_k p ; l_1 + l_2 - k i l_1 a l_2 b & = & v^*_k p l_1 + l_2 - k i [ l_1 a l_2 b ] n_l_1 + l_2 - k i |n_l_1 a |n_l_2 b , + k^+_l_1 + l_2 - k i l_1 a l_2 b ; l_1^ + l_2^ - k i^ l_1^ a^ l_2^ b^ & = & _ l_1 + l_2 - k i l_1 a l_2 a ; l_1^ + l_2^ - k i^ l_1^ a^ l_2^ b^ ( _ l_1 a + _ l_2 b - _ l_1 + l_2 - k i ) + & & n_l_1 + l_2 - k i |n_l_1 a |n_l_2 b , + c^+_l_1 + l_2 - k i l_1 a l_2 b ; l_1^ + l_2^ - k i^ l_1^ a^ l_2^ b^ & = & 0 .    the equations for  @xmath118 , @xmath119 , and @xmath120 are formally identical to eq . ( [ eq : k - co - adc2 ] ) apart from the change of labels as for the adaptation of eq . ( [ eq : adc - form - expanded ] ) to  @xmath121 and the occupation numbers in eqs . ( [ eq : k - co - adc2-u ] ) and ( [ eq : k - co - adc2-k ] ) which are in this case  @xmath122 . expression  ( [ eq : k - co - adc2 ] ) is derived in ref .   utilizing a gaussian basis set expansion of the bloch orbitals .    the spectral representation  ( [ eq:2h1p - spectral ] ) of the dynamic self - energy is obtained from the adc  form  ( [ eq : adc - form ] ) by solving the hermitian eigenvalue problem @xmath123 the vector notation in eqs . ( [ eq : adc - form ] ) and ( [ eq : adc - inversion ] ) is a convenient means to sum over all intermediate crystal momenta and band indices . the diagonal matrix  @xmath124 contains the eigenvalues of the secular matrix  @xmath98 , and @xmath125  denotes the eigenvector matrix . the eigenvalues are the pole positions of the dynamic self - energy , i.e. , @xmath126 , while the dyson amplitudes in eq . ( [ eq:2h1p - spectral ] ) are obtained via @xmath127 where the @xmath16-th column vector  @xmath128 of  @xmath125 is used here together with the adjoint of  @xmath129 . note that , unlike co - adc(2 ) , in co - adc(3 ) ( and higher orders ) the eigenvectors  @xmath128 couple several @xmath83- , @xmath84- or higher excited configurations via the modified coupling amplitudes  ( [ eq : k - co - adc2-c ] ) which are nonzero in this case . we derive in the following an expansion of the dynamic self - energy in terms of generalized wannier orbitals ; such orbitals are localized within unit cells and allow one to apply cutoff criteria inside a cell and between a cell and other cells . this affords one the ability to apply fine - grained configuration selection which is independent of the actual choice of the unit cell . in particular , it enables the treatment of crystals with large unit cells , not amenable to bloch - orbital - based approaches [ sec . [ sec : onion ] ] . @xcite    the generalized wannier transformation  @xcite and its inverse which mediate between bloch and wannier orbitals are    [ eq : wanniertrafo ] w_r ( r ) & = & _ k _ p=1^k u_p(k ) ^-ik r _ k p(r ) , + _ k p(r ) & = & _ r _ = 1^k u^*_p(k ) ^ ik r w_r ( r ) ,    where @xmath130  denotes the number of bands involved in the transformation and @xmath131  is a suitable unitary matrix as generated , e.g. , by the localization procedures in refs .  . both bloch and wannier orbitals are normalized to unity upon integration over the entire crystal which consists of @xmath10  unit cells . generally , we denote the wannier orbital indices that result from occupied bands with  @xmath132 , @xmath133 , @xmath134 , @xmath135 ,  @xmath3 , while the wannier orbital indices of unoccupied bands are denoted by  @xmath136 , @xmath137 , @xmath138 , @xmath139 ,  @xmath3 . indices of wannier orbitals from occupied or unoccupied bands are referred to by  @xmath140 , @xmath141 , @xmath142 , @xmath143 ,  @xmath3 .    for an evaluation of the one - particle green s function in terms of the feynman - dyson perturbation series , it is essential that there is no mixing between occupied and virtual wannier orbitals ; otherwise wick s theorem can no longer be applied . @xcite to fulfill this requirement , the wannier transformation  ( [ eq : wannier ] ) is applied to the occupied and virtual bloch orbitals separately , yielding two independent unitary matrices , one for the occupied bands and one for the virtual bands , respectively . hence , @xmath131 in eq . ( [ eq : wanniertrafo ] ) is block - diagonal . on the one hand , we want to evaluate the co - adc equations entirely in terms of wannier orbitals , in the so - called wannier or local representation ;  @xcite on the other hand , band structures are defined with respect to the crystal momentum quantum number . in what follows three different transformation schemes to switch between the local representation and the crystal momentum representation are applied to the co - adc formalism and their physical and methodological implications are discussed . the dynamic self - energy  @xmath144 depends on two external bloch orbitals @xmath145 and @xmath146 . hence , carrying out the inverse wannier transformation  ( [ eq : invwannier ] ) for the external orbitals in eq . ( [ eq : adc - form - expanded ] ) leads to @xmath147 for the left - hand side of eq . ( [ eq : adc - form - expanded ] ) and a similar expression for the right - hand side where the internal summations  @xmath112 , @xmath113 , @xmath114 , and @xmath115 still run over the brillouin zone .    in order to obtain a representation of the dynamic self - energy  @xmath148 entirely in terms of wannier orbitals , it is important to note that , although matrix elements such as  @xmath149 in eq . ( [ eq : adc - form - expanded ] ) only depend on three independent crystal momenta , they actually describe quantities which depend on four bloch orbitals , namely ,  @xmath23 , @xmath150 , @xmath151 , and @xmath152 and thus a fourfold inverse wannier transformation  ( [ eq : invwannier ] ) has to be applied . this transformation of eq . ( [ eq : adc - form - expanded ] ) to the wannier representation yields @xmath153 exploiting the unitarity of  @xmath154 and the basic relations  @xcite    _ r ^ik r & = & n_0 _ k , 0 , + _ k ^ik r & = & n_0 _ r , 0 .    the check accent on the quantities  @xmath155 , @xmath156 , and @xmath157 in eq . ( [ eq : mwannier ] ) indicates that the sums run over lattice vectors and wannier orbital indices rather than crystal momenta and band indices as in eq . ( [ eq : adc - form - expanded ] ) . the prime on the summation symbol indicates that  @xmath158 and @xmath159 must hold to avoid double counting of contributions . the corresponding equation for  @xmath160 is formally identical to eq . ( [ eq : mwannier ] ) apart from the changed occupation numbers  @xmath161 , @xmath162 , @xmath163 , @xmath164 , @xmath165 , and @xmath166 . one obviously hermitian form for the dynamic self - energy emerging from eqs . ( [ eq : sigmawannier ] ) and ( [ eq : mwannier ] ) is given by    [ eq : mcluster ] m_pq^(k , ) & = & _ , u_p(k ) u^*_q(k ) _ ^(k)^ + & & ( - ^ - ^)^-1 _ ^(k ) , + _ ^(k ) & = & _ r e^i k r _ the vector notation in eq . ( [ eq : mcluster ] ) combines the six internal summations over the intermediate lattice vectors and wannier orbital indices  @xmath167 , @xmath168 , @xmath169 , @xmath170 , @xmath171 , and @xmath172 in eq . ( [ eq : mwannier ] ) . we refer to eq . ( [ eq : mcluster ] ) as the _ supercell form _ of co - adc . the matrix  @xmath173 in eq . ( [ eq : mclusterass ] ) does not explicitly dependend on  @xmath24 which implies that its eigenvalues comprise all pole positions , i.e. , pole positions for each @xmath24-point . consequently , the resulting eigenvalues also do not explicitly depend on  @xmath24 . multiplying with  @xmath174 from the left projects on the desired set of eigenvectors . hence , the translational symmetry among the @xmath10  unit cells of the crystal is not exploited to reduce the size of the configuration space which is unfavorable in conjunction with the cutoff criteria discussed in sec . [ sec : onion ] . in practice , ( [ eq : mcluster ] ) is applied to a molecular cluster formed by @xmath10  unit cells of a crystal for which born von krmn boundary conditions are enforced . as a consequence , only  @xmath10 discrete crystal momenta are in the brillouin zone . the above given transformation from the wannier representation to the crystal momentum representation is simple and robust and can be used in conjunction with almost any electron correlation method , but it suffers from finite - size effects due to the , usually , small number of unit cells which are considered in practice ; see , e.g. , ref .  .    the translational symmetry of the dynamic self - energy  ( [ eq : mwannier ] ) can be exploited by applying a translation by  @xmath175 which removes the explicit dependence of  @xmath176 on two external lattice vectors such that it only depends on the difference of the lattice vectors according to @xmath177 this facilitates to remove the lattice summation over  @xmath178 and the prefactor  @xmath179 in eq . ( [ eq : sigmawannier ] ) which results in ( dropping the prime on  @xmath180 ) @xmath181 to continue , we give up the idea of an adc  form for  @xmath182 , where the modified coupling amplitudes on the left and on the right are related by hermitian conjugation , as is the case in eq . ( [ eq : mcluster ] ) , and arrive at    [ eq : m_k ] m_pq^(k , ) & = & _ , u_p(k ) u^*_q(k ) _ 0 ^ , + & & ( - ^ - ^)^-1 _ ^(k ) + _ ^(k ) & = & _ r ^i k r _ r ^ ,    where a slightly changed definition of the transformed modified coupling amplitude  @xmath183 is employed which is indicated by a tilde accent on the quantities  @xmath184 , @xmath185 , and @xmath186 . because of its asymmetric nature , we denote eq . ( [ eq : m_k ] ) as the _ semitransformed form_.    by exploiting the translational symmetry  ( [ eq : transdyn ] ) , we remove the summation over all translationally equivalent octuples consisting of the two external wannier orbitals  @xmath187 and @xmath188 of the dynamic self - energy  @xmath189 and two @xmath83- or @xmath84-configurations . we arrive at a summation over  @xmath190 , @xmath191 , and @xmath141 in conjunction with the remaining six internal lattice vectors and wannier orbital indices . this formulation thus provides a far better starting point for reducing the number of configurations by means of the cutoff criteria of sec . [ sec : onion ] than the supercell form . using the unsymmetric adc  form  ( [ eq : m_k ] ) implies that the pole positions still do not dependent explicitly on  @xmath24 . moreover , the modified coupling amplitudes  ( [ eq : m_k_dyson_amp ] ) are , in contrast to eq . ( [ eq : mclusteradd ] ) , no longer related by hermitian conjugation ; yet , for sure , the dynamic self - energy itself stays hermitian .    this way of exploiting translational symmetry in conjunction with the wannier orbitals has already been applied successfully before to devise _ ab initio _ electron correlation methods for crystals like the local hamiltonian approach  @xcite or the green s function based method of albrecht  _ et al._.  @xcite    the internal structure of the dynamic self - energy is not regarded in the two formulas  ( [ eq : mcluster ] ) and ( [ eq : m_k ] ) . expression  ( [ eq : m_k ] ) results from eqs . ( [ eq : sigmawannier ] ) and ( [ eq : mwannier ] ) by exploiting the translational symmetry of  @xmath116 in eq . ( [ eq : transdyn ] ) and utilizing its exclusive dependence on two external wannier orbitals . however , the translational symmetry of the matrix  @xmath192 in eq . ( [ eq : m_k ] ) can also be exploited additionally which already was harnessed to derive eq . ( [ eq : adc - form - expanded ] ) from its original molecular orbital formulation . we obtain    [ eq : m_k_hermitian ] m_pq^+(k , ) & = & _ , u_p(k ) u^*_q(k ) _ g_1 , g_2 , g_1^ , g_2^ _ , , , ^ , ^ , ^ |u _ ; g_1 g_2 ^+ * ( k ) + & & ( - ^+(k ) - ^+(k))^-1_g_1 g_2 ; g_1^ g_2^ ^ ^ ^ + & & |u _ ; g_1^ g_2^ ^ ^ ^^+(k ) , +    ( ^+(k ) + ^+(k))_g_1 g_2 ; g_1^ g_2^ ^ ^ ^ & = & _ r ^ik r ( ^+ + ^+)_0 g_1 g_2 ; r ^ g_1^ + r ^ g_2^ + r ^ ,    where the identities  @xmath193 , @xmath194 , and @xmath195 hold . the corresponding relations for  @xmath196 are formally identical to eq . ( [ eq : m_k_hermitian ] ) apart from the changed occupation numbers  @xmath197 , @xmath162 , @xmath163 , @xmath198 , @xmath165 , and @xmath199 . alternatively , we can derive eq . ( [ eq : m_k_hermitian ] ) by inserting the inverse wannier transformation  ( [ eq : invwannier ] ) into the co - adc equations for bloch orbitals  ( [ eq : adc - form - expanded ] ) which already are fully adapted to translational symmetry . we refer to eq . ( [ eq : m_k_hermitian ] ) as the _ fully translational symmetry adapted form_.    the modified coupling amplitudes  ( [ eq : m_k_h_u ] ) are constructed by considering only wannier orbitals relative to an origin cell ( or reference cell ) where the hole with index  @xmath132 is assumed to reside . the external orbital index  @xmath141 of a modified coupling amplitude can be viewed to represent a bloch orbital  @xmath200 . this orbital interacts with the @xmath83-configurations that are pinned with one lattice vector to the origin cell and extend with the two remaining lattice vectors  @xmath201 and @xmath202 over up to two different unit cells . alternatively , one can use the translational symmetry of the  @xmath203 to arrive at @xmath204 now the external orbital is  @xmath205 which is independent of  @xmath178 and consequently eq .  ( [ eq : m_k_h_u_trans ] ) is interpreted to describe the interaction of a wannier orbital in the origin cell with @xmath83-configurations centered around all lattice vectors  @xmath178 which are combined to give an intermediate @xmath83-configuration with a total crystal momentum  @xmath206 . the matrix  @xmath207 in eq . ( [ eq : m_k_h_kc ] ) describes the coupling among the @xmath83-configurations . one of the lattice vectors of the first triple of wannier orbitals is pinned to the origin cell , and the two other lattice vectors  @xmath201 and @xmath202 are offsets to it . the remaining three lattice vectors belong to a wannier orbital in an arbitrary cell  @xmath178 and two further wannier orbitals in the cells  @xmath208 and @xmath209 relative to the former . obviously the entire eq . ( [ eq : m_k_h_kc ] ) can be interpreted to describe the interaction of a @xmath83-configuration of wannier orbitals centered around the origin cell with another @xmath83-configuration with crystal momentum  @xmath24 . due to the full exploitation of translational symmetry , also of the intermediate @xmath83-configurations , the fully translational symmetry adapted form yields the smallest configuration space in conjunction with the cutoff criteria of sec .  [ sec : onion ] .    in order to give explicit co - adc expressions for the matrices showing up in eq . ( [ eq : m_k_hermitian ] ) , we transform the hamiltonian  ( [ eq : bloch_hamilton ] ) with the help of formula  ( [ eq : invwannier ] ) to the wannier representation . a similar decomposition as in eq . ( [ eq : bloch_hamilton ] ) is chosen for a subsequent application of perturbation theory : @xmath210 with  @xmath211 being the diagonal elements of the fock matrix  @xmath212 . here @xmath213 } \ , n_{\vec r \ , \tau } + ( 1 - \delta_{\vec r_1 \ , \varrho \ ; \vec r_2 \ , \sigma } ) \ , f_{\vec r_1 \ , \varrho \ ; \vec r_2 \ , \sigma}$ ] comprises the negative of the hartree - fock potential and also the off - diagonal elements of the fock matrix . note that  @xmath212 is block - diagonal and thus the wannier orbitals are still separated into occupied and virtual wannier orbitals . using the inverse generalized wannier transformation  ( [ eq : invwannier])as is done throughout this article  corresponds to the use of the full fock matrix as zeroth - order hamiltonian . in addition , a perturbative expansion of the off - diagonal elements in the fock matrix is applied in succession of the inverse wannier transformation . this conduct is equivalent to employing the partitioning of the fock matrix introduced in eq . ( [ eq : wannier_hamilton ] ) . the one - electron interaction is mediated by the off - diagonal fock matrix elements , and the two - electron interaction is represented by the two - electron integrals ; let the one - electron interaction be treated using feynman diagrams up to @xmath19-th  order and account for diagrams which exclusively involve two - electron interactions of highest order  @xmath16 . this separation facilitates to take higher - order terms into account for the slowlier convergent series ; the resulting scheme is denoted by co - adc(@xmath19,@xmath16 ) . the co - adc(2,2 ) approximation for wannier orbitals is obtained similarly to the case of the co - adc(2 ) approximation for bloch orbitals  ( [ eq : k - co - adc2 ] ) and reads    [ eq : r - co - adc2 ] |u^+_r ; g_1 g_2 g_3 & = & v^*_r g_1 [ g_2 g_3 ] n_g_1     g_3^ ^ ( _ g_2 + _ g_3 p_3 - _ g_1 ) n_g_1 |n_g_2 |n_g_3 , +     the equations for  @xmath214 , @xmath215 , and @xmath216 are formally identical to eq . ( [ eq : r - co - adc2 ] ) upon modifying the labels according to the adaptation of eq . ( [ eq : m_k_hermitian ] ) to  @xmath89 and replacing the occupation numbers in eqs . ( [ eq : m_r_h_u ] ) and ( [ eq : m_r_h_k ] ) by  @xmath217 . up to second order , the perturbative treatment of the off - diagonal elements of the fock matrix does not introduce any new diagrams for the dynamic self - energy and thus @xmath218 remains diagonal . @xcite however , in a formulation which is invariant under unitary transformations in the occupied and virtual space , this is no longer the case because the full fock matrix needs to be chosen as the zeroth - order hamiltonian in place of  @xmath219 then . the static self - energy has not been determined so far . it is represented by all those diagrams in the diagrammatic expansion of the self - energy where the external points of the diagrams correspond to equal times . @xcite in a strict second - order treatment of electron correlations , using bloch orbitals  ( [ eq : bloch_hamilton ] ) , the static self - energy is zero because the first static self - energy diagrams arise the earliest in third order . moreover , it turns out that the diagrammatic series for the static self - energy does not converge reasonably in many cases . @xcite however , a self - consistent solution is possible  @xcite utilizing  @xcite @xmath220 } \\      & & { } \times \biggl [ \frac{1}{2\pi{{\rm i } } } \oint g_{sr}(\vec k^{\prime } ,        \omega ) { \>\textrm{d}}\omega \biggr ] \ ; ,    \end{array}\ ] ] here given in terms of bloch orbitals . the contour integration in eq . ( [ eq : statics ] ) runs along the real axis and closes in the upper complex @xmath77-plane , hence enclosing only the poles of the advanced green s function  @xmath221 in eq . ( [ eq : p - decomp ] ) .    the self - consistent solution of eq . ( [ eq : statics ] ) is computationally expensive . yet a stable and efficient dyson expansion method  ( dem ) for determining the static self - energy is devised in ref .  . there , the first two terms of the dyson expansion  ( [ eq : dysonexp ] ) are inserted into eq . ( [ eq : statics ] ) . the term  @xmath222 in eq . ( [ eq : statics ] ) cancels the result of the contour integration over the free green s function  @xmath223 . carrying out the contour integration over the product of the two free green s functions and the static self - energy yields    [ eq : statick ] _ pq^(k ) & = & _ k^ _ r , s a_k p k q , k^ s k^ r _ sr^(k^ ) + b_pq(k ) , + a_k p k q , k^ s k^ r & = & v_k p k^ r [ k q k^ s ] ( _ k^ r - _ k^ s)^-1 , + b_pq(k ) & = & _ k^ _ r , s v_k p k^ r [ k q k^ s ] q_rs(k^ ) , + q_rs(k^ ) & = & g^0_ss(k^ , ) m_sr(k^ , ) g^0_rr(k^ , ) .    the  @xmath224 can now be determined by inserting the spectral representation of the dynamic self - energy  ( [ eq:2h1p - spectral ] ) into eq . ( [ eq : kqqp ] ) . carrying out the contour integration yields    [ eq : qkl_mlk - k ] q_rs(k^ ) & = & q^+_rs(k^ ) + q^-_rs(k^ ) , + q^+_rs(k^ ) & = & _ n \{n+1 } m_s^+,(n)(k^ ) m_r^+,(n)*(k^ ) , + q^-_rs(k^ ) & = & _ n \{n-1 } m_s^-,(n)(k^ ) m_r^-,(n)*(k^ ) .    in order to evaluate eq .  ( [ eq : qkl_mlk - k ] ) , from now on , we resort to the adc  form of the dynamic self - energy in terms of bloch orbitals  ( [ eq : adc - form - expanded ] ) . by affixing a bar accent to the quantities  @xmath117 , @xmath93 , and @xmath94 in all following equations of this subsection , they become formally identical to the ones obtained with the fully translational symmetry - adapted form  ( [ eq : m_k_hermitian ] ) .    the direct usage of eq . ( [ eq : qkl_mlk - k ] ) is in practice very time consuming because the  @xmath225 matrix has to be diagonalized fully  ( [ eq : adc - inversion ] ) to obtain the pole positions  @xmath226 of the spectral representation  ( [ eq:2h1p - spectral ] ) . for the co - adc(2 ) approximation of the dynamic self - energy in crystal momentum representation  ( [ eq : k - co - adc2 ] ) , however , @xmath227  remains diagonal . this leads to  @xmath228 and @xmath229 where the @xmath16-th  eigenvector of  @xmath230 is unity for a single @xmath83- or @xmath84-configuration , respectively , and zero otherwise . in general , the  @xmath224 are determined efficiently from eq . ( [ eq : qkl_mlk - k ] ) by the inversion method or by a single - vector lanczos diagonalization to circumvent a full diagonalization of  @xmath225 . @xcite    the inversion method has been found to be more efficient  @xcite than the single - vector lanczos diagonalization , so we concentrate on the former . we define auxiliary vectors @xmath231 as the solution of the inhomogeneous systems of linear equations    @xmath232    inserting the eigenvector matrix  @xmath233 from eq . ( [ eq : adc - inversion ] ) into eq . ( [ eq : vr ] ) , one obtains , by using formula  ( [ eq : moddysontraf ] ) , @xmath234 which reveals the usefulness of the new vectors and allows us to rewrite eq . ( [ eq : qkl_mlk - k ] ) as follows :    [ eq : qrs_vr ] q_rs^+(k^ ) & = & -v_s^+(k^ ) v_r^+(k^ ) n_k^ r n _ k^ s + ( _ k^ r - _ k^ s)^-1 + & & , + q_rs^-(k^ ) & = & v_s^-(k^ ) v_r^-(k^ ) |n_k^ r | n_k^ s - ( _ k^ r - _ k^ s)^-1 + & & . we have reduced the problem of determining the  @xmath235 to the problem of determining the  @xmath231 . inspecting eqs . ( [ eq : vr ] ) and ( [ eq : qrs_vr ] ) , we note that we have to solve for the  @xmath236 in the large @xmath83-configuration space but only for the usually small set of all occupied orbitals . conversely , the  @xmath237 have to be calculated in the small @xmath84-configuration space but for the usually large set of all virtual orbitals . the diagonal parts  @xmath238 in eq . ( [ eq : vr ] ) are , at least , of a magnitude of twice the hartree - fock band gap , implying usually a diagonal dominance of the full matrices in eq . ( [ eq : vr ] ) . therefore , a solution by jacobi iterations  @xcite is suggested : @xmath239 which turns out to converge rapidly .    the inhomogeneous linear system of equations  ( [ eq : kstaticimplicit ] ) can now be solved for  @xmath240 by a matrix inversion  @xcite : @xmath241 where  @xmath242 , @xmath243 , and @xmath244 are composed by numerating the compound indices  @xmath245 and @xmath246 in  @xmath240 , @xmath247 , and @xmath248 by integer numbers  @xmath35 and @xmath36 , respectively . inspecting eq . ( [ eq : ka ] ) , we note that @xmath249 couples only to the components of @xmath250 where the @xmath251  indices denote a combination of particle - hole- or hole - particle - orbitals . as all quantities in eq . ( [ eq : getsigmak ] ) are hermitian , we can further restrict the band indices to  @xmath252 . as soon as the number of @xmath24-points becomes large , the solution of eq . ( [ eq : getsigmak ] ) gets cumbersome as the dimension of the system of linear equations is given by the number of entries of the upper triangle ( including the diagonal entries ) of the static self - energy matrix times the number of @xmath24-points . thus an iterative linear equations solver has to be employed . @xcite      the static self - energy can also be determined directly in wannier representation . having expressed the static self - energy in terms of wannier orbitals , we obtain the static self - energy in crystal momentum representation by @xmath253 to find an approximation for  @xmath254 , we insert the inverse wannier transformation  ( [ eq : invwannier ] ) and its hermitian conjugate into eq . ( [ eq : statick ] ) and arrive at    [ eq : getsigma ] _ ^(r ) & = & _ r^ _ , _ ^(r^ ) + & & + b_(r ) , + b_(r ) & = & _ r^ _ , q _ ( r^ ) . note that the translational symmetry of  @xmath255 could only be utilized once to remove one lattice summation such that two lattice summations ( one over @xmath180 and one over @xmath256 ) show up in eq . ( [ eq : getsigma ] ) while only one summation over  @xmath257 is necessary in eq . ( [ eq : statick ] ) . there is no first - order contribution to the static self - energy in eq . ( [ eq : getsigma ] ) because the static self - energy used here evolves from an inverse wannier transformation of the self - energy obtained in mller - plesset partition  @xcite and a subsequent perturbative expansion with respect to the off - diagonal elements of the fock matrix .    for the remaining equations  ( [ eq : ka ] ) and ( [ eq : kqqp ] ) , we assume the partition  ( [ eq : wannier_hamilton ] ) of the hamiltonian in wannier representation , to obtain    [ eq : staticcoeff ] a_0 r , g g + r^ & = & v_0 g + r^ [ r g ] + & & ( _ g + r^ - _ g ) ^-1 + & & ( n_g + r^ |n_g - |n_g + r^ n_g ) , + q_(r^ ) & & q_r^ 0 = g_^0(0 , ) + & & m_(r^ , ) g_^0(0 , ) ,    with the free green s function in wannier representation given by @xmath258 \ ; . \end{array}\ ] ]    to determine the  @xmath259 , we introduce the spectral representation of the dynamic self - energy in wannier representation [ the analog of eq . ( [ eq:2h1p - spectral ] ) ] @xmath260 which is obtained from the adc  form of the dynamic self - energy in wannier representation  ( [ eq : mwannier ] ) by letting  @xmath261 and diagonalizing  @xmath173 analogously to eq . ( [ eq : adc - inversion ] ) . since translational symmetry can not be exploited in eq . ( [ eq : mwannier ] ) as  @xmath180in contrast to  @xmath24is not a good quantum number , this leads to the redundancies that have already been mentioned in conjunction with the supercell form of the dynamic self - energy in eq . ( [ eq : mcluster ] ) ; namely , the full configuration space of @xmath83- and @xmath84-configurations from the @xmath10  unit cells of the crystal is necessary here . however , the cutoff criterion of sec . [ sec : onion ] , nevertheless , can be applied . the  @xmath259 can now be determined by inserting eq . ( [ eq:2h1p - spectral - r ] ) into eq . ( [ eq : qtauupsilon ] ) . carrying out the contour integration yields    [ eq : qkl_mlk ] q_(r^ ) & = & q_^+(r^ ) + q_^-(r^ ) , + q^+_(r^ ) & = & _ n \{n+1 } m_^+,(n)(0 ) m_^+,(n)*(r^ ) , + q^-_(r^ ) & = & _ n \{n-1 } m_^-,(n)(0 ) m_^-,(n)*(r^ ) . we proceed as before by defining the auxiliary vectors @xmath262 setting  @xmath263 , and rewrite eq . ( [ eq : qkl_mlk ] ) :    [ eq : qkl - r ] q_^+(r^ ) & = & - v_^+(0 ) v_^+(r^ ) n_r^ n_0 + ( _ r^ - _ 0 ) ^-1 + & & , + q_^-(r^ ) & = & v_^-(0 ) v_^-(r^ ) |n_r^ |n_0 - ( _ r^ - _ 0 ) ^-1 + & & . the algorithm for jacobi iterations now reads @xmath264 and the inhomogeneous linear system of equations  ( [ eq : staticinhomo ] ) can be solved for  @xmath254 by a matrix inversion  @xcite : @xmath265 where a check accent is affixed to indicate that wannier orbitals rather than bloch orbitals are used . we employ a mapping of the compound indices  @xmath266 and @xmath267 to the integer numbers  @xmath35 and @xmath36 , respectively , for the quantities  @xmath268 , @xmath269 , and @xmath270 . note that @xmath271 is not hermitian . hence , the full static self - energy matrix has to be included in  @xmath272 explicitly for several lattice vectors  @xmath178 . however , the property @xmath273 holds which can be utilized to reduce the number of equations in formula  ( [ eq : getsigmar ] ) . nevertheless , in many cases , the linear system of equations needs to be solved iteratively . having determined approximations for the static and the dynamic self - energy in terms of wannier orbitals in eqs . ( [ eq : sigmakr ] ) and ( [ eq : m_k_hermitian ] ) , respectively , we can finally determine the positions and strengths of the poles of the one - particle green s function from eq . ( [ eq : solvedysonformal ] ) . [ the pole positions and strengths of the green s function , using bloch orbitals to represent the self - energy  ( [ eq : statick ] ) and ( [ eq : adc - form - expanded ] ) , can be obtained following a nearly identical line of argument . ] to this end , the small hermitian eigenvalue problem @xmath274 with @xmath275  being the diagonal matrix of bloch orbital energies and @xmath276  denoting eigenvectors , has to be solved self - consistently , i.e. , such that the energy  @xmath77 entering the dynamic self - energy  @xmath277 is identical to the resulting eigenvalue  @xmath77 .    to avoid the unitary matrices  @xmath131 and @xmath278 in eqs . ( [ eq : m_k_h_m ] ) and ( [ eq : sigmakr ] ) arising in the following expressions , we multiply eq . ( [ eq : selfconsistent ] ) from the left with  @xmath278 and insert  @xmath279 between the right parenthesis and the eigenvector  @xmath280 on the left - hand side of the formula . defining the translational symmetry - adapted fock matrix  @xmath281 , one arrives at @xmath282 where the new quantities  @xmath283 , @xmath284 , @xmath285 , and @xmath286 are related to the old ones  @xmath287 , @xmath240 , @xmath288 , and @xmath289 by  @xmath290 , @xmath291 , @xmath292 , and @xmath293 .    inserting the adc  form of the dynamic self - energy  ( [ eq : m_k_h_m ] ) into eq . ( [ eq : fixedomega ] ) and defining @xmath294 yields the following form of the eigenvalue problem  ( [ eq : fixedomega ] ) : @xmath295 rewriting eq . ( [ eq : fixeddynamic ] ) as @xmath296 it becomes evident that eqs . ( [ eq : unfixomega ] ) and ( [ eq : evsub ] ) can be recast as a joint hermitian eigenvalue problem    @xmath297    & { \hbox{\boldmath{$b$}\unboldmath}}(\vec k ) = \left (    \begin{array}{ccc }      { \hbox{\boldmath{$\bar f$}\unboldmath}}(\vec k ) + { \hbox{\boldmath{$\bar \sigma$}\unboldmath}}^{\infty}(\vec k )      & { \hbox{\boldmath{$\bar u$}\unboldmath}}^{+\dagger}(\vec k ) & { \hbox{\boldmath{$\bar u$}\unboldmath}}^{-\dagger}(\vec k ) \\      { \hbox{\boldmath{$\bar u$}\unboldmath}}^+(\vec k ) & { \hbox{\boldmath{$\bar k$}\unboldmath}}^+(\vec k ) + { \hbox{\boldmath{$\bar c$}\unboldmath}}^+(\vec k ) & { \hbox{\boldmath{$0$}\unboldmath } } \\ { \hbox{\boldmath{$\bar u$}\unboldmath}}^-(\vec k ) & { \hbox{\boldmath{$0$}\unboldmath } } & { \hbox{\boldmath{$\bar k$}\unboldmath}}^-(\vec k ) + { \hbox{\boldmath{$\bar c$}\unboldmath}}^-(\vec k )    \end{array }    \right ) \ ; . & \end{aligned}\ ] ]    the new matrix  @xmath298 is called the _ band structure matrix_. it has to be diagonalized for several @xmath24-points yielding eigenvalues  @xmath299 and eigenvectors  @xmath300 . formula  ( [ eq : solvepadc ] ) is similar to the result for molecules . @xcite one can envisage that the diagonalization of the subblocks  @xmath301 in  @xmath298 is carried out first by solving eq . ( [ eq : adc - inversion ] ) ; afterwards , the full eigenvalue problem  ( [ eq : solvepadc ] ) is treated with all subblocks modified accordingly . @xcite therefore , via eq . ( [ eq : adcsuminfinity ] ) , infinitely many proper self - energy diagrams are summed before the result is put in the dyson equation  ( [ eq : dyson ] ) to sum all improper self - energy diagrams that derive from those contained in the adc  form . the spectral representation of the one - particle green s function  ( [ eq : p - spectrum ] ) reads , in terms of the eigenpairs of the band structure matrix  ( [ eq : solvepadc ] ) , @xmath302 where @xmath303  denotes the transition amplitude of the @xmath16-th state . it is given by either the  @xmath304 or the  @xmath305 in eq . ( [ eq : transamp ] ) whether it belongs to an @xmath306 or an @xmath307  particle state .    in some cases , the external wannier orbitals of the dynamic self - energy are restricted to a single cell of the crystal , i.e. , @xmath308 , the so - called one - lattice - site approximation for the dynamic self - energy . note , however , that the @xmath83- and @xmath84-configurations in neighboring unit cells are not entirely excluded . the @xmath24-dependence , in this approximation , is solely due to the one - particle matrix  @xmath309 . the static self - energy can also be evaluated in one - lattice - site approximation . then the @xmath24-dependence is exclusively due to  @xmath310 and this approximation to the co - adc equations becomes similar to the dynamical mean - field theory . all equations derived so far contain certain infinite lattice sums that run over the entire crystal forming the configuration space . to obtain a meaningful approximation of the one - particle green s function , the convergence of those lattice sums must be granted . sun and bartlett  @xcite study the convergence properties of the lattice sums in several quantum chemical correlation methods for ground states and exited states of crystals , based on the mller - plesset partition of the hamiltonian  ( [ eq : bloch_hamilton ] ) . following the reasoning of ref .  , the lattice sums in co - adc in crystal momentum representation converge . moreover , the equations in local representation also converge as they evolve from the former by means of an inverse wannier transformation  ( [ eq : invwannier ] ) . although the convergence of lattice sums is granted principally , we have to devise an algorithm for their proper truncation ;  @xcite i.e. , a dynamical building of the configuration space is required which means to meet a chosen accuracy in the lattice sums for a given crystal . the appropriate truncation of lattice sums is an essential ingredient of all _ ab initio _ methods for crystals and corresponds to a configuration selection procedure which ensures that the configuration space is sufficiently large for the desired accuracy of the band structure but still sufficiently small to be tractable on present - day computers . in molecular physics , configuration selection has been introduced in the context of the configuration interaction method  @xcite which suffers from an exponentially growing configuration space with respect to the number of atoms in molecules . the techniques discussed in ref . are not applicable here , as they are designed with the configuration interaction method in mind , but there are analogs for crystals which are similar in spirit .    the incremental scheme for ground states of crystals  @xcite is based on wannier orbitals . the total configuration space of a crystal is partitioned into certain subsets which are used in a subsequent calculation of the correlation energy to yield the so - called energy increments . frequently , explicit manual configuration selection is applied . this allows the partition of the configuration space to be chosen following chemical intuition . yet the situation is much more cumbersome for band structures because the number of distinct excited states grows with the number of correlated electrons , i.e. , the number of bloch orbitals and wannier orbitals which are considered in the sums in eqs . ( [ eq : adc - form - expanded ] ) and ( [ eq : m_k_hermitian ] ) , respectively . in the local hamiltonian approach , the number of states to be calculated is restricted to a fixed number which is treated in an incremental way . @xcite albrecht  _ et al . _ @xcite ensure the convergence of an incremental series for the self - energy for  @xmath311 , a criterion which is independent of the number of states actually described . we devise here a configuration selection procedure that is perfectly adapted to the structure of the co - adc scheme . for the equations in terms of bloch orbitals  ( [ eq : statick ] ) and ( [ eq : adc - form - expanded ] ) , @xmath24-points are chosen harnessing born von krmn boundary conditions which lead to a net of equidistant @xmath24-points . @xcite in this case , configuration selection means choosing a sufficiently dense net in conjunction with a cutoff criterion which is similar to the one for co - adc in terms of wannier orbitals  ( [ eq : sigmakr ] ) and ( [ eq : m_k_hermitian ] ) discussed in the next paragraph . yet beforehand , one should interpret the impact of the number of @xmath24-points used for brillouin zone integration . using a given net of @xmath10 @xmath24-points to carry out a wannier transformation  ( [ eq : wannier ] ) , co - adc in terms of bloch orbitals can be immediately represented and analyzed in terms of the supercell form  ( [ eq : mcluster ] ) . however , one should emphasize that the redundancies of this representation  which is utilized here for interpretation only  are not present in the bloch orbital formulation of co - adc . the supercell form affords a tight - binding interpretation : within the supercell , the number of distinct neighbors to the origin cell is restricted by the volume of the supercell . distinct interaction terms among the cells , i.e. , the fock matrix and two - electron integrals , are thus also resticted and the number of @xmath24-points employed can thus be understood to imply the number of nearest - neighbor cells treated in the interaction terms . these arguments generalize the tight - binding arguments given , e.g. , in ref .  , in conjunction with a single brillouin zone integration of a crystal - momentum - dependent function to the multi - dimensional case . configuration selection in crystal momentum representation is not intuitive and does not allow a fine - grained selection of configurations within unit cells and between an origin cell and its neighbor cells . for crystals with large unit cells this is a significant restriction . to obtain a suitable cutoff criterion for co - adc in terms of wannier orbitals  ( [ eq : sigmakr ] ) and ( [ eq : m_k_hermitian ] ) , we evaluate the second - order diagram of the retarded dynamic self - energy  @xcite and examine the summand therein : @xmath312 } \ , v_{\vec r \ , \sigma \ ; \vec g_1 \ , \kappa \ ; [ \vec g_2 \ , \alpha    \ ; \vec g_3 \ , \beta]}^ * \over    \omega - \varepsilon_{\vec g_2 \ , \alpha } - \varepsilon_{\vec g_3 \ , \beta }    + \varepsilon_{\vec g_1 \ , \kappa } } \ ; n_{\vec g_1 \ , \kappa } \ ;    \bar n_{\vec g_2 \ , \alpha } \ , \bar n_{\vec g_3 \ , \beta } \ ; .\ ] ] the most delocalized occupied orbitals are the valence orbitals . therefore , the occupied valence orbitals couple to the most important @xmath83-configurations hence , @xmath77 can be assumed to be of the order of the band gap . the denominator in eq . ( [ eq : selfsecond ] ) is regarded to be constant to a reasonable degree . it is thus considered as being part of the cutoff threshold . therefore , if  @xmath313 } \ , v_{\vec r \ , \sigma \ ; \vec g_1 \ , \kappa \ ; [ \vec g_2 \ , \alpha \ ; \vec g_3 \ , \beta]}^ * \ ; n_{\vec g_1 \ , \kappa } \ ; \bar n_{\vec g_2 \ , \alpha } \ , \bar n_{\vec g_3 \ , \beta}$ ] is above a certain cutoff threshold for some combination of  @xmath314 and @xmath315 , we include the corresponding @xmath83-test - configuration  @xmath316 in the configuration space . if the product is small for all external wannier orbitals of the self - energy , @xmath317 is neglected completely by setting all two - electron integrals containing this configuration exactly to zero . by this , only a finite range of the residual coulomb interaction , which is not treated in hartree - fock approximation , is explored . the same arguments hold for the selection of @xmath84-configurations for the evaluation of the advanced dynamic self - energy , if eq . ( [ eq : selfsecond ] ) is adapted in the same way as eq . ( [ eq : r - co - adc2 ] ) is adapted to  @xmath214 , @xmath215 , and @xmath216 . convergence of the configuration space , i.e. , the number of unit cells contributing configurations to be considered in the cutoff criterion  ( [ eq : selfsecond ] ) , is expected to be sufficiently quick . van der waals dispersion interactions which are present both in the ground state of the @xmath9  particle system and the states of the @xmath59  particle system are nearly equal and thus almost cancel in band structure calculations . however , the effect of the extra coulomb charge that arises in @xmath59  particle states in contrast to the @xmath9  particle ground state is long range . the potential of an extra point charge falls off like  @xmath318 , where @xmath319  is the dielectric constant of the solid . this causes an overall polarization of the crystal which can be accounted for by adjusting the chemical potential  @xmath65 of the added or removed electron accordingly by using a continuum approximation . @xcite this procedure is not specific to the co - adc theory advocated here to treat the short - range interactions , but is required generally . having introduced a convenient method to truncate the configuration space for co - adc , we have to elucidate the physical impact of configuration selection on band structures . to this end , we consider the fully translational symmetry adapted form  ( [ eq : m_k_hermitian ] ) . the modified coupling amplitudes  ( [ eq : m_k_h_u ] ) in the band structure matrix  ( [ eq : solvepadc ] ) of a crystal carry two independent crystal lattice vectors  @xmath320 and @xmath321 . if @xmath320 and @xmath321 are sufficiently far away from the origin cell , all two - electron matrix elements  ( [ eq : m_r_h_u ] ) which contribute to the modified coupling amplitude  @xmath322 fall for all  @xmath141 below the configuration selection threshold and are therefore zero . because the modified interaction matrices  ( [ eq : m_r_h_c ] ) vanish in the co - adc(2,2 ) approximation and eq . ( [ eq : m_k_h_kc ] ) thus denotes a diagonal matrix , the eigenvalue associated with the configuration  @xmath323 is simply given by the diagonal element  @xmath324 , with an eigenvector being unity on position  @xmath325 and zero elsewhere . consequently , such an eigenstate of  @xmath298 is removed by deleting the column and the row  @xmath326 in  @xmath298 because the transition amplitude  ( [ eq : transamp ] ) vanishes and thus it does not contribute to the spectral representration of the one - particle green s function  ( [ eq : spectral_co - adc ] ) . a further consequence of configuration selection involves the lattice summations over  @xmath178 in eqs . ( [ eq : m_k_h_u ] ) and ( [ eq : m_k_h_kc ] ) . in the first place , we consider a crystal with a macroscopic lattice constant which consists of  @xmath10 unit cells . as a result , only fock and two - electron matrix elements that involve wannier orbitals from the origin cell are selected . therefore , in eqs . ( [ eq : m_k_h_u ] ) and ( [ eq : m_k_h_kc ] ) , the lattice sums run over only the single term for  @xmath327 with  @xmath328 and the resulting band structure matrix  ( [ eq : solvepadc ] ) is consequently independent of  @xmath24 . as the brillouin zone contains @xmath10  crystal momenta , we obtain a @xmath10-fold degenerate spectrum similarly to simple tight - binding models . @xcite in typical crystals , interactions with neighboring unit cells are important . yet a sufficiently large supercell which consist of @xmath329  unit cells can be chosen such that interactions between supercells are negligible leading to a @xmath330  fold degenerate spectrum . in other words , a cutoff threshold , to be used for the selection of fock matrix elements and @xmath83- and @xmath84-configurations , implies a definition of degeneracy of the physical states in our model of the crystal .    in order to selectively diagonalize the band structure matrix , one uses an iterative eigenvalue solver , e.g. , a block - lanczos algorithm ,  @xcite that is capable of exploiting the sparsity of  @xmath298 . the most expensive step of the block - lanczos algorithm is a matrix times vector product  denote the number of rows of the square matrix  @xmath249 . then , specifically for the co - adc(2 ) approximation in terms of bloch orbitals  ( [ eq : adc - form - expanded ] ) , ( [ eq : k - co - adc2 ] ) , and ( [ eq : solvepadc ] ) , a single matrix times vector operation requires  @xmath331 = o(n_{\mathbf f(\vec k ) } \ ,    n_{\mathbf k^+(\vec k)})$ ] flops . by affixing bar accents to the matrices showing up in the expression for  @xmath332 , we immediately obtain the corresponding formula for the number of flops required in the co - adc(2,2 ) approximation for wannier orbitals  ( [ eq : m_k_hermitian ] ) , ( [ eq : r - co - adc2 ] ) , and ( [ eq : solvepadc ] ) . ] between  @xmath298 and lanczos vectors which determines the overall performance of the eigenvalue solver . @xcite to investigate the scaling of a selective computation of eigenpairs , we consider a supercell consisting of two unit cells of an original lattice . let us assume a crystal with a macroscopic lattice constant to investigate the asymptotic scaling behavior of the problem . then , the configuration selection method of the previous paragraphs selects only configurations local to the individual unit cells . hence , the resulting configuration space scales linearly with the system size , i.e. , the number of atoms per unit cell . this is a necessary condition for the selective diagonalization of the band structure matrix  ( [ eq : solvepadc ] ) to scale linearly as well . the total fock matrix of the supercell also is block - diagonal . similarly , the modified coupling amplitudes  ( [ eq : m_k_h_u ] ) only contain two nonzero blocks which describe the coupling of @xmath83- and @xmath84-configurations to the fock matrix of each of the two constituting unit cells of the supercell such that the total band structure matrix decomposes into two subproblems that can be solved independently . a single matrix times vector product for the supercell requires twice as many floating point operations than are needed for a matrix times vector product for one of the two unit cells in the supercell . hence , the computation of matrix times vector products scales linearly with the system size . ) with a given accuracy , is assumed to be the same for a single unit cell and a supercell consisting of two unit cells . ] yet doubling the system size also usually means that we are interested in twice as many excited states . hence , the overall effort to determine all excited states of a crystal in a given energy range scales quadratically . it is a _ quadratically _ scaling problem , where linear scaling can only be achieved by an _ a priori _ restriction to a few excited states of the system .    for typical crystals , the fock matrix of a supercell is not block - diagonal which implies another doubling of the number of floating point operations upon doubling the system size . yet this factor cancels as the number of @xmath24-points needed for a given accuracy of the integration over the brillouin zone of the supercell is halved because the volume of its brillouin zone is half the volume of the brillouin zone which corresponds to the original crystal lattice . lithium fluoride crystallizes in a face - centered - cubic  ( fcc ) rocksalt structure described by the space group  @xmath333 as shown in fig .  [ fig : lifgeom ] . the crystal lattice has a two - atomic basis li(0,0,0 ) and f(@xmath334 ) , given in units of the lattice constant is reported therein which is consistent with early and recent experimental values . in order to allow comparison , we adopt the previously used value of  @xmath335 . although the deviation is small , it has a quite noticeable influence on the fundamental band gap which decreased by  @xmath336 when changing from  @xmath337 to @xmath338 in a recent density functional theory study at the lda  level . @xcite ] @xmath339 . to investigate a lif  crystal , we have chosen the simplest possible _ ab initio _ description . the main purpose here is to demonstrate the feasibility of the co - adc formulas for a three - dimensional crystalline solid . thus , only a minimal gaussian sto-6 g basis set  @xcite is used which describes the lithium atom by a @xmath340 $ ]  contraction and the fluorine atom by  @xmath341 $ ] . each shell of the sto-6 g basis set is constructed by fitting six gaussian functions to each occupied orbital of the isolated atoms . @xcite the @xmath342  shell of a lithium atom is rather diffuse due to the single outer valence electron . as the lithium atoms are ionized in lif crystals , we remove the two most diffuse gaussian functions in the li@xmath7  contraction because they do not reflect the physical situation of a compact li@xmath343  ion , and arrive at a @xmath344 $ ]  contraction scheme . the hartree - fock equations are solved self - consistently by means of the wannier program  @xcite which directly yields wannier orbitals . a finite cluster of unit cells is utilized as support for the wannier orbitals of the origin cell which , in our case , consist of up to third - nearest - neighbor cells ( 43  unit cells alltogether ) . as the program determines only occupied wannier orbitals , crystal projected atomic orbitals  ( crystal paos )  @xcite@xmath345 have been devised on the basis of the projected atomic orbitals introduced by saeb and pulay  @xcite and also used by hampel and werner . @xcite the resulting virtual functions are subsequently lwdin orthonormalized . @xcite for the co - adc(2,2 ) method  ( [ eq : m_k_hermitian ] ) , ( [ eq : r - co - adc2 ] ) , and ( [ eq : solvepadc ] ) , the fock matrix and the two - electron integrals in wannier representation , entering the hamiltonian  ( [ eq : wannier_hamilton ] ) , need to be obtained by a transformation of the corresponding quantities in terms of gaussian basis functions . this transformation is presently the bottleneck in practical computations since the wannier program  @xcite carries it out for clusters only and neither uses translational nor point group symmetry . at present , this permits only the use of a minimal basis set for a three - dimensional lif  crystal . the subsequent correlation calculations with the co - adc program  @xcite are very efficient and as rapid as typical molecular adc calculations . .convergence of the top of the valence bands  @xmath346 , the bottom of the conduction band  @xmath347 , the fundamental band gap  @xmath348 , and the bandwidth of the f@xmath0  valence band complex  @xmath349 of a lif  crystal . `` cells ''  designates the number of unit cells taken into account in co - adc(2,2 ) calculations of the quasiparticle band structure where zero refers to the hartree - fock result . unity denotes the inclusion of the @xmath83- and @xmath84-configurations from the origin cell only in correlation calculations . 13 and 19 indicate the additional inclusion of configurations from nearest- and next - nearest - neighbor cells , respectively . for a lif crystal , @xmath346 and @xmath347 are both situated at the @xmath350  point . all data are given in electronvolts . [ cols="^,^,^,^,^",options="header " , ]     the hartree - fock band structure that results from this minimal basis is shown in fig . [ fig : lif_19cell ] . all band energies are given with respect to the top of the valence bands which is located at the @xmath350  point and is set to zero . below zero there is a complex of three f@xmath0  valence bands shown as a closeup in fig . [ fig : lif_19cell_valence ] . they originate from the three f@xmath0  valence energy levels of the lif molecule around @xmath351 where the lowest ionization potential of the molecule is doubly degenerate . @xcite upon crystallization , this degeneracy is lifted due to interactions with neighboring atoms and , except for the high symmetry lines  @xmath350l and @xmath350x , one observes three distinct f@xmath0  valence bands . as a minimal basis set is utilized , there is only a single conduction band of mainly li@xmath7  character . a slight admixture of f@xmath0  character to this band is observed in the w  k  panel by comparing with the band structure resulting from a larger basis set , e.g. , the one of ref .  , where a conduction band of f@xmath0  character mixes with the li@xmath7  band . nevertheless , the single conduction band represents the energetically lower edge of the conduction - band complex well . lif has a direct band gap ; i.e. , the maximum of the f@xmath0  valence bands and the minimum of the li@xmath7  conduction band are located at the same crystal momentum , the @xmath350  point , here .    - and @xmath84-configurations involving wannier orbitals in the origin cell and the nearest- and next - nearest - neighbor cells . ]    in order to investigate electron correlations , we start from a single unit cell , to form the @xmath83- and @xmath84-configurations entering eq . ( [ eq : m_k_hermitian ] ) for  @xmath116 , respectively . successively , configurations which contain wannier orbitals in nearest- and next - nearest - neighbor cells are considered . this leads to three distinct quasiparticle band structures whose convergence with respect to the number of unit cells where electron correlations are regarded explicitly is monitored for the following key quantities : the top of the valence bands , the bottom of the conduction band , and the band gap and the width of the f@xmath0  valence - band complex [ tab . [ tab : cablif ] ] . the results are compared with the plain hartree - fock data ( with configurations out of `` 0 ''  unit cells ) . taking into account electron correlations of a single unit cell causes a slight upwards shift of the valence bands by  @xmath352 while the conduction band essentially remains unchanged . upon inclusion of configurations in the nearest - neighbor cells , the band structure changes drastically , leading to a significant reduction of the band gap by  @xmath353 . accounting for configurations which extend to next - nearest - neighbor cells has three times less impact on the band gap . therewith , it can clearly be seen that the major contributions of electron correlations have been covered and that the effect of configurations which extend to wannier orbitals beyond next - nearest - neighbor cells will yield much less significant changes of the quantities than those discernible in tab . [ tab : cablif ] . as pointed out in sec . [ sec : onion ] , the effect of the more distant unit cells can be determined using a continuum approximation . @xcite the observed correlation corrections in tab . [ tab : cablif ] with a growing number of neighboring unit cells result primarily from the formation of a quasiparticle and the polarization of the surrounding by the coulomb charge of the added particle or hole . @xcite note that , additionally , the beneficial basis set extension  ( bse )  @xcite effect is noticeable in solids which leads to an improved description of electron correlations with an increasing number of unit cells taken into account .    . ]    the quasiparticle band structure of  lif is displayed in fig . [ fig : lif_19cell ] for the f@xmath0  valence bands and the li@xmath7  conduction band . the former bands are additionally displayed on an enlarged scale in fig . [ fig : lif_19cell_valence ] . the data are taken from the computation with configurations from alltogether 19  unit cells . valence and conduction bands do not shift by the same amount , upwards and downwards , respectively . instead , the former bands are much more strongly influenced by electron correlations than the latter bands . this is due to the fact that the occupied wannier orbitals are spatially less extended than the virtual wannier orbitals . for electrons in compact orbitals , the hartree - fock approximation becomes progressively worse due to an insufficient description of the coulomb hole . @xcite although the quasiparticle bands are considerably shifted with respect to the corresponding hartree - fock bands , they essentially keep their width . we also observe this effect for  ( hf)@xmath354  chains . @xcite this is in contrast to the significant reduction of bandwidths observed for quasiparticle band structures of covalently bonded polymers like _ trans_-polyacetylene  @xcite and covalently bonded crystals like diamond ,  @xcite silicon  @xcite , and germanium . @xcite    the _ ab initio _ description of a lif  crystal using a minimal basis set provides valuable insights . however , we can not expect full quantitative agreement neither at the hartree - fock nor at the correlation level . the accuracy of our approach can be determined by comparing with experimental and theoretical data from related studies . the fundamental band gap is accessible by photoelectron spectroscopy and has been measured to be  @xmath355 . another experimentally accessible quantity , the width of the f@xmath0  valence - band complex , was determined to lie in the range  @xcite of  @xmath356@xmath357 and thus is inconclusive , unfortunately . with our sto-6g - like basis set , we find a hartree - fock band gap of  @xmath358 [ tab . [ tab : cablif ] ] which deviates considerably from the value of  @xmath359 communicated by kunz  @xcite or the value of  @xmath360 reported by albrecht  @xcite who both use larger basis sets . for the width of the f@xmath0  valence - band complex , we find  @xmath361 [ tab . [ tab : cablif ] ] in the hartree - fock approximation which is appreciably larger than the value of  @xmath362 from kunz  @xcite and the value of  @xmath363 from albrecht . @xcite    including electron correlations , the quasiparticle band gap reduces to  @xmath364 [ tab . [ tab : cablif ] ] and thus falls short by  @xmath365 compared to the experimental value . albrecht  @xcite finds @xmath366 , and calculations of kunz  @xcite yield a theoretical band gap of  @xmath367 , in good agreement with the recent measurements of  @xmath355 . for the f@xmath0  valence band complex , we obtain a width of  @xmath368 [ tab . [ tab : cablif ] ] which is nearly unchanged compared to the hartree - fock value , a fact that has also been found in previous studies of bulk  lif . @xcite consequently , our value of the width remains much higher than the value of  @xmath369 from kunz  @xcite or the value of  @xmath370 from albrecht . @xcite    the above results show that the co - adc theory can be applied without problems to compute the energy bands of three - dimensional crystalline solids . what remains to be improved considerably is the transformation of the fock matrix and the two - electron integrals from the gaussian basis set representation to a representation in terms of wannier orbitals . in this article , a general local orbital _ ab initio _ green s function method is devised for band structures . the band structure is given by the pole positions of the one - particle green s function in terms of bloch orbitals which in turn is evaluated efficiently utilizing dyson s equation . to approximate the self - energy up to @xmath16-th order , we devise a crystal orbital formulation of the well - established algebraic diagrammatic construction  ( adc )  @xcite which is termed crystal orbital adc  ( co - adc ) . @xcite the pole search of the one - particle green s function is recast into a hermitian eigenvalue problem which is a numerically stable and efficient formulation , permitting us to explore strong correlations that occur for the energetically lower lying bands of crystals . @xcite wannier orbitals are used in internal summations in the self - energy and allow us to exploit the fact that electron correlations are predominantly local . the lattice summations in the co - adc equations must be truncated to render the problem tractable . to this end , we devise a configuration selection procedure which can equally well be used in conjunction with the calculation of ionization potentials and electron affinities of large molecules to speed up computations . configuration selection is shown to lead to a definition of degeneracy among the states of a crystal . the derivation of the local orbital co - adc theory sets out from the equations in terms of bloch orbitals . afterwards internal indices are transformed to a wannier representation . we consider this line of argument to be compelling due to the close analogy of the equations in crystal momentum representation to the equations of molecular physics . alternatively , the derivation can be conducted by starting from co - adc in terms of wannier orbitals . a transformation of the equations to crystal momentum representation is then carried out by utilizing the wannier transformation  ( [ eq : wannier ] ) . the co - adc computer program  @xcite has been developed for computations based on the co - adc(2,2 ) approximation . it requires the fock matrix and the two - electron integrals in local representation which are obtained by an integral transformation using the wannier orbitals from a hartree - fock calculation with the _ ab initio _ program wannier . @xcite at present , this transformation poses a bottleneck because the current implementation has been designed for molecular clusters and does not exploit any symmetries . therefore , only a minimal basis set is used to study a lithium fluoride crystal in order to demonstrate that the method is working . a fundamental band gap of  @xmath371 and a bandwidth of the f@xmath0  valence - band complex of  @xmath368 is obtained . these values are to compare with the experimental data of  @xmath355 and @xmath356@xmath357 and the theoretical results of  @xmath372 , @xmath367 , and @xmath373 , @xmath374 for the band gap and the bandwidth , respectively . increasing the size of the basis set will improve the co - adc(2,2 ) results . yet it is demonstrated that the co - adc theory is well suited for tackling the problem of calculating energy bands of solids with controlled approximations . we would like to point out that the block of @xmath83-configurations in the band structure matrix  ( [ eq : solvepadc ] ) is much larger than the block of @xmath84-configurations . in order to reduce the computational effort , the advanced part and the retarded part of the one - particle green s function can be independently evaluated diagrammatically . this leads to a different , so - called non - dyson co - adc scheme , which decouples the computation of occupied and virtual bands completely  @xcite and results in two independent hermitian matrix eigenvalue problems to be solved . we would like to thank jochen schirmer for fruitful discussions and critically reading the manuscript . is indebted to thomas sommerfeld for providing his complex - symmetric block - lanczos algorithm and helpful advice .
Abstract:

Answer: using the green s function formalism , an _ ab initio _ theory for band structures of crystals is derived starting from the hartree - fock approximation . it is based on the algebraic diagrammatic construction scheme for the self - energy which is formulated for crystal orbitals  ( co - adc ) . in this approach , the poles of the green s function are determined by solving a suitable hermitian eigenvalue problem . the method is not only applicable to the outer valence and conduction bands , it is also stable for inner valence bands where strong electron correlations are effective . the key to the proposed scheme is to evaluate the self - energy in terms of wannier orbitals before transforming it to a crystal momentum representation . exploiting the fact that electron correlations are mainly local , one can truncate the lattice summations by an appropriate configuration selection scheme . this yields a flat configuration space ; i.e. , its size scales only linearly with the number of atoms per unit cell for large systems and , under certain conditions , the computational effort to determine band structures also scales linearly . as a first application of the new formalism , a lithium fluoride crystal has been chosen . a minimal basis set description is studied , and a satisfactory agreement with previous theoretical and experimental results for the fundamental band gap and the width of the f@xmath0  valence band complex is obtained .