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: cloud radio access networks ( c - rans ) have drawn considerable attention for their potential to sustain the explosive traffic demand in wireless communications . unlike traditional cellular networks , a c - ran splits the low - cost and light - weighted remote radio heads ( rrhs ) from the baseband processing units ( bbus ) , and merges the latter into a data center . the rrhs and bbus are connected by a low - latency , high - bandwidth fiber network . the special architecture of c - ran allows full - scale rrh coordination , which enables flexible interference management , dynamic resource allocation , and collaborative radio technology . this consequently leads to significant capacity enhancement . the full - scale coordination , however , also introduces a severe complexity issue . the state - of - the - art c - ran technology is able to support thousands of rrhs @xcite . full - scale rrh coordination over such a large network involves prohibitively high computational complexity . for example , the linear minimum mean square error ( mmse ) detector requires cubic complexity in the network size ( in terms of the number of rrhs ) , or equivalently a quadratic complexity normalized by the number of rrhs@xcite . this implies that the detection complexity quickly becomes unaffordable as the network size grows . as such , a main challenge of c - ran is to design _ scalable _ coordination algorithms , where _ scalable _ means : 1 ) the performance is near the optimum performance of full - scale rrh coordination , 2 ) the normalized computational complexity per rrh does not grow with the network size , or equivalently , the total computational complexity grows linearly with the network size .    in a c - ran , users and rrhs are scattered over a large area . due to the attenuation effect in the propagation of electromagnetic waves , an rrh usually receives strong signals from only a small number of nearby users . intuitively , ignoring the signals from far - away users in general does not cause much performance loss . as shown in our previous work @xcite , with a distance - threshold - based channel sparsification approach , a vast majority of signals over the transmission links can be ignored if we can tolerate a small degradation in the signal - to - noise - plus - interference ratio ( sinr ) . as such , each rrh only needs to serve its nearby users whose distances to the rrh are below a certain threshold . based on the sparsified channel matrix , @xcite proposes an algorithm that greatly reduces the computational complexity of mmse detection from @xmath0 to @xmath1 , where @xmath2 is the total number of rrhs and @xmath3 $ ] is a constant depending on the computing implementations . however , the algorithm is not perfectly scalable yet , in the sense that the complexity grows faster than linear with the network size .    in this paper , we are interested in designing a perfectly scalable algorithm for joint signal detection in the uplink of c - ran . with channel sparsification @xcite , a c - ran system can be represented by a bipartite random geometric graph , as shown in fig . [ fig : sfg ] . here , rrhs and users are treated as vertices / nodes , and an edge connects an rrh and a user if the distance between them does not exceed the threshold . then , signal detection in a c - ran is converted to a statistical inference problem over a bipartite random geometric graph , where an inference problem is to estimate the signals from unobserved nodes ( i.e. , user nodes ) , conditional on signals from observed nodes ( i.e. , rrh nodes ) . message passing ( a.k.a . belief propagation ) is well known as a low - complexity iterative solution to inference over random sparse graphs @xcite . in a message - passing algorithm , messages are exchanged between nodes with edge connection . thus , the complexity of message passing is proportional to the number of edges in the network . in a c - ran with channel sparsification , the number of messages per rrh is proportional to the number of nearby users in its neighborhood , which does not scale with the network size . thus , the total complexity per iteration of message passing in c - rans is linear in the network size . unfortunately , the convergence of message passing over a bipartite random geometric graph is not guaranteed . the reason is that the random geometric graph is locally dense and always contains loops . it is well - known that message passing is not guaranteed to converge when the graph is loopy @xcite . indeed , numerical simulations indicate a non - trivial probability that the message - passing algorithm for c - rans does not converge .      in this paper , we propose a novel algorithm , referred to as randomized gaussian message passing ( rgmp ) , to address the convergence issue of message passing over a bipartite random geometric graph . unlike conventional message passing with synchronous message updating , the rgmp algorithm updates messages asynchronously with a random sequential order . to the best of our knowledge , this is the first work to introduce randomness to the update schedule of message passing algorithms . intuitively , when messages are exchanged among nodes of a loopy graph , errors may accumulate along loops , which eventually leads to the divergence of the algorithm @xcite . updating messages sequentially in a random order effectively weakens the effect of loops and thus improves the convergence . the randomness of the message update schedule also simplifies the analysis , and allows to derive a necessary and sufficient condition for the expected convergence of the rgmp algorithm . we show by both analysis and numerical results that the rgmp algorithm converges with a much higher probability than conventional message passing . indeed , we have never observed a single divergence in our simulations when the network size is moderately large ( i.e. , when the network has more than five rrhs ) . moreover , our numerical results indicate that the number of iterations of rgmp does not increase with the network size . this implies that the total computational complexity is linear with the network size , and thus the algorithm is perfectly scalable .      in the previous work @xcite , we proposed a dynamic clustering algorithm to reduce the computational complexity of the mmse detector . the complexity of the algorithm is reduced from cubic to no more than quadratic in the number of rrhs . in @xcite , shi _ et al . _ presented a two - stage approach to solve large - scale convex optimization problems for dense wireless cooperative networks , such as c - rans . matrix stuffing and alternating direction method of multipliers ( admm ) were used to speed up the computation . particularly , it was shown in @xcite that the expected output of randomly permuted admm converges to the unique solution of the optimal linear detector . however , we will see that the admm algorithm converges much more slowly than the proposed rgmp algorithm when applied to large networks like c - rans . in addition , there are other iterative solvers of the linear mmse detector , such as the preconditioned conjugate gradient ( pcg ) method @xcite and the generalized approximate message passing ( gamp ) algorithm @xcite . we show by numerical simulation that the number of iterations needed for convergence in these algorithms are roughly linear in the network size . this translates to quadratic computational complexity in total . therefore , all these algorithms are not scalable . the rest of the paper is organized as follows . in section ii , we describe the system model . in section iii , we introduce a gaussian message - passing algorithm with channel sparsification for signal detection in c - rans with linear complexity per iteration , and then discuss the convergence issue . in section iv , we propose the rgmp algorithm to address the convergence issue of gaussian message passing . in section v , the convergence condition of the rgmp algorithm are analysed . in section vi , simulation results are demonstrated to compare the rgmp algorithm with other existing algorithms . conclusions and future works are discussed in section vii . in this paper we consider the uplink transmission of a c - ran with @xmath2 single - antenna rrhs and @xmath4 single - antenna users . suppose that both the rrhs and the users are randomly located over an area . let @xmath5 be the signal transmitted by user @xmath6 , and @xmath7 be the received signal at rrh @xmath8 . denote @xmath9^t$ ] and @xmath10^t$ ] . then , the received signal vector @xmath11 at the rrhs is @xmath12 where @xmath13 denotes the channel matrix , with the @xmath14-th entry @xmath15 being the channel coefficient between the @xmath6-th user and the @xmath8-th rrh ; @xmath16 is the transmission power allocated to each user ; and @xmath17 is a noise vector received by the rrhs . the transmitted signals are assumed to have zero mean and unit variance , i.e. , @xmath18=\mathbf{0}$ ] and @xmath19=\mathbf{i}$ ] . we further assume @xmath20 , where @xmath21 is the i.i.d . rayleigh fading coefficient with zero mean and unit variance , @xmath22 is the distance between the @xmath8-th rrh and the @xmath6-th user , and @xmath23 is the path loss exponent . here , @xmath24 is the path loss from the @xmath6-th user to the @xmath8-th rrh .    in this paper , we employ linear mmse detection to estimate the transmitte signal vector @xmath25 , with the decision statistics given by @xmath26 in the above , the inversion of the @xmath27 matrix @xmath28 requires computational complexity of @xmath0 . this complexity is prohibitively high for a large - scale c - ran with hundreds and thousands of rrhs , thus posing a serious scalability problem . in what follows , we endeavour to develop a scalable algorithm to detect @xmath25 with complexity @xmath29 under the assumption that @xmath4 grows at the same rate as @xmath2 . in other words , the average computational complexity per rrh ( or per unit network size ) does not scale with @xmath2 . in this section , we first describe the channel sparsification approach introduced by the authors in @xcite to model a c - ran as a bipartite random geometric graph . then , we apply the gaussian message - passing algorithm over bipartite random geometric graphs for signal detection . we borrow the channel sparsification approach in our recent work @xcite to sparsify the channel matrix , as described below . the entries of @xmath30 are discarded based on the distances between rrhs and users . specifically , the @xmath31-th entry in the resulting sparisified channel matrix @xmath32 is given by @xmath33 where @xmath34 is a distance threshold . given the sparsified channel matrix @xmath32 , the received signal can be represented as @xmath35 where @xmath36 the mmse estimator of @xmath25 is approximated by @xmath37 with @xmath38 + n_0 $ ] for arbitrary rrh @xmath8 . as proven in @xcite , the channel matrix can be sparsified without considerably compromising the sinr . the reason is that as the rrhs and users are uniformly distributed over a large area , an rrh can only receive reasonably strong signals from a small number of nearby users , and vice versa . therefore , the majority of the elements of @xmath30 are relatively small in magnitude , and ignoring them in signal detection leads to marginal loss in the overall system performance . indeed , according to @xcite , when @xmath2 scales in the same order as @xmath4 , the distance threshold @xmath34 does not increase with the network size to achieve a certain sinr performance . thus , in this paper , we assume that @xmath34 is a predetermined constant regardless of the network size . this implies that the average number of users connecting to an rrh does not scale with the network size . channel sparsification simplifies the signal detection in a c - ran to an inference problem over a bipartite random geometric graph ( see fig . [ fig : sfg ] ) . in the bipartite random geometric graph , rrhs and users in a c - ran are referred to rrh nodes and user nodes respectively , and edge connections exist only between rrh nodes and user nodes . more specifically , an rrh node is connected to a user node only if the distance between them falls below the threshold @xmath34 , and the weight over such an edge is the channel coefficient from the corresponding user to the corresponding rrh . suppose that the entries in @xmath25 follow an independent complex gaussian distribution . does not follow a gaussian distribution , the message - passing algorithm presented in this paper gives an approximation of the linear mmse estimation @xcite . ] then , @xmath39 and @xmath25 are jointly gaussian , and therefore the mmse detector in ( [ eqn : x ] ) is also the maximum _ a posteriori _ probability ( map ) detector that maximizes the _ a posteriori _ probability @xmath40 . that is , @xmath41 the probability density function @xmath40 can be factorized as @xmath42 recall that we sparsify the channel matrix by using the channel sparsification approach given in @xcite . based on ( [ eqn : app_x ] ) , the factorization of @xmath40 is approximated as @xmath43 where @xmath44 contains all @xmath45 with @xmath46 and @xmath47 is the set of user indices with @xmath48 . we now transfer the bipartite random geometric graph to a factor graph with the factorization in ( [ eqn : pdf ] ) . as illustrated in fig . [ fig : sfg_2 ] , a factor graph is also a bipartite graph comprising two types of nodes , namely , variable nodes ( denoted by circles ) and check nodes ( denoted by squares ) , together with edges connecting these two types of nodes . the relation between the factorization ( [ eqn : pdf ] ) and its associated factor graph is as follows . a check node @xmath49 is connected to a variable node @xmath5 by an edge when there is an edge connecting the @xmath8-th rrh node and the @xmath6-th user node in the corresponding random geometric graph ( i.e. , @xmath48 ) , or equivalently , when the function @xmath49 takes @xmath5 as input . @xmath32 , @xmath39 @xmath50 for all @xmath6 initial @xmath51 , @xmath52 for all @xmath53 set @xmath54 for all @xmath55 such that @xmath56 , compute @xmath57 @xmath58 @xmath59 @xmath60 compute @xmath61 @xmath62    we are now ready to introduce the gaussian message - passing algorithm for signal detection . the algorithm will be implemented in the centralized data center . the messages , namely , the marginals of @xmath63 and @xmath64 , are exchanged along the edges . in this paper , both @xmath63 and @xmath64 are gaussian distributed , and therefore the messages are gaussian probability density functions and can be completely characterised by mean and variance . denote by @xmath65 and @xmath66 the mean and variance sent from check node @xmath49 to variable node @xmath5 at iteration @xmath67 , respectively , and denote by @xmath68 and @xmath69 the mean and variance sent from variable node @xmath5 to check node @xmath49 at iteration @xmath67 , respectively . the detailed steps of message passing are presented in algorithm 1 . we refer to this algorithm as gaussian message passing ( gmp ) , as all the messages involved are gaussian marginals . note that each rrh only serves users located in a circle with a constant radius @xmath34 . thus , the average number of messages to be exchanged and computed at each node does not scale with the network size . therefore , the complexity per iteration of the gmp algorithm is linear in the number of rrhs and users .    in spite of its linear complexity per iteration , the gmp algorithm is not guaranteed to converge on the factor graphs induced by c - rans . it is known that the gaussian message - passing algorithm always converges to the optimal solution on a tree - type factor graph @xcite . it is also known that , if a factor graph is random and sparse enough , the corresponding message - passing algorithm converges asymptotically as the network size grows to infinity @xcite . however , the factor graph for a bipartite random geometric graph induced from a c - ran is locally dense and far from being a tree . this is due to the fact that every rrh needs to simultaneously serve multiple nearby users . for example , @xmath70 in fig . [ fig : sfg ] form a loop of length 4 . indeed , we observe in simulations that the gmp algorithm diverges in c - ran with a non - trivial probability . even worse , the probability of divergence grows with the network size , as illustrated later in fig . [ fig : prob_con ] . we focus on improving the convergence performance of gmp in the rest of the paper . the gmp algorithm for a c - ran with channel sparsification can be simply extended to the case without channel sparsification by setting the distance threshold to infinity . however , this leads to an increase of the computational complexity per iteration . we see that in each iteration of algorithm 1 , messages need to be updated on every edge of the factor graph . from channel randomness , the entries of @xmath30 are non - zero with probability one . thus , in the factor graph without channel sparsification , every rrh check node @xmath71 is connected to all user variable nodes @xmath72 . this implies that the total number of edges in the factor graph is @xmath73 , implying that the complexity of the gmp algorithm is @xmath74 per iteration , which is unaffordable for a large - scale c - ran . @xmath32 , @xmath39 @xmath50 for all @xmath6 initialize @xmath51 , @xmath52 for all @xmath53 . set @xmath54 . pick a permutation @xmath75 of @xmath76 uniformly at random . for @xmath77 , and @xmath78 , compute @xmath79 @xmath80 @xmath81 @xmath82 compute @xmath83 @xmath84      in this section , we propose the rgmp algorithm to address the convergence issue of gmp . the main novelty of the rgmp algorithm is on the scheduling strategy for message updating . the conventional gmp algorithm employs synchronous message passing , i.e. , messages are updated in parallel . as aforementioned , synchronous message passing does not work well in c - rans due to local loops in the factor graph . to address this issue , we propose the rgmp algorithm with random asynchronous scheduling , i.e. , messages are updated sequentially in a randomly permuted order .    the rgmp algorithm is described as follows . define @xmath85 as @xmath86 at each iteration , we draw a permutation @xmath75 of @xmath87 uniformly from @xmath85 , and update the messages at the user variable node side in the order of @xmath75 . for example , at the @xmath67-th iteration , when @xmath88 and @xmath89 , we first update all the messages on the edges connecting the user variable node @xmath90 . then , the messages on the edges connecting the user variable node @xmath91 are updated . finally , messages related to variable node @xmath92 are updated . the rgmp algorithm is given in algorithm 2 . in this subsection , we use a toy example to illustrate the difference between our proposed rgmp algorithm and synchronous / asynchronous gmp . consider the following randomly generated channel matrix @xmath93,\ ] ] and let the transmit snr ( i.e. , @xmath94 ) be @xmath95db . the corresponding received signal @xmath39 is @xmath96^t.\ ] ] for fairness of comparison , we do not conduct channel sparsification in this example . that is , the distance threshold is set to infinity . [ fig : err_mp ] plots the relative error versus the number of iterations for the rgmp algorithm and the gmp algorithm with different message update strategies , i.e. , synchronous update and asynchronous update with different fixed orders , @xmath97 and @xmath98 . the relative error is defined as @xmath99 , where @xmath100 is the estimation of the transmitted signal after the @xmath67-th iteration . we see that the synchronous gmp algorithm and the asynchronous one with order @xmath98 diverge , but the asynchronous gmp with order @xmath97 and the proposed rgmp algorithm converge .    the examples in fig . [ fig : err_mpc ] and [ fig : err_mpd ] show that convergence of asynchronous gmp heavily depends on the update order . unfortunately , there is no systematic way to derive a fixed update order that guarantees convergence . in general , finding such an update order is difficult , especially in large networks . this issue is avoided in the proposed rgmp algorithm by randomizing the update schedule instead of fixing one . indeed , the randomization significantly weakens the loopy effect of the graph , and thus convergence is almost ensured in rgmp .    in fig . [ fig : prob_con ] , we plot the empirical probability of convergence against the network size , where users and rrhs are uniformly located in a circular network area with user density @xmath101 and rrh density @xmath102 . the distance threshold @xmath34 is @xmath103 m . for each simulated point in fig . [ fig : prob_con ] , both gmp and rgmp are run for over @xmath104 times that are randomized over both rrh / user location and channel fading . for gmp , the convergence probability decreases when the network size becomes large . in contrast , no divergence has been observed for the rgmp algorithm throughout our simulations . in this section , we establish a necessary and sufficient condition for the expected convergence of the proposed rgmp algorithm . for self - containedness , we start with existing results on the analysis of the convergence condition for conventional gmp . the factor graph of a c - ran contains loops with high probability . the convergence of gmp on a loopy factor graph has been previously studied in @xcite , with the main result summarized below .    from algorithm 1 , we see that the evolution of the variances @xmath105 is independent of the means @xmath106 , @xmath107 and the received signal @xmath39 . substituting ( [ eqn : v2 ] ) into ( [ eqn : v1 ] ) , we obtain @xmath108 denote ( [ eqn : evo_v ] ) in a vector form as @xmath109 where @xmath110 is the evolution function determined by ( [ eqn : evo_v ] ) , and @xmath111 is a vector consisting of @xmath112 for all @xmath8 and @xmath6 with @xmath113 . note that @xmath110 is a standard function , the definition of which is given below . a function @xmath114 is standard if for all @xmath115 the following properties are satisfied .    * _ positivity : _ @xmath116 . * _ monotonicity : _ if @xmath117 , then @xmath118 . * _ scalability : _ for all @xmath119 , @xmath120 .    from theorem 5.1 of @xcite , the sequence of @xmath111 always converges to a unique fixed point if the initial point @xmath121 and the evolution function is standard . we now consider the convergence of means . a vector of means , @xmath122 , is constructed with its @xmath123-th entry being @xmath124 the recursion of the means is given by ( [ eqn : m1 ] ) and ( [ eqn : m2 ] ) . as the variances always converge , the evolution of the means can be written as follows : @xmath125 where @xmath126 is an @xmath127 vector with its @xmath123-th entry being @xmath128 and @xmath129 is an @xmath130 matrix with the @xmath131-th entry being @xmath132 with @xmath133 and @xmath134 . then , a necessary and sufficient condition for the convergence of ( [ eqn : evo_m ] ) is given in theorem 5.2 , @xcite . that is , in algorithm 1 , the sequence of @xmath122 converges to a unique fixed point if and only if the spectral radius @xmath135 .      in this subsection , we first show that the message variances always converge in the rgmp algorithm . then , we focus on the convergence condition of the means in rgmp . recall that the evolution function of the variances in ( [ eqn : evo_v ] ) is a standard function . as proven in @xcite , if the evolution function of a synchronous algorithm is standard , then the corresponding asynchronous algorithm converges . based on that , we obtain the following theorem .    in the rgmp algorithm , the sequence of @xmath112 always converges to the same unique fixed point as in algorithm 1 if the initial point @xmath121 .    with theorem 1 , it suffices to focus on the convergence condition of the means in the rgmp algorithm . we first derive the evolution function of the means . denote the permutation at the @xmath67-th iteration as @xmath136 . combining ( [ eqn : rm1 ] ) and ( [ eqn : rm2 ] ) , we obtain the evolution of means @xmath122 as @xmath137 where @xmath138 is an @xmath139 subvector of @xmath122 with the @xmath8-th entry being @xmath140 and @xmath141 is an @xmath139 subvector of @xmath142 with the @xmath8-th entry being @xmath143 @xmath144 is the @xmath27 evolution matrix from user @xmath145 to user @xmath146 with the @xmath147-th entry being @xmath148 @xmath144 is the @xmath149-th submatrix of @xmath150 . more specifically , @xmath151        \mathbf{\omega}_{2,1 } & \mathbf{0 } & \mathbf{\omega}_{2,3 } & \cdots & \mathbf{\omega}_{2,k}\\[0.6em ]      \vdots     & \ddots & \ddots & \ddots&\vdots\\[0.6em ]         \mathbf{\omega}_{k-1,1 } & \cdots & \mathbf{\omega}_{k-1,k-2 } & \mathbf{0}&\mathbf{\omega}_{k-1,k}\\[0.6em ]         \mathbf{\omega}_{k,1 } & \mathbf{\omega}_{k,2}&\cdots&\mathbf{\omega}_{k-1,k}&\mathbf{0 }       \end{bmatrix}.\ ] ]    we can further rewrite the equation ( [ eqn : random_mp ] ) as @xmath152 where @xmath153_{i , j}\in \mathbb{c}^{nk\times nk}$ ] with its @xmath149-th submatrix being @xmath154 and @xmath155_{i , j}\in \mathbb{c}^{nk\times nk}$ ] with its @xmath149-th submatrix being @xmath156 based on the definition of @xmath157 , the determinant of @xmath157 is always @xmath158 or @xmath159 . it implies that @xmath157 is nonsingular . then , multiplying both sides of ( [ eqn : evo_m_r ] ) by @xmath160 , we obtain @xmath161 consequently , we obtain the following condition for the convergence of the rgmp algorithm .    for a given sequence of permutations @xmath162 , the rgmp algorithm converges to the fixed point @xmath163 if and only if @xmath164 , where @xmath136 is the permutation at the @xmath67-th iteration .    for an arbitrary permutation @xmath75 , the fixed point of @xmath165 is given by @xmath166 substituting @xmath167 into ( [ eqn : update_3 ] ) , we obtain @xmath168 . clearly , @xmath169 is independent of the choice of permutation @xmath75 . define @xmath170 . then , @xmath171 , for any iteration number @xmath67 . by recursion , we obtain @xmath172 . therefore , @xmath173 provided @xmath174 as @xmath175 . this concludes the proof . proposition 2 discusses the convergence for a given sequence of permutations @xmath162 . to quantify the average performance of rgmp over random permutations , we consider expected convergence in the following , where the expectation is taken over all possible permutations . let the expected output be @xmath176,\ ] ] where @xmath177 is the set of the permutations after iteration @xmath67 . we are now ready to present a necessary and sufficient condition for the convergence of @xmath178 . the expected output @xmath179 $ ] converges to the unique point @xmath163 if and only if the spectral radius @xmath180 , where @xmath181=\frac{1}{k!}\sum_{\sigma \in \sigma}(\mathbf l_{\sigma}^{-1}\mathbf r_{\sigma } ) . \label{eqn : condition } \ ] ]    denote @xmath182 $ ] . based on the definition of @xmath183 , we obtain @xmath184\\ & = & & \mathrm{e}_{\sigma_t}\left[\mathrm{e}_{\xi_{t-1}}\left[\mathbf l_{\sigma_t}^{-1}\mathbf r_{\sigma_t}\mathbf{m}^{(t)}+\mathbf{l}_{\sigma_t}^{-1}\mathbf{z}\right]\right]\\ & = & & \mathbf{\lambda}\bm \phi^{(t)}+\mathbf{a}\mathbf{z}. \end{aligned}\ ] ]    from theorem 5.3 in @xcite , the sequence of @xmath178 converges to the fixed point @xmath185 if and only if the spectral radius @xmath180 . then , it suffices to show that @xmath186 . note that @xmath187 . substituting @xmath188 into ( [ eqn : condition ] ) , we obtain @xmath189\\ & = & & \mathrm{e}_{\sigma}[\mathbf l_{\sigma}^{-1}(\mathbf{\omega}-\mathbf{i})+\mathbf{i}]\\ & = & & \mathbf{a}(\mathbf{\omega}-\mathbf{i})+\mathbf{i}. \end{aligned}\label{eqn : proofa}\ ] ] recall that @xmath180 , and thus @xmath190 is nonsingular . together with @xmath191 from ( [ eqn : proofa ] ) , we see that both @xmath192 and @xmath193 are nonsingular . hence , @xmath194 which concludes the proof .    from theorem 1 , both gmp and rgmp have guaranteed convergence for variances . however , they have different conditions to ensure the convergence of means : gmp requires @xmath195 while rgmp requires @xmath196 . since both @xmath150 and @xmath197 highly depend on the network geometry , it is difficult to theoretically compare these two convergence conditions . to shed light on the difference , we plot the cumulative distribution function ( cdf ) of the spectral radius of @xmath150 and @xmath197 in fig .  [ fig : cdf_rho ] . we assume that users and rrhs are randomly located in a circular network area with radius @xmath198 . the user density is @xmath199 and the rrh density is @xmath200 . the path loss exponent is @xmath201 , and the average transmit snr at the user side equals to @xmath95db . that is @xmath202db . we see that @xmath203 is more likely to take small values than @xmath204 . this implies that rgmp converges with a higher probability than gmp . indeed , we have run over @xmath205 times for each setting , and @xmath196 for all the cases .      in this subsection , we present some special cases to show that the convergence condition in theorem 3 is less stringent than that of the synchronous gmp . when @xmath150 is a hermitian matrix , the condition of expected convergence in theorem 2 can be further simplified as given in the following corollary . the expected output @xmath179 $ ] converges to the unique point @xmath163 if @xmath150 is a hermitian matrix , and all the eigenvalues of @xmath150 is no greater than @xmath158 .    if @xmath150 is hermitian and all the eigenvalues of @xmath150 is no greater than @xmath158 , the matrix @xmath193 is a positive - semidefinite matrix and can be rewritten as @xmath206 . as proven by lemma 2 in @xcite , all eigenvalues of @xmath207 , or equivalently @xmath208 , lie in @xmath209 . substituting @xmath188 into ( [ eqn : condition ] ) , we obtain @xmath210 as all the eigenvalues of @xmath208 lie in @xmath209 , @xmath196 holds all the time .    obviously , the condition that all the eigenvalues of @xmath150 is no greater than one is much weaker than @xmath135 . it implies that in special cases where @xmath150 is hermitian , @xmath196 is a much less stringent condition than @xmath135 , and thus our rgmp algorithm converges with a higher probability than gaussian message passing . we now consider the c - ran with two users . we show that if the gmp algorithm converges , the expected output @xmath211 of rgmp always converges , as proven in the following corollary . when the number of users @xmath212 , the expected output @xmath211 converges to the unique point if and only if the spectral radius @xmath213 .    from theorem 2 , it suffices to show that @xmath214 . note that there are only two different permutations when @xmath212 . that is , @xmath215 . for simplicity of notation , let @xmath216 and @xmath217 . then , @xmath218 . noting @xmath219 for all @xmath220 , we obtain @xmath221 by using taylor series expansion . similarly , @xmath222 . then @xmath223    denote by @xmath224 the eigenvalues of @xmath150 . then , the corresponding eigenvalues of @xmath225 are @xmath226 . then , the spectral radius of @xmath150 and @xmath225 are @xmath227 and @xmath228 respectively . thus , @xmath213 always holds if @xmath195 . therefore , the condition for expected convergence of the rgmp algorithm is less stringent than that of the synchronous message passing . in this section , we compare the performance of rgmp with other existing algorithms . unless specified otherwise , we assume that both users and rrhs are uniformly at random located in a circular network area with user density @xmath199 and rrh density @xmath200 . the path loss exponent is @xmath201 , and the average transmit snr at the user side equals to @xmath229db . that is @xmath230db . moreover , the stopping criteria is @xmath231 , where @xmath232 is the relative error after the @xmath67-th iteration . in particular , @xmath233 , with @xmath100 being the estimated transmitted signal after @xmath67 iteration .      in this subsection , we compare the converge speed of rgmp with other algorithms including admm @xcite , gamp @xcite , and pcg @xcite . for fairness , the channel sparsification approach with distance threshold @xmath234 m is adopted in all algorithms . in this way , all the algorithms have a linear per - iteration computational complexity with the network size . thus , we only focus on the convergence speed of these algorithms .    in fig . [ fig : err_admm ] , the relative error @xmath232 is plotted against the number of iterations for @xmath235 and @xmath236 . we see that both rgmp and pcg converge very fast . for example , the relative error of rgmp reduces to @xmath237 within @xmath238 iterations . however , the performance of the admm algorithm is unsatisfactory . around @xmath239 iterations are needed for the admm algorithm to reduce the relative error to @xmath240 . in fact , from simulation results not presented here , admm requires over @xmath241 iterations on average to reduce the relative error to @xmath237 for the network configuration in fig . [ fig : err_admm ] . therefore , even though the computational complexity per iteration of admm is linear in the number of rrhs and the expected convergence is guaranteed @xcite , it is still impractical to adopt the admm algorithm in c - ran due to the extremely slow convergence .    in fig . [ fig : rate_size ] , we plot the convergence speed of the rgmp algorithm against the network size , where the convergence speed is measured by the critical number of iterations to achieve @xmath242 . due to the extremely slow convergence speed of admm as shown in fig . [ fig : err_admm ] , we ignore admm and only plot the convergence speed of pcg and gamp for comparison . we observe that the number of iterations needed by both pcg and gamp grow roughly linearly with the network size . in contrast , the convergence speed of rgmp is constant with the network size . note that the computational complexity per iteration of pcg / gamp / rgmp is linear in the network size . thus , the total computational complexity of the proposed rgmp algorithm is linear in the network size , while that of pcg and gamp grows quadratically with the network size .      in this subsection , we compare the performance of the rgmp algorithm with a disjoint clustering algorithm . the disjoint clustering algorithm divides the whole network into disjoint square clusters with area @xmath243 , and mmse detection is done independently in each disjoint cluster . channel sparsification is also applied in the disjoint clustering algorithm . in fig . [ fig : disjoint ] , we plot the mean squared error ( mse ) against the distance threshold , where mse refers to @xmath244 $ ] . the network area is 200km@xmath245 . thus , the numbers of rrhs and users are @xmath246 and @xmath247 , respectively . we see that the gap between the rgmp algorithm and the disjoint clustering algorithm is very large . for example , when the distance threshold @xmath34 is @xmath248 m , the mse of the rgmp algorithm is less than @xmath249 , which is only half of the mse of the disjoint clustering algorithm with cluster area @xmath250km@xmath245 . in this paper , we proposed the rgmp algorithm for scalable uplink signal detection in c - rans . with channel sparsification , signal detection in a c - ran was converted to an inference problem over a bipartite random geometric graph . a random message - update schedule was employed to address the convergence issue of gmp over a bipartite random geometric graph . we analysed the convergence condition of the proposed rgmp algorithm and showed that the convergence condition of rgmp is much less stringent than that of gmp . moreover , numerical results demonstrated that rgmp exhibits much faster convergence than existing algorithms , such as pcg , gamp , and admm . the work in this paper sheds light on the design of message - passing algorithms on general loopy graphs , which has been a challenging topic in the field for years . future work can be done in a number of interesting directions , including the extension to the downlink processing of c - rans and the design of signal detectors with limited fronthaul capacity . y.  shi , j.  zhang , b.  odonoghue , and k.  b.  letaief , `` large - scale convex optimization for dense wireless cooperative networks , '' _ ieee trans . signal process . 18 , pp . 4729 - 4743 , sept . 2015 . r.  barrett , m.  berry , t.  chan , j.  demmel , j.  donato , j.  dongarra , v.  eijkhout , r.  pozo , c.  romine , and h.  van  der  vorst . _ templates for the solution of linear systems : building blocks for iterative methods . _ philadelphia , pa : siam , 1994        t.  j.  richardson , m.  a.  shokrollahi , and r.  l.  urbanke , `` the capacity of low - density parity - check codes under message - passing decoding , '' _ ieee trans . inf . theory _ , 599 - 618 , february 2001 .
Abstract:

Answer: cloud radio access network ( c - ran ) is a promising architecture for unprecedented capacity enhancement in next - generation wireless networks thanks to the centralization and virtualization of base station processing . however , centralized signal processing in c - rans involves high computational complexity that quickly becomes unaffordable when the network grows to a huge size . this paper endeavours to design a _ scalable _ uplink signal detection algorithm , in the sense that the complexity per unit network area remains constant when the network size grows . to this end , we formulate the signal detection in c - ran as an inference problem over a bipartite random geometric graph . by passing messages among neighboring nodes , message passing ( a.k.a . belief propagation ) provides an efficient way to solve the inference problem over a sparse graph . however , the traditional message - passing algorithm is not guaranteed to converge , because the corresponding bipartite random geometric graph is locally dense and contains many short loops . as a major contribution of this paper , we propose a randomized gaussian message passing ( rgmp ) algorithm to improve the convergence . instead of exchanging messages simultaneously or in a fixed order , we propose to exchange messages asynchronously in a random order . numerical results show that the proposed rgmp algorithm has significantly better convergence performance than conventional message passing . the randomness of the message update schedule also simplifies the analysis , and allows the derivation of the convergence conditions for the rgmp algorithm . cloud - ran ; signal processing ; message passing ; belief propagation


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: much of the computational effort in solving the kohn - sham equations of density functional theory ( dft)@xcite is spent in iteratively updating the screened effective potential . this procedure requires the calculation of the electronic density , obtained over a sum of all occupied kohn - sham eigenstates . the required number of states ( bands ) therefore scales with the number of atoms . for the calculation of optical or transport properties , only the bands close to the band gap are required , which represents , for a large number of atoms , only a small fraction of the total number of bands required for the calculation of the total energy in a self - consistent formalism . an approach that bypasses the calculation of the entire spectrum and focuses on states close to the band gap is the semi - empirical pseudopotential approach.@xcite in this approach the pseudopotential is derived from a dft calculation in the local density approximation ( lda ) and subsequently empirically corrected to reproduce experimental _ target _ properties , such as the optical band gap . these semi - emipirical pseudpotentials ( seps ) were derived for cdse and inp and their accuracy demonstrated.@xcite in this work , we follow a similar procedure for the derivation of nitride seps . former treatment of nitrides were at the level of empirical pseudopotentials ( as opposed to seps ) which assume an empirical analytic functional dependence ( a sum of gaussians ) instead of being obtained from dft - lda and include only a local component , while the seps are angular momentum dependent non - local ( or semi - local ) potentials.@xcite these empirical pseudopotentials have a parametric strain dependence that is important and more intricate in the case of nitrides than for the case of iii - v semiconductors.@xcite our derivation of seps for nitrides leads to pseudopotentials that are slightly more expensive computationally ( they have an angular momentum dependent semi - local component and a slightly larger energy cut - off ) but more accurate and transferable . there is no need for an explicit empirical parametric strain dependence . the advantages of seps compared to the local empirical pseudopotentials have been discussed earlier , generally to zinc - blende semiconductors.@xcite in this work we have extended the sep formulation to wurtzite semiconductors , and more specifically to gan .    the second development that was necessary to properly address nitride nanostructures was the problem of the passivation . in earlier work @xcite the dangling bonds of nanostructures were passivated by empirical gaussian potentials . the three parameters describing the passivant , namely the width and height of the gaussian and its distance to the passivating atom , were varied until the surface states , that fill the gaps of most unpassivated semiconductors , disappears . this fitting procedure was originally done manually and subsequently automatized using a minimization procedure.@xcite even in the case of the exhaustive numerical search of the configuration space , some materials could not be properly passivated . this fact is probably due to the lack of flexibility in the gaussian function used . in this paper we describe our approach to obtain passivant potentials . it is based on the dft - lda calculation of a passivated and an unpassivated ( bare ) slab . from the difference of both calculations we extract an effective passivant potential in real space . we can accurately fit the numerical results to a yukawa - type potential . once the bulk seps for a given material are found , together with its corresponding passivants pseudopotentials , the study of the electronic structure and related properties of nanostructures can be performed . we have chosen as an example to apply the sep method to gan nanowires,@xcite a promising nanostructure in the field of solar cells and optoelectronic devices due to the wide band gap of gan ( 3.5 ev ) , that together with indium nitride , a narrow band gap semiconductor ( 0.67 ev ) can cover the full visible solar spectrum by a proper alloying.@xcite more specifically , we explore the effects of confinement in the valence and conduction band states close to the nanowire band gap and we discuss the influence of confinement on the optical properties . first , we solve the kohn - sham equations for the bulk wurtzite structure within the lda using standard _ ab initio _ norm - conserving nonlocal pseudopotentials.@xcite    @xmath0 \psi_i(\bm{r } ) = \varepsilon_i\psi_i(\bm{r } ) . \label{kseq}\ ] ]    the kohn - sham potential @xmath1 includes the local and nonlocal ionic pseudopotentials and the electron - electron interaction :    @xmath2    the ensuing local part of the self - consistent screened effective potential is defined as : @xmath3 where @xmath4 is the local part of the non - screened ionic pseudopotential and the term @xmath5 includes the hartree and exchange - correlation contributions . a more detailed discussion concerning the nonlocal part of the pseudopotential can be found elsewhere.@xcite the potential @xmath6 is a periodic function and can be expanded in a fourier series : @xmath7 where @xmath8 is the set of the reciprocal lattice vectors of the corresponding crystal structure with unit cell volume @xmath9 . we now express the periodic potential @xmath6 as a sum over atom - centered potentials @xmath10 : @xmath11 where we use capital letters for crystal potentials and lower case letters for the corresponding atomic potentials . the vectors @xmath12 are the basis vectors of the crystal unit cell and @xmath13 are the atomic positions of the atoms @xmath14 in the unit cell corresponding to @xmath12 . the goal now is to obtain an expression for @xmath15 that we could use in other environments , as give for instance in nanostructures . these _ screened _ atomic pseudopotentials @xmath15 contain , in addition to the ionic part , in an average way , the contribution related to the electron - electron interaction , expressed by @xmath16 . however , the procedure of obtaining these atomic potentials from the total local potential must be clearly specified and checked for consistency . in the following , we explain how the construction of the screened atomic pseudopotential is carried out . first , by combining eqs .   and , we obtain the following connection between the unknown @xmath17 and the known self - consistent fourier coefficients @xmath18 , for every value of @xmath19 : @xmath20 with @xmath21 where the integral is over the atomic volume @xmath22 , as a consequence of eq .   and after ref .  . the lattice parameters are taken as @xmath23  , and @xmath24 . the wurtzite unit cell has four atoms , with the primitive vectors @xmath25 and the atom positions of the two anions @xmath26 and two cations @xmath27 @xmath28 \bm{\tau}_{a_2}&=+\frac{1}{6}\bm{a}_1-\frac{1}{6}\bm{a}_2+\frac{1}{16}\bm{a}_3\longrightarrow -\bm{\tau}_{c},\\[2pt ] \bm{\tau}_{c_1}&=-\frac{1}{6}\bm{a}_1+\frac{1}{6}\bm{a}_2-\frac{1}{16}\bm{a}_3\longrightarrow \bm{\tau}_{c } ,   \\[2pt ] \bm{\tau}_{c_2}&=+\frac{1}{6}\bm{a}_1-\frac{1}{6}\bm{a}_2+\frac{7}{16}\bm{a}_3\longrightarrow   -\bm{\tau}_{a},\\ \end{aligned}\ ] ] we then obtain for @xmath18 : @xmath29 , \end{aligned}\ ] ] with the symmetric ( @xmath30 ) and antisymmetric ( @xmath31 ) pseudopotentials @xmath32 we obtain @xmath33\\ \\ & + & i\frac{1}{\omega_c}v_-(\bm{g})\left[\sin\left(\bm{\tau}_c\cdot\bm{g}\right)-\sin\left(\bm{\tau}_a\cdot\bm{g}\right)\right].\\ \end{aligned}\ ] ]    since @xmath34 have inversion symmetry , @xmath35 are real , and they can be written in terms of the real and imaginary parts of @xmath18    @xmath36    the potential @xmath37 is calculated from the fourier transform of the lda crystal potential , @xmath38 . once @xmath35 are obtained , the reciprocal - space atomic pseudopotentials @xmath34 are obtained from eq .  ( 11 ) . through this procedure we can obtain @xmath39 only at discrete values @xmath8 . to palliate this problem the lda potential of the bulk system is solved for a set of unit cell volumes , at different lattice constants , around the optimized value , as shown below . so far the procedure to obtain @xmath40 is exact . however , for an efficient implementation it is convenient to assume that the screened pseudopotentials are well represented by their spherically averaged counterparts @xmath41 . at this stage we have a spherically symmetric screened atomic lda potential @xmath42 .     and @xmath43 for gan as defined in the main text . we have obtained five sets of data points by varying the lattice constant @xmath44 as indicated in the legends . the @xmath45 ratio is kept constant at the ideal value . the lines represent a fit to the date points using cubic splines.,width=7 ]    figures  [ vplusminus](a ) and [ vplusminus](b ) show the lda results ( symbols ) for @xmath30 and @xmath31 . in order to have a dense grid of @xmath19 points the calculations are performed for five different lattice constants @xmath44 . the data points from different unit cell volumes fall nearly over the same line , with some slight deviations for the data points around @xmath46  bohr@xmath47 in @xmath30 . this shows that the spherical approximation is justified . concerning the pseudopotentials at @xmath48 , @xmath49 , their values are arbitrary , and the variation of @xmath50 produces a overall shift of the band structure . the value of @xmath50 is set to the work function , that can be obtained from experimental data@xcite or obtained from _ ab initio _ calculations . in nanostructures , as superlattices or quantum wells , @xmath50 determines the band alignment and its value has to be set carefully.@xcite in regions of small @xmath51 ( @xmath52 bohr@xmath47 ) , the changes in the pseudopotentials @xmath53 and @xmath43 do not substantially modify the band structure around the @xmath54 point . in all the cases , we have fitted the data points with cubic splines . in regions of large @xmath51 ( @xmath55 bohr@xmath47 ) , small oscillations persist , requiring a careful fitting .        in fig . [ bulkpot ] we plot the corresponding gallium and nitrogen semiempirical pseudopotentials . concerning the general behavior of the pseudopotentials , the oscillations persist until @xmath56  bohr@xmath47 , which determines the energy cut - off in the subsequent calculations ( with a energy cutoff of 30 ry ) . the last step is the solution of the kohn - sham equation  , where the kohn - sham potential @xmath1 is replaced by the semiempirical pseudopotentials calculated above , and the nonlocal term @xmath57 . this equation needs to be solved only once , since there is no need to achieve self - consistency . we proceed by comparing the band structures calculated using the sep method and the lda . as the spin - orbit interaction only produces minor splittings in gan , we have neglected this interaction for the overall band structure comparison , but it will be introduced later on , once the validity of the sep method has been demonstrated . the effects on the band structures caused by a reduction of the energy cut - off will also be analyzed . in fig . [ sepbulk ] we plot the following gan band structures : the lda results for @xmath58 ry ( black solid lines ) , the sep results with @xmath59 ry , ( green dashed lines ) and the band gap - corrected sep results ( red dots ) .    a comparisons of lda and sep calculations with the same energy cut - off shows a nearly perfect agreement for all the bands and throughout the brillouin zone and is not shown . the reduction of the enegy cut - off from 90 ry to 30 ry lead to a rigid displacement of the conduction bands towards lower energy , reducing the band gap approximately by 1 ev . the valence band states most affected by the reduction of @xmath60 are those far from the top of the valence band , at around @xmath61 ev . the topmost valence band states remain almost unaltered , and belong to the @xmath62 representation , in agreement with the _ ab initio _ band structure . we conclude that the lda and the sep band structures are in good agreement for the bands close to the band gap , i.e. , the most important part of the band structure in , e.g. , photoluminescence experiments .     and @xmath63 indicates the path in the brillouin zone . the lda band - gap is 2.7 ev , 1.6 for unshifted sep and 3.5 for shifted sep . , width=7 ]    another important issue in _ ab initio _ calculations of semiconductors is the underestimation of the band gap undermining the lda.@xcite such underestimation is enhanced , in this special case , if the energy cut - off is reduced . in the sep - method one can correct this problem easily by multiplying the local potentials by a gaussian function or renormalizing the nonlocal pseudopotentials.@xcite in this work we have chosen the second possibility where we act on the non - local angular momentum @xmath64-channel . figure  [ sepbulk ] shows the band structure of the band gap corrected seps . the band gap opening does not influence significantly the valence band , although a slight modification is appreciated in the vicinity of the @xmath54 point , which produces a change in the curvature , as expected from the decreased coupling with the valence bands . the conduction band effective mass for the corrected seps is now @xmath65 , in good agreement with the reported experimental value of @xmath66 ( being @xmath67 the free - electron mass).@xcite    for energy cut - offs lower than 30 ry , the changes in the band structure exceed the admissible margin . the wave functions of the top of the valence band and of the bottom of the conduction band have the same symmetry and shape than the lda wave functions . as noted above , the atoms at the free surfaces of a nanostructure have dangling bonds . the existence of these dangling bonds generates electronic states localized at the surface , with energies usually within the band gap of the semiconductor.@xcite these states can have energies close to the conduction and valence band edges , and can induce an artificial modification of the conduction and valence band states , and they have to be eliminated.@xcite in the literature two methods are used to overcome this difficulty . the first possibility is to achieve quantum confinement by embedding the structure into a virtual material with a high band gap @xcite . the main disadvantage of using this virtual material is the associated increase in the number of atoms . another way consists of attaching an atom , usually hydrogen , to the dangling bonds ( the so - called passivation ) , which displaces the energies of the surfaces states far from the band gap.@xcite however , the determination of the hydrogen pseudopotentials is not a trivial task , and it can depends on the material and of the surface orientation . the simplest way is to model the passivant pseudopotential by an analytical function , such as a gaussian . the method to optimize the form of the gaussian function can be a simple trial and error procedure , or a more sophisticated approach , such as a genetic algorithm , as developed in ref . however , the empirical parametrized gaussian pseudopotentials involve iterative calculations of the electronic structure , which complicates the computational task , and does not always guarantee to achieve an adequate passivation . in this work , we propose to derive the passivant pseudopotentials from lda calculations of a small free - standing structure . in the following we describe the necessary steps and illustrate them by an example for the @xmath68 surface of gan . _ step 1 : ab initio calculation of the self - consistent potential . _ the self - consistent potential is calculated for two free - standing @xmath68 gan layers . a bare layer ( unsaturated dangling bonds ) and a passivated layer ( a hydrogen attached to each dangling bond ) . specifically , the gan layer is approximately 15 angstrom thick ( 24 atoms ) , and the vacuum layer is about 10 angstrom on both sides to prevent coupling between layers . the input hydrogen norm - conserving potential is taken to be the same for the hydrogen atoms bound to the gallium and the nitrogen atom . in the case of materials with metallic surfaces ( e.g. cdse ) , a satisfactory passivation may require hydrogen potentials with fractional charges . for the atomic positions , we fix the nitrogen and gallium atoms to their bulk ideal positions and relax only the distance between the hydrogen atoms ( usually two atoms on each side of the slab ) and the corresponding passivated atom.@xcite    the band structures of the bare slab ( without geometry optimization ) and passivated @xmath69 $ ] slab ( with optimization of the distance of hydrogen atoms with respect to gan ) are shown in figs . [ abo_layer](a ) and ( b ) . note that the bands are rather flat along the @xmath70 direction , which confirms the decoupling between adjacent layers . the energy dispersion along the @xmath71 line of the brillouin zone is shown , exhibiting the same behavior as in the case of a quantum well.@xcite we observe that conduction and valence bands are very similar in both cases . however , the bare layer has a band gap larger by 0.37 ev . the main difference is obviously the presence of surface states , located within the band gap . these states are grouped in two energy ranges . the states @xmath72 ( index assigned in increasing order of energy ) , are closer to the valence band edge and have a flatter dispersion . the other group , formed by the states @xmath73 and @xmath74 , located in the middle of the band gap , have a more pronounced dispersion . on the contrary , the band structure of the passivated layer has no state within the band gap . this analysis is further clarified by examining the charge density in real space . figures  [ abo_layer](c ) and ( d ) show the plane - averaged density corresponding to some relevant states of the bare and passivated @xmath68-layers , together with the plane - averaged self - consistent potential profile , along the @xmath69 $ ] direction . in the case of the bare layer , the mid - gap surface states @xmath73 and @xmath74 present different profiles having a strong hybridization with the layer atoms , penetrating appreciably inside the layer . for the state @xmath75 the penetration into vacuum is significant for both , the passivate and the bare layers , which can be explained by the small electronic effective mass . in the case of the bare layer , the state @xmath75 develops a maximum at the surface of the layer , which is the signature of a hybridization with a surface state . on the other hand , the charge density of the valence band states seems to vanish at the surface without any sign of hybridization , in either the bare or the passivated layer . it is worth to mention that for larger passivated nanostructures the charge density at the border will be negligible . we can also conclude that the dangling bonds affect mainly the states of the conduction band . _ step 2 : determination of the screened hydrogen pseudopotential . _ we calculate the effective pseudopotential for the atom positions of the unpassivated slab using our semi - empirical pseudopotentials for ga and n , @xmath76 . we use the same fft grid as for the self - consistent potential for the passivated slab calculated in the previous step via lda , @xmath77 . the difference between the two sets of data @xmath78 represent the modification of the potential in response to the presence of the passivated surface . this includes the effects of the self - consistent charge redistribution around the surface ga and n atoms , as well as the additional hydrogen potential . the potential @xmath79 is therefore a potential that represents the entire modification of the total potential trough the introduction of a passivated surface . accordingly , we obtain a different passivant  hydrogen " potential for a hydrogen atom attached to a ga atom and for a hydrogen atom attached to a n atom . to extract a spherically averaged atomic quantity from @xmath79 we calculate the center of potential ( akin a center of mass ) @xmath80 of a set of @xmath81 points enclosed within a sphere around the hydrogen atomic positions @xmath82 : @xmath83 if the potential would be spherical and centered around @xmath82 then @xmath84 and @xmath82 would coincide . in our case , the geometric potential center is slightly shifted from the atomic position and we solve the equation self - consistently until @xmath85 . for the value of @xmath86 we use a value of 1.0 a.u ..    our passivant pseudo - hydrogen pseudopotential is given by @xmath87 and is shown in fig .  [ fig : scattered ] , where we plot the discrete set of fft grid points as black circles . the position of the hydrogen atom @xmath82 has been determined from eq . ( [ eq : rh ] ) . the smooth curves in fig .  [ fig : scattered ] represent our fits to the data points , where we used three cubic spline functions for @xmath88 . in this region , the data points lie mainly on a curve , which shows that the potential is close to spherically symmetric . for larger radii , the data points show significant scattering along the enegy axis , which is expected , since the potential at the surface will have some non - spherical character . we obtain good results when the data points are fitted by a yukawa potential : @xmath89    for the n ( ga ) passivant potential we determine @xmath90 = 1.864344 ( 2.012838 ) and @xmath91 = 3.7635 ( 4.060362 ) with @xmath92 = 0.80 a.u . for a potential in hartree units .    ) for a 15  thick gan film . the three different fits to the data points represent three different values of @xmath93 , @xmath94 and @xmath95 a. u. ( where the cubic splines are connected to the yukawa potential ) . we obtain the best results with @xmath92 = 0.80 a.u.,width=7 ]    we use a 1d fourier transformation to obtain the passivant potential in reciprocal space : @xmath96 with @xmath97    in figure  [ hydro ] we show the pseudo - hydrogen passivant pseudopotentials in reciprocal space for the potential attached to ga and to n. both curves are simple and smooth , and deviate from each other only for small values of @xmath98 , i.e. in their long range response , as expected . _ step 3 : consistency checks on wave functions . _ the accuracy of the passivation procedure has been checked by applying the sepm to the same layer used for the lda calculations . figure  [ sepm_layer ] shows the potential and the charge densities of one conduction @xmath75 and two valence band @xmath99 states . the calculation are done with an energy cut - off of 30 ry and a non - local part of the pseudopotential renormalized to the experimental bulk band gap . all spurious states are removed from the gap and the states are well localized inside the structure . the states @xmath100 and @xmath101 have the same profile than their counterparts lda states . the small discrepancies can be attributed to the small dimensions of the layer , and gradually disappear for larger nanostructures.@xcite     and two valence band states , @xmath99 of a passivated gan slab calculated with the sepm . lower panel : atomic semiempirical pseudopotential across the slab . the quantities are averaged in successive planes perpendicular to the surface.,width=283 ] to illustrate the performance of the sepm in the study of nitride nanostructures , we present here calculations of the electronic structure of free - standing gan nws . the nws are assumed to be oriented along the @xmath102-axis and infinite in length . we assume a hexagonal cross section , with @xmath68 lateral faces . this is the most usual shape adopted by gan nws grown by plasma - assisted molecular beam epitaxy technique.@xcite this nw faceting has also been supported by theoretical calculations,@xcite showing a lower formation energy for @xmath68 than @xmath103 surfaces . the electronic states are obtained by solving @xmath104 \psi_{k , n}(\bm{r } ) = \varepsilon_n(k)\psi_{k , n}(\bm{r } ) \quad , \ ] ] where the wire potential @xmath105 is constructed as a superposition of atomic pseudopotentials .    in real free - standing nws a surface reconstruction is expected that would minimizes the total energy of the system . however , in our calculations we will assume for simplicity that the atomic arrangement corresponds everywhere to the perfect wurtzite crystal structure . when building the potential @xmath105 we have passivated the facets using the potential described in section [ surface ] , which removes states from the gap in all the cases investigated . equation ( [ vnw ] ) is solved by imposing artificial periodic boundary conditions on the wire surrounded by @xmath106 layers of vacuum . this transforms eq . ( [ vnw ] ) into a bloch - periodic band - structure problem with a large simulation cell , solved here by expanding @xmath107 in plane waves . we use a sufficiently large wire - wire separation @xmath106 , so that the solutions become independent of @xmath106 . a sketch of the modeled nw is shown in fig . [ sketch ] . in the following , the term _ size _ refers to the largest lateral dimension , as shown in in fig . [ sketch ] and is labeled as @xmath108 . we have considered a range of sizes @xmath108 from 1.5 to 6.5 nm , which corresponds to a number of atoms varying between 100 to 1500 . the matrix diagonalization corresponding to the problem eq . ( [ vnw ] ) can be solved using the folded spectrum method.@xcite the computational times , within the sepm , are significantly shorter than they would be in a corresponding self - consistent _ ab initio _ calculation . $ ] direction . the lateral surfaces of the nw are @xmath109 , and the corresponding surface atoms are passivated by pseudo - hydrogen atoms.,width=7 ]    here , we will focus on the analysis of the zone center ( @xmath110 ) nw electronic states . these states are naturally classified into valence band ( vb , @xmath111 ) and conduction band ( cb , @xmath102 ) nw states , and in each case they will be indexed in ascending order , as @xmath112 and @xmath113 , according to their separation in energy from the respective @xmath114 bulk band edges @xmath115 and @xmath116 . any of the considered states exhibits kramer s double degeneracy . when analyzing the corresponding charge densities ( here meaning , i.e. , the ( vertically - averaged ) squared wave functions ) , an envelope pattern is recognized in the cross - section , which corresponds very closely to the radial pattern of the bessel functions @xmath117 . we find convenient to incorporate this symmetry information to the labeling of the states by adding the notation @xmath64 ( meaning : _ approximately the same symmetry as @xmath118 _ ) , @xmath119 ( _ approximately the same symmetry as @xmath120 _ ) , etc ...    in order to further characterize the symmetry of the nw states @xmath121 , we can calculate their projections on the bulk band states , as explained in ref .  . in particular , in the analysis presented below , we have considered the projections @xmath122 where @xmath123 are the near - gap @xmath54 ( @xmath124 ) bulk states labeled schematically by their symmetries : @xmath125 ( bottom of the cb ) and @xmath126 , @xmath127 , @xmath128.@xcite    ) in the text , whereas the red dots represents the confinement energy of the lowest cb nw . the solid line refers to the right scale and it represents the projection of the @xmath75 wave function onto the cb edge bulk state . the charge density of the cb state @xmath129 for @xmath130 nm is shown as an inset . ] figure  [ cond ] shows the confinement energies ( i.e. , the energies referred to the bulk cb edge energy ) of the cb nw states as a function of the nw size . it is also shows the projection of the lowest energy state @xmath75 onto the cb bulk state @xmath131 , revealing that @xmath75 has between 65 % and 97 % conduction band character . additionally , the charge density of @xmath75 exhibits a clear @xmath64-type envelope for all the sizes explored here , in agreement with expectations from the ema . on the other hand , we have fitted the energy of the @xmath129 state to the function @xmath132 and have found the values @xmath133 and @xmath134 . the fitted value of the parameter @xmath135 differs significantly from the prediction of the single - band ema ( @xmath136 ) , which is represented by a dotted line in fig . we conclude that , even for the case of the almost isolated conduction band , the confinement imposed by the nanowire geometry , while leading to an overall symmetry of the cb states similar to that predicted by the ema , leads to a remarkably different energy dispersion as a function of nw size . , is shown as an inset . the solid line connects the size dispersion of the hybrid @xmath137-@xmath138 state with @xmath119-type envelope symmetry , which is discussed in the main text . , width=283 ]        we now turn to the analysis of the vb electronic structure of the nw . figure  [ vale ] shows an overall view of the confinement effects in the vb by displaying the highest energy levels for every nw size . in order to give further insight into the complex evolution of the vb states with sizes , we have represented in fig . [ evo ] the charge densities of the first four vb states for various nw sizes . a first examination reveals that the wave functions are well confined inside the nws , and only for some states in the smallest nw show a slight state localization on surface atoms . in particular , the analysis of the charge density of the highest vb state @xmath100 shows that its envelope symmetry changes suddenly from @xmath119-type for the largest nw to @xmath64-type for smaller nws . for the size 2.6 nm and below , the first three vb states show a @xmath64-type envelope , whereas the vb state with @xmath119-type envelope is moved to the fourth place in the energy spectrum . .projections ( @xmath139 , @xmath140 , @xmath141 ) of the upper vb nw states @xmath142 ( @xmath143 ) , onto the bulk vb states at @xmath54 , ( @xmath137 , @xmath138 , @xmath144 ) . [ cols="<,^,^,^,^",options="header " , ]     the above study must be complemented with the analysis of the projections of the nw vb states onto the bulk states @xmath137 , @xmath145 and @xmath146 , at the @xmath54 point , which will give us additional information about the symmetry of the states . in table . [ tabla ] we have summarized the values of the corresponding projections ( @xmath139 , @xmath140 and @xmath141 ) , for the same set of states as in fig . [ evo ] . for the largest nanowire shown here , the composition of the topmost state , which has envelope @xmath119 symmetry , is almost equally divided between the bulk states @xmath137 and @xmath138 . the next two vb states , @xmath101 and @xmath147 , of @xmath64-symmetry , clearly show a dominant contribution from bulk bands @xmath137 ( @xmath148 ) and @xmath138 ( @xmath149 ) , respectively . neither of these states has a significant contribution from the @xmath144-band . when the nanowire size is reduced to 3.7 nm , we observe that the two highest vb states have now dominant contributions from @xmath137- and @xmath138-bands , with @xmath64 symmetry , whereas the mixed @xmath137-@xmath138 state becomes now the third state , @xmath150 . if one further reduces the nw size to 2.6 nm and below , the mixed @xmath137-@xmath138 state now becomes @xmath151 , whereas the first three states show a high degree of @xmath137-@xmath138-@xmath144 mixing . thus , our detailed analysis of the nw wave functions as a function of size has allowed us to identify a state with @xmath119-type envelope symmetry which is essentially an equal - weight mixture of bulk bands @xmath137 and @xmath138 . in fig . [ vale ] we have traced a line connecting these states . the change in the symmetry of the highest vb state as a function of size is just one of the consequences of the nontrivial interplay between symmetry mixing , spin - orbit coupling and confinement effects on the nw valence band electronic structure . this is in contrast to the simpler situation of the cb electronic structure . this unconventional trend of the nanowire electronic structure has also been reported by calculations using the tight - binding method.@xcite moreover , studies in other wurtzite systems , as the one made in zno nanocrystals by s. baskoutas and g. bester @xcite shows also that the highest state of the valence band is characterized by a @xmath119-type envelope . we have followed the idea @xcite to use the screened effective potential from a dft bulk calculation to extract a spherically averaged atomic pseudopotential . we have derived the analytic connection between the screened dft effective potential and the atomic quantities in the case of wurtzite bulk unit cells . this allowed us to derive atomic semi - empirical pseudopotentials for nitride materials . the pseudopotentials are non - local , in the same way as dft norm - conserving pseudopotentials , and reproduce the dft results for the eigenvalues of bulk up to a few mev . our potentials are not aimed at the calculation of total energies but at the calculations of the eigenvalues and eigenfunctions around the band gap , which gives some leverage on the reduction of the energy cut - off , compared to dft calculations . the advantage of the method compared to the original empirical pseudopotential approach ( see ref . and references therein ) resides in the direct link to ab - intio calculations and the ensuing rather simple fitting procedure of a curve through a dense set of data points . the transferability of the potentials are good for the situations studied here but are expected to degrade when charge transfer significantly deviates from the bulk situation . this represents a general limitation of the method that we , however , palliate for the case of a free surface by deriving effective passivant potentials ( see next paragraph ) . the defective band gap inherited from dft is correct by a slight adjustment of the non - local channel of the potential which has only a very marginal effect on the wave functions , as shown previously @xcite . the derived method should be applicable to other materials . eventually , the method can be built on the basis of gw calculations . however , full self - consistent gw calculations are still prohibitive in the applications to extended systems and only work for small molecules.@xcite    in the second part of this work , we have established a simple and direct methodology to obtain the passivating pseudo - hydrogen pseudopotentials , essential for the treatment of free surfaces in nanowires or nanoparticles . we thereby use the difference between the screened effective dft potentials of a passivated structure and the semiempirical pseudopotential of the corresponding unpassivated structure . the derivation is done in real space and we find that a cubic spline interpolation up to a cut - off radius @xmath92 , connected to a yukawa potential leads to satisfactory results .    the suitability of the semi - empirical pseudopotential method and of our passivation procedure has been explored in gan nanowires . due to the high confinement of nw sizes below 10 nm , the use of accurate atomic potentials is unavoidable , as continuous methods are expected to fail in the description of such systems . we have found that , while the conduction band states follow a behavior compatible with the ema predictions , the valence band states exhibit a complex evolution as a function of the nw size . we find drastic variations in the character of the wave functions with the nw size , which will influence optical properties such as the polarization of the photoluminescence . a. molina - snchez thanks the deutscher akademischer austausch dienst ( daad ) for funding the research stay at the max - planck institute fr festkrperforschung in stuttgart ( germany ) , and the group of gabriel bester for its kind hospitality . all dft calculations are done with abinit . x. gonze , b. amadon , p .- anglade , j .- beuken , f. bottin , p. boulanger , f. bruneval , d. caliste , r. caracas , m. cote , t. deutsch , l. genovese , ph . ghosez , m. giantomassi , s. goedecker , d.r . hamann , p. hermet , f. jollet , g. jomard , s. leroux , m. mancini , s. mazevet , m.j.t . oliveira , g. onida , y. pouillon , t. rangel , g .- rignanese , d. sangalli , r. shaltaf , m. torrent , m.j . verstraete , g. zerah , j.w . zwanziger , computer phys . commun . * 180 * , 2582 ( 2009 ) .
Abstract:

Answer: we present a semiempirical pseudopotential method based on screened atomic pseudopotentials and derived from _ ab initio _ calculations . this approach is motivated by the demand for pseudopotentials able to address nanostructures , where _ ab initio _ methods are both too costly and insufficiently accurate at the level of the local - density approximation , while mesoscopic effective - mass approaches are inapplicable due to the small size of the structures along , at least , one dimension . in this work we improve the traditional pseudopotential method by a two - step process : first , we invert a set of self - consistently determined screened _ ab initio _ potentials in wurtzite gan for a range of unit cell volumes , thus determining spherically - symmetric and structurally averaged atomic potentials . second , we adjust the potentials to reproduce observed excitation energies . we find that the adjustment represents a reasonably small perturbation over the potential , so that the ensuing potential still reproduces the original wave functions , while the excitation energies are significantly improved . we furthermore deal with the passivation of the dangling bonds of free surfaces which is relevant for the study of nanowires and colloidal nanoparticles . we present a methodology to derive passivant pseudopotentials from _ ab initio _ calculations . we apply our pseudopotential approach to the exploration of the confinement effects on the electronic structure of gan nanowires .