Instruction: Write an abstract for the following acadamic paper.

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

Paper: much experimental attention has recently been given to the topic of rabi - oscillations @xcite in microwave irradiated josephson systems at very low temperatures @xcite . a common motivation for the studies is the proposed macroscopic quantum behavior outlined in ref . @xcite , which has opened the way for potential applications of josephson technology in quantum information processing @xcite . this is proposed by operating a josephson junction in its zero - voltage state at temperatures below the quantum transition temperature , while manipulating the possible energy states with the application of commensurate microwave frequencies . the significance of rabi - oscillations is many fold , and includes _ i _ ) a direct connection to a well known quantum mechanical concept in perturbed atomic systems , and _ ii _ ) a method by which the control of quantum states of a system can be evaluated . thus , it is of significant importance to understand the nature of these observations in order for us to evaluate how to interpret the josephson system under investigation . recent classical analyses of josephson junctions , perturbed by microwaves and low temperature thermal fluctuations , have revealed that key signatures of microwave - induced multi - peaked switching distributions @xcite , used to display quantum mechanical features of josephson junctions , have direct classical analogs , which may obscure the interpretation of the observations @xcite . inspired by the work on slowly modulated transients to phase - locking in josephson systems by lomdahl and samuelsen @xcite , it has most recently been demonstrated @xcite that also rabi - oscillations have a classical analogue in microwave perturbed josephson junctions , providing a system response very similar to reported observations under the same conditions . we therefore denote this classical phenomenon rabi - type oscillations . it is the aim of this paper to complement and expand on the analysis provided in ref . @xcite in order to extract more detailed information from the classical model as it relates to observations of rabi - oscillations . we further demonstrate the validity of the analysis through direct quantitative comparisons with the observations of ref . @xcite , in which enough information is provided to conduct direct quantitative comparisons with our theory . specifically , we present an analysis based on direct perturbations in the dynamical variables . this approach results in closed form results that include dissipaion , which is of direct importance for experimental measurements of , e.g , coherence time . we will start section ii by providing the system equations that govern the classical model as it relates to specific reported measurements of rabi - oscillations . next , in section iii , we present a new , more complete , perturbation analysis of the slow transient modulations to phase - locking , and provide specific analytical results for ( rabi ) oscillation frequency and attenuation as a function of the system parameters . the theory is numerically validated in section iv . direct observations of rabi - type oscillations in the classical model are produced in section v through numerical simulations of a system parameterized to mimic reported experiments , and close agreement is found between these simulations , the developed theory , and the reported experiments . we conclude the paper in section vi . a normalized classical equation for a josephson junction can be written @xcite @xmath0 where @xmath1 is the difference between the phases of the quantum mechanical wave functions defining the junction , @xmath2 represents the dc bias current , and @xmath3 and @xmath4 represent microwave current amplitude and frequency , respectively . all currents are normalized to the critical current @xmath5 , and time is measured in units of the inverse plasma frequency @xmath6 , where @xmath7 , @xmath8 being the capacitance of the junction and @xmath9 the flux quantum . the phase of the applied microwave field is @xmath10 . tunneling of quasiparticles is represented by the dissipative term , where @xmath11 is given by the shunt resistance @xmath12 , and the accompanying thermal fluctuations are defined by the dissipation - fluctuation relationship @xcite @xmath13 @xmath14 being the temperature and @xmath15 is the characteristic josephson energy . a current pulse for probing the state of the system is represented by @xmath16 .    from the above normalizations , we define the normalized energy from the time independent perturbation terms in equation  ( [ eq : eq_1 ] ) @xmath17 the undamped bias current dependent plasma frequency is given by @xmath18{1-\eta^2 } \ ; . \label{eq : eq_6}\end{aligned}\ ] ] within this model , we will make analysis and conduct simulations according to experiments reported in ref . system parameters , such as @xmath2 , @xmath19 , @xmath4 , @xmath3 , and @xmath16 are matched as closely as possible to reported values , and rabi - type oscillations are observed statistically through simulations of the distribution of probe field induced switching from the zero - voltage state as a function of applied microwave field @xmath20 . microwave and probe perturbations have the form sketched in figure  [ fig : fig1 ] . we assume the phase @xmath10 of the microwave field is random for each switching event . notice the difference in probe field between this presentation ( along with ref . @xcite ) compared to refs . @xcite , where the probe was a microwave field with frequency slightly lower than @xmath4 . following the analysis of refs . @xcite , we will use the monochromatic ansatz for a solution of the unperturbed equation  ( [ eq : eq_1 ] ) ( @xmath21 ) @xmath22 where @xmath23 is the average phase , @xmath24 is an oscillation amplitude ( which is not necessarily small ) , and @xmath25 is the natural frequency . inserting this ansatz into equation  ( [ eq : eq_1 ] ) , we obtain the anharmonic relationships between the parameters @xmath26 where @xmath27 is the @xmath28th order bessel function of the first kind @xcite . notice that the ansatz dictates @xmath29 as a condition for existence of a zero - voltage state . also notice that @xmath30 and that @xmath31 for @xmath32 . the average system energy of the monochromatic ansatz is found by inserting equations  ( [ eq : eq_7 ] ) and ( [ eq : eq_8 ] ) into ( [ eq : eq_4 ] ) , then average over the period given by ( [ eq : eq_9 ] ) , @xmath33 where the minimum energy @xmath34 for a given @xmath2 is @xmath35 the applicability of this very simple monochromatic ansatz ( [ eq : eq_7 ] ) is illustrated in figure  [ fig : fig2 ] . here we have shown the average phase @xmath23 and oscillation frequency @xmath36 as a function of oscillation amplitude @xmath24 for two different values of dc bias @xmath2 . data are obtained from direct numerical simulations of equation  ( [ eq : eq_1 ] ) ( markers ) and from equations  ( [ eq : eq_8])-([eq : eq_29 ] ) ( curves ) . the agreement between the results of the ansatz and the simulations is very close , except at the highest amplitudes ( when @xmath37 ) , where higher harmonics of the dynamics contribute significant components to the dynamics . we use this ansatz to develop a perturbation theory for the system response to the application of microwaves and damping .    for @xmath38 we now separate the parameters of the ansatz into steady state and small time dependent deviations ; @xmath39 , @xmath40 , and @xmath41 , where @xmath42 , @xmath43 , and @xmath44 . inserting equation  ( [ eq : eq_7 ] ) into ( [ eq : eq_1 ] ) for @xmath45 and @xmath46 ( @xmath47 being heavyside s step function ) , yields the steady state phase - locked relationships @xcite @xmath48 and the linearized expressions for the small deviations from steady state : @xmath49\,\delta{\!a } = 2j_1(\bar{a})\sin\bar{\varphi}_0\,\delta\varphi_0 \nonumber \\ & & \label{eq : eq_13 } \\ & & \bar{a}\ddot{\delta{\theta}}+\alpha(\bar{a}\dot{\delta{\theta}}+ \omega_s\delta{\!a})+2\omega_s\dot{\delta{\!a}}+\left[\left(j_0(\bar{a})+ j_2(\bar{a}\right)\cos\bar{\varphi}_0-\omega_s^2\right]\bar{a}\,\delta{\theta } = 0 \ ; . \label{eq : eq_14}\end{aligned}\ ] ] equation  ( [ eq : eq_12 ] ) represents the small amplitude , slowly varying terms in ( [ eq : eq_1 ] ) , while ( [ eq : eq_13 ] ) and ( [ eq : eq_14 ] ) represent small amplitude terms oscillating with @xmath4 . evaluating a solution to equations  ( [ eq : eq_12])-([eq : eq_14 ] ) is not simple , but making the assumption that all three variables oscillate slowly with the same frequency @xmath50 and decay with the same attenuation @xmath51 , we can write the simple relationships : @xmath52 , @xmath53 , and @xmath54 , with @xmath55 . this results in the following polynomial @xmath56\omega^4 + i\alpha\left[\frac{5}{2}\alpha^2 - 2a_4\right]\omega^3   + \left[\frac{15}{16}\alpha^4-\frac{3}{2}\alpha^2a_4+a_2\right]\omega^2 \nonumber \\ & & -i\alpha\left[\omega_s^2\alpha^2+\gamma_1\right]\omega = \alpha^2\omega_s^2j_0(\bar{a})\cos\bar{\varphi}_0+\gamma_2 \ ; , \label{eq : eq_15}\end{aligned}\ ] ] where the parameters @xmath57 and @xmath58 are given in equations  ( [ eq : eq_17])-([eq : eq_21 ] ) below . the polynomial can be simplified considerably by inserting @xmath55 , and realizing that @xmath59 is the complex part of the roots ( for all roots ) with a non - zero real part @xmath50 ; i.e. , for underdamped modulations . notice that dissipative systems ( @xmath60 ) will naturally imply that @xmath61 for a nonzero threshold value of the microwave amplitude @xmath62 @xcite . with @xmath59 the ( real ) modulation frequency @xmath50 is determined by the following real polynomial : @xmath63 where @xmath64 \ ; . \label{eq : eq_21}\end{aligned}\ ] ] the steady - state resonance frequency @xmath65 is given by equation  ( [ eq : eq_9 ] ) as a function of @xmath66 . generally speaking , the three roots @xmath67 are real and positive . one is located near @xmath68 , one near @xmath69 , and one much smaller than @xmath69 . it is the latter we are interested in when studying transients and modulations to phase - locking , and it can be very well approximated by either of the following approximations , since the normalized frequency is small : @xmath70    a further simplification to the coefficients in the polynomial can be produced by noticing that the damping parameter @xmath19 in relevant josephson experiments is very small , allowing for omission of several terms of high order in @xmath19 in equations  ( [ eq : eq_17])-([eq : eq_19 ] ) . a specific limit of the above solution can be derived for small oscillation ( @xmath24 ) and microwave ( @xmath71 ) amplitudes for @xmath72 and @xmath73 @xmath74 the resonant choice @xmath75 is consistent with the experimental studies @xcite . in fact , measurements are usually conducted for @xmath4 very close to the linear resonance frequency ( see , e.g. , ref . @xcite ) . however , if @xmath76 , then @xmath77 , and the relationship between @xmath67 and @xmath78 becomes linear for small values of @xmath24 and @xmath71 . further , if @xmath79 , @xmath24 may become a multi - valued function of @xmath71 , corresponding to several resonant energy states of the system for the same parameter values @xcite . the solution to equation  ( [ eq : eq_16 ] ) , which can be given in closed form as a solution to a third order polynomial , is a complementary , more convenient , approach to the one presented in ref . it provides the important additional information regarding the attenuation @xmath51 of the transient modulations and it provides explicit expressions for the modulation frequency as a function of all system parameters . the two different perturbation methods agree qualitatively , and are quantitatively similar . the main difference between the two approaches is that we have here considered phase @xmath80 and amplitude @xmath24 as harmonically varying variables , instead of using the energy balance approach to transients @xcite , which treats the total energy of the system as the linearized dynamical variable . notice that the linearized equations  ( [ eq : eq_12])-([eq : eq_14 ] ) provide information about frequency , attenuation , and internal phase - relationships , but not about the specific magnitude of the deviations . however , the magnitude can be quantitatively estimated from the following reasoning . for a simulation , in which a josephson junction is described by equation  ( [ eq : eq_1 ] ) with @xmath45 ( or at least very low temperature ) and @xmath46 , the system is at rest @xmath81 for @xmath82 ( @xmath83 ) . the onset of the microwave field at @xmath84 will therefore , within the harmonic approximation , produce an envelope function of oscillation given by @xcite @xmath85 the two other modulated variables , @xmath86 and @xmath87 , can then be found from the coefficients @xmath88 the modulated system energy ( for @xmath89 ) can be calculated from equation  ( [ eq : eq_28 ] ) as @xmath90 where @xmath91 is the steady state energy of the phase - locked state in equations  ( [ eq : eq_10 ] ) and ( [ eq : eq_11 ] ) , and @xmath92 is the transient modulation . the average phase @xmath23 can be calculated either as a function of @xmath93 as given by equations  ( [ eq : eq_8 ] ) and ( [ eq : eq_25 ] ) , or it can be extracted from the theory through equation  ( [ eq : eq_26 ] ) . the two results are very similar , but since the latter is a first order approximation in the perturbation theory , it will have a second order error at @xmath84 . using equation  ( [ eq : eq_8 ] ) , @xmath94 will satisfy the true value both at @xmath84 and for @xmath95 . numerical validation of the expressions for the transient modulation can be directly acquired from simulating equation  ( [ eq : eq_1 ] ) for @xmath45 and @xmath96 for different values of @xmath97 and @xmath19 . simulations are conducted by choosing @xmath98 and @xmath99 , then initiating the system in the static state @xmath100 . at @xmath84 , the microwave field switches on , and we measure the phase @xmath101i , the average phase @xmath102 , and the system energy @xmath103 as given by equation  ( [ eq : eq_4 ] ) . figure  [ fig : fig3 ] shows a typical comparison between the analysis of the previous section and direct numerical simulations . we have chosen system parameters @xmath104 , @xmath105 , @xmath106 , and @xmath107{1-\eta^2}=0.6527147 $ ] ( inspired by ref . we can clearly see that the transient in phase and amplitude is a very slow modulation , and that this modulation provides a slow transient oscillation in the system energy . notice that the phase - locked frequency is _ much _ higher than the depicted modulation , @xmath108 . we have shown results of the theory on the left ( figures  [ fig : fig3]a and b ) and the corresponding simulation results on the right ( figures  [ fig : fig3]c and d ) . the quantitative agreement between the results is obvious . we do see a slightly larger predicted oscillation amplitude compared to the simulated . we also notice a slightly higher predicted modulation frequency than what is observed in the simulations . despite these slight discrepancies , which are usually about 5% and no larger than about 10% , we consider this remarkable agreement given the simple monochromatic trial function and the linearization in the perturbation treatment . figure  [ fig : fig4 ] shows the details of the perturbation variables as calculated from the theory , and they reveal that the actual oscillation frequency @xmath109 is only insignificantly different from @xmath4 . direct comparisons between the predicted modulation @xmath110 and the simulated values as a function of microwave amplitude for different values of @xmath2 and @xmath19 are summarized in figures  [ fig : fig5 ] and [ fig : fig6 ] . comparisons for @xmath50 are provided in figure  [ fig : fig5 ] and the attenuation @xmath51 is validated in figure  [ fig : fig6 ] . the general agreement is very good for all parameters tested . we observe the simple relationship @xmath111 for all oscillating solutions and the general trend of the polynomial ( equation ( [ eq : eq_15 ] ) or ( [ eq : eq_16 ] ) ) is to provide a modulation frequency slightly larger than what is observed in simulations . this is quantitatively similar to the values obtained from the energy balance perturbation theory developed in ref . @xcite , where the true modulation frequency was consistently slightly under - estimated . we finally observe the broad applicability of the asymptotic value of the modulation frequency given in equation  ( [ eq : eq_24 ] ) . as was argued in ref . @xcite , the rabi - oscillations observed in , e.g. , refs . @xcite may be closely related to the classical transient modulations described above . the numerical simulations in ref . @xcite were conducted in close agreement with the procedures described in ref . @xcite , by applying a small temperature @xmath112 ( @xmath113 ) to a josephson junction , which is perturbed by a microwave pulse starting at @xmath84 with a frequency @xmath107{1-\eta^2}$ ] . varying either the microwave amplitude or duration may yield rabi - type oscillations in switching probability when a subsequent microwave pulse with a smaller frequency @xmath114 is applied to probe the energy state of the junction as it is left by the first microwave pulse . the frequency of the second ( probe ) field was chosen to match an anharmonic amplitude that could excite the junction beyond the energy saddle point and lead to escape from the well . the measurements were conducted for varying microwave amplitude and fixed microwave duration . the experiments described in ref . @xcite were conducted slightly differently . first , the josephson system was a small - inductance @xmath115 interferometer , which we can conveniently approximate with a single josephson junction , since the dynamics for small @xmath115 is known to be well represented by a single degree of freedom @xcite . second , the probe pulse was in this case a `` dc '' pulse as sketched in figure  [ fig : fig1 ] , and the measurements were presented for fixed microwave amplitude and variable duration . we will here provide simulations of the classical system described in equation  ( [ eq : eq_1 ] ) and parameterized by information provided in ref . characteristic current and frequency are @xmath116 and @xmath117 , which lead to a normalized temperature of @xmath118 ( @xmath119 ) , and a normalized microwave frequency @xmath120 , which is close to the resonance such that @xmath121 . we have estimated the dissipation parameter @xmath19 from the reported decay of rabi - oscillations , @xmath104 . the normalized applied microwave amplitude @xmath71 is varied between 0 and @xmath122 , and the normalized duration is in the range @xmath123 $ ] . simulations are conducted with a probe pulse @xmath16 as shown in figure  [ fig : fig1 ] , followed by a short time in which we determine whether or not the junction has switched from the zero - voltage state . for every specific set of parameters , this type of simulation is conducted 25,000 times , each with a randomly chosen value of microwave phase @xmath10 , and the switching probability is recorded , before a new microwave duration is chosen . typical results are displayed in figure  [ fig : fig7]a , where the switching probability @xmath124 , reported as population occupancy of excited quantum state in ref . @xcite , is shown as a function of microwave duration . we clearly observe the rabi - type oscillations of this classical system with a distinct frequency . moreover , the oscillation frequency is in near perfect agreement with the reported comparable figure of the experimental measurements as well as the theoretical value @xmath50 provided by the analysis in the above section . figure  [ fig : fig7]b shows the simulated system energy as an average of many thermal realizations , each with a randomly chosen @xmath10 , of the trajectory ( at @xmath125 ) . the different choices of @xmath126 $ ] yield slightly different trajectories of @xmath127 , @xmath128 , and @xmath129 , which in turn make possible switching a function of @xmath10 . the importance of this ensemble average can be viewed in figure  [ fig : fig7]c , where a single trajectory of energy is shown for the applied temperature of @xmath119 . after a short transient , the thermal effects are dominating the individual trajectory , exciting the modulation frequency @xmath50 at random phase . only after an appropriate average over trajectories do we observe the smoothly decaying envelope of the modulated energy curve in figure  [ fig : fig7]b , which in turn is very similar to the ( single ) zero - temperature trajectory seen in figure  [ fig : fig7]d . the close connection between the transient energy modulation and the rabi - type modulation is obvious , and the close relationship between classical numerical simulations , classical theory , and the experimental measurements is further emphasized in figure  [ fig : fig8 ] , where the experimental data ( open markers ) of ref . @xcite are shown alongside theory ( solid curve ) and simulations ( closed markers ) . the three - way close agreement indicates that much about these measurements can be understood from classical theory of driven , damped josephson junctions . we have presented a perturbation analysis of the classical , nonlinear model describing the microwave - perturbed josephson junction . results show direct quantitative analogy between experimentally reported rabi - oscillations and the classical transients to phase - locking , and results are presented as specific functions of experimental system parameters . the analysis in this paper is both a complement and an extension to the analysis of ref . @xcite , where an energy balance perturbation approach was applied to produce similar connection between reports of rabi - oscillations and the nonlinear classical josephson system . the extension provides direct closed form solutions of a ( rabi - type ) modulation frequency as well as attenuation , and the resulting frequency agrees closely with the previously obtained result . this consistency lends credibility to the value of a classical interpretation of the reported rabi - oscillations , since the quantum mechanical and classical intuition turns out to be very similar ( as is also the case in laser physics @xcite ) . the macroscopic quantum picture interprets the observed oscillations as a result of a microwave induced temporal variation in the population probability of two , or more , intrinsic quantum mechanical energy levels in the josephson washboard potential . this slow variation is then indirectly observed through the associated variation in tunneling probability , switching from the zero - voltage state as a result of the application of the probe pulse . the stochastic nature of the system is due to the absorption of microwave photons as well as the subsequent tunneling probability of the measurement . the classical picture presents a system with a microwave induced temporally modulated energy . the slow variation is indirectly observed through the associated variation in escape from the energy well when the probe is applied , with high probability for passing the energetic saddle point during times of elevated energy and relatively little escape probability during times of depressed energy content . the stochastic nature of this system is due to the random phase of the microwave signal as well as the thermal fluctuations . additional analogies between intrinsic quantum mechanical energy levels and the multi - valued resonances of the classical nonlinear system briefly outlined above ( as well as in ref . @xcite ) further bridges the gap between what can be expected from the two interpretations of the microwave induced measurements . since the experimental measurements are concerned only with detecting the escape event , and since the classical and quantum mechanical interpretations seem to provide the same signatures for that escape , we are faced with a fundamental ambiguity of how to read this information . potential applications of josephson technology for quantum information processing will therefore benefit from the development of new unambiguous measurements , which must present signatures of macroscopic quantum behavior that can not be explained classically . however , we submit that , given the very close quantitative agreement between our classical theory and the experimental measurements , the intuition and the closed form expressions presented in this paper can be directly applied to guide future experiments and microwave manipulations of josephson systems . this work was supported in part by the uc  davis center for digital security under afosr grant fa9550 - 04 - 1 - 0171 and in part by miur ( italy ) cofin04 . n. grnbech - jensen , m. g. castellano , f. chiarello , m. cirillo , c. cosmelli , v. merlo , r. russo , and g. torrioli , in  quantum computing : solid state systems  , edited by b. ruggierio , p. delsing , c. granata , y. paskin , and p. silvestrini ( kluwer academic and springer publishers , new york , in press ) , cond - mat/0412692 . a.  barone and g.  patern , _ physics and applications of the josephson effect _ ( wiley , new york , 1982 ) ; t. van duzer and c. w. turner , _ principles of superconductive devices and circuits _ , 2nd ed.(prentice - hall , new york , 1998 ) .
Abstract:

Answer: we present a classical analysis of the transient response of josephson junctions perturbed by microwaves and thermal fluctuations . the results include a specific low frequency modulation in phase and amplitude behavior of a junction in its zero - voltage state . this transient modulation frequency is linked directly to an observed variation in the probability for the system to switch to its non - zero voltage state . complementing previous work on linking classical analysis to the experimental observations of rabi - oscillations , this expanded perturbation method also provides closed form analytical results for attenuation of the modulations and the rabi - type oscillation frequency . results of perturbation analysis are compared directly ( and quantitatively ) to numerical simulations of the classical model as well as published experimental data , suggesting that transients to phase - locking are closely related to the observed oscillations .


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: general relativity describes very well gravity at the classical level , but a quantum theory of gravity requires the introduction of a more general framework . one of the most promising candidates is presently string theory . this theory is believed to change drastically the short - range behavior of classical gravity , but also some of its global properties can be modified , such as black hole thermodynamics . for example , the study of the solutions of effective low - energy theory has shown that , due to the presence of non - minimal couplings , non - trivial scalar hair can arise @xcite , in contrast with classical general relativity , where no - hair theorems @xcite rule out this possibility . the effects of string theory on gravitational physics are usually investigated by means of effective field theory actions , obtained through a perturbative expansion in the string tension @xmath0 . at the tree level , the effective action of the heterotic ( but also other types of ) string contains a coupling of the dilaton with gravity via the term . the solutions of this model have been extensively studied in the literature , both in a perturbative @xcite and numerical @xcite approach . it turns out that the model admits solutions with non - trivial dilatonic hair . the scalar charge is not an independent parameter , but is a function of the mass of the black hole , and is therefore an example of secondary hair @xcite . also the thermodynamics is different from that of black holes . in particular , it was shown that the theory predicts a lower bound on the mass of black holes @xcite , which corresponds to the state of highest ( but finite ) temperature and lowest entropy . the configuration of minimal mass should be identified with the ground state of the hawking evaporation process .    in order to build realistic models , one should however take into account that in string theory other scalar fields are present in the spectrum in addition to the dilaton , as for example the moduli , which originate from the compactification of the higher - dimensional spacetime . these also couple to gravity through one - loop effects . at leading order the coupling term is proportional to the logarithm of a dedekind @xmath1-function of the moduli , which multiplies the term @xcite . the effect of the non - minimal coupling of the moduli to gravity in a cosmological context has been studied in several papers @xcite and it has been shown that in some cases it may lead to models without initial singularities , but to our knowledge no investigation has been devoted till now to its implications on physics . on the other hand , it is well known that in effective string actions the electromagnetic field exhibits a non - minimal coupling to the dilaton and the moduli similar to that of the term @xcite . the solutions have been thoroughly studied in this case : if one neglects the moduli , one obtains exact magnetically charged solutions , with secondary scalar hair , the scalar charge being a function of the mass and the magnetic charge @xcite . if one instead takes into account also one modulus , the general solution can no longer be written in analytic form , but it can nevertheless be shown to depend on three parameters @xcite : thus in this case a new independent parameter arises , besides the mass and the charge , and one may speak of a primary scalar hair .    in this paper , we investigate if a similar phenomenon can occur in the purely gravitational sector . since we are mainly interested in showing the existence of primary scalar hair in the scalar - gravity sector , we consider a simplified model with a dilaton and a unique modulus which couples exponentially with the term . we study it both in a perturbative and numerical setting using the techniques developed in ref . @xcite and @xcite , respectively . we find that in fact the qualitative features are similar to the case of maxwell coupling . we also find that an upper limit must be imposed on the scalar charges for given mass , in order to avoid naked singularities . this is reminiscent of the extremality bounds in the multiscalar einstein - maxwell case @xcite . the structure of our paper is the following . in section 2 we present the perturbative solution and discuss its thermodynamical properties . in section 3 we describe the numerical solution and the occurrence of an upper bound for the scalar charges . section 4 contains a discussion and the main conclusions . the bosonic sector of the effective action for the heterotic string in absence of yang - mills and axionic fields , is given at leading order in @xmath2 by i_eff=116d^4 x where @xmath3 , @xmath4 is a coupling constant of order unity , @xmath5 is the dilaton , @xmath6 is a modulus , whose coupling with the term , @xmath7 has been taken for simplicity to be of exponential form . the equations can be written as [ fe ] & & g_=t_^()+t_^( ) , + & & ^2=2,^2=2 , where + & & -8 + & & + 2_m_n - g_mn_m_n , and an alougous expression for @xmath8 . the total energy - momentum is conserved but , as noticed in @xcite for the case of a single scalar , its time component , corresponding to the total energy , is not positive definite , due to the contribution of the term , leaving room for the possibility of a violation of the no - hair conjecture . we look for solutions , with scalars @xmath9 , @xmath10 . a generic metric can be written as [ metric ] -(r)dt^2+^2(r)(r)dr^2+r^2(r)d^2 . instead of substituting the ansatz ( [ metric ] ) into the field equations ( [ fe ] ) , it is easier to substitute it directly into the action , and then vary with respect to the fields @xmath5 , @xmath6 , @xmath11 , @xmath12 and @xmath13 . the action reads :     + & & -2(1-^r^2 ) ]    the can then be written in the form : [ feq ] +   +   + & & - 2v^(1-^r^2)- r^2 , +   +   + & & = 4v(^-3r)+ 4v^-3r-2r , where @xmath14 and a prime denotes derivative with respect to @xmath15 . only four of the previous equations are independent .    in order to find approximate solutions to the field equations , we adopt the approach of ref . @xcite and expand the fields around the background constituted by the schwarzschild metric with vanishing scalar fields , which is of course a solution for @xmath16 . our expansion is in the parameter @xmath17 , @xmath18 being the mass of the solution . since @xmath19 is believed to be of order unity in planck units , the expansion is valid for large @xmath18 , in the region where @xmath20 , i. e. for @xmath21 . for macroscopic black holes ( @xmath22 ) this condition is always satisfied , except in a neighborhood of the singularity , well inside the horizon ( region of strong curvature ) . in particular , the approximation is valid for the discussion of the asymptotic properties of the fields , and the questions concerning the scalar hair . near the physical singularity , however , the higher order corrections to the effective string lagrangian become important and the perturbation theory is no longer reliable . we shall however discuss this regime using numerical techniques in the next section . at this point , it must be observed that the ansatz ( [ metric ] ) for the metric is too general and still leaves the possibility of a choice of gauge . in order to perform the perturbative calculations , the most convenient choice @xcite is to impose @xmath23 , i.e. [ ghs ] -(r)dt^2 + 1(r)dr^2+r^2(r)d^2 . this gauge was also used for finding exact charged solutions in effective string theory @xcite . we expand the fields as follows : [ exp ] & = & _ 0(1+_1+^2_2+^3_3 +  ) , + r & = & r+_1+^2_2+^3_3 +  , + & = & _ 1+^2_2+^3_3 +  , + & = &  _ 0+_1 + ^2_2 + ^3_3 +  , where @xmath24 . we have normalized @xmath5 such that @xmath25 at infinity . this is always possible , by rescaling the coupling constant @xmath0 ( this means that our expansion is actually in @xmath26 ) . however , it is not possible to rescale independently also @xmath6 , and hence we take @xmath27 at infinity . the parameters @xmath4 and @xmath28 will always appear in the combination @xmath29 . substituting the expansion ( [ exp ] ) into the ( [ feq ] ) , one obtains at first order ]=24m^2r^4,]=24m^2zr^4 . with the previous , requiring regularity at the horizon @xmath30 , the scalar fields are uniquely determined at first order : = z = -1 m ( 1r+mr^2 + 4m^23r^3 ) , the equations for the metric fields are given at the same order by  = 0,]=-2mr^2 . we impose the boundary conditions that @xmath31 const , @xmath32 at infinity . we are still free to choose the boundary conditions at @xmath30 . changing the at @xmath30 yields a reparametrization of the solutions , but no change in their physical properties : in particular , the relations between the physical quantities , like mass , temperature and entropy , are independent of the parametrization . the most convenient choice is to require that the @xmath33 and @xmath34 are regular at @xmath30 . this is equivalent to fix the location of the horizon at @xmath30 . with these boundary conditions , @xmath35 , @xmath36 . this choice greatly simplifies the higher order calculations . we can now evaluate the second order corrections . with the stated boundary conditions , the equations for the second order perturbations give +   +   +   + & & + 8mr-2mr^2(+z)-16mr-3mr^3(+z),whose solution is & & _ 2=_2z^2=-1m^3 ( 7360r+73m60r^2 + 73m^245r^3 + 73m^330r^4 + 112m^475r^5 + 8m^59r^6 ) , + & & _ 2=-1 + z^2m^2 ( 12r+2m3r^2 + 7m^23r^3 + 16m^35r^4 + 24m^45r^5 ) , + & & _ 2=- 1 + z^2m^3 ( 16r+m3r^2 + 4m^23r^3 - 14m^33r^4 - 136m^415r^5 -272m^515r^6).up to this order , @xmath5 and @xmath6 are proportional in this gauge . however , in order to clarify the structure of the solutions , it is useful to go to the next order , even if such corrections would of course be modified by taking into account terms of order @xmath37 in the action . a long but straightforward calculation gives _ 3=-73 ( 1 + z^3)+,_3=-73 ( 1 + z^3)+,_3=- 16480 + 3969z^2+,_3 = - z ( 3969 + 16480z^2 ) + .from these results appears that the metric functions are expanded in terms of @xmath38 , while the scalar fields depend in a more involved way from the parameters . moreover , the functional dependences of the dilaton and the modulus on @xmath15 are different if @xmath39 . the perturbative solutions have the following properties : a horizon is present at @xmath30 , while a singularity is located at the zero of @xmath13 ; the evaluation of this zero is however outside the range of validity of our approximation . it can be expected nevertheless that for small values of the mass , or great values of @xmath40 , the zero can occur for @xmath41 , leading to the presence of naked singularities . this will be confirmed by the numerical results of next section . the mass @xmath42 of the black hole can be deduced from the asymptotic behavior of the metric function @xmath11 and is given by [ mass ] m = m(1 + 112 ^ 2(1+z^2)m^4 + 73360 ^ 3(1+z^3)m^6 ) . its value is greater than that of the with equal radius . analogously , from the asymptotic behaviour of @xmath5 and @xmath6 one can deduce the scalar charges @xmath43 and @xmath44 , which , in terms of the mass @xmath42 , turn out to be & & d_=m(1 + 7360m^2 + 17110 + 4599z^27560 ^ 2m^4 ) , + & & d_=zm(1 + 7360zm^2 + 4599 + 17110z^27560 ^ 2m^4 ) . it is clear that , contrary to the case in which only one scalar field is present @xcite , the scalar charges are no longer function only of the mass of the black hole , but depend also on another parameter , which we identify with @xmath40 . hence , in analogy with the dilaton - modulus gravity non - minimally coupled to the electromagnetic field @xcite , also in this case a primary scalar hair is present in the solution . this gives an example of primary scalar hair in pure gravity models . we notice that , at leading order , @xmath45 . the temperature of the can be defined as usual as the inverse periodicity of the time coordinate which renders regular the euclidean section of the metric . this is given by @xcite [ temp ] t=14dg_00dr _ hor , which yields , at order @xmath46 : t=18m(1+^2(2m)+^3_3(2m)).taking into account ( [ mass ] ) , a straightforward calculation leads to t=18m(1 + 73120 ^ 2(1+z^2)m^4 + 125117560 ^ 3(1+z^3)m^6).the temperature is higher than that of a of equal mass , but is still a monotonically decreasing function of the mass . the entropy @xmath47 can be defined by means of the euclidean formalism @xcite as [ entropy ] s =- i_e , where @xmath48 is the inverse temperature and @xmath49 the euclidean action [ eaction ] i_e&=&-116_m[-2()^2 - 2()^2 +  ( + ) ) ] dv + & & -18_m(k - k_0)dwhere @xmath50 is the exterior curvature . a lenghty calculation gives s=4m^2(1 + 73120 ^ 2(1+z^2)m^4 + 1251115120 ^ 3(1+z^3)m^6 ) .    in the range of validity of our approximation , the thermodynamical quantities of course do not differ much from their background values , except that now they depend on the further parameter @xmath40 . they behave differently only for @xmath51 , where however the approximation breaks down . it is interesting to notice that all the thermodynamical quantities are expanded in terms of @xmath52 . as discussed previously , the perturbative solution is not valid in the whole domain of definition of the solution . since we are not able to find the exact analytical solution of the field equations , we use the numerical approach described in ref . for this calculation it is more convenient to choose a gauge in which the metric function @xmath13 is identified with the radial coordinate @xmath15 , i.e. @xmath53 where @xmath54 , @xmath55 . comparison with ( [ metric ] ) yields @xmath56 . of course , the physical quantities do not depend on the choice of gauge . however , in order to make the comparison of results easier , we give the perturbative expansions in the new coordinates : +   +   + where use has been made of the condition @xmath57 , and the functions @xmath33 , @xmath58 @xmath59 and @xmath60 are those obtained above , evaluated at @xmath61 . in particular , the metric function are , up to order @xmath37 , & = & 1 - 2mr-^2(1+z^2)m^4(m6r+5m^33r^3 -6m^4r^4 + 74m^515r^5 + 32m^65r^6 + 688m^715r^7 ) + & = & 1-^2(1+z^2)m^4(m^22r^2 + 4m^33r^3 + 7m^4r^4 + 64m^55r^5 + 24m^6r^6 )    in the parametrization ( [ e2.1 ] ) the einstein - lagrange equations can be written in a matrix form ( we set all the string coupling constants to be equal to one for simplicity ) @xmath62 where @xmath63 and the entries of the matrices @xmath64 and @xmath65 are @xmath66    the last ( constraint ) equation is @xmath67 \nonumber \\ & + & 4 \biggl [ \lambda_\phi e^{-2 \phi } \phi '      + \lambda_\sigma q e^{-2 q \sigma } \sigma ' \biggr ]      \delta ' ( \sigma^2 - 3 \delta ) = 0 \nonumber\end{aligned}\ ] ]    the numerical integration is performed as follows : we start from a neighborood of the horizon and integrate towards infinity . the mass and charges of the solution are then evaluated from the asymptotic behaviour of the metric functions and the integration is performed again backwards . more technical details on the numerical procedure can be found in @xcite . the behaviour of the solution outside the horizon agrees with that obtained by perturbative methods in the previous section ( see figures 1 - 4 ) and is similar to that of the single - scalar solutions discussed in @xcite . when @xmath68 vanishes , we recover the solutions described in ref . in particular , one must impose a lower bound on the mass in order to avoid the occurrence of naked singularities . for non - zero value of @xmath68 the situation changes significantly . taking @xmath42 and @xmath43 fixed , for small values of @xmath68 the behavior of the solution does not differ much from the case of a single scalar field , except that the position of the inner black hole singularity slowly moves up . this situation is shown by the solid lines in figs.1 - 4 , where the dependence of the functions @xmath69 , @xmath70 , @xmath71 and @xmath72 against the radial coordinate @xmath15 is plotted . when @xmath68 reaches a critical value @xmath73 the positions of the singularity @xmath74 and the horizon @xmath75 coincide . the dependence of @xmath73 against black hole mass is represented on figure 5 . when @xmath76 a naked singularity appears . this situation is shown by dashed lines in figs . 1 - 4 . this naked singularity is a continuation of the black hole inner singularity and has the same nature . of course , since the system is symmetric for the exchange of @xmath5 and @xmath6 , an analogous behavior is expected when @xmath5 is varied and @xmath6 held fixed . the dependence of @xmath73 from the black hole mass @xmath42 is approximately linear . hence , according to the cosmic censorship conjecture , the mass of the gives an upper limit for the modulus / dilatonic field charge . = 0.7=0.4    = 0.7=0.4    = 0.7=0.4    = 0.7=0.4    from a mathematical point of view the appearance of this singularity is a consequence of the vanishing of the second factor ( in brackets ) of the determinant @xmath77 of the system ( [ syst1 ] ) @xmath78 , \qquad \mbox{where } \nonumber \\ a & = & 64 r^2 \sigma^5 \biggl [          4 r^2 \sigma^2 e^{-4\phi } \phi'^2       +   8 r^2 \sigma^2 e^{-2\phi } \phi ' e^{-2\sigma } \sigma '       +   4 r^2 \sigma^2 e^{-4\sigma } \sigma'^2 \nonumber \\     & - &   \sigma^2 e^{-4\phi }       -   \sigma^2 e^{-4\sigma }       +   12 e^{-6\phi } \phi ' \delta '       +   12 e^{-4\phi } e^{-2\sigma } \sigma ' \delta ' \nonumber \\     & + & 12 e^{-2\phi } \phi ' e^{-4\sigma } \delta '       +   12 e^{-6\sigma } \sigma ' \delta ' \biggr ] \nonumber \\ b & = &   32 r^2 \sigma^6 \biggl [           r^3 \sigma^3 e^{-2\phi } \phi '        +   r^3 \sigma^3 e^{-2\sigma } \sigma '        +   r^3 \sigma^2 e^{-2\phi } \phi '        +   r^3 \sigma^2 e^{-2\sigma } \sigma ' \nonumber \\      & + & 2 r \sigma e^{-4\phi } \delta '        +   2 r \sigma e^{-4\sigma } \delta '        +   2 r e^{-4\phi } \delta '        +   2 r e^{-4\sigma } \delta '        +   4 \sigma^3 e^{-4\phi }   \nonumber \\      & + & 4 \sigma^3 e^{-4\sigma }        -   16 \sigma e^{-6\phi } \phi ' \delta '        -   16 \sigma e^{-4\phi } e^{-2\sigma } \sigma ' \delta ' \nonumber \\      & - & 16 \sigma e^{-2\phi } \phi ' e^{-4\sigma } \delta '        -   16 \sigma e^{-6\sigma } \sigma ' \delta'\biggr ] \nonumber \\ c & = &   4 r^2 \sigma^8 \biggl [            r^4 \sigma^2         -   16 r \sigma e^{-4\phi } \delta '         -   16 r \sigma e^{-4\sigma } \delta '         -   16 r e^{-4\phi } \delta ' \nonumber \\       & - & 16 r e^{-4\sigma } \delta '         -   16 \sigma^3 e^{-4\phi }         -   16 \sigma^3 e^{-4\sigma }         -   64 \sigma e^{-6\phi } \phi ' \delta ' \nonumber \\       & - & 64 \sigma e^{-4\phi } e^{-2\sigma } \sigma ' \delta '         -   64 \sigma e^{-2\phi } \phi ' e^{-4\sigma } \delta '         -   64 \sigma e^{-6\sigma } \sigma ' \delta ' \biggr ] \nonumber\end{aligned}\ ] ] a naked singularity occurs when the factor in the bracket of @xmath77 vanishes before the metric function @xmath69 . figure 6 represents the dependence of the position of the zero of @xmath77 against the parameter @xmath40 . = 0.7=0.4    = 0.7=0.4    the thermodynamical parameters can evaluated numerically and compared with the perturbative results . this is interesting especially in order tounderstand their behaviour for small mass , where the perturbative approach fails . the temperature can be obtained from ( [ temp ] ) in a straightforward way . for the calculation of the entropy , eq . ( [ entropy ] ) has been used . in particular , the euclidian action @xmath49 was evaluated by adding its definition as an additional equation to the main system ( [ syst1 ] ) . the numerical evaluation of the black hole temperature @xmath79 and entropy @xmath47 agrees with the perturbative results for great @xmath42 , see the table 1 . a similar agreement holds for the entropy . .the dependence of the black hole temperature @xmath79 against the mass @xmath42 and the parameter @xmath80 for the numerical and perturbative solutions . the agreement is better for great @xmath42 or small @xmath40 , in accordance with the limit of validity of the perturbative approach . [ cols="<,<,<,<",options="header " , ]     finally , we notice that the thermodynamical parameters stay finite in the extremal case . the thermodynamics is therefore analogous to that studied in absence of modulus fields @xcite , except that now one has one further independent parameter ( the scalar charge ) . we have shown perturbatively the occurrence of primary scalar hair in black hole solutions of models with more than one scalar field non - minimally coupled to gravity via the term . this result has been checked numerically . from the numerical calculations also follows that naked singularities can appear for small values of the mass ( as in pure dilaton - models ) , or for large values of the scalar charges . this is a novel feature of the model under study , and can be compared with a similar phenomenon occurring in multi - scalar einstein - maxwell models @xcite . in that case some analytical relations for the extremality condition of the black holes can be obtained , while in our case this seems not to be possible . we can however conjecture that a relation of the same kind exists also in our case . s. m. wishes to thank sternberg astronomical institute for kind hospitality during the last stages of this work . this work was partially supported by a coordinate research project of the university of cagliari and by `` universities of russia : fundamental investigations '' program via grant no . 99 d. garfinkle , g.t . horowitz , and a. strominger , _ phys . * d 43 * , 3140 ( 1991 ) . p. kanti , n.e . mavromatos , j. rizos , k. tamvakis and e. winstanley , _ phys . * d 54 * , 5049 ( 1996 ) ; p. kanti and k. tamvakis , _ phys . * b 392 * , 30 ( 1997 ) ; t. torii , h. yajima , and k. maeda , _ phys . * d 55 * , 739 ( 1997 ) . r. easther and k. maeda , _ phys . _ * d 54 * , 7252 ( 1996 ) ; p. kanti , j. rizos , and k. tamvakis , _ phys . rev . _ * d 59 * , 083512 ( 1999 ) ; s. alexeyev , a. toporensky , v. ustiansky , _ class . * 17 * , 2243 ( 2000 ) .
Abstract:

Answer: we show that the black hole solutions of the effective string theory action , where one - loop effects that couple the moduli to gravity via a term are taken into account , admit primary scalar hair . the requirement of absence of naked singularities imposes an upper bound on the scalar charges . # 1#2 c i u


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: hadron ( baryons and mesons ) spectroscopy opens the gate for the development of the fundamental theory of the strong interaction : quantum chromodynamics ( qcd ) . however the non - perturbative complication of low energy qcd makes it impossible to calculate the hadron structure analytically from qcd directly . the unique color structure of the known hadrons makes the construction of quark models very efficient ( fig.1 ) . a variety of quark models employing two - body interactions give a good description of hadron properties . ( for baryon , the three - body interaction can be well approximated by two - body interaction . ) however the unique color structure also limits our understanding of the properties of other color structures available in qcd . in order to understand low energy qcd , to study system with more quarks is indispensable . hadron - hadron scattering provides a window on the nature of other color structures , but it is not enough because the colorless meson exchange model and chiral perturbation theory both describe the low energy hadron - hadron scattering data well . qcd does not rule out the existence of glueballs , quark - gluon hybrids , multi - quark states , _ etc . _ based on the present understanding . multi - quark systems are important samples for providing information on low energy qcd interaction , especially for complex color structures ( fig.2 ) . = 3.5 in    fig.1 the color structure of baryons and mesons . since 1977 , interest in multi - quark systems has continued although there have been rising and falling . the search for the h - particle ( 6-quark system ) , which was predicted in 1977 by r. jaffe @xcite with mit bag model , has not succeeded for almost 30 years . in 1993 , the dibaryon state @xmath2 appeared unexpectedly @xcite , observed in double charge exchange reactions but disappeared nine years later . other six quark states , : @xmath3 @xcite , @xmath4 @xcite and @xmath5 @xcite _ etc . _ , had been proposed but none of them has been established experimentally . in 2003 , the pentaquark state , @xmath6 , aroused a new enthusiasm in multi - quark systems . more than ten groups claimed that they observed the @xmath6 signal , but almost the same number of groups did not @xcite . its appearance had raised great trouble for theory of hadron spectroscopy , almost none of the models constrained by hadron properties and hadron - hadron scattering can account for the @xmath6 @xcite . in studying pentaquark @xmath6 , various color and spatial structures of pentaquark had been proposed : color singlet hadron molecules [ k(@xmath7)n(@xmath8 ) ] , color anti - triplet diquarks [ @xmath9 @xcite , diquark - triquark [ @xmath10 @xcite , quark methane [ @xmath11 @xcite , _ etc_. however recent high statistic experiments did not confirm the @xmath6 signal @xcite . today , about three years later , pentaquark @xmath6 seems to be about to disappear also . after the pentaquark @xmath6 ( a member of anti - decuplet ) , other states in 27-plet and 35-plet @xcite , had been proposed . in order to understand the nucleon spin structure within constituent quark model the five - quark component is necessary @xcite . recently to explain the strange magnetic moment of proton , five - quark component was introduced in the nucleon @xcite . tetra - quark states are re - interested both experimentally and theoretically because of new discoveries since 2003 @xcite . = 4.5 in    fig.2 the color structure of multi - quark states . fig.2 shows possible color structures for 5- and 6-quark systems based on the lattice qcd calculation @xcite and general color confinement idea . clearly they have more possible color structures than 2- and 3-quark hadrons . due to color confinement , only color singlet combinations of quarks can be separated . colorful clusters will be confined in a genuine multi - quark system . it will transit to color singlet sub - systems through color flux rearrangement first then decay . this will induces resonance similar to compound nucleus formation but due to color confinement and could be called color confinement resonance . to include these intermediate hidden color configurations in the quark model , a multi - channel coupling calculation with multi - body interactions is required . this is quite involved @xcite and finding a method to make the calculation tractable is an important element in the study of multi - quark systems with constituent quark models . moreover up to now we almost do nt have any idea about the transition interaction between different color configurations . to develop model to include effects of these hidden color configurations is expected .    the group theory classification of 4 , 5 and 6 quark states had been published @xcite . the fractional parentage ( fp ) expansion technique has been proven as a powerful method for few - body problem . to full play the powerful group theoretic method for quark model calculation one needs not only the fractional parentage expansion coefficients of these multi quark states but also a relation between various quark model states ( hereafter called physical bases ) and the group theoretical classification states ( hereafter called symmetry bases ) . such a method has been developed and successfully applied in systematic search of dibaryon @xcite , where the physical bases were transformed to symmetry bases first , then the 6-body matrix elements calculation of hamiltonian ( with two body interaction ) on the symmetry bases was done which can be reduced to four - body overlap and two - body matrix elements calculation by means of fractional parentage expansion . at last the matrix elements on the symmetry bases were transformed back to physical bases . the main content of this paper is to provide the transformation coefficients between physical bases and symmetry bases of five - quark systems to facilitate the calculation of many - body hamiltonian matrix elements . the physical bases discussed in this paper is the jaffe - wilczek diquark model ones . another useful physical bases is the meson - baryon bases which will be given in a company paper . to illustrate the application of this group theory method for five - quark systems three quark models were employed for pentaquark calculation . they are the naive quark model , i.e. , the glashow - isgur model @xcite ; the salamanca chiral quark model @xcite and the quark delocalization color screening model ( qdcsm ) developed by our group @xcite . the calculation of the related fractional parentage coefficients will be mentioned but the needed results have been published elsewhere @xcite .    in section ii , the physical bases and symmetry bases are introduced and the transformation between them is derived . the fractional parentage technique applied to calculate the matrix elements on the symmetry bases is also explained in this section . section iii explains three quark models we used . the results of the systematic calculation of pentaquark in the @xmath12 3-flavor world are given in section iv . the last section gives the summary . the physical bases are constructed as follows , first the wave function of each quark cluster was constructed based on the group chain classification @xmath13\right\},\ ] ] then the quark cluster wave functions of the system was coupled to definite color , spin and isospin quantum numbers by clebsch - gordan ( cg ) coefficients of color @xmath14 , spin @xmath15 and isospin @xmath16 group , and finally anti - symmetrized . = 2.5 in    fig.3 the jaffe - wilczek configuration of pentaquark . the symmetry basis of 4-quark system is just the group chain ( eq.([group - chain ] ) ) classification basis ( @xmath38 ) , @xmath39w_{\nu_4 } \\ \left [ c_4 \right ] w_{c_4 } [ \mu_4][f_4]y_4i_4m_{i_4}j_4m_{j_4 } \end{array } \right\rangle,\ ] ] where @xmath40 stands for @xmath41 . similar to eq.([phys_q5 ] ) , the symmetry bases for 5-quark is @xmath42\bar{i}\bar{j}}(\bar{q}_5 ) \right]^{[c]ij}_{wm_im_j}. \label{symm_q5}\ ] ] the cluster bases and symmetry bases for 5-quark can be transformed to each other @xcite . @xmath43[c_4][\mu_4]}_{[\tilde{\nu}_2][c_2][\mu_2 ] , [ \tilde{\nu}'_2][c'_2][\mu'_2 ] } c^{[\mu_4][f_4][j_4]}_{[\mu_2][f_2][j_2],[\mu'_2][f'_2][j'_2 ] } c^{[f_4]y_4i_4}_{[f_2]y_2i_2,[f'_2]y'_2i'_2 } \phi_{\alpha k}(q^4\bar{q } ) , \label{trans}\end{aligned}\ ] ] @xmath44s are the isoscalar factors of @xmath1 , which can be obtained from the book @xcite . all the transformation coefficients are listed in the appendix . a physical 5-quark state with quantum number @xmath45 is expressed as a channel coupling wave function @xmath46 the channel coupling coefficient @xmath47 is determined by the diagonalization of the 5-quark hamiltonian as usual . the calculation of hamiltonian matrix elements in the cluster bases is tedious and it can be replaced by the matrix elements in the symmetry bases by the transformation eq.([trans ] ) , @xmath48 in the symmetry bases , the matrix elements @xmath49 can be calculated by the well known fp expansion method . because there is an antiquark , we have to use different fp expansion for @xmath50 interaction and @xmath7 interaction . for @xmath50 interaction , @xmath51 @xmath52[\nu_4][\tilde{\nu}_4]}_{[1 ^ 2][\nu_1][\tilde{\nu}_1],[1 ^ 2][\nu_2][\tilde{\nu}_2 ] } c^{[\tilde{\nu}_4][c_4][\mu_4]}_{[\tilde{\nu}_1][c_1][\mu_1],[\tilde{\nu}_2][c_2][\mu_2 ] } c^{[\mu_4][f_4]j_4}_{[\mu_1][f_1]j_1,[\mu_2][f_2]j_2 } c^{[f_4]y_4i_4}_{[f_1]y_1i_1,[f_2]y_2i_2 }     \nonumber\\ & & c^{[1 ^ 4][\nu_4][\tilde{\nu}_4]}_{[1 ^ 2][\nu'_1][\tilde{\nu}'_1],[1 ^ 2][\nu'_2][\tilde{\nu}'_2 ] } c^{[\tilde{\nu}_4][c_4][\mu_4]}_{[\tilde{\nu}'_1][c'_1][\mu'_1],[\tilde{\nu}'_2][c'_2][\mu'_2 ] } c^{[\mu_4][f_4]j_4}_{[\mu'_1][f'_1]j'_1,[\mu'_2][f'_2]j'_2 } c^{[f_4]y_4i_4}_{[f'_1]y'_1i'_1,[f'_2]y'_2i'_2 }      \nonumber\\ & & c^{[\nu_4]w_{x_4}}_{[\nu_1]w_{x_1},[\nu_2]w_{x_2 } } c^{[\nu_4]w_{x'_4}}_{[\nu'_1]w_{x'_1},[\nu'_2]w_{x'_2 } } \langle \alpha_1 k_1 | \alpha^{\prime}_1 k^{\prime}_1 \rangle \langle \alpha_2 k_2 |h_{34 } | \alpha^{\prime}_2 k^{\prime}_2 \rangle .\end{aligned}\ ] ] here @xmath44s are the @xmath1 isoscalar factors , @xmath53 is the two quark overlap and @xmath54 represents two - body operator for the second pair . the interacting pair number is @xmath55 . for @xmath7 interaction , @xmath56 is used , @xmath57[\nu_4][\tilde{\nu}_4]}_{[1 ^ 3][\nu_3][\tilde{\nu}_3],[1][\nu_1][\tilde{\nu}_1 ] } c^{[\tilde{\nu}_4][c_4][\mu_4]}_{[\tilde{\nu}_3][c_3][\mu_3],[\tilde{\nu}_1][c_1][\mu_1 ] } c^{[\mu_4][f_4]j_4}_{[\mu_3][f_3]j_3,[\mu_1][f_1]j_1 } c^{[f_4]y_4i_4}_{[f_3]y_3i_3,[f_1]y_1i_1 }     \nonumber\\ & & c^{[1 ^ 4][\nu_4][\tilde{\nu}_4]}_{[1 ^ 3][\nu'_3][\tilde{\nu}'_3],[1][\nu'_1][\tilde{\nu}'_1 ] } c^{[\tilde{\nu}_4][c_4][\mu_4]}_{[\tilde{\nu}'_3][c'_3][\mu'_3],[\tilde{\nu}'_1][c'_1][\mu'_1 ] } c^{[\mu_4][f_4]j_4}_{[\mu'_3][f'_3]j'_3,[\mu'_1][f'_1]j'_1 } c^{[f_4]y_4i_4}_{[f'_3]y'_3i'_3,[f'_1]y'_1i'_1 }      \nonumber\\ & & u(c_3c_1cc_{\bar{1}};c_4c_2 ) u(i_3i_1ii_{\bar{1}};i_4i_2)u(j_3j_1jj_{\bar{1}};j_4j_2 )   \nonumber\\ & & u(c'_3c'_1cc'_{\bar{1}};c_4c'_2 ) u(i'_3i'_1ii'_{\bar{1}};i_4i'_2)u(j'_3j'_1jj'_{\bar{1}};j_4j'_2 )   \nonumber\\ & & c^{[\nu_4]w_{x_4}}_{[\nu_3]w_{x_3},[\nu_1]w_{x_1 } } c^{[\nu_4]w_{x'_4}}_{[\nu'_3]w_{x'_3},[\nu'_1]w_{x'_1 } } \langle \alpha_3 k_3 | \alpha^{\prime}_3 k^{\prime}_3 \rangle \langle \alpha_2 k_2 |h_{45 } | \alpha^{\prime}_2 k^{\prime}_2 \rangle .\end{aligned}\ ] ] here @xmath58s are racah coefficients , @xmath59 is the three quark overlap and @xmath60 is the two body matrix element and @xmath61 represents quark - antiquark operator . the interacting pair number is @xmath62 . the calculation of every factor in above equations can be found in ref.@xcite . at last the matrix elements of 5-body hamiltonian can be obtained , @xmath63 the matrix elements of five body hamiltonian on the physical bases can be obtained from the matrix elements on the symmetry bases and the transformation coefficients . the hamiltonian of naive quark cluster model is @xcite , @xmath64    @xmath65    @xmath66 , \label{oge}\ ] ]    @xmath67    all of the symbols retain their original meaning as in the ref.(@xcite ) . the single particle orbital wavefunction in the naive quark model is as follows : @xmath68 @xmath69 @xmath70 where @xmath71 and @xmath72 are the coordinates of the reference center of each quark cluster . we choose the salamanca model as representative of this class of models . details of the model can be found in ref.@xcite . here we display only the hamiltonian , @xmath73      @xmath75 and @xmath76 are pseudo scalar and scalar meson exchange potentials . the color confinement and one - gluon - exchange potentials and the single particle orbital wavefunctions are the same as those in the naive quark model . the quark delocalization , color screening model ( qdcsm ) is an extension of the naive quark cluster model and was developed with the aim of addressing multi - quark systems @xcite . first of all , a quark - delocalization similar to the percolation of electrons in atoms is introduced to take into account the contribution of orbital excitation or the mutual distortion of hadrons in the interaction region . second , a different parametrization of the confinement interaction is assumed for the quark pairs in different states . the parametrization is an effort to account for the qcd interactions corresponding to various hidden color configurations in the multi - quark system which have not been modeled in two body interaction model . the main advantage of qdcsm is that it allows the multi - quark system to choose its most favorable configuration through its own dynamics . this is accomplished by varying the energy of the system with respect to the delocalization parameter , which is a tentative approach to take into account the self - consistency of the quark and gluon distributions in the course of hadron interaction process . this model reproduces the existing baryon - baryon interaction data well @xcite ( bound - state deuteron as well as @xmath77 , @xmath78 and @xmath79 scattering ) . it is therefore interesting to apply the qdcsm to study the pentaquark system . some generalizations are needed here : the quark can delocalize among clusters and the color confinement between quarks in different clusters is screened as before , but now the clusters may be colorful as well as colorless . we admit there should be difference of qcd vacuum between colorless hadrons and colorful ones and so color screening should be different but in order to avoid new parameters we employ the original one tentatively . the quark delocalization is realized by replacing the single particle orbital wavefunctions ( @xmath80 and @xmath81 by @xmath82 and @xmath83 ) : @xmath84 @xmath85 here , @xmath86 , @xmath87 and @xmath88 are variational parameters determined by the dynamics of the multi - quark system rather than adjustable ( fitting ) parameters . the color screeening is realized by re - parameterizing the color confinement interaction as follows : @xmath89 details of the model can be found in ref.@xcite . the adiabatic approximation is used here to do a systematic study . for each given separations @xmath17 and @xmath18 , the energy of 5-quark state is calculated . ( for qdcsm , the energy for given @xmath17 and @xmath18 is obtained by varying the energy with the delocalization parameters . ) if there exists minimum energy at finite separations @xmath90 and @xmath91 , then the energy @xmath92 is taken as the mass of the state . the calculated transformation coefficients are listed in table a1 of appendix . ( the index of diquark cluster is given in table a2 . ) in this table , all the channels are included , so it can be used to expand the physical bases in terms of the symmetry bases , and vice versa . this table can be used for any quark model hamiltonian but restricted to jw @xmath0 cluster configuration . for meson - baryon cluster configuration one needs other transformation coefficients which will be given in a company paper . all possible states within the @xmath93 , @xmath94 , @xmath95 three - flavor world have been calculated . both single - channel and channel coupling calculations have been carried out with three quark models : the naive quark model , the chiral quark model and the qdcsm . the results are listed in table 1 . because the parity is a good quantum number , so the results for positive and negative parity states are listed separately . to save space , only the lowest single channel ( sc ) and channel coupling ( cc ) results are given . all the states with @xmath96 but only several states with @xmath97 are given . because most states with @xmath97 have quite similar features to the ones with @xmath96 part . some general features are listed below . \(1 ) the parity of the lowest channel is negative in all three quark models , which is different from jaffe - wilczek s estimation . the diquark with orbital , color , spin , flavor symmetry : [ 2 ] , [ 11 ] , [ 11 ] , [ 11 ] does have the lowest energy under the color - magnetic interaction , however the pauli principle excludes the @xmath17-wave orbital motion between two such diquarks . the @xmath98-wave excitation and the  residue " interaction between two diquarks make the energy of jw state higher than the state where two diquarks have the symmetry : @xmath99\times [ 2]$ ] ( orbital ) , @xmath100\times [ 2]$ ] ( color ) , @xmath100\times [ 2]$ ] ( spin ) and @xmath100\times [ 11]$ ] ( flavor ) , where @xmath17-wave orbital is permitted . \(2 ) generally there exist effective attractions for both positive and negative parity states , thus forming resonance in both parity states is possible . this is due to the hidden color structure of jw configuration as discussed in the introduction part . however in most cases , the attraction is not enough to make the energy lower than the corresponding threshold , i.e. , the sum of corresponding baryon and meson mass . therefore there will be no narrow resonances except there is special mechanism to prevent the decay . fig.4 gives the effective potential for the channel - coupling result of @xmath102 in qdcsm . the other two models have similar effective potentials as qdcsm . other states have similar effective potentials as that of @xmath102 shown in fig.4 . in qdcsm , the mass of @xmath102 is 1786 mev which is too high to match the experimentally claimed 1540 mev . other two models give even higher mass of @xmath6 . for @xmath103 , the qdcsm mass is about 1884 mev , a little higher than the claimed 1862 mev . other two models mass of @xmath103 are much higher ( 1974 mev for chiral quark model and 2145 for naive quark model ) . \(3 ) the spectroscopy of qdcsm and chiral quark model are quite similar except a mass shift . this result is quite consistent with our former work on @xmath77 interaction @xcite , where the two models with quite different mechanism of intermediate range attraction give similar results of deuteron properties and @xmath77 scattering phase shifts , and shows the @xmath104 meson effect can be replaced by qdcsm mechanism for pentaquark systems as well . however the naive qauark model results are quite different . the naive quark model is quite possible not a realistic one for multi - quark system even though it is a good model for single hadrons because it can not reproduce the intermediate range @xmath77 attraction and the deuteron properties . \(4 ) in qdcsm , the majority of the lowest energy states have the triangle pattern : the separations between the two diquarks are @xmath105 fm and the separation between the anti - quark and the center of the two diquarks is @xmath106 fm except few positive parity states . while , in the chiral quark model and naive quark model , the lowest energy states always have the linear pattern : the separations between the two diquarks are @xmath107 fm and the separation between the anti - quark to the diquark center is 0.1fm . we do nt have any experimental data to check which one is more realistic , lattice qcd calculation might provide useful information . r. l. jaffe , phys . lett . * 38 * , 195 ( 1977 ) . r. bilger , h. a. clement and m. g. schepkin , phys . * 71 * , 42 ( 1993 ) . t. goldman , k. maltman , g.j . stephenson , k.e . schmidt and f. wang , phys . rev . * c39 * , 1889 ( 1989 ) ; j.l . ping , f. wang and t. goldman , nucl . a688 * , 871 ( 2001 ) ; j.l . ping , h.r . pang , f. wang and t. goldman , phys . rev . * c65 * , 044003 ( 2002 ) . t. goldman , k. maltman , g. j. stephenson , k.e . schmidt and f. wang , phys . lett . * 59 * , 627 ( 1987 ) ; h.r . pang , j.l.ping , f. wang , t. goldman and e.g. zhao , phys . rev . * c69 * , 065207 ( 2004 ) ; h.r . pang , j.l.ping , l.z . chen , f. wang and t. goldman , phys . rev . * c70 * , 035201 ( 2004 ) . kopeliovich , b. schwesinger and b. e. stern , phys . lett . * b242 * , 145 ( 1990 ) ; q.b . li , p. n. shen z.y . zhang and y.w . yu , nucl . * a683 * , 487 ( 2001 ) ; h.r . pang , j.l . ping , f. wang and t. goldman , phys . rev . * c66 * , 025201 ( 2002 ) . t. nakano , d. s. ahn , j. k. ahn @xmath108 @xmath109 , phys . 91 * , 012002 ( 2003 ) . k. hicks , prog . part . * 55 * , 647 ( 2005 ) and references therein ; v. d. burkert , int . phys . * a21 * , 1764 ( 2006 ) . m. oka , prog . . phys . * 112 * , 1 ( 2004 ) ; k. goeke , _ et al . . phys . * 55 * , 350 ( 2005 ) and references therein ; s.l . zhu , int . phys . * a20*,1548 ( 2005 ) and references therein . ping , d. qing , f. wang and t. goldman , phys . b602 * , 197 ( 2004 ) . r. jaffe and f. wilczek , phys . * 91 * , 232003(2003 ) . m. karliner , h.j . lipkin , phys . lett . * b575 * , 249 ( 2003 ) . x. c. song and s. l. zhu , mod.phys.lett . * a19 * , 2791 ( 2004 ) . clas collaboration , phys . . lett . * 97 * , 032001 ( 2006 ) ; * 96 * , 042001 ( 2006 ) ; * 96 * , 212001 ( 2006 ) ; * 97 * , 102001 ( 2006 ) ; phys . * d74 * , 032001 ( 2006 ) ; j. m. link _ et al_. , [ focus collaboration ] , phys . lett . * b639 * , 604 ( 2006 ) ; cosy - tof collaboration , hep - ex/0612048 . b. wu and b.q . ma , phys . rev . * d69 * , 077501 ( 2004 ) ; q.h . zhou and b.q . j. * c49 * , 897 ( 2007 ) . b. wu and b.q . ma , phys . rev . * d70 * , 094042 ( 2004 ) . r. bijker , _ et al _ , eur . j. * a22 * , 319 ( 2004 ) . d. qing , x.s . chen and f. wang , phys . c57 * , r31 ( 1998 ) ; phys . rev . * d58 * , 114032 ( 1998 ) . zou and d.o . riska , phys . * 95 * , 072001 ( 2005 ) . b. aubert , _ et al _ , phys . lett . * 90 * , 242001 ( 2003 ) ; s.k . choi , _ et al _ , phys . lett . * 91 * , 262001 ( 2003 ) ; a.v . evdokimov , _ et al _ ( selex collaboration ) phys . 93 * , 242001 ( 2004 ) ; b. aubert , _ et al _ , phys . lett . * 95 * , 142001 ( 2005 ) ; r.l . jaffe , phys . report * 409*,1 ( 2005 ) . f. okiharu , h. suganuma and t.t . takahashi , phys . rev . * d72 * , 014505 ( 2005 ) . paris and v.r . pandharipande , phys . rev . * c62 * , 015201 ( 2000 ) e. santopinto and g. galata , phys . rev . * c75 * , 045206 ( 2007 ) . f. wang , j. l. ping , and t. goldman , phys . rev . * c51 * , 1648 ( 1995 ) . f. wang , j.l . ping , g.h . teng and t. goldman , phys . * c51 * ( 1995 ) 3411 . t. goldman , _ et al _ , mod . a13 * , 59 ( 1998 ) ; h.r . pang , j.l . ping , f. wang and t. goldman , phys . * c65 * , 014003 ( 2001 ) . a. de rujula , h. georgi and s.l . glashow , phys . * d12 * , 147 ( 1975 ) ; n. isgur and g. karl , phys d18 * , 4187 ( 1978 ) ; * d19 * , 2653 ( 1979 ) ; * d20 * , 1191 ( 1979 ) . a. valcarce , h. garcilazo , f. fernandez , and p. gonzalez , rept . * 68 * ( 2005 ) 965 . f. wang , g.h . wu and l.j . teng and t. goldman , phys . lett . * 69 * , 2901 ( 1992 ) . chen , j.l . ping and f. wang , _ group representation theory for physicists _ , ( world sci . , singapore , 2002 ) . chen , x.b . wu and m.j . gao , _ tables of the su(mn)@xmath110su(m)@xmath111su(n ) coefficients of fractional parentage _ ( world sci . , singapore , 1991 ) . m. harvey , nucl . a352 * , 301 ( 1981 ) ; * 481 * 834 ( 1988 ) . j.q . chen , y.j . shi , d.h . feng and m. vallieres , nucl . * a393 * , 122 ( 1983 ) . g. h. wu , j. l. ping , l. j. teng , f. wang and t. goldman , nucl . a673 * , 279 ( 2000 ) ; j. l. ping , f. wang and t. goldman , nucl . phys . * a657 * , 95 ( 1999 ) . l. z. chen , h. r. pang , h. x. huang , j. l. ping and f. wang , phys . * c76 * , 014001 ( 2007 ) . b. wu and b.q . ma , phys . rev . * d70 * , 097503 ( 2004 ) .
Abstract:

Answer: group theoretic method for the systematic study of multi - quark states is developed . the calculation of matrix elements of many - body hamiltonian is simplified by transforming the physical bases ( quark cluster bases ) to symmetry bases ( group chain classified bases ) , where the fractional parentage expansion method can be used . five - quark system is taken as an example in this study . the jaffe - wilczek @xmath0 configuration is chosen as one of examples to construct the physical bases and the transformation coefficients between physical bases and symmetry ones are shown to be related to the @xmath1 isoscalar factors and a complete transformation coefficients table is given . the needed isoscalar factors and fractional parentage coefficients had been calculated with our new group representation theory and published before . three quark models : the naive glashow - isgur , salamanca and quark delocalization color screening , are used to show the general applicability of the new multi - quark calculation method and general results of constituent quark models for five - quark states are given .