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: synchronization of chaotic systems has been an active research area in recent years  @xcite . now the notion of synchronization of chaotic systems is extended far beyond complete synchronization between identical systems@xcite , and various kinds of nonlinear synchronization have been proposed  @xcite . a significant extension is _ generalized synchronization _ ( gs ) @xcite , which is defined by a time - independent nonlinear functional relation @xmath0 between the states @xmath1 and @xmath2 of two systems @xmath3 and @xmath4 . experimental detection and characterization of gs from observed data is a challenging problem , especially in biology ; e.g. , for study on nonlinear interdependence observed in binding of different features in cognitive process  @xcite and epilepsies in the brain  @xcite . in unidirectionally coupled systems , a way to detect gs is to make an identical copy @xmath5 of the response system @xmath4 driven by the common signal from the driver system @xmath3 , then investigate whether orbits of both @xmath4 and @xmath5 coincide after transient  @xcite . however , it is very difficult to prepare an identical copy of the response system , and this approach does not give any information on the structure of the synchronization manifold @xmath6 in the state space . on the other hand , several indices have been proposed  @xcite , to quantify nonlinear dependence between @xmath3 and @xmath4 based on a local relation between them . while the usefulness of these indices were shown in some examples  @xcite , those approaches are still ad - hoc , lacking systematic methods for unifying the local relation of different regions in the state space . recently , _ kernel methods _ have been attracting much attention in the machine learning community as powerful tools for analyzing data with high nonlinearity  @xcite . in this paper , we propose a novel approach for analyzing gs based on _ kernel canonical correlation analysis _ ( kernel cca )  @xcite which is a version of the kernel methods  @xcite . we will demonstrate that the kernel cca provides a suitable index for measuring the nonlinear interdependence , and gives global coordinates characterizing the synchronization manifold @xmath7 of gs . note that the latter has not been addressed in conventional studies on the analysis of gs . the proposed method does not require explicit knowledge on the underlying equations behind time series . although we restrict our attention to the analysis of numerical experiments here , it is also applicable to experimental data . the present paper is organized as follows . in section  2 , we introduce a formulation of the kernel cca . sections  3 and 4 provide several successful applications of the kernel cca for treating gs . in section  5 , it is demonstrated that the nonstationary change of the coupling strength between two chaotic systems can be captured by the proposed approach . in section  6 , we discuss an issue on the choice of kernel parameters . our conclusion is given in section  7 . canonical correlation analysis ( cca ) is a multivariate analysis method to find a pair of vectors that maximize the correlation coefficient between projections of signals onto them  @xcite . cca is useful for detecting a linear relation between a pair of multidimensional data sets , but can not deal with the case like gs where there is a high nonlinear relation between two data sets . extending the ability of cca for analyzing data with nonlinearity , kernelization of cca has been proposed by several researchers  @xcite . an intuitive formulation of the kernel cca is as follows . for a pair of multivariate variables @xmath8 with @xmath9 and @xmath10 , the kernel cca seeks a pair of nonlinear scalar functions @xmath11 and @xmath12 such that an estimator of the correlation coefficient @xmath13 between @xmath14 and @xmath15 is maximized . methods based on the expression  ( [ rho_f ] ) of the generalized correlation coefficient have a long history  @xcite . however the use of the kernel approach gives a notable progress in the subject , with a remarkably simple and flexible way of treating  ( [ rho_f ] ) . the essence of the kernel methods is an assumption that nonlinear functions @xmath16 and @xmath17 are well approximated by linear combinations of _ kernels _ on data points @xmath18 , @xmath19 where @xmath20 denotes a training data set of @xmath8 . examples of kernels are @xmath21 , @xmath22 , and @xmath23  @xcite . a rich representation of nonlinear functions is allowed by suitably chosen weight coefficients @xmath24 , @xmath25 and kernel parameters . at first sight , the dependence of the expressions for @xmath16 and @xmath17 on the given data @xmath20 brings some ad - hoc nature to the method . the assumption is , however , equivalent to the assumption that @xmath16 and @xmath17 are minimum norm functions in a reproducing kernel hilbert space ( rkhs ) . it is proved by invoking _ representer theorem _ in the theory of rkhs  @xcite . substituting ( [ fg ] ) into ( [ rho_f ] ) , and replacing covariance @xmath26 and variances @xmath27 with empirical averages over the data set @xmath20 , the kernel cca reduces to the maximization of @xmath28 where @xmath29 , @xmath30 , and @xmath31 are the _ gram matrices _ @xmath32 and @xmath33 defined with the given sample . the gram matrix contains relevant topological information of data points in the state space . for simplicity , we tentatively assume that the averages of @xmath34 and @xmath35 are zero . later we will show that subtraction of the averages from data can be done in an implicit way and does not cause any technical difficulty . a naive generalization of cca leads to the maximization of @xmath36 subject to @xmath37 . the corresponding lagrangian is @xmath38 where @xmath39 and @xmath40 are lagrange multipliers . taking derivatives of @xmath41 with respect to @xmath42 and @xmath43 , and @xmath44 from the constraint , the kernel cca is formulated as the following generalized eigenvalue problem : @xmath45 there is , however , a difficulty when the gram matrices @xmath46 and @xmath47 are invertible . by multiplying both sides of ( [ kcca_1 ] ) by @xmath48 from the left hand side , we have @xmath49 substituting ( [ invert_case_beta ] ) into ( [ invert_case_alpha ] ) , we obtain @xmath50 which gives a trivial solution @xmath51 . such a situation is not uncommon , because the gram matrix is invertible when a gaussian kernel is used and data points are distinct each other . thus the naive kernelization ( [ kcca_1 ] ) of cca does not provide useful information on the nonlinear dependence .    to overcome this problem we introduce small regularization terms in the denominator of the right hand side of ( [ rho_2 ] ) as @xmath52 where @xmath53 and @xmath54 are quadratic norms of @xmath16 and @xmath17 on the rkhs defined as @xmath55 and @xmath56 is its parameter . note that although the parameter @xmath57 is required for a nontrivial result , the precise value of @xmath56 is not important and @xmath17 to the optimal ones is proved when @xmath16 and @xmath17 belong to a rkhs , under the condition that @xmath58 with @xmath59 . fukumizu k , bach f r , and gretton a 2005 consistency of kernel canonical correlation analysis , _ ism research memorandum _ no.942 ] this insensitivity to the value of @xmath56 will be confirmed by a numerical experiment in the next section . as well as ( [ kcca_1 ] ) , the maximization of the numerator @xmath60 in ( [ rho_3 ] ) subject to @xmath61 reduces to the following generalized eigenvalue problem @xmath62 the first eigenvalue of ( [ kcca_2 ] ) gives the maximal value @xmath63 of @xmath64 in ( [ rho_3 ] ) . @xmath63 is called the _ canonical correlation coefficient _ , and the variables @xmath65 and @xmath66 transformed by @xmath16 and @xmath17 are called the _ canonical variates _ of the kernel cca . so far , the averages of @xmath34 and @xmath35 are set to zero . in fact , it is not necessary to subtract the averages explicitly , because the subtraction of the means is equivalent to the replacement of the gram matrices @xmath67 with @xmath68 where @xmath69 is a @xmath70-dimensional vector such that each component equals to the unity  @xcite . in the examples of this paper , the data @xmath71 and @xmath72 are normalized within an unit interval @xmath73 $ ] before applying the kernel cca . the kernel methods such as the kernel cca have been mainly applied to problems of pattern recognition  @xcite and bioinformatics  @xcite . in this paper , we propose that the kernel cca can also be a powerful tool for analyzing nonlinear dynamics . et al . _ @xcite proposed a method for analyzing dynamical system by the maximization the expression  ( [ rho_f ] ) . it is based on the _ alternating conditional expectation _ ( ace ) algorithm  @xcite , where the transformations @xmath16 and @xmath17 in  ( [ rho_f ] ) are restricted to a linear combination of univariate functions such as @xmath74 while the kernel cca allows any nonlinear function @xmath75 of @xmath76 variables in a rkhs . this difference makes the proposed method much more powerful than the one based on ace , especially in the cases where nonlinear dependence between variables @xmath77 is essential . unlike ace , the kernel cca does not need any iterative procedure and it is enough to solve a generalized eigenvalue problem by standard numerical algebra , only once for each set of data and parameters .    in this section , we study two unidirectionally coupled hnon maps  @xcite as an illustrative example and demonstrate the ability of the proposed approach . +      two unidirectionally coupled hnon maps are described by the following difference equations : @xmath78 where @xmath79 and @xmath80 are state variables at time @xmath81 , and @xmath82 denotes the coupling strength between two systems @xmath3 and @xmath4 . figures  [ revise_figure1 ] ( a ) and ( b ) show the projections of the strange attractors of ( [ henon_1 ] ) and ( [ henon_2 ] ) with @xmath83 and @xmath84 onto the @xmath85 plane , respectively . when the coupling strength @xmath82 is large , a complicated driver - response relation with high nonlinearity is formed between two different chaotic dynamics . it can be identified in figure  [ revise_figure1 ] ( b ) with a visual inspection , but it is not easy to represent it by naive statistical tools such as correlation coefficients .     plane with @xmath86 in ( a ) and @xmath87 in ( b ) . for both of ( a ) and ( b ) , an orbit with length @xmath88 is used for plotting . ( c , d ) scatter plots of the canonical variates of the kernel cca with @xmath86 in ( c ) and @xmath89 in ( d).,width=453 ]    we apply the kernel cca to the cases shown in figures  [ revise_figure1 ] ( a ) and ( b ) . here , we employ a gaussian kernel with @xmath90 and @xmath56 is set to @xmath91 . we prepare an orbit with length @xmath92 as training data . as already remarked , each of variables @xmath93 , @xmath94 , @xmath95 , and @xmath96 is normalized within the unit interval @xmath73 $ ] . figures  [ revise_figure1 ] ( c ) and ( d ) illustrate results of the kernel cca . figure  [ revise_figure1 ] ( d ) shows that gs is clearly identified as a cloud of points along the diagonal on the plane of the canonical variates @xmath97 and @xmath98 . on the other hand , when @xmath3 and @xmath4 are independent , the correlation between the canonical variates is very weak as shown in figure  [ revise_figure1 ] ( c ) . figure  [ revise_figure2 ] ( a ) shows the dependence of the canonical correlation coefficient @xmath99 on the coupling strength @xmath82 . the rapid increase of @xmath99 against @xmath82 in an interval @xmath100 is indicated . we will show that the kernel cca can be used with time - delay embedding scheme for time series . in time - delay embedding scheme , a pair of @xmath101-dimensional vectors @xmath102 and @xmath103 , where @xmath104 denotes the delay time , is used as a sample of training data . figure  [ revise_figure2 ] ( b ) shows the graph of @xmath99 vs. @xmath82 for two different values of the embedding dimension @xmath101 . shapes of the graphs shown in figure  [ revise_figure2 ] ( b ) does not change significantly compared to those shown in figure  [ revise_figure2 ] ( a ) . this result indicates that the kernel cca works for data obtained by using the time - delay embedding scheme .     of the kernel cca as a function of the coupling strength @xmath82 for three different values of @xmath105 with @xmath106 and @xmath92 . inset an enlargement of ( a ) . ( b ) the index @xmath99 vs. @xmath82 for the embedding dimension @xmath107 and @xmath108 with @xmath109 , and @xmath92 . the delay time @xmath104 for embedding is set to @xmath110 . ( c ) the time average of the synchronization error between the response system and its identical copy as a function of @xmath82 . ( d ) the maximal lyapunov exponent @xmath111 of the response system as a function of @xmath82.,width=453 ]      in @xcite , an identification of gs based on the occurrence of the complete synchronization between the response system and its identical copy is proposed . for the coupled hnon maps ( [ henon_1 ] ) and ( [ henon_2 ] ) , we investigate the synchronization errors @xmath112 between the system @xmath4 and its identical copy @xmath5 with variables @xmath113 . the result is shown in figure  [ revise_figure2 ] ( c ) . here , the index @xmath112 between @xmath4 and @xmath5 is measured by the average of @xmath114 time steps and plotted as a function of @xmath82 . the index @xmath112 decreases with increasing @xmath82 and becomes zero at @xmath115 . our results shown in figures  [ revise_figure2 ] ( a ) and ( b ) are consistent with this result . we also investigate the maximal lyapunov exponent @xmath111 of the state @xmath116 . the result is shown in figure  [ revise_figure2 ] ( d ) . the index @xmath111 is determined from the first eigenvalue of the product of the jacobian matrix of ( [ henon_2 ] ) with respect to the variables @xmath2 , and an orbit with length @xmath114 is used for its numerical evaluation . there is an interval where @xmath111 changes nonmonotonically against @xmath82 . such nonmonotonic change is also observed in the graph of @xmath99 vs. @xmath82 as shown in the inset of figure  [ revise_figure2 ] ( a ) . these results tell us that @xmath99 defined by the kernel cca can characterize subtle change as well as global tendency of gs . we consider how the performance of the kernel cca is influenced by the introduction of the observational noise and the change of the size of training data . first we consider the influence of noise . in numerical simulation , for each variable of ( [ henon_1 ] ) and ( [ henon_2 ] ) normalized within the interval @xmath73 $ ] , gaussian random numbers with the mean zero and the standard deviation @xmath17 are added as the observational noise . results are shown in figure  [ revise_figure3 ] . for a low noise level ( @xmath117 ) , the graph of @xmath99 vs. @xmath82 is almost same as that for the noise free case . for higher noise levels ( @xmath118 ) , there is a moderate decrease of @xmath99 in the whole interval of @xmath82 , however , the global tendency of the graph of @xmath99 vs. @xmath82 is not lost . the proposed approach is fairly robust against noise except for extremely high noise level ( @xmath119 ) .     as a function of @xmath82 with @xmath120 , and @xmath121 . the average and the standard deviation of @xmath99 over @xmath122 realizations are plotted as symbols and vertical bars , respectively.,width=453 ]    second , we study the influence of the size of training data . figure  [ revise_figure4 ] ( a ) shows the average of @xmath99 over @xmath122 realizations as a function of @xmath82 for several different values of the size @xmath70 . with the exception of @xmath123 , the graph of @xmath99 vs. @xmath82 does not depend much on @xmath70 . the result shows that the kernel cca works even with relatively small size of training data . for the cases of @xmath123 and @xmath89 , the graphs of @xmath99 vs. @xmath70 are shown in figure  [ revise_figure4 ] ( b ) . here , the vertical bars denote the corresponding standard deviation . the average of @xmath99 for @xmath123 increases monotonically with decreasing @xmath70 whereas the one for @xmath89 does not almost change against @xmath70 . the result with @xmath124 suggests that a proper choice of the kernel parameter for a given data is important for obtaining correct results . this issue is discussed in section 6 . on @xmath82 for three different values of the size @xmath70 of training data with @xmath120 . ( b ) dependence of the index @xmath99 on @xmath70 for @xmath124 and @xmath89 with @xmath120 . the average of @xmath99 over @xmath122 realizations are plotted as symbols in ( a ) , and in addition to the average , the standard deviation is also plotted as vertical bars in ( b ) . , width=453 ]      in order to investigate the ability of the kernel cca to nonlinear dependence between two systems ( [ henon_1 ] ) and ( [ henon_2 ] ) , we use a method of surrogation@xcite . multivariate surrogate data are generated as follows : first , the fourier transform of the time series is calculated for each of variables , then the _ common _ random numbers are added to the phase variables , and finally the inverse fourier transform is applied . the resulting multivariate time series have the same power spectra and cross spectra as those of the original time series . by changing random numbers added to the phases , an arbitrary number of different time series which preserve the linear properties of the original is obtained . see papers @xcite for technical details .    in numerical simulation , 19 realizations of the surrogate data for the time series @xmath125 of ( [ henon_1 ] ) and ( [ henon_2 ] ) are prepared by using the tisean package@xcite . we take @xmath107 and @xmath126 for time - delay embedding .    in figure  [ revise_figure5 ] ( a ) , the index @xmath99 defined by the kernel cca for the original data and that for the surrogate data as functions of @xmath82 are shown . for the surrogate data , the average over 19 realizations is plotted as a function of @xmath82 , and the corresponding maximal and minimal values are also shown as the both edges of vertical bars . for both of the original and the surrogate data , the index @xmath99 increases with increasing @xmath82 . except for @xmath123 and @xmath127 , however , the value of @xmath99 for the original data is significantly higher than that for the surrogate data . for larger values of @xmath82 , @xmath99 for the surrogate data increases monotonically , and @xmath128 as @xmath129 . this coincides with the fact that the attractor of the systems ( [ henon_1 ] ) and ( [ henon_2 ] ) is located around the plane @xmath130 and @xmath131 and the relation between two systems becomes almost linear one . the results suggest that the kernel cca is sensitive to nonlinearity of the dependence between two systems . we also investigate the performance of the linear cca in the same way and results are shown in figure  [ revise_figure5 ] ( b ) . in this case , the difference between the maximal canonical correlation coefficients @xmath132 of the linear cca for the original data and the one for the surrogate data is not significant for any value of @xmath82 . this indicates that the linear cca can detect only the linear dependence between two systems .     of the kernel cca ( a ) and @xmath132 of the linear cca ( b ) as functions of @xmath82 for the original data and the surrogate data . for the surrogate data , the average , and the maximal and minival values of @xmath99 over @xmath133 realizations are plotted as symbols , and both edges of vertical bars , respectively . parameters are @xmath134 , and @xmath121 in ( a ) and @xmath135 and @xmath121 in ( b).,width=453 ]      as mentioned in the preceding section , the regularization terms are required for nontrivial results . figure  [ revise_figure6 ] shows the dependence of @xmath99 on the regularization parameter @xmath56 for three different values of @xmath105 . although @xmath136 for too small @xmath56 , the value of @xmath99 decreases gradually and does not depend on the precise value of @xmath56 . as a function of the regularization parameter @xmath56 for three different values of @xmath105 . @xmath137 and @xmath138.,width=377 ] in order to illustrate the capability of the kernel cca in more complicated situations , we add the following three examples .      first , we consider gs in a lorenz system driven by a rssler system@xcite : @xmath139 the coupling term is introduced in the first equation of ( [ lore ] ) , and @xmath82 is its strength . we confirm that there is a sharp transition of gs at @xmath140 by investigating the long time average of the synchronization error between the system ( [ lore ] ) and its identical copy driven by the common signal @xmath141 of the system ( [ ross2 ] ) . there coexist two attractors in the state space after the transition of gs and we choose one of them . figures  [ revise_figure7 ] ( a )  ( c ) show the projections of the strange attractor of the system ( [ ross2 ] ) and ( [ lore ] ) with @xmath142 onto the planes of @xmath85 , @xmath143 , and @xmath144 respectively . for each of pairs shown in figures  [ revise_figure7 ] ( a )  ( c ) , the time series is embedded as points in a @xmath101-dimensional state space , and we apply the kernel cca . we set the embedding dimension @xmath108 and the delay time @xmath145 . the size of training data is @xmath121 . results are shown in figure  [ revise_figure8 ] ( a ) . all indices @xmath99 shown in figure  [ revise_figure8 ] ( a ) take large values for @xmath146 . the value @xmath147 agrees with the transition point of gs . in figure  [ revise_figure8 ] ( b ) , results of the linear cca applied to the same data are also shown . the indices @xmath132 of the linear cca for the cases of @xmath148 vs. @xmath149 and @xmath150 vs. @xmath151 take large values after the gs transition as well as the indices @xmath99 obtained by the kernel cca . for the case of @xmath152 vs. @xmath153 , however , @xmath132 of the linear cca is smaller than @xmath99 of the kernel cca . this result suggests that the kernel cca outperforms the linear cca when the relation between two observed time series has high nonlinearity .     onto the planes of @xmath85 in ( a ) , @xmath154 in ( b ) , and @xmath144 in ( c ) . , width=453 ]     of the kernel cca as a function of the coupling strength @xmath82 for different pairs of variables with @xmath155 and @xmath106 . ( b ) the canonical correlation coefficient @xmath132 of the linear cca as a function of the coupling strength @xmath82 for different pairs of variables with @xmath156 , and @xmath145 . the results for the case where the original state variables @xmath157 and @xmath158 are used as training data are also shown.,width=453 ]      second , we analyze the following fitzhugh - nagumo ( fhn ) neuron model modulated by chaotic dynamics of the rssler model : @xmath159 where the term @xmath160 defines a chaotic inputs to the neuron , and @xmath161 is a parameter that controls the dominant time scale of the rssler dynamics . the system composed of ( [ ross ] ) and ( [ fhn ] ) has been investigated from the viewpoint of the problem whether the information on the input signal can be decoded from the output interspike intervals ( isis ) generated by a neuron or not@xcite . here , we focus on the relation between the input chaotic stimulus and the output isis from the viewpoint of gs between two oscillators with different dynamics . we define the @xmath162-th isi as @xmath163 where @xmath164 is onset time of the @xmath162-th spike defined as the time when the variable @xmath95 makes upward crossing over some fixed threshold @xmath165 . the value of @xmath165 is set to @xmath166 here . we also define the chaotic stimulus associated with the @xmath162-th isi as @xmath167 , which is the value of @xmath93 in ( [ ross ] ) at @xmath164 . by using the delay embedding scheme , we transform the time series @xmath168 and @xmath169 into the state points @xmath170 and @xmath171 in @xmath172 and @xmath173-dimensional state spaces , respectively . we set @xmath174 here and apply the kernel cca to these data sets . in figures [ revise_figure9 ] ( a1)(a3 ) , significant nonlinear dependences between the input chaos and the output isi are observed in scatter plots of @xmath175 . figures  [ revise_figure9 ] ( b1)(b3 ) also show scatter plots of the canonical variates of the kernel cca associated with figures  [ revise_figure9 ] ( a1)(a3 ) . it is easy to see that the correlation between canonical variates @xmath97 and @xmath98 shown in figures  [ revise_figure9 ] ( b1)(b3 ) correspond to the complexity of nonlinear dependence between @xmath176 and @xmath177 shown in figures  [ revise_figure9 ] ( a1)(a3 ) . the relation between @xmath176 and @xmath178 changes according to the value of the control parameter @xmath161 . in figure  [ revise_figure9 ] ( c ) , the canonical correlation coefficient @xmath99 is plotted as a function of @xmath161 , which visualizes the change of the input - output relation between two systems of ( [ ross ] ) and ( [ fhn ] ) . the index @xmath99 changes nonmonotonically with the increase of @xmath161 , and there is a regime around @xmath179 where the value of @xmath99 is large . in addition to gs , chaotic phase synchronization ( cps )  @xcite occurs between two systems in this regime  @xcite . the increase of @xmath99 in this regime can be attributed to the occurrence of cps .     and @xmath178 of ( [ ross ] ) and ( [ fhn ] ) with @xmath180 in ( a1 ) , @xmath181 in ( a2 ) , and @xmath182 in ( a3 ) . ( b1  b3 ) scatter plots of canonical variates @xmath97 and @xmath98 of the kernel cca with @xmath183 in ( b1 ) , @xmath181 in ( b2 ) , and @xmath182 in ( b3 ) . @xmath184 , @xmath185 , and @xmath92 . ( c ) the canonical correlation coefficient @xmath99 of the kernel cca as a function of @xmath161 with @xmath184 , @xmath185 , and @xmath92 . , width=453 ]      the notion of gs is not resricted to unidirectinally coupled systems . in  @xcite , the occurrence of gs for bidirectionally coupled systems is also discussed . as an example of gs in bidirectinally coupled systems , we consider the following coupled lorenz - rssler systems  @xcite : @xmath186 a mutual interaction between a lorenz system ( [ bi_lore ] ) and a rssler system ( [ bi_ross ] ) is introduced as diffusion terms in the first equations of both ( [ bi_lore ] ) and ( [ bi_ross ] ) , and @xmath82 is its strength . in  @xcite , zheng et al . try to define gs in bidirectionally coupled systems by considering identical copies @xmath187 and @xmath5 of the systems @xmath3 and @xmath4 . using the approach in @xcite , two transitions of gs are found in the systems ( [ bi_lore ] ) and ( [ bi_ross ] ) . at @xmath188 , the rssler system @xmath4 is entrained by the lorenz system @xmath3 in the sense that orbits of the systems @xmath4 and @xmath5 completely coincide with each other by receiving the common signal of the system @xmath3 . with the further increase of @xmath82 , the system @xmath3 is also entrained by the system @xmath4 at @xmath189 which means that the complete synchronization between @xmath3 and @xmath187 also occurs . with @xmath190 in ( a1 ) , @xmath191 in ( b1 ) @xmath192 in ( c1 ) . ( a2 , b2 , c2 ) scatter plots of the canonical covariates of the kernel cca with @xmath190 in ( a1 ) , @xmath191 in ( b1 ) @xmath192 in ( c1 ) . parameters are @xmath193 , and @xmath92.,width=453 ]    we apply the kernel cca to the systems ( [ bi_lore ] ) and ( [ bi_ross ] ) and results are shown in figures  [ amend_revise1 ] and [ amend_revise2 ] . figure  [ amend_revise1 ] ( a1 ) , ( b1 ) , and ( c1 ) show the projections of attractors of the systems ( [ bi_lore ] ) and ( [ bi_ross ] ) with @xmath194 , @xmath195 , and @xmath196 onto the @xmath85 planes , respectively . figures  [ amend_revise1 ] ( a2 ) , ( b2 ) and ( c2 ) also show the corresponding results of the kernel cca . here , we use a gaussian kernel with @xmath197 . the state variables @xmath198 where @xmath199 are used as the training data set . it is observed that the state of gs is clearly identified as the high linear correlation between canonical variates @xmath97 and @xmath98 of the kernel cca . application of the kernel cca to the bidirectionally coupled systems is straightforward while the approach of @xcite requires rather subtle procedures . figure  [ amend_revise2 ] shows the index @xmath63 as a function of the coupling strength @xmath82 . there is a rapid increase of @xmath63 against @xmath82 in an interval @xmath200 . this result is consistent with the first transition of gs defined in @xcite . the second transntion of gs at @xmath189 defined in @xcite is not seen in the graph of @xmath63 vs. @xmath82 . a reason is that the value of @xmath63 already becomes nearly one at @xmath201 . another reason is that by definition , the proposed approach based on the kernel cca is insensitive to the directionality of synchronization . it will be interesting to study modifications of the approach which can deal with the directionality of synchronization .     of the kernel cca as a function of the coupling strength @xmath82 with @xmath193 , and @xmath92.,width=453 ] so far , we have focused on the dependence of the first eigenvalue @xmath99 on the coupling strength @xmath82 . we turn our attention to the eigenvector @xmath202 in ( [ kcca_2 ] ) and investigate changes of the structure of the dynamical system with the kernel cca . as an illustration , let us consider the problem of extracting nonstationary changes of the coupling strength @xmath82 from time series data generated by the coupled hnon maps . an example of such changes is shown in figure  [ revise_figure10 ] ( b ) . when the value of @xmath82 varies from @xmath203 , an orbit leaves from the synchronization manifold @xmath7 with @xmath203 , and is attracted again to it when the value of @xmath82 is restored to @xmath203 . as shown in figure  [ revise_figure10 ] ( a ) , the time series of the original variables @xmath85 does not tell us whether the orbit is lying on @xmath7 or not at a given time @xmath81 . by using the kernel cca , a deviation of the orbit from @xmath7 will be detected as a large value of @xmath204 , where @xmath205 and @xmath206 are nonlinear transformations determined by the first eigenvector @xmath207 . numerical experiments are performed with @xmath208 in two different conditions , and results are shown in figures  [ revise_figure10 ] ( c ) and ( d ) , respectively . figure  [ revise_figure10 ] ( c ) shows the case where @xmath205 and @xmath206 are estimated from an orbit with @xmath209 , which is prepared separately from the orbit to be analyzed . in figure  [ revise_figure10 ] ( c ) , time intervals where the values of @xmath82 are different from @xmath203 is clearly detected as successive bursts in the time series of @xmath204 . in figure  [ revise_figure10 ] ( d ) , we show the case where the orbit to be analyzed is also used for estimating @xmath205 and @xmath206 . we see that even when estimation of nonlinear transformations is affected by  noise " from the points not lying on @xmath7 , successive bursts are still observed in the time series of @xmath204 except for the one in a time interval @xmath210 .     and @xmath95 of the coupled hnon maps when the value of @xmath82 changes temporally as shown in ( b ) . ( c , d ) time series of the difference of the canonical variates of the kernel cca with @xmath211 and @xmath185.,width=453 ] in the preceding section , we set the value of the kernel parameter such as the width @xmath105 of a gaussian kernel in an ad - hoc manner . if @xmath105 is too small , the nonlinear transformations @xmath205 and @xmath206 in ( [ fg ] ) can not interpolate between data points of the training sample . contrary , if @xmath105 is too large , ( [ fg ] ) can not represent a highly nonlinear structure such as the synchronization manifold @xmath7 of gs . thus a proper choice of the kernel parameter is crucial to obtain the good performance of the method . in this section , we discuss how the kernel parameter can be suitably chosen from a given data set . a naive way of choosing @xmath105 is to find the value of @xmath105 that maximizes @xmath63 . the dotted line with open circles in figure  [ revise_figure11 ] shows the graph of an index @xmath63 as a function of @xmath105 when two hnon maps ( [ henon_1 ] ) and ( [ henon_2 ] ) are uncoupled ( @xmath123 ) . here , an orbit of ( [ henon_1 ] ) and ( [ henon_2 ] ) with length @xmath121 is used as a set of training data , and the average and the standard deviation over @xmath122 realizations are plotted as symbols and vertical bars , respectively . the index @xmath63 increases monotonically with the decrease of @xmath105 , and @xmath212 is attained in the limit of @xmath213 , while there is no interaction between two systems @xmath3 and @xmath4 . this monotonic tendency of @xmath63 does not determine an optimal value of the kernel parameter @xmath105 . as a way to overcome this difficulty , the following procedure is proposed . first , we set aside the data @xmath214 for assessing the performance of the kernel cca ( we call this  test " data ) separately from the data used for training the kernel cca . then , we estimate the nonlinear transformations @xmath205 and @xmath206 from the training data , and calculate the correlation coefficient @xmath215 between the variables @xmath216 and @xmath217 defined as @xmath218 where @xmath219 is the empirical average over @xmath220 samples . this strategy for assessing the performance of the estimated model with new data is regarded as a version of _ cross - validation _ ( cv )  @xcite . first we check that spurious detection of synchronization can be avoided by the procedure based on cv . the solid line with filled squares in figure  [ revise_figure11 ] shows the dependence of an index @xmath215 on @xmath105 with @xmath123 . when there is no interaction between two systems , the value of @xmath215 is nearly equal to zero for any @xmath105 . corresponding results using time - delay embedding are shown in figure  [ revise_figure12 ] . in figure [ revise_figure12 ] , the dependencies of @xmath63 and @xmath215 on @xmath105 with @xmath124 are plotted for several different values of the embedding dimension @xmath101 . again , the value of @xmath215 is almost equal to zero for any @xmath105 and @xmath101 , whereas the value of @xmath63 as a function of @xmath105 increases monotonically with the increase of @xmath101 . this result suggests that the procedure based on cv effectively erases spurious detection even when data is embedded in a high - dimensional state space .    ) . the dotted line with open circles shows @xmath63 vs. @xmath105 , and the solid line with filled squares shows @xmath215 vs @xmath105 . we prepare 20 data sets with size @xmath121 for training the kernel cca , and a test data with size @xmath221 for cross - validation . the average over @xmath122 realizations are plotted as symbols . the standard deviation is also plotted as vertical bars , but the ones with small values of errors are hidden by symbols . @xmath56 is set to @xmath222 . , width=453 ]     for four different values of the embedding dimension @xmath101 . we take the delay time @xmath126 . the condition for numerical experiments is the same as figure  [ revise_figure11].,width=453 ]    next we show that the cross validation procedure is useful for choosing optimal values of @xmath105 . figures  [ revise_figure13 ] shows the values of @xmath63 and @xmath215 as functions of @xmath105 for the coupled hnon maps ( [ henon_1 ] ) and ( [ henon_2 ] ) with @xmath87 . the condition for the numerical experiment is the same as @xmath124 . for all of graphs , the value of @xmath215 with the dotted line takes its maximum at a nonzero value of @xmath105 , whereas @xmath63 increases monotonically with the decrease of @xmath105 . this result indicates that we can choose this value as an optimal width of a gaussian kernel .    . the embedding dimension @xmath223 in ( a ) , @xmath107 in ( b ) , @xmath108 in ( c ) , and @xmath224 in ( d ) , and the delay time @xmath126 . the condition for numerical experiments is the same as figure  [ revise_figure11].,width=453 ] in conclusion , we have proposed a new approach for analyzing gs in a unified framework of a kernel method . we have tested the proposed approach by applying it to several examples exhibiting gs , and demonstrated that the canonical correlation coefficient of the kernel cca is a suitable index for the characterization of gs . in addition , it has been shown that nonstationary changes of the coupling are detected from the time series by the difference between canonical variates of the kernel cca . it has been also discussed how the parameter of the kernel function can be suitably chosen from data by the procedure of cross - validation . our experiments show that a method based on cv gives promising results in optimizing the parameter @xmath105 . the cross - validation procedure is also useful to circumvent spurious detection of gs by overfitting . the approach based on the kernel cca provides not only an index for measuring nonlinear interdependence . it also provides global nonlinear coordinates , and these coordinates allow a representation of the interaction between the dynamical systems under investigation . note that the linear cca also provides global coordinates , but it can not discriminate between linear and nonlinear relation . in this respect , it goes beyond conventional methods of analyzing gs . our attempts open a new possibility of the kernel methods for analyzing complex dynamics observed in nonlinear systems . we thank s.  akaho for stimulating discussions on the kernel methods , and k.  fukumizu for useful comments and informing us his results prior to publication . we also thank y.  hirata for useful suggestions on surrogate data analysis . 99 pikovsky a , rosenblum m and kurths j 2002 _ synchronization , a universal concept in nonlinear sciences _ ( cambridge : cambridge university press ) fujisaka h and yamada t 1983 stability theory of synchronized motion in coupled - oscillator system _ prog . phys . _ * 69 * 3247                                      melzer t , reiter m , and bischof h 2001 nonlinear feature extraction using generalized canonical correlation analysis in _ proceedings of the international conference on artficial neural networks ( icann ) _ 353360 ( london : springer - verlag )
Abstract:

Answer: a unified framework for analyzing generalized synchronization in coupled chaotic systems from data is proposed . the key of the proposed approach is the use of the _ kernel methods _ recently developed in the field of machine learning . several successful applications are presented , which show the capability of the kernel - based approach for detecting generalized synchronization , and dynamical change of the coupling strength between two chaotic systems can be captured by the proposed approach . it is also discussed how the kernel parameter is suitably chosen from data .


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: resonant inelastic x - ray scattering ( rixs ) is a technique that matures rapidly due to the recent increase in brilliance of the new generation synchrotron x - ray sources , where high flux photon beams with energies that are tunable to resonant edges are now becoming widely available  @xcite . the probability for x - rays to be scattered from a solid state system can be enhanced by orders of magnitude when the energy of the incoming photons is in the vicinity of an electronic eigenmode a resonant edge of the system . rixs experiments are performed on e.g. the k - edges of transition metal ions , where the frequency of the x - rays is tuned to match the energy of an atomic @xmath1-@xmath2 transition , which is around 5 - 10 kev  @xcite . at this resonant energy a @xmath1 electron from the inner atomic core is excited into an empty @xmath2 state , see fig . 1 . it is a well - known fact that the @xmath1 core - hole that is created has an ultrashort lifetime , of the order of femtoseconds . the reason is that the core - hole has a very high energy and is prone to decay via all sorts of radiative and non - radiative processes , severely cutting down the efficiency of rixs . in theoretical treatments of rixs this life - time effect is normally introduced as a core - hole broadening and disregarded from that point on .    in a previous study  @xcite , however , we have shown that from the theory perspective there is a great advantage of the very short lifetime of the core - hole . the ultrashort lifetime implies that for the other electrons in the system particularly for the slow ones that are close to the fermi - energy the core - hole potential is almost an instantaneous delta - function in time . this allows for a systematic expansion of the scattering cross section in terms of the lifetime , for which we present a detailed derivation and various generalizations in this paper . we shall see that the most important consequence of the ultrashort core - hole lifetime is that for indirect rixs the effective scattering cross section is proportional to the charge structurefactor @xmath0 and the longitudinal spin structure factor that is associated with it . the indirect rixs process is shown schematically in fig . 1 . in transition metal systems the photo - electron is promoted from a @xmath1 core - orbital to empty @xmath2 states that are far ( 10 - 20 ev ) above the fermi - level . so the x - rays do not cause direct transitions of the @xmath1 electron into the lowest @xmath3-like conduction bands of the system . still rixs is sensitive to excitations of electrons near the fermi - level . the coulomb potential of @xmath1 core - hole causes e.g. very low energy electron - hole excitations in the valence / conduction band : the core - hole potential is screened by the valence electrons . when the excited @xmath2-electron recombines with the @xmath1 core - hole and the outgoing photon is emitted , the system can therefore be left behind in an excited final state . experimentally the momentum @xmath4 and energy @xmath5 of the elementary excitation is determined from the difference in energy and momentum between incoming and outgoing photons . since the excitations are caused by the core - hole , we refer to this scattering mechanism as _ indirect _ resonant inelastic x - ray scattering ( rixs ) .    at present energy resolutions of about 100 mev can be reached . in the near feature it seems experimentally feasible for rixs to become sensitive to the low energy excitations of the solid , where excitation energies are of the order of room temperature . recently it has been shown that also magnetic excitations , magnons , can be measured in rixs  @xcite . other interesting low - lying electronic excitations that potentially can be probes by rixs are , for example , collective features such as plasmons , orbitons , excitons , but also single - particle - like continua related to the band structure . rixs provides a new tool to study these elementary excitations .    for the interpretation of spectroscopic data , it is very important to express the scattering cross section for a technique in terms of physical correlation functions . in this paper , we derive in detail the dynamical correlation function that is measured in indirect resonant inelastic x - ray scattering . for local core - hole potentials and ultrashort lifetimes , the dynamical correlation function turns out to be a linear combination of the charge density and longitudinal spin density response function . we show that for a single band system the actual linear combination that is measured depends on the energy of the incoming photons and we determine the precise energy dependence of its coefficients . a sum - rule is derived and we generalize these results to the case of finite temperature and for multiband systems . the kramers - heisenberg formula  @xcite for the resonant x - ray scattering cross section at finite temperature is @xmath6 where @xmath7 and @xmath8 denote the final and initial state of the system , respectively . the sum is over all final states and the brackets denote the statistical average over initial states @xmath8 for a temperature @xmath9 . the momentum and energy of the incoming / outgoing photons is @xmath10 and @xmath11 and the loss energy @xmath12 is equal to the energy difference between the final and initial state @xmath13 . in the following we will take the groundstate energy of our system as reference energy : @xmath14 . the scattering amplitude @xmath15 is given by @xmath16 where @xmath17 the resonant energy , @xmath18 denotes the intermediate states and @xmath19 the ( dimensionless ) dipole operator that describes the excitation from initial to intermediate state and the de - excitation from intermediate to final state . the energy of the incoming x - rays with respect to the resonant energy is @xmath20 ( this energy can thus either be negative or positive : @xmath21 ) and @xmath22 is the energy of intermediate state @xmath23 with respect to the resonance energy .    in the intermediate state a core - hole and a photo - excited electron are present . when we take the coulomb interaction between the intermediate state core - hole and the valence band electrons into account , we obtain a finite inelastic scattering amplitude . in that case there is a non - zero probability that an electron - hole excitation is present in the final state , see fig.1 . the intermediate state , however , is not a steady state . the highly energetic @xmath1 core - hole quickly decays e.g. via auger processes and the core - hole life - time is very short . the heisenberg time - energy uncertainty relationship then implies that the core - hole energy has an appreciable uncertainty . this uncertainty appears in the formalism above as the core - hole energy broadening @xmath24 which is proportional to the inverse core - hole life - time , which is of the order of electron volts as the lifetime is ultrashort , of the order of femtoseconds . note that the life - time broadening only appears in the intermediate states and not in the final or initial states as these both have very long life times . this implies that the core - hole broadening does not present an intrinsic limit to the experimental resolution of rixs : the loss energy @xmath5 is completely determined by kinematics .    when the incoming energy of the x - rays is equal to a resonant energy of the system @xmath25 and we see from eqs . ( [ eq : kramers],[eq : amplitude ] ) that the resonant enhancement of the x - ray scattering cross section is @xmath26 , which is @xmath27 for a transition metal k - edge  @xcite .    in a resonant scattering process , the measured system is generally strongly perturbed . formally this is clear from the kramers - heisenberg formula  ( [ eq : kramers ] ) , in which both the energy and the wavefunction of the intermediate state where a potentially strongly perturbing core - hole is present appear . this is in contrast with canonical optical / electron energy loss experiments , where the probing photon / electron presents a weak perturbation to the system that is to be measured .    to calculate rixs amplitudes , one possibility is to numerically evaluate the kramers - heisenberg expression . to do so , all initial , intermediate and final state energies and wavefunctions need to be known exactly , so that in practice a direct evaluation is only possible for systems that , for example , consist of a small cluster of atoms  @xcite . in this paper , however , we show that under the appropriate conditions we can integrate out the intermediate states from the kramers - heisenberg expression . after doing so , we can directly relate rixs amplitudes to linear charge and spin response functions of the unperturbed system . for non - resonant scattering , one is familiar with the situation that the scattering intensity is proportional to a linear response function , but for a resonant scattering experiment this is a quite unexpected result . let us proceed by formally expanding the scattering amplitude in a power series @xmath28 where we introduced the matrix elements @xmath29 the formal radius of convergence of this power series is given by @xmath30 , so that the series is obviously convergent when the incoming x - ray energy is e.g far enough below the resonance , i.e. when @xmath31 . but also at resonance , when @xmath32 the series is convergent for intermediate energies that are smaller than the core - hole broadening @xmath24 . thus this expansion is controlled for ultrashort core - hole lifetimes , which implies that @xmath24 is large . in the following we will be performing re - summations of this series . we denote the denominator of the expansion parameter @xmath33 by the complex number @xmath34 , so that @xmath35 where @xmath36 is the hamiltonian in the intermediate state . we thus obtain the following series expansion for the resonant cross section @xmath37 we will first calculate the resonant x - ray cross section at zero temperature in the case where the valence and conduction electrons are effectively described by a single band of spinless fermions : spin , and orbital degrees of freedom of the valence electron system are suppressed . physically this situation can be realized in a fully saturated ferromagnet . the final and initial states of the system are determined by a hamiltonian @xmath38 that describes the electrons around the fermi - level . the generic form of the full many - body hamiltonian is @xmath39 where @xmath40 and @xmath41 denote lattice sites with lattice vectors @xmath42 and @xmath43 . note that the sum is over each pair @xmath44 once , with @xmath44 ranging from 1 to @xmath45 , where @xmath45 is the number of sites in the system . the hopping amplitudes of the valence electrons are denoted by @xmath46 and the @xmath47-operators annihilate / create such electrons . the coulomb interaction between valence electrons is @xmath48 , as the coulomb interaction only depends on the relative distance between two particles . the intermediate states are eigenstates of the hamiltonian @xmath49 , where @xmath50 accounts for the coulomb coupling between the intermediate state core - hole and the valence electrons : @xmath51 where @xmath52 creates a core - hole on site @xmath40 . we assume that the core - hole is fully localized and has no dispersion . we will see shortly that this leads to major simplifications in the theoretical treatment of indirect rixs . the core - hole  valence electron interaction is attractive : @xmath53 . the dipole operators are given by @xmath54 where @xmath55 creates a photo - excited electron in a @xmath2 state and @xmath56 denotes the hermitian conjugate of both terms .      in order to calculate the cross section , we need to evaluate the operator @xmath57 in equation ( [ eq : mterms ] ) . a direct evaluation of this operator is complicated by the fact that @xmath58 \neq 0 $ ] . we therefore proceed by approximating @xmath59 with a series that contains the leading terms to the scattering cross section for both strong and weak core - hole potentials , if the lifetime is short . after that we will do a full re - summation of that series . this approximation is central to the results in this paper . expanding @xmath60 gives a series with @xmath61 terms : @xmath62 using @xmath63 , this series reduces to @xmath64 \hat{o}{\left|\textsc{i}\right>}.\end{aligned}\ ] ] using in addition that @xmath65 , we find @xmath66 \hat{o}{\left|\textsc{i}\right>}.      \label{eq : exp_lhs}\end{aligned}\ ] ]    for strong core - hole potentials , the leading term of @xmath59 is @xmath67 . corrections to this term contain at least one factor of @xmath38 and are therefore smaller by a factor of at least @xmath68 . for weak core - hole potentials , the term @xmath69 vanishes because @xmath70=0 $ ] . the leading term for this limit therefore is @xmath71 . correction terms contain at least two factors of @xmath50 , which make them at least a factor of @xmath72 smaller . now we consider the approximation @xmath73 it can be seen that the leading order terms for both strong @xmath74 and weak @xmath75 core - hole potentials are included in the sum . the other terms are included only for mathematical convenience lateron ; they can be neglected if we consider either limit . note that the @xmath76 term in eq . ( [ eq : h_expand ] ) is @xmath77 , so that it can be removed from the sum . after performing the same manipulations as above , we obtain @xmath78 \hat{o } { \left|\textsc{i}\right>}.      \label{eq : exp_rhs}\end{aligned}\ ] ] comparing eqs . ( [ eq : exp_lhs ] ) and ( [ eq : exp_rhs ] ) , it can be seen that the approximation ( [ eq : h_expand ] ) is exact in the limit of both strong and weak core - hole potentials . we can also represent the series expansion and its approximation graphically ( fig . [ fig : expand ] ) . when we expand @xmath79 , where @xmath80 and @xmath81 are non - commuting operators , each term in the series corresponds to a graph on the grid of fig . [ fig : expand].1 . each graph occurs only once and can be constructed by starting at the lower left corner of the grid and moving either to the right , representing an @xmath80 , or up , representing a @xmath81 . at the next vertex a new move ( right or up ) we perform this procedure @xmath82 times and in this way we can obtain @xmath61 distinct graphs , each corresponding to a term in the expansion of @xmath79 . for example moving @xmath82 times to the right represents the term @xmath83 and moving @xmath82 times up corresponds to @xmath84 , see fig . [ fig : expand].2 and  [ fig : expand].3 . all other terms in the series can be constructed by moving up and right a different number of times and in different order . as we consider a fixed value of @xmath82 ( @xmath85 in fig . [ fig : expand ] ) , all graphs must end on the diagonal of the triangle that forms the grid . in the series for @xmath86 ( @xmath87 and @xmath88 ) we have the simplification that terms ending with @xmath38 acting on the groundstate give zero . these terms can thus be removed from the expansion . the graphs for this expansion now live on a reduced grid where the horizontal grid - lines at the diagonal of the triangle are absent , see fig . [ fig : expand].5 : these represent all terms ending on @xmath80 .    , where @xmath89 and @xmath90 are two non - commuting operators . in this example @xmath85 . ]    in fig . [ fig : expand ] we also represent the approximate series of the r.h.s . of eq . ( [ eq : h_expand ] ) . graphically this sum corresponds to the set of graphs on the reduced grid of fig . [ fig : expand].5 , with either one kink ( fig . [ fig : expand].6 , [ fig : expand].7 , [ fig : expand].8 , [ fig : expand].9 ) or without kinks ( fig . [ fig : expand].2 , [ fig : expand].3 ) . thus , in our approximation in eq . ( [ eq : h_expand ] ) of the exact series for @xmath91 we neglect all graphs with two or more kinks ( fig . [ fig : expand].4 , [ fig : expand].10 , [ fig : expand].11 , [ fig : expand].12 ) . in the limit of either very @xmath80 or very large @xmath81 , the graphs that we neglect correspond to sub - leading order corrections . when @xmath80 is largest then the leading terms are , first , graph [ fig : expand].2 , which is however zero because it ends on @xmath80 . the leading term is therefore of the order @xmath92 and shown in graph [ fig : expand].6 . other higher order terms are shown in the graphs [ fig : expand].7 , [ fig : expand].8 , [ fig : expand].10 and [ fig : expand].11 . the last two graphs are neglected in our approximate expansion . in case @xmath81 is dominating , the leading term is @xmath93 , graph [ fig : expand].3 , and next to leading is graph [ fig : expand].9 , with @xmath94 . the highest order terms that are neglected in our approximate series are of the type shown in graph [ fig : expand].12 .      in order to obtain @xmath95 and from there the scattering amplitude @xmath15 and finally the scattering cross - section , in eq.([eq : exp_rhs ] ) we need to evaluate expressions of the kind @xmath96 in the initial state no core - hole is present : just one core - hole is created by the dipole operator . we therefore have that @xmath97 . inserting this in equation  ( [ eq : h_c^n , n-1 ] ) , we obtain @xmath98 and by recurrence @xmath99 ^n | \textsc{i } \rangle \label{eq : h_c^n}\end{aligned}\ ] ]    let us for the moment consider the strong core - hole potential limit and keep in the expansion eq.([eq : h_expand ] ) only the term @xmath100 . inserting the results above in eq.([eq : mterms ] ) , we find that @xmath101^l   | \textsc{i } \rangle , \label{eq : m_strong}\end{aligned}\ ] ] where the transfered momentum @xmath102 . the first important observation is that the term @xmath103 does not contribute to the inelastic x - ray scattering intensity because @xmath104 , which only contributes to the elastic scattering intensity at @xmath105 and other multiples of the reciprocal lattice vectors . from inspection of equation ( [ eq : m_l ] ) we see immediately that the @xmath103 term actually vanishes irrespective of the strength of the core - hole potential . this is of relevance when we consider the scattering cross section in the so - called `` fast - collision approximation ''  @xcite . this approximation corresponds to the limit where the core - hole life time broadening is the largest energy scale in system ( @xmath106 or , equivalently , @xmath107 \rightarrow -\infty$ ] ) . in this limit only the @xmath103 term contributes to the indirect rixs amplitude and the resonant inelastic signal vanishes . in any theoretical treatment of indirect resonant scattering one therefore needs to go beyond the fast - collision approximation . physically this vanishing of spectral weight is ultimately due to an interference effect . if we study a process in which we start from the initial state and reach a certain final state , we need to consider all different possible paths for this excitation  de - excitation process . when the core - hole broadening is very large we can reach the final state via any intermediate state and in order to obtain the scattering amplitude we thus add up coherently the contributions of all intermediate states . we then obtain @xmath108 . when the set of intermediate states that we sum over is complete ( which by definition is the case when @xmath109 ) , this leaves us with @xmath110 which is , because of the orthogonality of eigenstates , only non - zero when the initial and final state are equal hence only when the scattering is elastic . the second observation is that @xmath95 is a @xmath61-particle correlation function . if we measure far away from resonance , where @xmath111| \gg 0 $ ] , the scattering cross section is dominated by the @xmath112 , two - particle , response function . when the incoming photon energy approaches the resonance , gradually the four , six , eight etc . particle response functions add more and more spectral weight to the inelastic scattering amplitude . generally these multi - particle response functions interfere . we will show , however , that in the local core - hole approximation the multi - particle correlation functions in expansion ( [ eq : h_expand ] ) collapse onto the dynamic two - particle ( charge - charge ) and four - particle ( spin - spin ) correlation function . in hard x - ray electron spectroscopies one often makes the approximation that the core - hole potential is local . this corresponds to the widely used anderson impurity approximation in the theoretical analysis of e.g. x - ray absorption and photo - emission , introduced in refs . this approximation is reasonable as the coulomb potential is certainly largest on the atom where the core - hole is located .    in the present case , moreover , we can consider the potential generated by both the localized core - hole and photo - excited electron at the same time . as this exciton is a neutral object , its monopole contribution to the potential vanishes for distances larger than the exciton radius . the multi - polar contributions that we are left with in this case are generally small and drop off quickly with distance . we insert a local core - hole potential @xmath113 in our equations and aim to re - sum the approximate series expansion in eq.([eq : h_expand ] ) for arbitrary values of the local core - hole potential . we find from eq.([eq : h_c^n ] ) that @xmath114^n | \textsc{i }   \rangle \label{eq : h_c^nloc}\end{aligned}\ ] ] using that for fermions @xmath115^n = c^{\tiny \dagger}_i c^{\phantom{\tiny \dagger}}_i$ ] , we obtain for our spinless fermions @xmath116 the sum over @xmath117 can easily be performed : @xmath118 and we obtain @xmath119 using that @xmath120 is the density operator , we have to perform the sum over @xmath82 in equation  ( [ eq : a_fi ] ) . the @xmath103 term is zero , as we discussed above , so that the scattering amplitude is @xmath121 using @xmath122 we finally find that the indirect resonant inelastic scattering amplitude for spinless fermions is @xmath123 where the resonant enhancement factor is @xmath124 ^{-1}$ ] and @xmath125 . for spinless fermions with a local core - hole potential the scattering cross section thus turns out to be the density response function a two - particle correlation function with a resonant prefactor @xmath126 that depends on the loss energy @xmath5 , the resonant energy @xmath127 , on the distance from resonance @xmath128 $ ] ) , on the core - hole potential @xmath129 and on the core - hole life time broadening @xmath130)$ ] . we see that the resonant enhancement is largest when the energy of the incoming photons is either equal to the core - hole potential ( @xmath131 ) or to the loss energy ( @xmath132 ) , which one could refer to as a ` final - state resonance ' .    the density response function is related to the dielectric function @xmath133 and the dynamic structure factor @xmath0  @xcite , so that we obtain for the resonant scattering cross section @xmath134 \nonumber \\ & \propto & |p_1(\omega)|^2   \ s({\bf q},\omega),\end{aligned}\ ] ] for a fixed value of the core - hole potential @xmath129 . @xmath135 is the fourier transform of the coulomb potential . for weak core - hole potentials the total scattering intensity is proportional to @xmath136 and for strong core - hole potentials , where @xmath137 , the scattering intensity at resonance ( @xmath138 ) is to first order independent of the strength of the core - hole potential . far away from the edge , however , where @xmath139 , the scattering intensity is again proportional to @xmath136 , just as for weak core - hole potentials . integrating @xmath140 over all incoming photon energies , we obtain the integrated inelastic intensity at fixed loss energy @xmath5 and momentum @xmath141 @xmath142 it seems that the resonant enhancement factor of the integrated intensity has a maximum when the loss energy is equal to the core - hole potential . however , the core - hole potential is attractive and therefore lower than zero , and the loss energy @xmath5 is by definition greater than zero . so the integrated intensity is maximal at energy loss @xmath143 . in this section , we generalize the previous calculation to the case of finite temperature . the starting point is as before @xmath144 where @xmath145 is the partition function and @xmath146 . equation ( [ eq : kramerst ] ) represents the statistical average over all the initial states @xmath147 , where now the more general relation @xmath148 holds . we expand the scattering amplitude @xmath149 , using again the ultrashort life time of the core - hole as in eq . ( [ eq : a_fi ] ) . we are left with the evaluation of the operator @xmath150 . we proceed by expanding it in the following way : @xmath151 where we neglected the term @xmath69 , as it will not contribute to the inelastic scattering cross section . this approximation reproduces the correct leading order terms , which represent the strong and weak coupling case , respectively . moreover , it is a generalization of ( [ eq : h_expand ] ) , that takes into account that the initial state is no longer the ground state so that @xmath152 . in our graphical representation , with respect to the @xmath153 case , it corresponds to retain all the additional terms , having more than one kink , that start and finish with a horizontal step . in doing this , we are neglecting again the sub - leading order terms @xmath154 .    after inserting expansion ( [ eq : h_expandt ] ) in the expression ( [ eq : mterms ] ) for @xmath155 , we finally have to evaluate @xmath156 in the local core - hole approximation , we can resum this approximate series expansion . by using the results of eqs . ( [ eq : h_c^nloc ] ) , we obtain for spinless fermions @xmath157 by performing the sums over @xmath18 and @xmath117 @xmath158 and after summing over @xmath82 , we finally obtain @xmath159 this equation clearly shows that one of the main effects of finite temperature is to modify the resonant enhancement factor , nevertheless preserving the same structure for the scattering amplitude .    at this point we observe that at resonance @xmath160 , which is of the order of electron volts and thus several orders of magnitude larger than @xmath161 , even at high temperature . this allows us to approximate the prefactor in eq . ( [ eq : a_spinlesst ] ) as @xmath162 at the lowest order in @xmath163 , the prefactor is not modified by @xmath9 at all , hence we conclude that the major modifications to the cross section are induced by thermal averaging of the correlation function . after integrating over all the incoming photon energies , we get the following approximate expression for the thermal average of the inelastic intensity at loss energy @xmath5 and momentum @xmath141    @xmath164    in this expression the temperature dependence is entirely due to the temperature dependence of @xmath0 . the pre - factor is in leading order temperature independent . at finite temperatures energy gain scattering can occur : the photon can gain an energy of the order of @xmath165 from the system , which corresponds to negative energy loss . we generalize the calculation above to the situation where the electrons have an additional spin degree of freedom . in the hamiltonians  ( [ eq : h0],[eq : hcore ] ) we now include a spin index @xmath166 ( with @xmath167 or @xmath168 ) to the annihilation and creation operators : @xmath169 and @xmath170 and sum over these indices , taking into account that the hopping part of the hamiltonian is diagonal in the spin variables . in order to re - sum the series in equation   ( [ eq : h_expand ] ) we now need to evaluate expansions of the number operators of the kind @xmath171 . using @xmath172 for @xmath173 , we obtain @xmath174 + 2 p_2 ( \omega ) \rho^{\uparrow \downarrow}_{\bf q}| \textsc{i } \rangle,\ ] ] with @xmath175 and @xmath176 . we see that in the case that each site can only be occupied by at most one valence electron , this equation immediately reduces to eq . ( [ eq : a_spinless ] ) with @xmath177 . the two terms in the scattering amplitude can also be written in terms of density and spin operators . using @xmath178 , we obtain @xmath179 , where we introduce the longitudinal spin density correlation function @xmath180 . in terms of these correlation functions the scattering amplitude for spinfull fermions is @xmath181 \langle \textsc{f}| { \bf s}_{\bf q}^2 | \textsc{i } \rangle    + p_2(\omega )   \langle \textsc{f}|\rho_{\bf q }   | \textsc{i } \rangle,\ ] ] clearly the contributions to the scattering rate from the dynamic longitudinal spin correlation function and the density correlation function need to be treated on equal footing as they interfere  @xcite . moreover , the spin and charge correlation functions have different resonant enhancements , see fig . [ fig : resonance ] . for instance when @xmath182=u$ ] , the scattering amplitude is dominated by @xmath126 and hence by the longitudinal spin response function . at incident energies where @xmath182=2u$ ] , on the other hand , @xmath183 is resonating so that the contributions to the inelastic scattering amplitude of charge and spin are approximately equal . ( black ) and @xmath184 ( blue ) at fixed loss energy @xmath5 as a function of incoming photon energy @xmath185 for different values of the local core - hole potential @xmath186 let us consider systems with more than one band and take as an explicit example a transition metal with a @xmath3 and a @xmath187 band . the coulomb attraction between the @xmath1 core - hole and an electron in the @xmath3 state ( @xmath188 ) is much larger than the interaction with a @xmath187 electron ( @xmath189 ) . neglecting spin degrees of freedom we would naively expect that the indirect rixs response in the two - band system is simply the sum of the responses of the two individual electronic systems , with possible interference between the two scattering channels : we expect the scattering amplitude to be equal to @xmath190 however , already from the calculation for the spinfull fermions we know that the situation should be more complicated , as in that case the full response function is not just the sum of the two response functions for spinless fermions . the point is that when both a @xmath3 and @xmath187 electron screen the core - hole , the intermediate state is at a lower energy ( at @xmath191 ) compared to the situation where only a single @xmath192/@xmath193 electron screens the core - hole ( with a resonance at @xmath194/@xmath189 , respectively . ) in the situation that both electrons screen the core - hole , the resonance therefore appears at a different incoming photon energy .    according to eq.([eq : h_c^n ] ) , we now need to evaluate expressions of the sort @xmath195 for @xmath196 . after using the binomial theorem and summing the resulting series , we obtain @xmath197,\end{aligned}\ ] ] which leads to a scattering amplitude @xmath198 \langle \textsc{f}| \rho^{ds}_{\bf q }     | \textsc{i } \rangle,\end{aligned}\ ] ] where @xmath199 . on the basis of the ultrashort life - time of the core - hole in the intermediate state we presented a series expansion of the indirect resonant inelastic x - ray scattering amplitude , which is asymptotically exact for both small and large local core - hole potentials . this algebraic series is also given in a graphical representation . by re - summing the terms in the series , we find the dynamical charge and spin correlation functions that are measured in rixs . the resonant pre - factor is only weakly temperature dependent . we have also derived a sum - rule for the total scattering intensity and considered rixs in both single and multi - band systems . on the basis of our results , the charge and spin structure factor of e.g. hubbard - like model hamiltonians can be used to be directly compared to the experimental rixs spectra . we thank michel van veenendaal and john p. hill for stimulating our interest in the theory of indirect resonant inelastic x - ray scattering , for intense discussions , for critical reading of the manuscript and for their generous hospitality . we thank stephane grenier , young - june kim , philip platzman and george sawatzky for fruitful discussions . we gratefully acknowledge support from the argonne national laboratory theory institute , brookhaven national laboratory ( de - ac02 - 98ch10996 ) , the dutch science foundation fom .    for a review see : a. kotani and s. shin , rev . mod . phys . * 73 * , 203 ( 2001 ) . hasan , e.d . isaacs , z .- x . shen , l.l . miller , k. tsutsui , t. tohyama and s. maekawa , science * 288 * , 1811 ( 2000 ) . hasan , p.a . montano , e.d . isaacs , z .- x.shen , h. eisaki , s.k . sinha , z. islam , n. motoyama and s. uchida , phys . lett * 88 * 177403 ( 2002 ) . kim , j.p . hill , c.a . burns , s. wakimoto , r.j . birgeneau , d. casa , t. gog and c.t . venkataraman , phys . * 89 * , 177003 ( 2002 ) . j.p . hill , c .- c . kao , w.a.l . caliebe , m. matsubara , a. kotani , j.l . peng and r.l . greene , phys . . lett . * 80 * , 4967 ( 1998 ) . k. hmlinen , j.p . hill , s. huotari , c .- c . kao , l.e . berman , a. kotani , t.id , j.l . peng and r.l . greene , phys . b * 61 * , 1836 ( 2000 ) . isaacs , p.m. platzman , p. metcalf and j.m . honig , phys . lett . * 76 * , 4211 ( 1996 ) . kao , w.a.l . caliebe , j.b . hastings and j .- m . gillet , phys . b * 54 * , 16361 ( 1996 ) . et al . _ , b * 67 * , 45108 ( 2003 ) . p. abbamonte _ et al . _ , lett . * 83 * , 860 ( 1999 ) . g. dring _ et al . _ , phys . rev . b * 70 * , 085115 ( 2004 ) . s. suga , s. imada , a. higashiya , a. shigemoto , s. kasai , m. sing , h. fujiwara , a. sekiyama , a. yamasaki , c. kim , t. nomura , j. igarashi , m. yabashi , and t. ishikawa , phys . b * 72 * , 081101 ( 2005 ) . t. nomura and j .- i . igarashi , phys . b * 71 * , 035110 ( 2005 ) . s. wakimoto , y.j . kim , h. kim , h. zhang , t. gog , r.j . birgeneau , phys . rev . b * 72 * 224508 ( 2005 ) . e. collart , a. shukla , j.p . rueff , p. leininger , h. ishii , i. jarrige , y.q . cai , s.w . cheong , g. dhalenne , phys . lett * 96 * 147004 ( 2006 ) . j.w.seo , k. yang and d.w . lee , phys . rev . b , * 73 * 161104 ( 2006 ) . j. van den brink and m. veenendaal , europhysics letters,*73 * 121 ( 2006 ) , cond - mat/0311446 . hill , g. blumberg , y .- j . kim , d. ellis , s. wakimoto , r.j . birgeneau , s. komiya , y. ando , d. casa , t. gog , unpublished . j. van den brink , preprint cond - mat/0510140v1 ( 2005 ) . h.a . kramers and w. heisenberg , z. phys . * 31 * , 681 ( 1925 ) . p.m. platzman and n. tzoar , phys . rev . * 182 * , 510 ( 1969 ) . klein in light scattering in solids i , ed . m. cardona , springer , berlin ( 1983 ) . m. blume , j. appl . * 57 * , 3615 ( 1985 ) . t. tohyama , h. onodera , k. tsutsui and s. maekawa , phys . lett . * 89 * , 257405 ( 2002 ) , k. tsutsui , t. tohyama and s. maekawa , phys . lett . * 91 * , 117001 ( 2003 ) . m. van veenendaal , p. carra and b.t . thole , phys . b. * 54 * , 16010 ( 1996 ) . o. gunnarsson and k. schnhammer , phys . rev . b * 28 * , 4315 ( 1983 ) . a. kotani , t. mizuta , t. jo and j.c . parlebas , solid state comm . * 53 * , 805 ( 1985 ) . j. zaanen , c. westra and g.a . sawatzky , phys . b * 33 * , 8060 ( 1986 ) . see , e.g. , g.d . mahan , _ many - particle physics _ , plenum press , new york , pp 426 ( 1990 ) . on a heuristic basis sometimes the charge density response function is considered for the calculation of rixs spectra , e.g. in t.p . devereaux , g.e.d . mccormack and j.k . freericks , phys . lett . * 90 * , 67402 ( 2003 ) and t.p . devereaux , g.e.d . mccormack and j.k . freericks , phys . b * 68 * , 75105 ( 2003 ) , h. kondo , s. ishihara and s. maekawa , phys . b. * 64 * , 14414 ( 2001 ) , and s. ishihara and s. maekawa , phys . b * 62 * , 2338 ( 2000 ) . note that based on qualitative arguments , it was anticipated that the rixs cross section should depend on a dynamic four - particle correlation function : p.m. platzman and e.d . isaacs , phys . b. * 57 * , 11107 ( 1998 ) .
Abstract:

Answer: in indirect resonant inelastic x - ray scattering ( rixs ) an intermediate state is created with a core - hole that has an ultrashort lifetime . the core - hole potential therefore acts as a femtosecond pulse on the valence electrons . we show that this fact can be exploited to integrate out the intermediate states from the expressions for the scattering cross section . by this we obtain an effective scattering cross section that only contains the initial and final scattering states . we derive in detail the effective cross section which turns out to be a resonant scattering factor times a linear combination of the charge response function @xmath0 and the dynamic longitudinal spin density correlation function . this result is asymptotically exact for both strong and weak local core - hole potentials and ultrashort lifetimes . the resonant scattering pre - factor is shown to be weakly temperature dependent . we also derive a sum - rule for the total scattering intensity and generalize the results to multi - band systems . one of the remarkable outcomes is that one can change the relative charge and spin contribution to the inelastic spectral weight by varying the incident photon energy .