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: understanding strong correlation effect away from equilibrium has been one of the most interesting yet challenging problems in condensed matter physics . a prominent realization of such phenomena is embodied in quantum transport through a nanostructure under finite bias voltage . to understand the interplay of the correlation effect and nonequilibrium nature , nonequilibrium version of the single impurity anderson model ( siam ) and its extensions have been serving and continue to do so as a central theoretical model . the siam is indeed considered to be one of the best studied strongly correlated model , and despite its apparent simplicity , it exhibits rich physics already in equilibrium , such as the coulomb blockade and the kondo physics that have been observed in experiments . equilibrium properties of the siam have been well understood thanks to concerted efforts of several theoretical approaches over years  by perturbative treatment , fermi liquid description , as well as exact results by the bethe ansatz method and numerical renormalization group ( nrg ) calculations  ( see for instance , . ) in contrast , the situation of the nonequilibrium siam is not so satisfactory . each of the above approaches has met some difficulty in treating nonequilibrium phenomena . a theoretical approach that can deal with the strong correlation effect in nonequilibrium is still called for .    notwithstanding , a number of analytical and numerical methods have been devised to investigate nonequilibrium stationary phenomena : nonequilibrium perturbation approaches  @xcite and its modifications  @xcite , the non - crossing approximation  @xcite , the functional renormalization group treatment  @xcite , quantum monte carlo calculations on the keldysh contour  @xcite , the iterative real - time path integral method  @xcite and so on . unfortunately those approaches fail to give a consistent picture concerning the finite bias effect on the dot spectral function , particularly regarding a possible splitting of the kondo resonance .    as for the equilibrium siam , the second - order perturbation regarding the coulomb interaction @xmath0 on the dot  @xcite is known to capture essential features of kondo physics and agrees qualitatively well with exact results obtained by the bethe ansatz and nrg approaches  @xcite . such good agreement seems to persist in nonequilibrium stationary state at finite bias voltage . for the two - terminal particle - hole ( ph ) symmetric siam , a recent study by @xcite showed that the nonequilibrium second order perturbation calculation of the spectral function agrees with that calculated by the diagrammatic quantum monte carlo simulation , excellently up to interaction strength @xmath1 ( where @xmath2 is the total relaxation rate due to leads ) ; pretty well even for @xmath3 at bias voltage @xmath4 . a typical magnitude of @xmath5 of a semiconductor quantum dot is roughly @xmath6 depending on the size and the configuration of the dot . therefore , there is a good chance of describing a realistic system within the validity of nonequilibrium perturbation approach . the great advantage of semiconductor dot systems is to allow us to control several physical parameters . those include changing gate voltage as well as configuring a more involved structure such as a multi - terminal dot  @xcite or an interferometer embedding a quantum dot . theoretical treatments often limit themselves to a system with the ph symmetry where the dot occupation number is fixed to be one half per spin . although assuming the ph symmetry makes sense and comes in handy in extracting the essence of the kondo resonance , we should bear in mind that such symmetry is not intrinsic and can be broken easily in realistic systems , by gate voltage , asymmetry of the coupling with the leads , or asymmetric drops of bias voltage  @xcite . the ph asymmetry commonly appears in a multi - terminal dot or in an interferometer embedding a quantum dot . it is also argued that the effect of the ph asymmetry might be responsible for the deviation observed in nonequilibrium transport experiments from the `` universal '' behavior of the ph symmetric siam  @xcite . to work on realistic systems , it is imperative to understand how the ph asymmetry affects nonequilibrium transport .    in this paper , we examine the second - order nonequilibrium perturbation regarding the coulomb interaction @xmath0 of the multi - terminal siam . the ph symmetry is not assumed , and miscellaneous types of asymmetry of couplings to the leads and/or voltage drops are incorporated as a generic multi - terminal configuration . our main focus is to provide solid analytical results of the behavior of the nonequilibrium self - energy and hence the dot spectral function for the full range of frequency and bias voltage , within the validity of the second - order perturbation theory of @xmath0 . the result encompasses fermi liquid behavior as well as incoherent non - fermi liquid contribution , showing analytically that increasing finite bias voltage leads to a crossover from the kondo resonance to the coulomb blockade behaviors . the present work contrasts preceding perturbative studies  @xcite whose evaluations relied on either numerical means or a small - parameter expansion of bias voltage and frequency . the only notable exception , to the author s knowledge , is a recent work by @xcite , which succeeded in evaluating analytically the second - order self - energy for the two - terminal ph symmetric dot . intending to apply such analysis to a wider range of realistic systems and examine the effect that the two - terminal ph symmetric siam can not capture , we extend their approach to a generic multi - terminal dot where the ph symmetry may not necessarily be present . an embarrassing drawback of using the nonequilibrium perturbation theory is that when one has it naively apply to the ph asymmetric siam , it may disrespect the preservation of the steady current . @xcite as a result , one then needs some current - preserving prescription , and different self - consistent schemes have been proposed and adopted .  @xcite as will be seen , the current preserving condition involves all the frequency range , not only of low frequency region that validates fermi liquid description but also of the incoherent non - fermi liquid part ( see eq .   below ) . therefore , an approximation based on the low energy physics , particularly the fermi liquid picture , should be used with care . the self - energy we will construct analytically is checked to satisfy the spectral sum rule at finite bias voltage , so that we regard it as giving a consistent description for the full range of frequency in nonequilibrium . by taking its advantage , we also demonstrate a self - consistent , current - preserving calculation of the nonequilibrium spectral function for a system where the ph symmetry is not present . the paper is organized as follows . in sec . [ sec : model ] , we introduce the multi - terminal siam in nonequilibrium . we review briefly how to obtain the exact current formula by clarifying the role of the current conservation at finite bias voltage . section  [ sec : self - energy ] presents analytical expression of the retarded self - energy for a general multi - terminal dot up to the second - order of the interaction strength . subsequently in sec . [ sec : discussion ] , we examine and discuss its various analytical behaviors including high- and low - bias voltage limits . sec . [ sec : spectral - fn ] is devoted to constructing a non - equilibrium spectral function using the self - energy obtained in the previous section . we focus our attention on two particular situations : ( 1 )  self - consistent , current - preserving evaluation of the nonequilibrium spectral function for the two - terminal ph asymmetric siam , and ( 2 )  multi - terminal effect of the ph symmetric siam . finally we conclude in sec . [ sec : conclusion ] . mathematical details leading to our main analytical result eq .   as well as other necessary material regarding dilogarithm is summarized in appendices . the model we consider is the single impurity anderson model connecting with multiple leads @xmath7 whose chemical potentials are sustained by @xmath8 . the total hamiltonian of the system consists of @xmath9 , where @xmath10 , @xmath11 and @xmath12 represent the dot hamiltonian with the coulomb interaction , the hopping term between the dot and the leads , and the hamiltonian of a noninteracting lead @xmath13 , respectively . they are specified by @xmath14 where @xmath15 is the dot electron number operator with spin @xmath16 and @xmath17 are electron operators at the lead @xmath13 . in the following , we consider the spin - independent transport case , but an extension to the spin - dependence situation such as in the presence of magnetic field or ferromagnetic leads is straightforward . when applying the wide - band limit , all the effect of the lead @xmath13 is encoded in terms of its chemical potential @xmath8 and relaxation rate @xmath18 , where @xmath19 is the density of states of the lead @xmath13 . the dot level @xmath20 controls the average occupation number on the dot ; it corresponds roughly to @xmath21 , @xmath22 , @xmath23 for @xmath24 , @xmath25 , and @xmath26 , respectively . the ph symmetry is realized when @xmath27 and @xmath28 ( see eqs .   and below ) . we here briefly summarize how the current through the dot is determined in a multi - terminal setting . special attention is paid to the role of the current conservation because it has been known that nonequilibrium perturbation calculation does not respect it in general . @xcite we illustrate how to ensure the current conservation by requiring minimally . the argument below is valid regardless of a specific approximation scheme chosen , whether nonequilibrium perturbation or any other approach . following the standard protocol of the keldysh formulation  @xcite , we start with writing the current @xmath29 flowing from the lead @xmath13 to the dot in terms of the dot s lesser green function @xmath30 and the retarded one @xmath31 : @xmath32,\end{aligned}\ ] ] where @xmath33 is the fermi distribution function at the lead @xmath13 . as the present model preserves the total spin as well as charge , the net spin current flowing to the dot should vanish in the steady state , which imposes the integral relation between @xmath30 and @xmath31 , @xmath34 = 0 . \label{eq : gless - gr}\end{aligned}\ ] ] here we have introduced the total relaxation rates @xmath35 and the effective fermi distribution @xmath36 weighted by the leads , @xmath37 when we ignore the energy dependence of the relaxation rates @xmath38 , we can recast eq .   into a more familiar form @xmath39 because @xmath40 is the definition of the exact dot occupation number . note the quantity @xmath41 is nothing but the exact dot spectral function out of equilibrium . we emphasize that eq . ( or equivalently eq .  ) is the minimum , exact requirement that ensures the current preservation . it constrains the exact @xmath42 and @xmath43 that depend on the interaction as well as bias voltage in a nontrivial way . one can accordingly eliminate @xmath44 in @xmath29 , to reach the landauer - buttiker type current formula at the lead @xmath13 , @xmath45 or , the current conservation allows us to write it as @xmath46,\end{aligned}\ ] ] where @xmath47 is the exact number of states with spin @xmath16 at finite temperature in general , defined by @xmath48 it tells us that differential conductance @xmath49 with fixing all other @xmath50 s is proportional to the nonequilibrium dot spectral function , provided changing @xmath8 does not affect the occupation number  @xcite . such situation is realized , for instance , when a probe lead couples weakly to the dot . the case of a noninteracting dot always satisfies the current preserving condition eq .   as @xmath51 holds for any @xmath52 ; the distribution function of dot electrons @xmath53 is equal to @xmath54 . this is not the case for an interacting dot , however . as for the interacting case , not so much can be said . we only see the special case with the two - terminal ph symmetric dot satisfy eq .   by choosing @xmath55 irrespective of interaction strength . except for this ph symmetric case , a general connection between @xmath42 and @xmath43 is not known so far . it is remarked that , based on the quasiparticle picture , a noninteracting relation @xmath56 is sometimes used to deduce an approximate form of @xmath42 out of @xmath43 for an interacting dot . such approximation is called the ng s ansatz  @xcite . although it might be be simple and handy , its validity is far from clear . we will not rely on such additional approximation below . it is also important to distinguish in eq . the electron occupation number @xmath57 from the quasiparticle occupation number @xmath58 , as the two quantities are different at finite bias voltage since the luttinger relation holds only in equilibrium  @xcite . contribution to the dot occupation number comes from all range of frequency , including the incoherent part . one sees fulfilling the spectral weight sum rule @xmath59 crucial to have the dot occupation number @xmath57 normalized correctly . in general , one needs to determine @xmath57 appropriately to satisfy eq .   as a function of interaction and chemical potentials of the leads . the applicability of quasi - particle approaches that ignores the incoherent part is unclear . in this section , we evaluate analytically the nonequilibrium retarded self - energy up to the second order of interaction strength @xmath0 for the multi - terminal siam . we first examine the contribution at the first - order and the role of current preservation . then we present the analytical result of the second - order self - energy in terms of dilogarithm .    following the standard treatment of the keldysh formulation ,  @xcite the nonequilibrium green function and the self - energy take a matrix structure @xmath60 satisfying symmetry relations @xmath61 and @xmath62 . the retarded green function is defined by @xmath63 ; the retarded self - energy , by @xmath64 .    to proceed evaluation , it is convenient to classify self - energy diagrams into the two types : the hartree - type diagram ( fig . [ fig : hartree - type ] ) that can be disconnected by cutting a single interaction line , and the rest which we call the correlation part and reassign the symbol @xmath65 to . the latter starts at the second - order . the resulting green function ( matrix ) takes a form of @xmath66^{-1},\end{aligned}\ ] ] where @xmath67 represents a pauli matrix of the keldysh structure , and @xmath68 refers to the exact occupation number of the dot electron with the opposite spin . accordingly , the corresponding retarded green function becomes @xmath69 where @xmath70 is the hartree level of the dot .    . the double line refers to the exact green function . ]      before starting evaluating the correlation part @xmath71 that starts contributing at the second order , it is worthwhile to examine the current preserving condition eq .   at the first order . at this order , it reduces to the self - consistent hartree - fock equation for the dot occupation number @xmath72 : @xmath73 . \label{eq : hf - self - consistency}\end{aligned}\ ] ] it shows how the two - terminal ph symmetric siam is special by choosing @xmath74 , @xmath75 and @xmath76 ; the second term of the right - hand side vanishes by having the solution @xmath77 even at finite bias voltage . it also indicates that the current preservation necessarily has the occupation number depend on asymmetry of the lead couplings as well as interaction strength for the ph asymmetric siam . indeed , for a small deviation from the ph symmetry and bias voltage , we see the hartree - fock occupation number behave as @xmath78 where @xmath79 is the average chemical potential weighted by leads , @xmath80 note @xmath79 vanishes when no bias voltage applies , as we incorporate the overall net offset by leads into @xmath20 .              following the standard perturbation treatment of the keldysh formulation , we see there is only one diagram contributing to @xmath81 at the second order ( fig .  [ fig : sigma2 ] ) after eliminating the hartree - type contribution . the contribution is written as @xmath82 where @xmath83 refer to to the unperturbed green functions ( including the hartree term ) , whose concrete expressions are found in appendix a. the polarization matrix @xmath84 is defined by @xmath85 ( fig . [ fig : pi ] ) .. ]    as was shown by the current formula in the previous section , we need only the dot spectral function to study quantum transport , hence @xmath71 suffices . therefore it is more advantageous to work on the rak - representation , where the polarization parts become    @xmath86 , \\      & \pi^{a}_{12 } = \frac{i\hbar}{2 } \left [ g^{a}_{12 } g^{k}_{21 } +        g^{k}_{12 } g^{r}_{21 }      \right ] , \\      & \pi^{k}_{12 } = \frac{i\hbar}{2 } \left [ g^{k}_{12 } g^{k}_{21 } +        g^{r}_{12 } g^{a}_{21 } + g^{a}_{12 } g^{r}_{21 } \right ] ,    \end{aligned}\ ] ]    [ eq : def - pi ]    and their fourier transformations are given in appendix b. accordingly , we can express the retarded self - energy @xmath71 as @xmath87,\end{aligned}\ ] ] where @xmath88 the above second - order expression of @xmath71 is standard , but it has so far been mainly used for numerical evaluation , quite often restricted for the two - terminal ph symmetric siam . we intend to evaluate eqs . ( [ eq : i1],[eq : i2 ] ) analytically for the generic multi - terminal siam , along the line employed in ref .  . delegating all the mathematical details to appendices c and d , we summarize our result of the analytical evaluation of @xmath71 as follows :    @xmath89 , \label{eq : sigmar}\end{aligned}\ ] ]    here , functions @xmath90 ( @xmath91 ) are found to be ( using @xmath92 ) ,    [ eq : xi123 ] @xmath93 \nonumber \\ & \qquad \qquad   + \sum_{a , b,\beta }        \tfrac{4\gamma_{a }          \gamma_{b}}{\gamma^{2 } } \big [ { \mathop{\mathrm{li}_{2}}\nolimits}(\tfrac{-\varepsilon+ \beta \zeta_{a\bar{\sigma } } } { \beta          \zeta_{b \bar{\sigma } } + i\gamma } ) + \tfrac{1}{4 } { \mathop{\mathrm{log}}\nolimits}^{2 }        ( \tfrac{-\zeta_{a\bar{\sigma}}+i\gamma}{\zeta_{b\bar{\sigma } } +          i\gamma } )        \big ] , \\ \xi_{2}(\varepsilon ) & =   6\pi^{2 }        - \sum_{a , b,\beta } \tfrac{4\gamma_{a } \gamma_{b}}{\gamma^{2 } } \big [        \lambda(\tfrac{\varepsilon -\zeta_{a\sigma } +          2i\gamma}{\beta \zeta_{b\bar{\sigma } } + i\gamma } ) +        \lambda(\tfrac{\varepsilon - \beta \zeta_{b\bar{\sigma } }          + 2i\gamma}{\zeta_{a\sigma } + i\gamma } ) + \tfrac{1}{2 } { \mathop{\mathrm{log}}\nolimits}^{2 }        ( \tfrac{\zeta_{a\sigma } + i\gamma}{\beta \zeta_{b \bar{\sigma } } +          i\gamma } ) \big ] \nonumber \\ & \qquad   \qquad - \sum_{a , b,\beta }        \tfrac{4\gamma_{a } \gamma_{b}}{\gamma^{2 } } \big [ \lambda        ( \tfrac{\varepsilon+ \beta \zeta_{a\bar{\sigma } } +          2i\gamma } { \beta \zeta_{b\bar{\sigma } } + i\gamma } ) +        \tfrac{1}{4 } { \mathop{\mathrm{log}}\nolimits}^{2 } ( \tfrac{-\zeta_{a\bar{\sigma } } +          i\gamma}{\zeta_{b\bar{\sigma } } + i\gamma } )        \big ] , \\ \xi_{3 } & = \left[\sum_{a }        \tfrac{2\gamma_{a}}{\gamma }        { \mathop{\mathrm{log}}\nolimits}\left ( \tfrac{-\zeta_{a\bar{\sigma } } +            i\gamma}{\zeta_{a\bar{\sigma } } + i\gamma}\right ) \right]^{2 } ,      \end{aligned}\ ] ] [ eq : xi-123 ]    where the summations over @xmath94 as well as terminals @xmath95 are understood . function @xmath96 is dilogarithm , whose definition as well as various useful properties are summarized in appendix c ; @xmath97 is defined by   is equivalent to that given in ref .  , but we prefer writing it in this form because its analyticity is more transparent . ] @xmath98 { \mathop{\mathrm{log}}\nolimits}z . \label{def : lambda}\end{aligned}\ ] ] analytical formula @xmath71 given in eqs .   and constitutes the main result of the present paper . consequently , the nonequilibrium spectral function of the multi - terminal siam is given analytically for full range of frequency and bias voltage , once one chooses @xmath57 to satisfy eq .  . the result also applies to a more involved structured system , such as an interferometer embedding a quantum dot , by simply replacing @xmath20 and @xmath38 to take account of those geometric effect . having obtained an explicit analytical form of the second - order self - energy @xmath99 at arbitrary frequency and bias voltage , we now examine its various limiting behaviors . most of those limiting behaviors have been known for the two - terminal ph symmetric siam , so it s assuring to reproduce those expressions in such a case . simultaneously , our results below provide multi - terminal , ph asymmetric extensions of those asymptotic results . we can reproduce the equilibrium result by setting all the chemical potentials equal , @xmath100 . then we immediately see @xmath101 and @xmath102 , \\ & \xi_{2 } = 8\cdot \big[\frac{3\pi^{2}}{4 } - 3\lambda(\tfrac{\varepsilon_{\sigma } + 2i\gamma}{i\gamma } )    \big],\end{aligned}\ ] ] where @xmath103 . as a result , the correlation part of the self - energy in equilibrium becomes @xmath104 . \label{eq : sigma - equilibrium}\end{aligned}\ ] ] the ph symmetric case in particular corresponds to @xmath105 . it reproduces the perturbation results by yamada and yoshida  @xcite up to the second - order of @xmath0 , when we expand the above for small @xmath52 . the ph symmetric result is indeed identical with the one obtained in ref .   for arbitrary frequency ( see also eq .   below ) . @xcite have evaluated analytically the self - energy and the spectral function for the two - terminal ph symmetric siam . in our notation , it corresponds to the case @xmath106 , and @xmath107 . when we parametrize the two chemical potentials by @xmath108 with @xmath109 in eqs .  , the self - energy can be written as @xmath110 , \label{eq : ph - symmetric - sigmar}\end{aligned}\ ] ] where @xmath111 , \\ & \xi_{2 } =   6\pi^{2 }   - 6\sum_{a , b } \big [      \lambda(\tfrac{\omega -a ev/2 + 2i\gamma}{b ev/2 +   i\gamma } )   + \frac{1}{4 } { \mathop{\mathrm{log}}\nolimits}^{2 } ( \tfrac{a ev/2 + i\gamma}{b ev/2 + i\gamma } ) \big].\end{aligned}\ ] ] the above result are identical with what ref .   obtained . we now employ the small - parameter expansion of @xmath71 around the half - filled equilibrium system . here we assume parameters @xmath112 and @xmath113 are much smaller than the total relaxation rate @xmath2 . the expansion of @xmath114 is found to contain the first- and second - order terms regarding @xmath115 and @xmath116 , while @xmath117 do only the second - order terms . therefore the result of the expansion up to the second - order of these small parameters is presented as @xmath118.\end{aligned}\ ] ] functions @xmath90 can be expanded straightforwardly by using the taylor expansion of dilogarithm in appendix c. they are found to behave @xmath119 , \\ & \xi_{2}(\varepsilon ) \approx   - 8 \ ,   \big [   \frac{3 ( -\varepsilon^{2 } + 2    \varepsilon\bar{\mu } + 2    \bar{\mu}^{2 } -     \overline{\mu^{2 } } - 2    \overline{\mu^{2}})}{4(i\gamma)^{2 } }    \big ] \\ & \xi_{3 } \approx \frac{16 \bar{\mu}^{2}}{(i\gamma)^{2}},\end{aligned}\ ] ] where @xmath79 is defined in eq . and we have introduced @xmath120 combining all of these , we reach the small bias / frequency behavior of the self - energy @xmath81 as @xmath121 \nonumber \\ & \qquad \qquad -   \frac{iu^{2}}{2\pi^{2}\gamma^{3 } } \big [   ( \omega - \bar{\mu})^{2 } + 3\ , ( \delta \mu)^{2 } \big ] . \label{eq : small - parameter - expansion}\end{aligned}\ ] ]    small bias expansion of @xmath122 for the two - terminal system has been discussed and determined by the argument based on the ward identity . @xcite the dependence appearing in eq .   fully conforms to it ( except for the presence of the bare interaction instead of the renormalized one ) . indeed , correspondence is made clear by noting the parameters @xmath79 and @xmath123 for the two - terminal case become @xmath124 the presence of linear term in @xmath52 and @xmath125 for the two - terminal ph asymmetric siam was also emphasized recently . @xcite      one expects naively that the limit of large bias voltage @xmath126 corresponds to the high temperature limit @xmath127 in equilibrium ; it was shown to be so for the two - terminal ph symmetric siam . @xcite we now show that the same applies to the multi - terminal siam where bias voltages of the leads are pairwisely large , _ i.e. _ , half of them are positively large , and the others are negatively large .    in the large bias voltage limit , all the arguments of dilogarithm functions appearing in eqs . become @xmath128 , where the values of dilogarithm are known ( see appendix c ) . accordingly , the pairwisely large bias limit of @xmath90 is found to be @xmath129 correspondingly , the retarded self - energy becomes @xmath130 it shows that the result of the multi - terminal siam is the same with that of the two - terminal ph symmetric siam except for a frequency shift . accordingly , the retarded green function @xmath131 in this limit is given by @xmath132 the form indicates that for sufficiently strong interaction @xmath133 , the dot spectral function has two peaks at @xmath134 with broadening @xmath135 , so the system is driven into the the coulomb blockade regime . on the other hand , for weak interaction @xmath136 , it has only one peak with two different values of broadening that reduce to @xmath2 and @xmath137 in the @xmath138 limit .    what is the role of the current preservation condition eq .   in this limit ? it just determines the dot occupation number explicitly . in fact , the condition becomes @xmath139 with @xmath140 , and @xmath57 is independent of the interaction strength because bias voltage sets the largest scale . one can evaluate the above integral exactly to have @xmath141 where @xmath142 . as a result , expanding it up to the second order of @xmath143 leads to @xmath144 the first term corresponds to the occupation number that one expects naturally from the effective distribution @xmath36 ; it corresponds , for instance , @xmath145 for the two - terminal dot with @xmath146 with @xmath147 . the second term is the deviation from it , which is proportional to the average of the inverse chemical potential weighted by the leads . we now turn our attention to the behavior of the nonequilibrium dot spectral function , using our analytical expression of the self - energy eqs . ( [ eq : sigmar],[eq : xi123 ] ) . below we particularly focus our attention on the two cases : the two - terminal ph asymmetric siam where current preservation has been an issue , and the multi - terminal ph symmetric siam where the role of multiple leads has been raising questions . in all of the calculations below , we have checked numerically the validity of the spectral weight sum rule at each configuration of bias voltages .      as was emphasized in sec . [ sec : current - formula ] , when a dot system does not retain the ph symmetry , the stationary current is not automatically conserved and one must impose the current preservation condition eq .   explicitly . as the right - hand side of eq .   also depends on the dot occupation number @xmath57 , this requires us to determine @xmath57 self - consistently by using the retarded green function in a certain approximation  the second - order perturbation theory in the present case .    ) at finite bias voltage @xmath148 . the interaction strength is chosen as @xmath149 . the dotted line represents the result of @xmath150 and @xmath151 . ] figure  [ fig : ph - symmetric ] shows the result of nonequilibrium dot spectral function of the two - terminal ph symmetric siam at bias voltage @xmath152 , and @xmath153 , which is essentially the same result with ref . ( of a different set of parameters ) . the occupation number is fixed to be @xmath154 in this case , so its self - consistent determination is unnecessary . the results were compared favorably with those obtained by diagrammatic quantum monte carlo calculations  @xcite ; a relatively good quantitative agreement was observed up to @xmath155 ( where the bethe ansatz kondo temperature @xmath156  @xcite while the estimated half width of the kondo resonance @xmath157 ) and bias voltage @xmath158 . applying bias voltage gradually suppresses the kondo resonance without splitting it , and the two peaks at @xmath159 are developed at larger bias voltages , which corresponds to the discussion in the previous section . figure  [ fig : ph - asymmetric ] shows the result of our self - consistent calculation of the nonequilibrium spectral function for the two - terminal ph asymmetric siam at ( a ) @xmath160 and ( b ) @xmath161 . a paramagnetic - type solution is assumed in determining @xmath57 . as in the ph symmetric siam , one sees increasing bias voltage not split but suppress the kondo resonance while it develops a peak around @xmath162 . the kondo resonance peak is suppressed more significantly at @xmath160 than at @xmath163 , because the kondo temperature of the former ( @xmath164 ; @xmath165 ) is smaller than that of the latter ( @xmath166 ; @xmath167 ) . an interesting feature of the ph asymmetric siam is that spectral weight of the kondo resonance seems shifting gradually toward @xmath168 with increasing bias voltage , without exhibiting a three - peak structure in the ph symmetric case . this suggests a strong spectral mixing between the kondo resonance and a coulomb peak at finite bias voltage . because of it , the interval of the two peaks at finite bias is observed as roughly @xmath169 and gets widened up to @xmath0 for larger @xmath170 . the bias dependence somehow looks similar to what was obtained by assuming equilibrium noninteracting effective distribution for @xmath57  @xcite ( which is hard to justify in our opinion ) , though we emphasize our present calculation only relies on the current - preservation condition without using any further assumption . it is remarked that the effect shown by bias voltage is quite reminiscent of finite temperature effect that was observed in the ph asymmetric siam in equilibrium  @xcite .    , and ( b ) @xmath171 . all the other parameters are the same with fig . [ fig : ph - symmetric ] . as an eye guide , the ph symmetric result of @xmath150 and @xmath151 is shown as a dotted line.,title="fig : " ] , and ( b ) @xmath171 . all the other parameters are the same with fig . [ fig : ph - symmetric ] . as an eye guide , the ph symmetric result of @xmath150 and @xmath151 is shown as a dotted line.,title="fig : " ]    . results of the interaction strengths @xmath172 , @xmath173 , and @xmath174 are shown , while dotted lines refer to the noninteracting case as an eye guide . ( a ) spectral function of the ph symmetric siam at @xmath175 . ( b ) spectral function of the ph asymmetric siam at @xmath176.,title="fig : " ] . results of the interaction strengths @xmath172 , @xmath173 , and @xmath174 are shown , while dotted lines refer to the noninteracting case as an eye guide . ( a ) spectral function of the ph symmetric siam at @xmath175 . ( b ) spectral function of the ph asymmetric siam at @xmath176.,title="fig : " ]    more insight can be gained by examining how the spectral structure depends on the interaction strength at finite bias voltage . figures  [ fig : ldos - u ] ( a , b ) show a structural crossover from a noninteracting resonant peak ( the dotted line ) to correlation peaks , for ( a ) the ph symmetric siam @xmath175 , and ( b ) the ph asymmetric siam @xmath177 . the ph symmetric siam shows introducing @xmath0 leads to developing the correlation two peaks as well as the kondo peak that is suppressed by finite bias voltage . in contrast , the bias voltage effect on the ph asymmetric siam is more involved , because the kondo resonance is apparently shifted and mixed with one of the correlation peaks , eventually showing the two - peak structure at @xmath178 in large bias voltage limit .      , and @xmath179 ) . other parameters are chosen as the same as in fig . [ fig : ph - symmetric ] . the dotted line correspond to the two - terminal noninteracting unbiased case , while the dashed line , to the two - terminal interacting unbiased case . ]    to examine finite bias affects further and see particularly how the presence of multi - terminals affects the nonequilibrium spectral function , we configure a special setup of the multi - terminal siam that preserves the ph symmetry : the dot is connected with @xmath180 identical terminals , with bias levels being distributed equidistantly between @xmath181 and @xmath182 , and each of relaxation rates is set to be @xmath183 . the latter ensures that the unbiased spectral function is the same , hence the kondo temperature . results of the nonequilibrium spectral function are shown in fig . [ fig : nldos - by - v ] . again , we confirm that no splitting of the kondo resonance is observed in this multi - terminal setting . one sees further that the increasing the number of terminals _ enhances _ the kondo resonance . this can be understood by weakening the bias suppression effect on the kondo resonance for a larger @xmath180 . more precisely , one may estimate the suppressing effect from small bias behavior eq .  . hence @xmath184 is a control parameter . in the present multi - terminal ph symmetric setting , the quantity @xmath185 is found to be @xmath186 therefore @xmath185 decreases with increasing @xmath180 , which results in weakening the suppression and enhancing the kondo resonance for a larger @xmath180 . ( @xmath187 ) . other parameters are the same as in fig . [ fig : nldos - by - v ] . ( b )  the effective fermi distribution @xmath188 at zero temperature for the ph symmetric multi - terminal dot ( @xmath189 ) . the inset shows relative locations of bias levels with fixed @xmath190 as a function of the number of terminals.,title="fig : " ]   ( @xmath187 ) . other parameters are the same as in fig . [ fig : nldos - by - v ] . ( b )  the effective fermi distribution @xmath188 at zero temperature for the ph symmetric multi - terminal dot ( @xmath189 ) . the inset shows relative locations of bias levels with fixed @xmath190 as a function of the number of terminals.,title="fig : " ]    the preceding argument also tells us that if the spectral function bears any multi - terminal signatures at all , they would be more conspicuous by examining it with fixing @xmath185 rather than @xmath125 . this is done in fig . [ fig : nldos - by - dmu ] ( a ) ; figure  [ fig : nldos - by - dmu ] ( b ) shows how the effective dot distribution @xmath188 and the relative locations of bias levels ( in the inset ) evolve for a fixed @xmath185 when @xmath180 increases . no multi - terminal signature in the nonequilibrium spectral function is seen in fig . [ fig : nldos - by - dmu ] ( a ) ; results of different @xmath180 actually collapse , not only around zero frequency but in the entire frequency range . it suggests that the suppression of the kondo resonance deeply correlates the development of coulomb peaks , and a mixing between those spectral weight is important . the parameter @xmath185 controls a crossover from the kondo resonance to the coulomb blockade structure . we may also understand the similarity between bias effect and temperature effect by the connection through the large @xmath180 limit of the effective fermi distribution @xmath188 , as shown in fig . [ fig : nldos - by - dmu ] ( b ) . although there is a consensus that bias voltage starts suppressing the kondo peak , and eventually destroys it with developing the two coulomb peaks when bias voltage is much larger than the kondo temperature , there is a controversy as to whether the kondo resonance peak will be split or not in the intermediate range of bias voltage . all the results obtained by the second - order perturbation consistently indicate that there is no split of the kondo resonance ; finite bias voltage starts to suppress the kondo resonance , and develops the two coulomb blockade peaks by shifting the spectral weight from the kondo resonance . we should mention that some other approximations draw a different conclusion . here we make a few remarks on apparent discrepancy seen in various theoretical approaches as well as experiments , as well as some speculation based on our results . typically , several approaches that rely on the infinite @xmath0 limit , notably non - crossing approximation , equation of motion method , and other approaches investigating the kondo hamiltonian , observed the splitting of the kondo resonance under finite bias voltage  @xcite . those results , however , have to be interpreted with great care , in our view . generally speaking , the spectral function obtained by those approaches does not obey the spectral weight sum rule : ignoring the doubly occupied state typically leads to the sum rule @xmath191  @xcite , rather than the correct value . therefore only half of the spectral weight can be accounted for in those methods . simultaneously , such ( false ) sum rule in conjugation with the bias - suppression of the kondo resonance can not hep but lead to a two - peak structure of the spectrum within the range of attention . splitting of the kondo resonance might be an artifact of the approximation . not fulfilling the correct sum rule , those approaches may not be able to distinguish whether finite bias will split the kondo resonance or simply suppresses it with developing the coulomb peaks . as for the two - terminal ph symmetric siam , fourth - order contribution regarding the coulomb interaction @xmath0 has been evaluated numerically  @xcite . the results seem unsettled , though . while @xcite suggested the fourth - order term may yield the splitting of the kondo resonance in the intermediate bias region @xmath192 for sufficiently large interaction @xmath193 , which the second - order calculation fails to report , another numerical study indicates that the spectral function remains qualitatively the same with the second - order result  @xcite . experimental situation is not so transparent , either . while the splitting of the kondo resonance was reported in a three - terminal conductance measurement in a quantum ring system  @xcite , a similar spitting observed in the differential conductivity was attributed to being caused by a spontaneous formation of ferromagnetic contacts , not purely to bias effect  @xcite . it is also pointed out that it has been recently recognized that the rashba spin - orbit coupling induces spin polarization nonmagnetically in a quantum ring system with a dot when applying finite bias voltage  @xcite ; hence such spin magnetization might possibly lead to the splitting of the kondo resonance . the kondo resonance is a manifestation of singlet formation between the dot and the lead electrons . one may naively think that when several chemical potentials are connected with the dot , such singlet formation would take place at each lead _ separately _ , causing multiple kondo resonances . the results of the multi - terminal ph symmetric siam presented in the previous section tempt us to speculate a different picture . let us suppose that ( almost ) the same dot distribution function @xmath194 is realized for a fixed @xmath185 with a different terminal number @xmath180 , as fig . [ fig : nldos - by - dmu ] ( a ) suggests . note the assumption is fully consistent to the ward identity for low bias , but it invalidates a quasiparticle ansatz @xmath195 that explicitly depends on @xmath180 . in the large @xmath180 limit with a fixed @xmath185 , the effective fermi distribution @xmath188 resembles the fermi distribution at finite temperature @xmath196 . accordingly bias voltage may well give effect similar to finite temperature . it is seen in low and large bias limit for a dot with or without the ph symmetry . it implies that a dot electron can not separately form a singlet with the lead at each chemical potential , because it needs to implicates states at different chemical potentials through coupling with other leads . our second - order perturbation results seem to support this view . in summary , we have evaluated analytically the second - order self - energy and green function for a generic multi - terminal single impurity anderson model in nonequilibrium . various limiting behaviors have been examined analytically . nonequilibrium spectral function that preserves the current is constructed and is checked to satisfy the spectral weight sum rule . the multi - terminal effect is examined for the ph - symmetric siam , particularly . within the validity of the present approach , it is shown that the kondo peak is not split due to bias voltage . it is found that most of the finite bias effect is similar to that of finite temperature in low- and high - bias limits with and without the ph symmetry . such nature could be understood by help of the ward identity and the connection through the @xmath197 terminal limit . the present analysis serves as a viable tool that can cover a wide range of experimental situations . although there is still a chance that high - order contributions might generate a new effect such as split kondo resonances in a limited range of parameters , it is believed that the second - order perturbation theory can capture the essence of the kondo physics in most of realistic situations . moreover , having a concrete analytical form that satisfies both the current conservation and the sum rule , the present work provides a good , solid , workable result that more sophisticated future treatment can base on . the author gratefully acknowledges a.  sunou for fruitful collaboration that delivered some preliminary results in this work . the author also appreciates r.  sakano and a.  oguri for helpful discussion at the early stage of the work . the work was partially supported by grant - in - aid for scientific research ( c ) no . 22540324 and 26400382 from mext , japan . [ sec : appendix ] we start with the nonequilibrium green function @xmath198 without the coulomb interaction on the dot . its keldysh structure is specified by @xmath199 where @xmath36 is the effective dot distribution defined in eq .  . we incorporate the hartree - type diagram into the unperturbed part by replacing @xmath20 to @xmath200 . note @xmath68 is the exact dot occupation that needs to be determined consistently later . after mapping to the rak representation , those components are given by @xmath201    \left [ g^{r}_{\sigma } ( \omega ) - g^{a}_{\sigma } ( \omega ) \right ] . \label{eq : gk0}\end{aligned}\ ] ] the function @xmath202 reduces to @xmath203 at zero temperature . taking the fourier transformation of eqs .  , using eq .  , and making further manipulations , we can rewrite @xmath204 and @xmath205 as @xmath206^ { * } ,   \label{eq : pir}\\   \pi^{k}(\varepsilon ) & =   2i\sum_{a , b } \frac{\gamma_{a } \gamma_{b}}{\gamma^{2 } }   \coth      \tfrac{\beta(\varepsilon-\mu_{ab})}{2 }    { \mathop{\mathrm{im}}\nolimits}\left [ \frac{\gamma      b_{ab}(\varepsilon)}{\pi \varepsilon ( \varepsilon+2i\gamma ) } \right ] , \label{eq : pik}\end{aligned}\ ] ] where @xmath207 , @xmath208 is the inverse temperature , and @xmath209 is given by @xmath210 \!\left [    g^{a}(\varepsilon ' ) - g^{r}(\varepsilon'+\varepsilon ) \right].\end{aligned}\ ] ] in the present work , we are interested in the zero temperature limit , for which @xmath211 becomes @xmath212 . the function @xmath213 in this limit is evaluated as ( with @xmath112 ) @xmath214 this corresponds to a multi - terminal extension of the result obtained for the two - terminal ph symmetric siam . to complete evaluating the remaining integral over @xmath215 of eqs . ( [ eq : i1],[eq : i2 ] ) , we take full advantage of various properties of dilogarithm function @xmath96 . a concrete integral formula we have utilized will be given in appendix d. for the sake of completeness , we here collect its definition and properties necessary to complete our evaluation . [ [ definition ] ] definition + + + + + + + + + +    dilogarithm @xmath96 with a complex argument @xmath216 is not so commonly found in literature . as it is a multi - valued function , we need to specify its branch structure properly . one way to define dilogarithm @xmath96 all over the complex plane consistently is to use the integral representation @xmath217 the multivaluedness of dilogarithm @xmath218 originates from logarithm in the integrand . here we designate the principal value of logarithm as @xmath219 , which is defined by @xmath220 according to eq .  , @xmath96 has a branch - cut just above the real axis of @xmath221 . accordingly , its values just above and below the real axis are different for @xmath221 : @xmath222 but @xmath223 . some special values are known analytically . we need @xmath224 , @xmath225 , @xmath226 and @xmath227 for evaluation later . [ [ sec : functional - relations ] ] functional relations + + + + + + + + + + + + + + + + + + + +    dilogarithm @xmath96 has interesting symmetric properties regarding its argument @xmath228 ; values at @xmath228 , @xmath229 , @xmath230 , @xmath231 , @xmath232 and @xmath233 are all connected with one another . those points are ones generated by symmetric operations @xmath234 and @xmath235 defined by @xmath236 and @xmath237 forms a group . other operations correspond to @xmath238 applying a series of integral transformations in eq . , one can connect the values of dilogarithm at these values with one another  @xcite . note , those functional relations are usually presented only for real arguments . extending them for complex variables needs examining its branch - cut structure carefully . by following and extending the derivations in ref   for complex @xmath216 , we prove that the following functional relations are valid for any complex variable @xmath228 . @xmath239^{2 } ,   \label{eq : dilog - func - rel - start}\\ & { \mathop{\mathrm{li}_{2}}\nolimits}(tz ) =      -{\mathop{\mathrm{li}_{2}}\nolimits}(z ) + \tfrac{\pi^{2}}{6 } - { \mathop{\mathrm{log}}\nolimits}(tz ) { \mathop{\mathrm{log}}\nolimits}z , \\ & { \mathop{\mathrm{li}_{2}}\nolimits}(tstz ) = - { \mathop{\mathrm{li}_{2}}\nolimits}(z ) -\tfrac{1}{2 }      { \mathop{\mathrm{log}}\nolimits}^{2}(stz )   \nonumber \\ & \qquad \qquad \qquad - [ { \mathop{\mathrm{log}}\nolimits}(tz ) + { \mathop{\mathrm{log}}\nolimits}(stz ) ] { \mathop{\mathrm{log}}\nolimits}z , \\ & { \mathop{\mathrm{li}_{2}}\nolimits}\big ( tsz \big ) = { \mathop{\mathrm{li}_{2}}\nolimits}(z ) -      \tfrac{\pi^{2}}{6 }   \nonumber \\ & \qquad \qquad \qquad -\tfrac{1}{2 } { \mathop{\mathrm{log}}\nolimits}^{2}(sz )   - { \mathop{\mathrm{log}}\nolimits}(sz ) { \mathop{\mathrm{log}}\nolimits}(tz ) , \\ & { \mathop{\mathrm{li}_{2}}\nolimits}\big(stz \big ) = { \mathop{\mathrm{li}_{2}}\nolimits}(z ) +      \tfrac{\pi^{2}}{6 }   \nonumber \\ & \qquad \qquad \qquad + \tfrac{1}{2 } { \mathop{\mathrm{log}}\nolimits}^{2}(tz ) + { \mathop{\mathrm{log}}\nolimits}(tz ) { \mathop{\mathrm{log}}\nolimits}(tstz ) . \label{eq : dilog - func - rel - end}\end{aligned}\ ] ] to our knowledge , the above form of extension of functional relations of dilogarithm has not been found in literature . [ [ the - taylor - expansion ] ] the taylor expansion + + + + + + + + + + + + + + + + + + + +    to examine various limiting behaviors , we need the taylor expansion of dilogarithm , which is derived straightforwardly from eq .  , @xmath240 the presence of @xmath241 reflects the branch - cut structure of @xmath96 . in particular , we utilize the following expansion in our analysis . @xmath242 here we derive and present the central integral formula for evaluating eqs . ( [ eq : i1],[eq : i2 ] ) . by performing a simple integral transformation in eq . , we have the integration @xmath243 where all the parameters @xmath244 as well as @xmath228 may be taken as complex numbers . combined with fractional decomposition , we see the following integral be evaluated in terms of dilogarithm : @xmath245 the remaining task to complete calculating @xmath71 in the form of eqs . ( [ eq : sigmar][eq : xi123 ] ) is to collect all the relevant formulas and organize them in a form that conforms to eq .  . to write concisely , we introduce the following notations @xmath247 where the hartree level @xmath248 is defined as before . we express the terms @xmath249 and @xmath250 defined in eqs . ( [ eq : i1],[eq : i2 ] ) as @xmath251 here the leads @xmath95 in @xmath249 as well as @xmath252 in @xmath250 carry spin @xmath253 , while @xmath13 in @xmath250 does spin @xmath16 . singularity on energy integration is prescribed by the principal values . equation enables us to perform and express the above integrals in terms of dilogarithm . the resulting expressions are still complicated , but we can simplify them further using functional relations of dilogarithm eqs  ( [ eq : dilog - func - rel - start][eq : dilog - func - rel - end ] ) . these require straightforward but rather laborious manipulations . in this way , we reach the final expression of @xmath81 of eq .  .
Abstract:

Answer: we study the nonequilibrium spectral function of the single - impurity anderson model connecting with multi - terminal leads . the full dependence on frequency and bias voltage of the nonequilibrium self - energy and spectral function is obtained analytically up to the second - order perturbation regarding the interaction strength @xmath0 . high and low bias voltage properties are analyzed for a generic multi - terminal dot , showing a crossover from the kondo resonance to the coulomb peaks with increasing bias voltage . for a dot where the particle - hole symmetry is not present , we construct a current - preserving evaluation of the nonequilibrium spectral function for arbitrary bias voltage . it is shown that finite bias voltage does not split the kondo resonance in this order , and no specific structure due to multiple leads emerges . overall bias dependence is quite similar to finite temperature effect for a dot with or without the particle - hole symmetry .


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: photoacoustic tomography is an emerging imaging technique that combines the good contrast of optical absorption with the resolution of ultrasound images ( see for instance @xcite ) . in experiments an object is irradiated by a short - pulsed laser beam . depending on the absorption properties of the material , some light energy is absorbed and converted into heat . this leads to a thermoelastic expansion , which causes a pressure rise , resulting in an ultrasonic wave , called photoacoustic signal . the signal is detected by an array of ultrasound transducers outside the object . using this signal the pressure distribution at the time of the laser excitation is reconstructed , offering a 3d image proportional to the amount of absorbed energy at each position . this is the imaging parameter of photoacoustics . common measurement setups rely on small ultrasound sensors , which are arranged _ uniformly _ along simple geometries , such as planes , spheres , or cylinders covering the specimen of interest ( see for instance @xcite below ) .    for the planar arrangement of point - like detectors there exist several approaches for reconstruction , including numerical algorithms based on filtered backprojection formulas and time - reversal algorithms ( see for instance @xcite ) . the suggested algorithm in the present work realizes a fourier inversion formula ( see below ) using the _ non - uniform fast fourier transform _ ( nufft ) . this method has been designed for evaluation of fourier transforms at non - equispaced points in frequency domain , or non - equispaced data points in spatial , respectively temporal domain . the prior is called ner - nufft ( non - equispaced range - non uniform fft ) , whereas the latter is called ned - nufft ( non - equispaced data - non uniform fft ) . both algorithms have been introduced in @xcite . both nufft methods haven proven to achieve high accuracy and simultaneously reach the computational efficiency of conventional fft computations on regular grids @xcite . this work investigates photoacoustic reconstructions from ultrasound signals recorded at _ non - equispaced _ positions on a planar surface . to the best of our knowledge , this is a novel research question in photoacoustics , where the use of regular grids is the common choice . for the reconstruction we propose a novel combination of ned- and ner - nufft , which we call nedner - nufft , based on the following considerations :    1 . the discretization of the analytic inversion formulas ( see ) contains evaluation at non - equidistant sample points in frequency domain . 2 .   in addition , and this comes from the motivation of this paper , we consider evaluation at non - uniform sampling points . the first issue can be solved by a ner - nufft implementation : for * 2d * photoacoustic inversion with * uniformly * placed sensors on a measurement * line * such an implementation has been considered in @xcite . this method was used for biological photoacoustic imaging in @xcite . in both papers the imaging could be realized in 2d because integrating line detectors @xcite have been used for data recording . in this paper , however , the focus is on 3d imaging , because measurements are taken with point sensors . experimentally , we show the applicability and superiority of the ( ned)ner - nufft reconstruction formula in three spatial dimensions , compared with standard interpolation fft reconstruction . to be precise , in this paper we conduct _ three dimensional _ imaging implemented using nedner - nufft , with ultrasound detectors aligned _ non - uniformly _ on a measurement plane . to easily assess the effect of a given arrangement , also 2d numerical simulations have been conducted , to support the argumentation . we quantitatively compare the results with other computational imaging methods : as a reference we use the @xmath0-wave toolbox @xcite with a standard fft implementation of the inversion algorithm . the nedner - nufft yields an improvement of the lateral and axial resolution ( the latter even by a factor two ) . the outline of this work is as follows :    in section [ sec : reconst ] we outline the basics of the fourier reconstruction approach by presenting the underlying photoacoustic model . we state the fourier domain reconstruction formula in a continuous setting . moreover , we figure out two options for its discretization . we point out the necessity of a fast and accurate algorithm for computing the occurring discrete fourier transforms with non - uniform sampling points .    in section [ sec : nufft ] we briefly explain the idea behind the nufft . we state the ner - nufft ( subsection [ subsec : ner_nufft ] ) and ned - nufft ( subsection [ subsec : ned_nufft ] ) formulas in the form we need it to realize the reconstruction on a non - equispaced grid .    in section [ sec : comp1_ner_nufft ] we discuss the 3d experimental setup . the ner - nufft is compared with conventional fft reconstruction . a test chart is used to quantify resolution improvements in comparison to the _ k - wave _ fft reconstruction with linear interpolation . in axial direction this improvement was about 170% while reducing the reconstruction time by roughly 35% .    in section [ sec : comp2_ned_nufft ] we then turn to the nedner - nufft in 2d with simulated data , in order to test different sensor arrangements in an easily controllable environment . an equiangular arrangement turns out to yield an improvement of over 40% compared to the best choice of equispaced sensor arrangement . furthermore , we use the insights gained from the 2d simulations to develop an equi - steradian sensor arrangement for our 3d measurements . we apply our nedner - nufft approach on these data and quantitatively compare the outcomes with reconstructions from equispaced data obtained by the ner - nufft approach . our results show a significant improvement of the already superior ner - nufft . let @xmath1 be an open domain in @xmath2 , and @xmath3 a @xmath4 dimensional hyperplane not intersecting @xmath5 . mathematically , photoacoustic imaging consists in solving the operator equation @xmath6=p|_{\gamma\times(0,\infty)}\,,\ ] ] where @xmath7 is a function with compact support in @xmath5 and @xmath8 $ ] is the trace on @xmath9 of the solution of the equation    @xmath10    in other words , the photoacoustic imaging problem consists in identifying the initial source @xmath7 from measurement data @xmath11 . an explicit inversion formula for @xmath12 in terms of the fourier transforms of @xmath7 and @xmath13 $ ] has been found in @xcite . let @xmath14 . assume without loss of generality ( by choice of proper basis ) that @xmath3 is the hyperplane described by @xmath15 . then the reconstruction reads as follows : @xmath16\left(\bm{k}\right)= \frac{2k_{y}}{\kappa\left(\bm{k}\right)}\mathbf{f}[\mathbf{q}f]\left(\bm k_{\bm x},\kappa\left(\bm{k}\right)\right).\end{aligned}\ ] ] where @xmath17 denotes the @xmath18-dimensional fourier transform : @xmath19\left(\bm{k}\right):=\frac1{(2\pi)^{n/2}}\int\limits_{\mathbb{r}^{d}}{{\mathrm{e}}}^{-{{\mathrm{i}}}\bm{k}\cdot(\bm{x},y)}f(\bm{x})\mathrm{d}\bm{x}\,,\end{aligned}\ ] ] and @xmath20 here , the variables @xmath21 are in @xmath22 , whereas @xmath23 .    for the numerical realization these three steps have to be realized in discrete form : we denote evaluations of a function @xmath24 at sampling points @xmath25 by @xmath26 for convenience , we will modify this notation in case of evaluations on an equispaced cartesian grid . we define the @xmath18-dimensional grid @xmath27 and assume our sampling points to be located on @xmath28 , where @xmath29 and write @xmath30 where @xmath31 resp . @xmath32 are the occurring step sizes .    in frequency domain , we have to sample symmetrically with respect to @xmath33 . therefore , we also introduce the interval @xmath34 since we will have to deal with evaluations that are partially in - grid , partially not necessarily in - grid , we will also use combinations of and . in this paper , we will make use of discretizations of the source function @xmath7 , the data function @xmath35 and their fourier transforms @xmath36 resp . @xmath37 . let in the following @xmath38 denote the @xmath18-dimensional discrete fourier transform with respect to space and time . by discretizing formula via riemann sums it follows @xmath39 where @xmath40 this is the formula from @xcite . note that we use the interval notation for the integer multi - indices for notational convenience . moreover , we also choose the length of the fourier transforms to be equal to @xmath41 in the first @xmath4 dimensions , respectively . this could be generalised without changes in practice . now , we assume to sample @xmath35 at @xmath42 , not necessarily uniform , points @xmath43 : then , @xmath44 the term @xmath45 represents the area of the detector surface around @xmath46 and has to fulfil @xmath47 . note that the original formula can be received from by choosing @xmath48 to contain all points on the grid @xmath49 . formula can be interpreted as follows : once we have computed the fourier transform of the data and evaluated the fourier transform at non - equidistant points with respect to the third coordinate , we obtain the ( standard , equispaced ) fourier coefficients of @xmath7 . the image can then be obtained by applying standard fft techniques . the straightforward evaluation of the sums on the right hand side of would lead to a computational complexity of order @xmath50 . usually this is improved by the use of fft methods , which have the drawback that they need both the data and evaluation grid to be equispaced in each coordinate . this means that if we want to compute efficiently , we have to interpolate both in domain- and frequency space . a simple way of doing that is by using polynomial interpolation . it is used for photoacoustic reconstruction purposes for instance in the _ k - wave _ toolbox for matlab @xcite . unfortunately , this kind of interpolation seems to be sub - optimal for fourier - interpolation with respect to both accuracy and computational costs @xcite    a regularized inverse k - space interpolation has already been shown to yield better reconstruction results @xcite . the superiority of applying the nufft , compared to linear interpolation , has been shown theoretically and computationally by @xcite . this section is devoted to the brief explanation of the theory and the applicability of the non - uniform fourier transform , where we explain both the ner - nufft ( subsection [ subsec : ner_nufft ] ) and the ned - nufft ( subsection [ subsec : ned_nufft ] ) in the form ( and spatial dimensions ) we utilise them afterwards .    the nedner - nufft algorithm used for implementing essentially ( up to scaling factors ) consists of the following steps :    1 . compute a @xmath4 dimensional ned - nufft in the @xmath51-coordinates due to our detector placement . 2 .   compute a one - dimensional ner - nufft in the @xmath33-coordinate as indicated by the reconstruction formula . 3 .   compute an equispaced @xmath18-dim inverse fft to obtain a @xmath18 dimensional picture of the initial pressure distribution .      with the ner - nufft ( non equispaced range  non - uniform fft ) it is possible to efficiently evaluate the discrete fourier transform at non - equispaced positions in frequency domain . to this end , we introduce the one dimensional discrete fourier transform , evaluated at non - equispaced grid points @xmath52 : @xmath53 in order to find an efficient algorithm for evaluation of , we use a window function @xmath54 , an oversampling factor @xmath55 and a parameter @xmath56 that satisfy : [ en : psi ]    1 .   @xmath54 is continuous inside some finite interval @xmath57 $ ] and has its support in this interval and 2 . @xmath54 is positive in the interval @xmath58 $ ] . then ( see @xcite ) we have the following representation for the fourier modes occurring in ( [ eq : nudft1 ] ) : @xmath59 by assumption , both @xmath54 and @xmath60 are concentrated around @xmath61 . so we approximate the sum over all @xmath62 by the sum over the @xmath63 integers @xmath0 that are closest to @xmath64 . by choosing @xmath65 and inserting in , we obtain @xmath66 here @xmath67 denotes the interpolation length and @xmath68 where @xmath69 is the nearest integer ( i.e. the nearest equispaced grid point ) to @xmath70 .    the choice of @xmath54 is made in accordance with the assumptions above , so we need @xmath54 to have compact support . furthermore , to make the approximation in ( [ eq : ner_nufft ] ) reasonable , its fourier transform @xmath71 needs to be concentrated as much as possible in @xmath72 $ ] . in practice , a common choice for @xmath54 is the kaiser - bessel function , which fulfils the needed conditions , and its fourier transform is analytically computable . a second major aim of the present work is to handle data measured at non - equispaced acquisition points @xmath73 in an efficient and accurate way . therefore we introduce the non - equispaced data , @xmath4 dimensional dft @xmath74 the theory for the ned - nufft is largely analogous to the ner - nufft @xcite as described in subsection [ subsec : ner_nufft ] . the representation is here used for each entry of @xmath75 and inserted ( with now setting @xmath76 ) into formula , which leads to @xmath77 where the entries in @xmath78 are the nearest integers to @xmath79 . here we have used the abbreviations @xmath80 for the needed evaluations of @xmath54 and @xmath71 . further remarks on the implementation of the ned- and ner - nufft , as well as a summery about the properties of the kaiser - bessel function and its fourier transform can be found in @xcite . before we turn to the evaluation of the algorithm we describe the photoacoustic setup . our device consists of a fp ( fabry prot ) polymer film sensor for interrogation @xcite . a @xmath81 pulsed laser source and a subsequent optical parametric oscillator ( opo ) provide optical pulses . these pulses have a very narrow bandwidth and can be tuned within the visible and near infrared spectrum . the optical pulses are then transmitted via an optical fibre . when the light is emitted it diverges and impinges upon a sample with homogeneous fluence , thus generating a photoacoustic signal . this signal is then recorded via the fp - sensor head . the sensor head consists of an approximately @xmath82 thick polymer ( parylene c ) which is sandwiched between two dichroic dielectric coatings . these dichroic mirrors have a noteworthy transmission characteristic . light from @xmath83 to @xmath84 can pass the mirrors largely unabated , whereas the reflectivity from @xmath85 to @xmath86 ( sensor interrogation band ) is about 95% @xcite . the incident photoacoustic wave produces a linear change in the optical thickness of the polymer film . a focused continuous wave laser , operating within the interrogation band , can now determine the change of thickness at the interrogation point via fp - interferometry . we choose two 2d targets for comparison , a star and a usaf ( us air force ) resolution test chart . both targets are made of glass with a vacuum - deposited durable chromium coating . the star target has 72 sectors on a pattern diameter of 5 mm , with an unresolved core diameter of @xmath87 . the targets are positioned in parallel to the sensor surface at a distance of about @xmath88 , and water is used as coupling medium between the target and the sensor . the interpolation length for the ner - nufft reconstruction is @xmath89 . the computational times are shown in tab . [ tab : ner - nufft comparison ] showing that the linearly interpolated fft is about 30 % slower than the ner - nufft . -axis of a star sample , reconstructed using fft and linear interpolation ( top left ) and using the nufft ( top right ) with a square plotted around the center . the intensity of the reconstructed image along the sidelines of the square is plotted , for both reconstructions . the purple line indicates the frequency of the star sample.[fig : starphantomreconstruction ] ]    a segment of the reconstructed star target is shown in fig . [ fig : starphantomreconstruction ] . the intensity is plotted against the sides of an imaginary square ( @xmath90 ) placed around the center of the star phantom . it is clearly visible , that the fft reconstruction is not able to represent the line pairs , when the density exceeds @xmath91 , corresponding to a resolution of @xmath87 , whereas they are still largely visible for the ner - nufft reconstruction .    [ cols="^,^,^,^",options="header " , ]     for a quantitative comparison of the resolution we use the usaf chart . it is recorded on an area of @xmath92 sensor points corresponding to @xmath93 with a grid spacing of @xmath94 and a time resolution of @xmath95 . as yet there is no standardized procedure to measure the resolution of a photoacoustic imaging system . we proceed similar to @xcite , by fitting a _ line spread function _ ( lsf ) and an _ edge spread function _ ( esf ) to the intensities of our reconstructed data . for the lsf to be meaningful , its source has to approximate a spatial delta function . this is the case in the @xmath96-axis , since the chrome coating of the usaf target is just about @xmath97 thick . -mip of an usaf chart reconstruction conducted with fft ( left ) and ner - nufft ( right ) . on the bottom images the data - sets used for the quantitative resolution are marked . the black square depicts the 49 @xmath98-coordinates used for the lsf fit along the @xmath96-axis for the axial resolution . the data for a single point is shown in the bottom right inlay . the white lines show the intensity fit for the lateral resolution . the bottom left inlay shows the data for a single white line . [ fig : usaf - chart ] ]    we fit a cauchy - lorentz distribution to the @xmath96-axis of our reconstructed data for 49 adjacent @xmath98-coordinates . their positions are marked as a black square within the white square , depicted in the bottom images of fig . [ fig : usaf - chart ] . the reconstruction in the @xmath96-axis for a single point is shown in the bottom right inlay of fig . [ fig : usaf - chart ] , for 8 points , covering a distance of @xmath99 in the @xmath96-direction , around the maximum intensity . a fit of the lorentz distribution is shown for both reconstruction methods . the fwhm ( _ full width , have maximum _ ) of the lorentz distribution , @xmath100 is the parameter @xmath101 . the output @xmath102 is the intensity in dependence from the @xmath96-axis , and @xmath103 and @xmath104 are fitting parameters . the average and standard deviation of @xmath101 for both of the 49 datasets are shown in table [ tab : ner - nufft comparison ] . the line spread function fwhm of the fft - reconstruction turns out to be more than twice as big as the one of the ner - nufft reconstruction .    for the lateral resolution , there is no target approximating a delta function , so we have to use the esf instead : @xmath105 the @xmath101 here is the fwhm of the associated lsf , and @xmath106 , @xmath107 and @xmath104 are fitting parameters . the esf requires a step function as source , of which our target provides plenty . we choose the long side of 12 bars , marked by white lines in the bottom images of fig . [ fig : usaf - chart ] , for this fit . the data for a particular line is shown in the bottom left inlay . we omitted all datasets , where only one point marked the transition from low to high intensity , rendering the edge fit unreliable and resulting in unrealistic improvements of our new method well over 100% . finally we averaged over 15 edges . the results are shown in table [ tab : ner - nufft comparison ] . while the deviation between different edges is quite high , the improvement of our ner - nufft reconstruction for the 15 evaluated edges ranged from @xmath108 to @xmath109 . the current setups allow data acquisition at just one single sensor point for each laser pulse excitation . since our laser is operating at @xmath81 data recording of a typical sample requires several minutes . reducing this acquisition time is a crucial step in advancing photoacoustics towards clinical and preclinical application . therefore , in this work we try to maximize the image quality for a given number of acquisition points . we are able to do this , because our newly implemented nedner - nufft is ideal for dealing with non - equispaced positioned sensors . this newly gained flexibility of sensor positioning offers many possibilities to enhance the image quality , compared against a rectangular grid . for instance a hexagonal grid was found to yield an efficiency of @xmath110 compared with @xmath111 for an exact reconstruction of a wave - number limited function @xcite . also any non - equispaced grids that may arrive from a specific experimental setup can be efficiently computed via the nedner - nufft approach .    here , we will use it to tackle the limited view problem . many papers deal with the limited view problem , when reconstructing images @xcite . our approach to deal with this problem is different . it takes into account that in many cases the limiting factor is the number of sensor points and the limited view largely a consequence of this limitation .    in our approach we therefore use a grid arrangement that is dense close to a center of interest and becomes sparser the further away the sampling points are located . we realize this by means of an equiangular , or equi - steradian sensor arrangement , where for a given point of interest each unit angle or steradian gets assigned one sensor point .    ]    to understand the limited view problem , it is helpful to define a detection region . according to @xcite , this is the region which is enclosed by the normal lines from the edges of the sensor . pressure waves always travel in the direction of the normal vector of the boundary of the expanding object . therefore certain boundaries are invisible to the detector as depicted in fig . [ fig : fieldofview ] .      for the equiangular sensor arrangement a point of interest is chosen . each line , connecting a sensor point with the point of interest , encloses a fixed angle to its adjacent line . in that sense we mimic a circular sensor array on a straight line . the position of the sensor points can be seen on top of the third image in fig . [ fig : tree ] . the obvious expansion of an equiangular projection to 3d is an equi - steradian projection . this problem is analogous to the problem of placing equispaced points on a 3d sphere and then projecting the points , from the center of the sphere , through the points , onto a 2d plane outside the sphere . we developed an algorithm for this problem , which is explained in detail in appendix [ app : equi - ster ] . our input variables are the grid size , the distance of the center of interest from the sensor plane and the desired number of acquisition points , which will be rounded to the next convenient value . a sensor arrangement with 1625 points on a @xmath112 grid is shown in the top left image in fig . [ fig : sensor_masks ] .      to determine the weighting term @xmath45 in eq . [ eq : disc_rec ] for 3d we introduce a function that describes the density of equidistant points per unit area @xmath113 . in our specific case , @xmath113 describes the density on a sphere around a center of interest . further we assume that @xmath113 is spherically symmetric and decreases quadratically with the distance from the center of interest @xmath114 : @xmath115 . we now define @xmath116 for a plane positioned at distance @xmath117 from the center of interest . in this case @xmath118 attenuates by a factor of @xmath119 , where @xmath120 is the angle of incidence . hence @xmath121 . by applying the spacing of the regular grid . this yields a weighting term of : @xmath122    analogously we can derive @xmath45 for 2d :    @xmath123    we applied a normalization after the reconstruction to all measurements .    for the application of this method to the fp setup it is noteworthy that there is a frequency dependency on sensitivity which itself depends on the angle of incidence , which has been extensively discussed in @xcite . the angle of incidence for our specific setup is @xmath124 . at this angle , the frequency components around @xmath125 get attenuated by more than @xmath126 . below @xmath127 the frequency response remains quite stable ( attenuation below @xmath128 for the measurement angles occurring in our setup . ) equispaced sensor arrangement , with a distance of 13 points between each sensor . the third image shows the nedner - nufft reconstruction with equiangular arranged sensor positions . the bottom image shows the same sensor arrangement , but all omitted sensor points are linearly interpolated and afterwards a ner - nufft reconstruction was conducted . [ fig : tree ] ]    a tree phantom , designed by brian hurshman and licensed under cc by 3.0 , is chosen for the 2 dimensional computational experiments on a grid with @xmath129 @xmath130 points . a forward simulation is conducted via _ k - wave _ @xcite . the forward simulation of the k - wave toolbox is based on a first order k - space model . a pml ( perfectly matched layer ) of 64 gridpoints is added , as well as @xmath131 of noise .    in fig . [ fig : tree ] our computational phantom is shown at the top . for each reconstruction a subset of 32 out of the 1024 possible sensor positions was chosen . in fig . [ fig : tree ] their positions are marked at the top of each reconstructed image . for the equispaced sensor arrangements , we let the distance between two adjacent sensor points sweep from 1 to 32 . the sensor points where always centered in the @xmath132-axis .    to compare the different reconstruction methods we used the correlation coefficient and the tenenbaum sharpness . these quality measures are explained in appendix [ app : qualitmeas ] . we applied the correlation coefficient only within the region of interest marked by the white circle in fig . [ fig : tree ] . the tenenbaum sharpness was calculated on the smallest rectangle , containing all pixels within the circle . the results are shown in fig . [ fig : tree_cc ] .    . the straight lines indicate the results for the equiangular projection.[fig : tree_cc ] ]    the tenenbaum sharpness for the equiangular sensor placement was 23001 , which is above all values for the equispaced arrangements . the correlation coefficient was 0.913 compared to 0.849 , for the best equispaced arrangement . in other words , the equiangular arrangement is 42.3 % closer to a full correlation , than any equispaced grid .    in fig . [ fig : tree ] the competing reconstructions are compared . while the crown of the tree is depicted quite well for the equispaced reconstruction , the trunk of the tree is barely visible . this is due to the limited view problem . when the equispaced interval increases , the tree becomes visible , but at the cost of the crown s quality . in the equiangular arrangement a trade off between these two effects is achieved . additionally the weighting term for the outmost sensors is 17 times the weighting term for the sensor point closest to the middle . this amplifies the occurrence of artefacts , particularly outside of our region of interest .    the bottom image in fig . [ fig : tree ] shows the equiangular sensor arrangement , reconstructed in a conventional manner . the missing sensor points are interpolated to an equispaced grid , and a ner - nufft reconstruction is applied afterwards . we conducted a linear interpolation from our subset to all 1024 sensor points . the correlation coefficient for this outcome was 0.7348 while the sharpness measure was 21474 . this outcome exemplifies the clear superiority of the nufft to conventional fft reconstruction when dealing with non - equispaced grids . ]    for a qualitative assessment of our new sensor arrangement we need a 3d phantom . we choose a yarn which we record on a rectangular grid with @xmath133 sensor points , with a grid spacing of @xmath134 and time sampling of @xmath135 hence an area of @xmath136 is covered . as coupling medium water is used , in which the yarn is fully immersed .    to determine the utility of non - equispaced grid sampling , we follow a certain routine . first we acquire a densely sampled dataset . then we use a very small subset of the initially collected sensor data , to test different sensor arrangements . therefore we can always use the complete reconstruction as our model standard together with the quality measurements explained in appendix [ app : qualitmeas ] . an upsampling factor of 2 was used for all reconstructions , hence the reconstructed image of the mip for the @xmath98 plane consists of @xmath137 pixels . the complete reconstruction with the ner - nufft took 154 seconds . maximum intensity projections ( mip ) of this full reconstruction for all axis are shown in fig . [ fig : golden - standard ] . the mip in the @xmath98 plane is our model standard , for comparison with the other reconstructions . sensor points is displayed . the intervall between two adjacent sensor points is @xmath138 for this configuration . the bottom images show mips of the ned - ner - nufft reconstruction only using the sensor points shown above . [ fig : sensor_masks ] ]    for the equi - steradian sensor mask we choose our center of interest right in the center of the @xmath98-mip where the little knot can be seen , @xmath139 off the sensor surface . the resulting sensor mask , including the reconstruction is shown in fig . [ fig : sensor_masks ] , it consists of @xmath140 sensor points ( or 3.18 % of the initial number of sensor points ) . the weighting term , accounting for sensor sparsity is 9.7 times higher for the outermost sensor point , than for centermost sensor points in the @xmath98-plane . the reconstruction for this arrangement took 134 seconds . we compared this arrangement to rectangular grids , which all had @xmath141 sensor points ( or 3.29 % of the initial number of sensor points ) but varying distances between two adjacent points . the grid with an interval between sensor points of 5 is shown on the top left in fig . [ fig : sensor_masks ] .    ) . in the inlay the correlation coefficient mask is shown for 20087 points which corresponds to all pixels with at least 4% of the maximum value of the model standard and is the value used in fig . [ fig : correlation - coefficient - maskdiameter ] . [ fig : correlation_coefficient_mask ] ]    in fig . [ fig : correlation_coefficient_mask ] the equi - steradian grid is compared with 2 equispaced grids . to get a more precise measure of the correlation between the reconstructed images and the model standard , we calculated the correlation coefficient only within a region of interest . the region of interest is firstly confined by a centered disc , whose boundary is shown as a dotted circle in the inlay in fig . [ fig : correlation_coefficient_mask ] . the @xmath132-axis shows the number of pixels of interest used to calculate the correlation coefficient within this disc . these pixels are increased according to the intensity of the corresponding pixels of the model s standard mip . [ fig : correlation_coefficient_mask ] demonstrates that the correlation coefficient for the equi - steradian arrangement , always remains better , within the depicted disc , when competing against the two strongest equispaced grids . corresponds to the straight line in fig . [ fig : correlation_coefficient_mask ] . the other straight lines indicate the side length of the square of a particular equispaced sensor point arrangement . [ fig : correlation - coefficient - maskdiameter ] ]    in fig . [ fig : correlation - coefficient - maskdiameter ] the 4 equispaced arrangements , with intervals reaching from 2 to 5 are compared to the equi - steradian grid . here the pixels of interest are set to a threshold of at least 4% of the maximum value , while the diameter is increased . the dotted straight lines indicate the side length of the square of a particular equispaced sensor point arrangement . as expected the correlation coefficient for this equispaced sensor arrangement starts to decline around that threshold . the correlation coefficient for the equi - steradian grid starts to fall behind the equispaced grid of interval 5 towards the end . this is expected , since the equi - steradian grid is only meant to give better results for a region of interest around the center . this is very clearly the case . between a diameter of 2 to 8 mm the correlation coefficient is on average 25.14 % closer to a value of @xmath142 than its strongest equispaced contender of interval 4 , and 30.87 % better than the interval 5 grid . at a diameter of @xmath143 the correlation coefficient for the equi - steradian grid was 0.960 compared to 0.947 for the interval 4 grid and 0.944 for the interval 5 grid . this is 24.8 % and 28.6% closer to full correlation .    ] [ fig : tenenbaum - sharpness ] shows the normalized tenenbaum sharpness . the tenenbaum sharpness , unlike the correlation coefficient , can not be calculated on non - adjacent grid points , therefore it has been calculated on the smallest rectangle , that contains all pixels of interest . the equi - steradian has the highest sharpness for most values , with the interval 4 grid being very slightly better around a diameter of 3 - 4 mm . there is a drop of the tenenbaum sharpness towards the end . we computationally implemented a 3d non - uniform fft photoacoustic image reconstruction , called ner - nufft ( non equispaced range - non uniform fft ) to efficiently deal with the non - equispaced fourier transform evaluations arising in the reconstruction formula . this method was compared with the _ k - wave _ implemented fft reconstruction , which uses a polynomial interpolation . the two reconstruction methods where compared using 2d targets . the lateral resolution showed an improvement of @xmath144 which is in good agreement with the illustrative results for the star target . the axial resolution showed an improvement of @xmath145 . the computation time was about @xmath146 less , for the ner - nufft than the linearly interpolated fft reconstruction . in conclusion the ner - nufft reconstruction proved to be unequivocally superior to conventional linear interpolation fft reconstruction methods . we further implemented the ned - ner - nufft ( non equispaced data - ner - nufft ) , which allowed us to efficiently reconstruct from data recorded at non - equispaced placed sensor points . this newly gained flexibility was used to tackle the limited view problem , by placing sensors more sparsely further away from the center of interest . we developed an equiangular sensor placement for 2d and an equi - steradian placement in 3d , which assigns one sensor point to each angle / steradian for a given center of interest . in the 2d computational simulations we showed that this arrangement significantly enhances image quality in comparison to regular grids . in 3d we conducted experiments , where a yarn phantom was recorded . the maximum intensity projection ( mip ) of the full reconstruction was compared to mips of reconstructions that only used about 3% of the original data . within our region of interest , the correlation of our image was 0.96 , which is 24.8% closer to full correlation than the best equispaced arrangement , reconstructed from slightly more sensor points . the sensor placement to tackle the limited view problem , combined with the ned - ner - nufft gives significantly better results for an object located at the center of a bigger sensor surface . this result was confirmed in the 2d simulation as well as for real data in 3d . this work is supported by the medical university of vienna , the european projects famos ( fp7 ict 317744 ) and fun oct ( fp7 health 201880 ) , macular vision research foundation ( mvrf , usa ) , austrian science fund ( fwf ) , project p26687-n25 ( interdisciplinary coupled physics imaging ) , and the christian doppler society ( christian doppler laboratory `` laser development and their application in medicine '' ) . [ 1]#10=0=0 0 by1pt#1    10 [ 1]#1 urlstyle [ 1]doi  # 1    p.  beard . biomedical photoacoustic imaging . , 1:602631 , 2011 . p.c . beard . two - dimensional ultrasound receive array using an angle - tuned fabry - perot polymer film sensor for transducer field characterization and transmission ultrasound imaging . , 52(6):10021012 , june 2005 . beard , f.  perennes , and t.n . transduction mechanisms of the fabry - perot polymer film sensing concept for wideband ultrasound detection . , 46(6):15751582 , november 1999 . p.  burgholzer , c.  hofer , g.  paltauf , m.  haltmeier , and o.  scherzer . thermoacoustic tomography with integrating area and line detectors . , 52(9):15771583 , september 2005 . cox and p.c . the frequency - dependent directivity of a planar fabry - perot polymer film ultrasound sensor . , 54(2):394404 , 2007 . w.  dan , c.  tao , x .- j . liu , and x .- d . influence of limited - view scanning on depth imaging of photoacoustic tomography . , 21(1):014301 , 2012 . k.  fourmont . non - equispaced fast fourier transforms with applications to tomography . , 9(5):431450 , 2003 .    f.  c.  a. groen , i.  t. young , and g.  ligthart . a comparison of different focus functions for use in autofocus algorithms . , 6(2):8191 , 1985 . m.  haltmeier , o.  scherzer , and g.  zangerl . a reconstruction algorithm for photoacoustic imaging based on the nonuniform fft . , 28(11):17271735 , november 2009 . m.  jaeger , s.  schpbach , a.  gertsch , m.  kitz , and m.  frenz . fourier reconstruction in optoacoustic imaging using truncated regularized inverse k - space interpolation . , 23:s51s63 , 2007 . p.  kuchment and l.  kunyansky . mathematics of thermoacoustic tomography . , 19:191224 , 2008 . g.  paltauf , r.  nuster , m.  haltmeier , and p.  burgholzer . photoacoustic tomography with integrating area and line detectors . in l.  v. wang , editor , _ photoacoustic imaging and spectroscopy _ , optical science and engineering , pages 251263 . crc press , boca raton , fl , 2009 . d.  p. petersen and d.  middleton . sampling and reconstruction of wave - number - limited functions in @xmath147-dimensional euclidean spaces . , 5:279323 , 1962 . o.  scherzer , editor . . springer , new york , 2011 . r.  schulze , g.  zangerl , m.  holotta , d.  meyer , f.  handle , r.  nuster , g.  paltauf , and o.  scherzer . on the use of frequency - domain reconstruction algorithms for photoacoustic imaging . , 16(8):086002 , 2011 . funded by the austrian science fund ( fwf ) within the fsp s105 - `` photoacoustic imaging '' .    c.  tao and x.  liu . reconstruction of high quality photoacoustic tomography with a limited - view scanning . , 18(3):27602766 , feb 2010 . b.  e. treeby and b.  t. cox . -wave : matlab toolbox for the simulation and reconstruction of photoacoustic wace fields . , 15:021314 , 2010 . b.  e. treeby , t.  k. varslot , e.  z. zhang , j.  g. laufer , and p.  c. beard . automatic sound speed selection in photoacoustic image reconstruction using an autofocus approach . , 16(9):0905010905013 , 2011 . l.  v. wang , editor . . optical science and engineering . crc press , boca raton , 2009 . m.  xu and l.  v. wang . exact frequency - domain reconstruction for thermoacoustic tomography  i : planar geometry . , 21:823828 , 2002 . m.  xu and l.  v. wang . time - domain reconstruction for thermoacoustic tomography in a spherical geometry . , 21(7):814822 , 2002 . m.  xu and l.  v. wang . analytic explanation of spatial resolution related to bandwidth and detector aperture size in thermoacoustic or photoacoustic reconstruction . , 67(5):056605105660515 , 2003 . m.  xu and l.  v. wang . universal back - projection algorithm for photoacoustic computed tomography . , 71(1 ) , 2005 . m.  xu and l.  v. wang . photoacoustic imaging in biomedicine . , 77(4 ) , 2006 . y.  xu , d.  feng , and l.  v. wang . exact frequency - domain reconstrcution for thermoacoustic tomography  i : planar geometry . , 21(7):823828 , 2002 . y.  xu , l.  v. wang , g.  ambartsoumian , and p.  kuchment . reconstructions in limited - view thermoacoustic tomography . , 31(4):724733 , 2004 . e.  zhang , j.  laufer , and p.  beard . backward - mode multiwavelength photoacoustic scanner using a planar fabry - perot polymer film ultrasound sensor for high - resolution three - dimensional imaging of biological tissues . , 47:561577 , 2008 . in our algorithm , the grid size and the distance of the center of interest from the sensor plane is defined . the number of sensor points @xmath147 will be rounded to the next convenient value . our point of interest is placed at @xmath148 , centered at a square @xmath98 grid . the point of interest is the center of a spherical coordinate system , with the polar angle @xmath149 at the @xmath96-axis towards the @xmath98-grid . first we determine the steradian @xmath150 of the spherical cap from the point of interest , that projects onto the acquistion point plane via @xmath151 this leads to a unit steradian @xmath152 with @xmath147 being the number of sensors one would like to record the signal with . the sphere cap is then subdivided into slices @xmath0 which satisfy the condition @xmath153 where @xmath154 encloses exactly one unit steradian @xmath155 and @xmath156 has to be a power of two , in order to guarantee some symmetry . the value of @xmath156 doubles , when @xmath157 , where @xmath158 is the chord length between two points on @xmath0 and @xmath159 is the distance to the closest point on @xmath160 . these values are chosen in order to guarantee we apply some restrictions to approximate equidistance between acquisition points on the sensor surface . the azimuthal angles for a slice @xmath0 are calulated according to : @xmath161 with @xmath162 , where @xmath163 stems from the former slice @xmath160 . the sensor points are now placed on the @xmath98-plane at the position indicated by the spherical angular coordinates : @xmath164 the correlation coefficient @xmath165 is a measure of the linear dependence between two images @xmath166 and @xmath167 . its range is @xmath168 $ ] and a correlation coefficient close to 1 indicates linear dependence @xcite . it is defined via the variance , @xmath169 of each image and the covariance , @xmath170 of the two images : @xmath171    we did not choose the widely used @xmath172 distance measure because it is a morphological distance measure , meaning it defines the distance between two images by the distance between their level sets . therefore two identical linearly dependent images can have a correlation coefficient of @xmath142 and still a huge @xmath172 distance . normalizing the images only mitigates this problem , because single , high intensity artifacts or a varying intensity over the whole image can greatly alter the @xmath172 distance .    while the correlation coefficient is a good measure for the overall similarity between two images it does not include any sharpness measure . hence blurred edges are punished very little , in comparison to slight variations of homogeneous areas . to address this shortcoming we chose a sharpness measure or focus function as a second quality criterion . sharpness measures are obtained from some measure of the high frequency content of an image @xcite . they have also been used to select the best sound speed in photoacoustic image reconstruction @xcite . out of the plethora of published focus functions we select the tenenbaum function , because of its robustness to noise :        like the @xmath175 norm and unlike the correlation coefficient the tenenbaum function is an extensive measure , meaning it increases with image dimensions . therefore we normalized it to @xmath176 , where @xmath147 is the number of elements in @xmath5 .
Abstract:

Answer: to obtain the initial pressure from the collected data on a planar sensor arrangement in photoacoustic tomography , there exists an exact analytic frequency domain reconstruction formula . an efficient realization of this formula needs to cope with the evaluation of the data s fourier transform on a non - equispaced mesh . in this paper , we use the non - uniform fast fourier transform to handle this issue and show its feasibility in 3d experiments . this is done in comparison to the standard approach that uses polynomial interpolation . moreover , we investigate the effect and the utility of flexible sensor location on the quality of photoacoustic image reconstruction . the computational realization is accomplished by the use of a multi - dimensional non - uniform fast fourier algorithm , where non - uniform data sampling is performed both in frequency and spatial domain . we show that with appropriate sampling the imaging quality can be significantly improved . reconstructions with synthetic and real data show the superiority of this method . * keywords : * image reconstruction , photoacoustics , non - uniform fft