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: spontaneous chiral symmetry breaking is one of the most important properties to understand the low - energy phenomena of qcd in the real world . this chiral symmetry is expected to be restored in qcd at several extreme conditions such as qcd with a large number of massless quarks , large @xmath0 qcd ( see , e.g. , refs . @xcite ) and qcd in hot and/or dense matter ( see , e.g. , ref . @xcite ) . in ref . @xcite , based on the infrared ( ir ) fixed point existing at a two - loop beta function for a large number of massless quarks ( @xmath1 )  @xcite , it was found through the improved ladder schwinger - dyson ( sd ) equation that chiral symmetry restoration takes place for @xmath0 such that @xmath2 , where @xmath3(@xmath4 for @xmath5 ) . then , in ref .  @xcite this chiral restoration at @xmath6 was further identified with `` conformal phase transition '' which was characterized by the essential singularity scaling . moreover , such chiral restoration is also observed by other various methods such as lattice simulation  @xcite , dispersion relation  @xcite , instanton calculus  @xcite , effective field theoretical approach  @xcite , etc ..    more attention has been paid to the property of the phase transition . especially , it is interesting to ask what are the light degrees of freedom near the phase transition point in the large @xmath0 qcd : for example , in the manifestation of the chiral symmetry restoration  la the linear sigma model , the scalar bound state becomes a chiral partner of the pseudoscalar bound state and becomes massless at the phase transition point . on the contrary , in the vector manifestation ( vm ) @xcite obtained by the effective field theoretical approach based on the hidden local symmetry model  @xcite , it is the vector bound state which becomes massless as a chiral partner of the pseudoscalar bound state . besides , from the viewpoint of the conformal phase transition  @xcite , it is natural to suppose that all the existing bound states become massless near the phase transition point when approached from the broken phase ( see ref . @xcite ) .    then , it is quite interesting to study which types of the bound states actually exist near the phase transition point , and investigate the critical behavior of their masses directly from qcd . such studies from the first principle will give us a clue to understand the nature of the chiral phase transition in large @xmath0 qcd .    a powerful tool to study the bound states of quark and anti - quark directly from qcd is the homogeneous bethe - salpeter ( hbs ) equation in the ( improved ) ladder approximation ( see , for reviews , refs . @xcite ) . when the mass of the quark is regarded as a constant , we can easily solve the hbs equation by using a so - called fictitious eigenvalue method  @xcite . however , for consistency with the chiral symmetry , the quark propagator in the hbs equation must be obtained by solving the sd equation with the same kernel as that used in the hbs equation  @xcite , and as a result , the quark mass becomes a certain momentum dependent function . then , in order to obtain the masses and the wave functions of the bound states , it is necessary to solve the hbs equation and the sd equation simultaneously . when we try to solve these two equations in real - life ( @xmath7 ) qcd , however , we encounter difficulties . first of all , for the consistency of the solution of the sd equation with qcd in a high energy region , we need to use the running coupling which obeys the evolution determined from qcd @xmath8-function in the high energy region ( see , for reviews , refs . since the running coupling diverges at some infrared scale , @xmath9 , we have to regularize the running coupling in the low energy region , for which there exist several ways ( see , e.g. , refs . even if we fix the infrared regularization in such a way that we can solve the sd equation on the real ( space - like ) axis , another problem arises when we try to solve the hbs equation . since the argument of the quark mass function in the hbs equations for the massive bound states becomes a complex quantity after the wick rotation has been made , we have to solve the sd equation on the complex plane , which requires an analytic continuation of the running coupling . several works such as in refs . @xcite proposed models of running couplings for a general complex variable which are consistent with perturbative qcd for large space - like momentum . however , they still have branch cuts on the complex plane , and it is a complicated task to obtain the solution of the sd equation for a general complex variable . one way to avoide a such complication is solving the inhomogeneous bs equation for vertex functions to obtain the current correlators in the space - like region which we can fit the mass of the relevant bound state to ( see , e.g. , refs . @xcite ) . another way might be replacing the entire running coupling with an ad hoc analytic function ( see , e.g. , ref . anyway , it is impossible to solve the sd equation on the complex plane without modeling the running coupling .    in this paper , we point out that the situation dramatically changes when we increase the number of massless quarks . when @xmath0 becomes larger than @xmath10 , the running coupling obtained from the renormalization group equation ( rge ) with two - loop approximation takes a finite value for all the range of the energy region due to the emergence of the non - trivial ir fixed point . then , we need no infrared regularization , and we do not have any ambiguities coming from the regularization scheme which do exist in the case of small @xmath0 . moreover , an explicit solution of the two - loop rge can be written in terms of the lambert @xmath11 function  @xcite , and when @xmath0 is close to @xmath6 the solution of the rge has no singularity on the complex plane except for the time - like axis  @xcite . consequently , we can solve the sd equation on the complex plane without introducing any models for the running coupling . based on these facts , we solve the hbs equations for the bound states of quark and anti - quark in large @xmath0 qcd with the improved ladder approximation in the landau gauge . the mass function for complex arguments needed in the hbs equation is obtained by solving the sd equation with the same kernel as that used in the hbs equation . we find the solution of the hbs equation in each of the scalar , vector , and axial - vector channels , which implies that the scalar , vector , and axial - vector bound states are actually formed near the phase transition point . our results show that the masses of the scalar , vector , and axial - vector bound states go to zero as the number of quarks @xmath0 approaches to its critical value @xmath6 where the chiral symmetry restoration takes place . this may suggest the existence of a new type of manifestation of chiral symmetry restoration in large @xmath0 qcd other than the linear sigma model like manifestation and a simple version of the vector manifestation proposed in ref . @xcite . this paper is organized as follows . in section  [ sec : sdeq ] we numerically solve the sd equation with an approximate form of the running coupling , and study the critical behavior of the nambu - goldstone boson decay constant . in section  [ sec : complex_sd ] we solve the sd equation for complex arguments . section  [ sec : hbseq ] is devoted to summarizing the numerical method for solving the hbs equation . section  [ sec : numerical_calculations ] is the main part of this paper . we first solve the hbs equation for the pseudoscalar bound state to show that the approximation adopted in the present analysis is consistent with the chiral symmetry . we next solve the hbs equation for the scalar , vector , and axial - vector bound states to obtain their masses . finally we give a summary and discussion in section  [ sec : discussion_and_conclusion ] . in appendix [ app : positronium ] we solve the hbs equation for the ortho - positronium with a constant electron mass to show the validity of the fictitious eigenvalue method . the bispinor bases for the bound states are listed in appendix  [ app : bispinor - bases ] . in appendix [ app : decay_const ] we calculate the coupling constants @xmath12 , @xmath13 , and @xmath14 of the vector , axial - vector , and scalar bound states to the vector current , axial - vector current , and scalar density . we briefly study numerical uncertainties in the present analysis in appendix  [ app : uncertainties ] . in this section we numerically solve the schwinger - dyson ( sd ) equation for the quark propagator with the improved ladder approximation in the landau gauge , and show the critical behaviors of the dynamical mass and the decay constant of the nambu - goldstone boson . we also show the behavior of the fermion - antifermion pair condensate @xmath15 near the phase transition point . schwinger - dyson ( sd ) equation is a powerful tool to study the dynamical generation of the fermion mass directly from qcd ( for reviews , see , e.g. , refs . @xcite ) . the sd equation for the full fermion propagator @xmath16 in the improved ladder approximation  @xcite is given by ( see fig .  [ fig : sdeq ] for a graphical expression ) @xmath17 where @xmath18 is the second casimir invariant , and @xmath19 is the running coupling . the landau gauge is adopted for the gauge boson propagator . the sd equation provides coupled equations for two functions @xmath20 and @xmath21 in the full fermion propagator @xmath22 . when we use a simple ansatz for the running coupling , @xmath23  @xcite , with @xmath24 being the euclidean momenta , we can carry out the angular integration and get @xmath25 in the landau gauge . then the sd equation becomes a self - consistent equation for the mass function @xmath26 . the resultant asymptotic behavior of the dynamical mass @xmath27 is shown to coincide with that obtained by the operator product expansion technique  @xcite . however , it was shown in ref . @xcite that the axial ward - takahashi identity is violated in the improved ladder approximation unless the gluon momentum is used as the argument of the running coupling as @xmath28 . in this choice we can not carry out the angle integration analytically since the running coupling depends on the angle factor @xmath29 . furthermore , we would need to introduce a nonlocal gauge fixing  @xcite to preserve the condition @xmath30 .    in ref . @xcite , however , it was shown that an angle averaged form @xmath31 gives a good approximation . then , in the present analysis we take the argument of the running coupling as @xmath32 after applying this angle approximation and carrying out the angular integration , we can show ( see , e.g. , refs . @xcite ) that @xmath20 always satisfies @xmath33 in the landau gauge . then the sd equation becomes @xmath34 where @xmath35 and @xmath36 . although the choice of arguments in eq . ( [ eq : approx_of_argument ] ) explicitly breaks the chiral symmetry as mentioned above , it will be shown later that the magnitude of the breaking is negligible .      in qcd with @xmath0 flavors of massless quarks , the renormalization group equation ( rge ) for the running coupling @xmath37 in the two - loop approximation is given by @xmath38 where @xmath39 from the above beta function we can easily see that , when @xmath40 and @xmath41 , i.e. , @xmath0 takes a value in the range of @xmath42 ( @xmath43 for @xmath44 ) , the theory is asymptotically free and the beta function has a zero , corresponding to an infrared stable fixed point  @xcite , at @xmath45    existence of the infrared fixed point implies that the running coupling takes a finite value even in the low energy region . actually , the solution of the two loop rge in eq . ( [ eq : rge_for_alpha ] ) can be explicitly written  @xcite in all the energy region as @xmath46^{-1 } , \label{lan w}\ ] ] where @xmath47 with @xmath48 is the lambert @xmath11 function , and @xmath49 is a renormalization group invariant scale defined by  @xcite @xmath50 . \label{eq : lambda}\ ] ] we note that , in the present analysis , we fix the value of @xmath49 to compare the theories with a different number of flavors , and that we have no adjustable parameters in the running coupling in eq . ( [ lan w ] ) , accordingly ( see discussion below ) . we show an example of @xmath51 for @xmath52 by the solid line in fig . [ fig : step ] .    ) ( dashed line ) for @xmath52 . , height=226 ]    the fact that the running coupling is expressed by a certain function as in eq . ( [ lan w ] ) implies that , in the case of large @xmath0 qcd , _ we do not need to introduce any infrared regularizations _ such as the ones adopted in refs . @xcite for studying real - life qcd with small @xmath0 in which the infrared regularization parameter must be chosen in such a way that the running coupling in the infrared region becomes larger than the critical value @xmath53 for realizing the dynamical chiral symmetry breaking  @xcite . the running coupling in large @xmath0 qcd takes a certain value in the infrared region for given @xmath0 , so that we can definitely determine , within the framework of the sd equation , whether or not the dynamical chiral symmetry breaking is realized . actually , the value of @xmath54 decreases monotonically with increasing @xmath0 , and the chiral symmetry restores when @xmath0 becomes large enough . in refs . @xcite , it was shown that the phase transition occurs at @xmath55 for @xmath44 ( corresponding to @xmath56 ) .    in order to reduce the task of numerical calculations in solving the hbs equation , we modify the shape of the running coupling . since the dynamics in the infrared region governs the chiral symmetry breaking , we adopt the following approximation for the running coupling  @xcite : @xmath57 in this approximation the coupling takes the constant value @xmath54 ( the value at the ir fixed point ) below the scale @xmath49 and entirely vanishes in the energy region above this scale . the dashed line in fig .  [ fig : step ] represents the approximated form of the running coupling for @xmath58 .      in this subsection we briefly explain how we solve the sd equation numerically . we first introduce the infrared ( ir ) cutoff @xmath59 and ultraviolet ( uv ) cutoff @xmath60 as @xmath61 then , we discretize the momentum variable @xmath62 and @xmath63 into @xmath64 points as @xmath65 \ , \ \ \     \big ( i = 0 , 1 , 2 , \cdots , ( n_{sd}-1 ) \big ) \ , \ ] ] where @xmath66 accordingly , the integration over @xmath63 is replaced with a summation as @xmath67 then , the sd equation in eq . ( [ eq : improved_ladder_sd ] ) with the running coupling in eq . ( [ eq : step ] ) is rewritten as @xmath68 this discretized version of the sd equation is solved by the recursion relation : @xmath69 starting from a suitable initial condition ( we choose @xmath70 ) , we update the mass function by the above recursion relation . then , we stop the iteration when the convergence condition @xmath71 ^ 2     \ < \   \varepsilon^2 \lambda^6\ ] ] is satisfied for sufficiently small @xmath72 , and regard this @xmath73 as a solution of eq . ( [ eq : discretized_largenf_sd_1 ] ) . in fig . [ fig : mass_function ] , we show the numerical solution for the mass function @xmath74 .    ) for @xmath52 . , height=226 ]    here , we took @xmath75 ( @xmath76 ) as an example and used the following parameters : @xmath77    now , let us study the critical behavior of the fermion mass as @xmath0 is varied . note that we can use @xmath54 instead of @xmath0 as an input parameter , because once we choose a value of @xmath0 , the value of @xmath54 is uniquely determined from eq . ( [ eq : alpha_ir ] ) . for example , @xmath78 implies @xmath79 and @xmath80 implies @xmath81 . we solve eq . ( [ eq : discretized_largenf_sd_1 ] ) for various values of @xmath54 and plot the values of @xmath82 in fig . [ fig : sigma_ir ] .     for several values of @xmath54 ( indicated by @xmath83 ) . solid line shows the function in eq . ( [ eq : fitting ] ) with the best fitted value @xmath84 . , height=264 ]    here , @xmath85 represents the dynamical mass defied by @xmath86 . it should be noticed that @xmath87 is defined in the space - like region which does not represent the pole mass of fermion . as we will show in section  [ sec : complex_sd ] , the present fermion propagator does not have any poles and then there are no pole masses of fermion . we compare this result with the analytic solution  @xcite : @xmath88 in the above form there is an ambiguity in the prefactor . then , we introduce the function @xmath89 and fit the value of the pre - factor @xmath90 by minimizing @xmath91 in the range of @xmath92 $ ] . the resultant best fitted value of @xmath90 is @xmath93 we plot the function in eq . ( [ eq : fitting ] ) with the best fitted value @xmath84 in fig . [ fig : sigma_ir ] ( solid line ) . this clearly shows that the @xmath54-dependence of the resultant @xmath94 from our numerical calculation is consistent with the analytic result : the dynamical mass function vanishes when @xmath54 reaches the critical value @xmath53 . noting that decreasing @xmath54 corresponds to increasing @xmath0 for fixed @xmath95 as we discussed in the previous subsection , we see that the chiral symmetry restoration actually occurs at @xmath96  @xcite . so far we have solved the sd equation and obtained the mass functions for the various values of @xmath54 . now we can calculate the pseudoscalar meson decay constant @xmath97 at each @xmath54 by using the pagels - stokar formula @xcite : @xmath98\ }         { \left(\   x \ + \ \sigma^2(x ) \ \right)^2 } .\ ] ] in fig . [ fig : f_pi ] , we plot the values of @xmath97 for several values of @xmath54 ( indicated by @xmath83 ) . to study the critical behavior of the pseudoscalar meson decay constant we use the function of the form in eq . ( [ eq : fitting ] ) and fit the value of @xmath90 by minimizing @xmath99 for @xmath92 $ ] . the resultant best fitted value of @xmath90 is @xmath100 we plot the fitting function with @xmath101 in fig . [ fig : f_pi ] ( dotted line ) . this shows that the results of the numerical calculations for @xmath97 are well fitted by the function of the form in eq . ( [ eq : fitting ] ) , and that the pseudoscalar meson decay constant has the same critical behavior as the mass function has . calculated from the pagels - stokar formula for several values of @xmath54 ( indicated by @xmath83 ) . the dotted line shows the function of the form in eq . ( [ eq : fitting ] ) with the best fitted value @xmath101 . , height=264 ]      in this subsection , we calculate the fermion - antifermion pair condensate @xmath102 in large @xmath0 qcd , and show that the system in the present analysis has the large anomalous dimension @xmath103 for the operator @xmath104 . we also show that the values of @xmath103 are not affected so much by the approximation for the running coupling used in the present analysis . the condensate @xmath102 is calculated from the following equation : @xmath105 where @xmath74 is the mass function obtained from the sd equation and @xmath106 represents uv cutoff introduced to regularize the uv divergence . in the improved ladder approximation , the high - energy behavior of the mass function is consistent with that derived using the operator product expansion ( ope ) . the chiral condensate calculated using the mass function was shown to obey the renormalization group equation derived with the ope ( see , e.g. , refs . then , as was adopted in refs.@xcite , we identify the condensate , which is calculated with uv cutoff @xmath106 , with that renormalized at the scale @xmath106 in qcd . ) , the integration in eq . ( [ eq : q - bar - q ] ) is effectively cut off at the scale of @xmath49 due to the truncation of the running coupling for any values of @xmath107 . ( see fig . [ fig : mass_function ] : @xmath108 for @xmath109 . ) ]    we expect that infrared dynamics in large @xmath0 qcd is similar to that of strong coupling qed or walking gauge theories  @xcite since the running coupling in large @xmath0 qcd is well approximated by the constant coupling ( see fig . [ fig : step ] )  @xcite . then , we also expect that the value of the anomalous dimension in large @xmath0 qcd becomes @xmath110 since the walking gauge theories have @xmath110  @xcite . when a considering system has the anomalous dimension @xmath103 , scaling properties of @xmath97 and @xmath111 with respect to @xmath54 near the critical point are expressed as follows  @xcite : @xmath112 @xmath113 where @xmath87 represents the dynamical fermion mass . these equations mean that the relation between @xmath102 and @xmath97 can be written as @xmath114 where @xmath115 is a certain positive constant . from this equation , we can express the anomalous dimension as @xmath116 where @xmath117 here , we note that @xmath118 approaches @xmath103 for @xmath119 since @xmath97 becomes small , i.e. , @xmath120 , near the critical point ( see fig .  [ fig : f_pi ] ) : @xmath121    in fig . [ fig : q - bar - q_over_fpi2 ] , we plot the values of @xmath118 for several values of @xmath54 as an estimation of the anomalous dimension .     for several values of @xmath54 . @xmath83 and @xmath122 represent the values of @xmath118 calculated from eq . ( [ eq : gamma_m ] ) with and without approximation for the running coupling . ( we call them @xmath123 and @xmath124 respectively . ) , height=264 ]    the data indicated by @xmath83 in fig .  [ fig : q - bar - q_over_fpi2 ] is obtained with the approximated running coupling ( dashed line in fig . [ fig : step ] ) in the sd equation . ( we call this kind of data @xmath125 . ) on the other hand , the data indicated by @xmath122 is the result from the calculation with the two - loop running coupling given in eq . ( [ lan w ] ) . ( we call this kind of data @xmath126.)[multiblock footnote omitted ] from these results , we conclude that large @xmath0 qcd with two - loop running coupling as well as with approximated running coupling actually possesses obtained by fitting to the data of @xmath102 , we find @xmath127 for the approximated running coupling , and @xmath128 for the two - loop running coupling . ] moreover , fig . [ fig : q - bar - q_over_fpi2 ] shows that the data of @xmath123 is in good agreement with that of @xmath124 , which implies that the approximation for the running coupling used in the present analysis works well . we also expect that the approximation does not affect the results so much when we calculate the hbs equations for the bound states . as we will discuss in section  [ sec : hbseq ] , we need the mass function for complex arguments when we solve the hbs equation for the massive bound state . in this section , we first introduce the sd equation for the complex argument following ref . @xcite ( see also ref . @xcite ) , and then solve it in the case of large @xmath0 qcd . the sd equation for the complex argument is expressed as  @xcite @xmath130                 \frac{\bar{g}^2(x+y)\ ,                 \sigma ( y)}{y + { \sigma^2 ( y ) } } \ , \label{eq : complex_sd}\ ] ] where @xmath131 is the integral path on the complex plane . here , we took the same argument of the running coupling as that in eq . ( [ eq : approx_of_argument ] ) , and carried out the angle integration . note that the integral path @xmath131 must be taken so as to avoid the branch cut appearing in the integral . we first study the structure of the running coupling appearing in the sd equation ( [ eq : complex_sd ] ) to clarify the branch cut . in the improved ladder approximation it is essential to use the running coupling determined from the @xmath8-function in the high energy ( space - like ) region for consistency with perturbative qcd . in qcd with small @xmath0 , however , the running coupling obtained from the perturbative @xmath8-function diverges at some infrared scale , @xmath9 . in the ordinary sd equation in the space - like region , the infrared singularity is avoided by introducing infrared regularization such as the so - called higashijima - miransky approximation  @xcite and its extension as in ref . however , since the running coupling in eq . ( [ eq : complex_sd ] ) is a complex function which has the complex argument , we need an extension with analyticity satisfied . several works such as in refs . @xcite proposed models of running coupling which are consistent with perturbative qcd in the high energy region . but they still have branch cuts on the complex plane , and it is a burdensome task to evade all the branch cuts by carefully selecting the integral path in eq . ( [ eq : complex_sd ] ) . one way to avoide such a complication might be to give up the consistency with perturbative qcd and use models of running coupling with analyticity such as the one used in ref .  @xcite . here we point out that the situation dramatically changes in the large @xmath0 qcd . in the case of large @xmath0 qcd , as we explained in subsection  [ sec : large_nf ] , the running coupling , as well as the two - loop @xmath8-function , is finite for any space - like momentum . this implies that we may be able to construct the running coupling by analytic continuation using the @xmath8-function . actually , an explicit solution of the two - loop renormalization group equation ( rge ) can be written in terms of the lambert @xmath11 function  @xcite . when @xmath0 is close to @xmath6 , the solution of the rge has no singularity on the complex plane except for the time - like axis  @xcite .    as a result , for general complex @xmath62 except on the time - like axis ( @xmath132 ) , we can take the integral path @xmath131 in such a way that it just avoids the branch cut coming from the angle integration . in fig . [ fig : integralpath ] we show the branch cut together with a simple choice of the integral path  @xcite .    ) . here , the branch cut appears from the four - dimensional angle integration . , height=226 ]    we stress again that the reason why we can take this simple integral path is that the running coupling has no singularity on the complex plane except for the time - like axis .    for solving the sd equation on the complex plane , we here study the explicit form of the running coupling . in fig . [ fig : analytic_cont ] we show the real part of the running coupling on the complex plane which is obtained by performing the analytic continuation from the running coupling on the real axis determined from the two - loop @xmath8-function .     on the complex plane obtained by the analytic continuation from the running coupling on the real axis ( we use the cauchy - rieman relation ) . the complex argument of @xmath133 is expressed as @xmath134 , where @xmath135 and @xmath136 are real . note that @xmath63 is in the space - like region for @xmath137 and in the time - like region for @xmath138 . , height=226 ]    this figure shows that @xmath139 , i.e. , @xmath140 , in the range of @xmath141 , and that @xmath142 in the range of @xmath143 . thus we take the following approximation for the running coupling on the complex plane : @xmath144 which is smoothly connected to the approximation adopted in eq . ( [ eq : step ] ) for the running coupling on the space - like axis . now , let us solve the sd equation ( [ eq : complex_sd ] ) to obtain the mass function for complex variable @xmath62 . along the integral path shown in fig . [ fig : integralpath ] , the variables @xmath62 and @xmath63 are parametrized as @xmath145 where @xmath146 , @xmath135 , and @xmath136 are real . then the sd equation ( [ eq : complex_sd ] ) is rewritten as @xmath147                 \frac{\bar{g}^2 ( x+y ) \ ,    \sigma ( y e^{i\theta})}{y + e^{-i\theta}\,\sigma^2 ( y e^{i\theta } ) } \ . \label{eq : complex_sd 2}\ ] ] from this we can easily see that the solution is expressed as @xmath148 where @xmath149 is real and satisfies the original sd equation on the real axis : @xmath150                      \frac{\bar{g}^2 ( x+y ) \ ,                      \bar\sigma ( y)}{y + { \bar\sigma ( y)}^2 }   \ . \label{eq : real sd}\ ] ] note that the fermion propagator @xmath151 does not have any poles since the kinetic part @xmath62 and the mass part @xmath152 have the same phases as @xmath153 [ see eq . ( [ eq : sigma asum ] ) ] and the mass function in the space - like region satisfies @xmath154 . we should note that the above solution in eq . ( [ eq : sigma asum ] ) is a double - valued function on the complex plane : the variable @xmath155 can be parametrized as @xmath156 for which the solution takes @xmath157 . this corresponds to the fact that the sd equation has two solutions : when @xmath74 is a solution , @xmath158 also satisfies the equation . when we choose the range of @xmath136 as @xmath159 $ ] , the branch cut emerges on the time - like axis . this choice is natural because the appearance of the branch cut in the time - like region seems consistent with the analytic structure of the running coupling . we will see that this branch cut does not matter in calculating the bound state masses . in this section we briefly review the homogeneous bethe - salpeter ( hbs ) equation for the bound states of quark and antiquark and show how to solve it numerically .      in this paper , we consider the scalar , pseudoscalar , vector , and axial - vector bound states of quark and antiquark , and we write these bound states as @xmath160 , @xmath161 , @xmath162 , and @xmath163 , respectively . here , @xmath164 represents the momentum of the bound states and @xmath165 represents the polarization vector satisfying @xmath166 and @xmath167 . now , we introduce the bethe - salpeter ( bs ) amplitudes @xmath168 for the bound states of quark and anti - quark as follows : @xmath169_{\alpha\beta},\ ] ] @xmath170_{\alpha\beta},\ ] ] @xmath171_{\alpha\beta},\ ] ] @xmath172_{\alpha \beta},\ ] ] where @xmath173 , @xmath174 is the generator of @xmath175 normalized as @xmath176 = 2 \delta_{ab}$ ] , and ( @xmath133 , @xmath8 ) , ( @xmath177 , @xmath178 ) , and ( @xmath179 , @xmath180 ) denote the spinor , flavor , and color indices , respectively .    we can expand the bs amplitude @xmath168 in terms of the bispinor bases @xmath181 and the invariant amplitudes @xmath182 as follows : @xmath183_{\alpha \beta }   = \ \sum_{i=1}^{4 } \left [ \gamma^i_{(s , p)}(p;q ) \right]_{\alpha \beta }    \chi_{(s , p)}^{i } ( p;q ) , \label{eq : chi - sp}\ ] ] @xmath184_{\alpha \beta }   = \ \sum_{i=1}^{8 } \left [ \gamma^i_{(v , a)}(p;q,\epsilon ) \right]_{\alpha \beta }    \chi_{(v , a)}^{i } ( p;q ) . \label{eq : chi - va}\ ] ] the bispinor bases can be determined from spin , parity , and charge conjugation properties of the bound states . the explicit forms of @xmath185 , @xmath186 , @xmath187 , and @xmath188 are summarized in appendix  [ app : bispinor - bases ] . we take the rest frame of the bound state as a frame of reference : @xmath189 where @xmath190 represents the bound state mass . after the wick rotation , we parametrize @xmath191 by the real variables @xmath192 and @xmath62 as @xmath193 consequently , the invariant amplitudes @xmath182 become functions in @xmath192 and @xmath62 : @xmath194 from the charge conjugation properties for the bs amplitude @xmath168 and the bispinor bases defined in appendix  [ app : bispinor - bases ] , the invariant amplitudes @xmath195 are shown to satisfy the following relation : @xmath196      the hbs equation is the self - consistent equation for the bs amplitude ( see , for a review , ref . @xcite ) , and it is expressed as ( see fig . [ fig : hbseq ] )        @xmath197    the kinetic part @xmath198 is given by @xmath199 where @xmath151 is the full fermion propagator @xmath200 , and the bs kernel @xmath201 in the improved ladder approximation is expressed as @xmath202 in the above expressions we used the tensor product notation @xmath203 and the inner product notation @xmath204    it should be noticed that the fermion propagators included in @xmath198 in eq . ( [ def t ] ) have complex - valued arguments after the wick rotation . the arguments of the mass functions appearing in two legs of the bs amplitude are expressed as @xmath205 in general , it is difficult to obtain mass functions for complex arguments . however , as we have shown in section  [ sec : complex_sd ] , it is easy to obtain them in the case of large @xmath0 qcd . we now modify eq . ( [ eq : hbseq ] ) so that we can solve it numerically . it is convenient to introduce the conjugate bispinor bases defined by @xmath206 multiplying these conjugate bispinor bases from the left , taking the trace of spinor indices , and summing over the polarizations , we rewrite eq . ( [ eq : hbseq ] ) into the following form : @xmath207 where the summation over the index @xmath180 is understood , and @xmath208 , \\    k_{ij}(u , x;v , y ) & = & \int_{-1}^{1 } d\cos\theta \ \sum_{\epsilon }       \frac{1}{3 }           \mbox{tr } \left [ \bar\gamma^i(p;q,\epsilon )      k(p , k ) \gamma^j(k;q,\epsilon ) \right ] \ , \end{aligned}\ ] ] with the real variables @xmath209 and @xmath63 introduced as @xmath210 and @xmath211 . here , we note that although the mass function @xmath74 has the branch cut on the time - like axis as mentioned in section [ sec : complex_sd ] , @xmath212 has no singularity and becomes a continuous function for all the range of @xmath192 and @xmath62 . as for @xmath213 , the branch cut of running coupling @xmath214 does not matter since its argument @xmath215 never takes a negative value .    using the property of @xmath182 in eq . ( [ eq : even_chi ] ) , we restrict the integration range as @xmath216 : @xmath217 \chi^j(v , y).\ ] ] then , all the variables @xmath192 , @xmath62 , @xmath209 , and @xmath63 can be treated as positive values .    to discretize the variables @xmath192 , @xmath62 , @xmath209 , and @xmath63 we introduce new variables @xmath218 , @xmath146 , @xmath219 , and @xmath135 as @xmath220 and set uv and ir cutoffs as @xmath221 , \ \     x , y \ \in \   [ \lambda_x,\lambda_x ] .\ ] ] we discretize the variables @xmath218 and @xmath219 into @xmath222 points evenly , and @xmath146 and @xmath135 into @xmath223 points . then , the original variables are labeled as @xmath224 } = \lambda                     \exp\left[\lambda_u/\lambda + d_u i_u \right ] , \ \ x_{[i_x ] } = \lambda                     \exp\left[\lambda_x/\lambda + d_x i_x \right ] , \nonumber\\    & & v_{[i_v ] } = \lambda                     \exp\left[\lambda_u/\lambda + d_u i_v \right ] , \ \         y_{[i_y ] } = \lambda                    \exp\left[\lambda_x/\lambda + d_x i_y \right ] , \nonumber\end{aligned}\ ] ] where @xmath225 and @xmath226 . the measures @xmath227 and @xmath228 are defined as @xmath229 as a result , the integration is converted into the summation : @xmath230 in order to avoid integrable singularities in the kernel @xmath231 at @xmath232 , we adopt the following four - splitting prescription  @xcite : @xmath233 , \end{aligned}\ ] ] where @xmath234 , \ \ y_\pm = \exp\left[y \pm \frac{d_x}{4}\right ] .\ ] ]      now that all the variables have become discrete and the original integral equation ( [ eq : hbseq ] ) turned into a linear algebraic one , we are able to deal it numerically . however , it is difficult to find the bound state mass @xmath190 and the corresponding bs amplitude @xmath168 directly since the hbs equation depends on @xmath190 nonlinearly . a way which enables us to solve the nonlinear eigenvalue problem is the _ fictitious eigenvalue method _  @xcite . in this method we introduce a fictitious eigenvalue @xmath235 and interpret the hbs equation  ( [ eq : hbseq ] ) as a linear eigenvalue equation for a given bound state mass @xmath190 : @xmath236 consequently , the hbs equation turns into an ordinary eigenvalue problem which we can solve by standard algebraic techniques . by adjusting an input mass @xmath190 such that an eigenvalue @xmath235 equals unity , we obtain the bound state mass and the corresponding bs amplitude as a solution of the original hbs equation  ( [ eq : hbseq ] ) . in appendix [ app : positronium ] , to show the validity of this method , we calculate the mass of the positronium using this method . in this section we show the results of our numerical analysis .      as discussed in subsection  [ sec : laddersd ] , the approximation to the argument of the running coupling in eq . ( [ eq : approx_of_argument ] ) breaks the chiral symmetry explicitly @xcite . so , before solving the hbs equation for the massive bound states , we solve that for the pseudoscalar bound state and see how much the chiral symmetry is explicitly broken by this approximation . the mass of the lowest - lying pseudoscalar bound state should become zero because it appears as a nambu - goldstone boson when the chiral symmetry is spontaneously broken . so , we substitute zero for the bound state mass and check whether the fictitious eigenvalue @xmath235 becomes unity . we use the following parameters for the calculations : @xmath237        =    [ \ -18.0,\ 0\ ] , \ \ [ \ \lambda_x , \ \lambda_x \ ]      =    [ \ -8.5,\ 0\ ] , \ ] ] @xmath238 we calculate the fictitious eigenvalues for several values of @xmath54 and show them in table  [ tab : pseudoscalar ] .    . fictitious eigenvalues obtained by solving eq . ( [ eq : fictitious hbs ] ) for the pseudoscalar bound state with zero mass used as an input . [ cols="^,^,^,^",options="header " , ]     we can see that @xmath239 is satisfied within @xmath240 uncertainty . this implies that our calculations actually reproduce the massless nambu - goldstone boson within the numerical error , and that the effect of explicit chiral symmetry breaking caused by the approximation for the running coupling is negligible .      in this subsection we show the results of the numerical calculations for the masses of the vector , axial - vector , and scalar bound states . for the uv and ir cutoffs we adjust the values of them in such a way that the dominant supports of the integrands of the decay constant in eq . ( [ eq : f_v ] ) and the normalization condition in eq . ( [ eq : normalization ] ) lie in the energy region between the uv and ir cutoffs . as an example , we show the integrands of the decay constant and the normalization condition for the vector bound states in appendix  [ app : uncertainties ] . from these figures , the dominant supports lie in the lower energy region for smaller value of @xmath54 . then , we use the following @xmath54-dependent uv and ir cutoffs for the vector and the axial - vector bound states : @xmath241        & = &    [ \ -12.0 + 22.0 \times ( \alpha_\ast - 1.0 )           \ , \ -1.0 + 35.0 \times ( \alpha_\ast - 1.0 ) \ ] ,   \label{eq : ir cutoffs } \\    \left[\ \lambda_x , \ \lambda_x \ \right ]      & = &    [ \ -5.0 + 22.0 \times ( \alpha_\ast - 1.0 )          \ , \ -2.0 + 20.0 \times ( \alpha_\ast - 1.0 ) \ ] . \label{eq : uv cutoffs}\end{aligned}\ ] ] for the scalar bound state , on the other hand , we use the following fixed uv and ir cutoffs : @xmath237        =    [ \ -18.0,\ 0\ ] , \ \ [ \ \lambda_x , \ \lambda_x \ ]      =    [ \ -10.0,\ 0\ ] .\ ] ] although the integrands of the normalization conditions are shown in appendix  [ app : uncertainties ] only for the vector bound states , we have checked that the dominant supports always lie within the energy region between uv and ir cutoffs for all kinds of bound states and for all values of @xmath54 . as for the numbers of the discretization , we use @xmath242 for the vector and the axial - vector bound states , and @xmath243 for the scalar bound states . in appendix [ app : uncertainties ] , we show that these numbers of discretization are large enough for the present analysis . we should stress that we actually found a solution for eq . ( [ eq : fictitious hbs ] ) reproducing @xmath244 for all the types of the bound states in the range of @xmath92 $ ] . this means that there do exist the vector , axial - vector , and scalar bound states near the phase transition point in the broken phase . ( or , equivalently , @xmath245 ) . this fact seems consistent with the property of the conformal phase transition which has no bound states in the symmetric phase  @xcite . ]    now , let us show the critical behavior of the masses of the existing bound states . in fig . [ fig : all_masses ] , we plot all the bound state masses calculated for several values of @xmath54 together with the pseudoscalar meson masses obtained in the previous subsection .    . , height=302 ]    this figure shows that the masses of the vector , axial - vector , and scalar bound states go to zero simultaneously as the coupling approaches its critical value ( or , equivalently , @xmath246 ) : @xmath247 next , to study the critical behavior of @xmath248 , @xmath249 , and @xmath250 we use the function of the form in eq . ( [ eq : fitting ] ) and fit the value of @xmath90 by minimizing @xmath251 the resultant best fitted values of @xmath90 for the scalar , vector , and axial - vector bound states are @xmath252 respectively . we also plot ( the square of ) the ratio of the bound state masses to @xmath97 for several values of @xmath54 in fig . [ fig : m_over_fpi ] .    , @xmath253 , @xmath254 and @xmath255 for several values of @xmath54 ( indicated by @xmath83 , @xmath122 , @xmath256 , and @xmath257 , respectively ) . dotted lines represent the values of @xmath258 , @xmath259 , and @xmath260 obtained from eqs . ( [ eq : d_mv ] ) and ( [ eq : d_fpi ] ) . , height=302 ]    the dotted lines plotted together with the data in this figure represent the values of the following ratios obtained from eqs . ( [ eq : d_mv ] ) and ( [ eq : d_fpi ] ) : @xmath261 this figure clearly shows that all the masses of the scalar , vector , and axial - vector bound states have the same scaling property as that of @xmath97 : @xmath262 we can also say that these masses have the same scaling property as that of @xmath94 since @xmath97 and @xmath94 have the same scaling property . their ratios are summarized as follows : @xmath263 one might think that the vector and axial - vector boundstates decay into a fermion and an antifermion since @xmath264 and @xmath265 . however , this does not happen : as we noticed above eq . ( [ eq : sigma_m ] ) , @xmath82 is not the pole mass but the dynamical mass defined in the space - like region . furthermore , as we have shown in section  [ sec : complex_sd ] , the fermion propagators do not have any poles in the entire complex plane including the time - like axis where the pole mass of fermion should be defined . thus the vector and axial - vector boundstates do not decay into a fermion and an antifermion . in this paper we first pointed out that , when we solve the schwinger - dyson ( sd ) equation in large @xmath0 qcd , we do not need to introduce any infrared regularizations for the running coupling since it takes a finite value for all the range of the energy region due to the existence of the ir fixed point . in the case of small @xmath0 , we have to regularize the infrared divergence of the running coupling , and different schemes of regularizations would give different results . furthermore , it is difficult to find the regularization which makes the analytic structure of the running coupling simple enough . on the contrary , the solution of the two - loop rge in large @xmath0 qcd is explicitly written in terms of the lambert @xmath11 function , and the running coupling does not have any singularities on the complex plane except for the time - like axis when @xmath0 is close to @xmath6 . this significant feature of the running coupling in large @xmath0 qcd enabled us to solve the sd equation on the complex plane . then , we solved the homogeneous bethe - salpeter ( hbs ) equations for the scalar , pseudoscalar , vector , and axial - vector bound states of quark and anti - quark in large @xmath0 qcd with the improved ladder approximation in the landau gauge . in the quark propagator included in the hbs equation , we used the quark mass function obtained from the sd equation with the same approximation , which is needed for the consistency with the chiral symmetry . we first showed that the hbs equation provides the massless pseudoscalar bound state in the broken phase which is identified with the nambu - goldstone boson associated with the spontaneous breaking of the chiral symmetry . next , we showed that there actually exist vector , axial - vector , and scalar bound states even near the phase transition point in the broken phase , and that their masses decreases as the number of massless quarks @xmath0 increases . at the critical point all the masses go to zero , showing the same scaling property as that of the pseudoscalar meson decay constant @xmath97 consistently with the picture expected from the conformal phase transition  @xcite . let us discuss the pattern of the chiral symmetry restoration by considering the representation of chiral @xmath266 of the low - lying mesons extending the analyses done in refs . @xcite . for @xmath267 the pseudoscalar meson denoted by @xmath268 and the longitudinal axial - vector meson denoted by @xmath269 are an admixture of @xmath270 and @xmath271 , since the chiral symmetry is spontaneously broken  @xcite @xmath272 where @xmath235 denotes the helicity in the collinear frame , and the experimental value of the mixing angle @xmath273 is given by approximately @xmath274 @xcite . on the other hand , the longitudinal vector meson denoted by @xmath275 belongs to pure @xmath276 and the scalar meson denoted by @xmath277 to pure @xmath278 : @xmath279    when the chiral symmetry is restored at the phase transition point , it is natural that the chiral representations coincide with the mass eigenstates : the representation mixing is dissolved . from eq . ( [ mix pi a ] ) one can easily see that there are two ways to express the representations in the wigner phase of the chiral symmetry : the conventional manifestation  la the linear sigma model ( called the gl manifestation in ref . @xcite ) corresponds to the limit @xmath280 in which @xmath268 is in the representation of pure @xmath281 of @xmath282 together with the scalar meson , both being the chiral partners : @xmath283 on the other hand , the vector manifestation ( vm )  @xcite corresponds to the limit @xmath284 in which the @xmath269 goes to a pure @xmath285 , now degenerate with the scalar meson in the same representation , but not with @xmath275 in @xmath286 : @xmath287 namely , the degenerate massless @xmath268 and ( longitudinal ) @xmath275 at the phase transition point are the chiral partners in the representation of @xmath286 . now , what does our result say on the chiral representation of low - lying mesons ? as can be seen from fig . [ fig : m_over_fpi ] , the resultant values of the masses obtained from the hbs equation roughly satisfies the following relation  @xcite : @xmath288 for all values of @xmath54 . this relation holds independently of the mixing angle @xmath273 given in eq . ( [ mix pi a ] ) when the low - lying mesons saturate the chiral algebra shown in ref .  @xcite . then , it is reasonable to discuss the chiral representation without worrying about the influence of the excited states of the bound states . by using the relation @xmath289 @xcite and the values of @xmath249 and @xmath248 obtained from the hbs equation in subsection  [ ssec : vb ] , the value of the mixing angle @xmath273 is roughly determined as @xmath290 this implies that @xmath268 and the longitudinal @xmath269 are still admixtures of the pure chiral representation even at the chiral restoration point : @xmath291 this may suggest the existence of a new type of manifestation of chiral symmetry restoration in large @xmath0 qcd which is neither of the gl manifestation nor the simple version of the vector manifestation ( vm ) . several comments are in order    in appendix  [ app : decay_const ] , we show the calculations of the coupling constants @xmath12 , @xmath13 , and @xmath14 of the vector , axial - vector , and scalar bound states to the vector current , axial - vector current , and scalar density . the results shows that they also have the same scaling properties as @xmath97 . these results indicate that all the dimension - full quantities determined by the infrared dynamics have the same scaling properties , as far as the ( improved ) ladder approximation is concerned . although the masses obtained from the hbs equation satisfy the condition ( [ sat cond ] ) needed for the saturation of the chiral algebra , the couplings @xmath97 , @xmath12 , and @xmath13 do not seem to satisfy the first weinberg s sum rule  @xcite : @xmath292 . we have not fully understood what this means for the pattern of chiral symmetry restoration . apparently , reducing the numerical uncertainty will help us to reach the final understanding .    in the present analysis we did not include the effect from the four - fermion interaction which is induced in the case of @xmath110 as was conjectured in strong coupling qed  @xcite and was demonstrated in walking gauge theories  @xcite . it is not clear at this moment whether or not the qualitative results in the present analysis will be changed when we include such an effect .    in the present analysis we stressed that the running coupling in large @xmath0 qcd determined from a two - loop @xmath8-function is expressed as the lambert @xmath11 function which enables us to solve the hbs and sd equations with mutual consistency near the critical point . apparently , this lambert @xmath11 function can be used as an infrared regularization to solve the hbs and sd equations with mutual consistency in the case of qcd with small @xmath0 . it will be very interesting to study meson masses near the chiral phase transition in hot and/or dense qcd by using such an infrared regularization . we think that it is important to clarify which effective field theory ( eft ) describes the new pattern of the chiral symmetry restoration expressed in eq . ( [ new pattern ] ) . especially , it is very interesting to see how the matching between the eft and the underlying qcd with large @xmath0 can be done to determine bare parameters of the lagrangian of the eft . this work was supported in part by the jsps grant - in - aid for the scientific research ( b)(2 ) 14340072 . the work of m.h . was supported in part by the brain pool program ( # 012 - 1 - 44 ) provided by the korean federation of science and technology societies and usdoe grant # de - fg02 - 88er40388 . m.h . would like to thank professor mannque rho , professor gerry brown , and professor dong - pil min for their hospitality during his stays at kias , at suny at stony brook , and at seoul national university where part of this work was done .    * appendices * in this appendix , to show the validity of the fictitious eigenvalue method for solving the hbs equation explained in subsection  [ ssec : fem ] , we calculate the mass of the ortho - positronium which is the vector bound state of the electron and the positron . the same analysis was done in ref . @xcite , and here we follow the analysis .    in the weak coupling limit the hbs equation for the ortho - positronium can be solved analytically , and the energy spectrum takes the following form  @xcite : @xmath293 where @xmath294 is the mass of the electron and the positron , and @xmath295 is the coupling constant of qed . we use following parameters in our calculation : @xmath296             =    [ \ -18.5,\ -2.9\ ] \ , \ \ [ \ \lambda_{x } , \ \lambda_{x } \ ]           =    [ \ -10.8,\ 2.2\ ] \ , \ ] ] @xmath297 ( we used the energy scale which satisfies the relation @xmath298 following ref . @xcite . ) in fig . [ fig : positronium_bssv ] , we show the resultant values of the fictitious eigenvalue @xmath235 for several values of the input parameter @xmath249 .        finding the point where the smallest @xmath235 becomes unity , we determine the value of the mass of the ground state @xmath299 as the solution of the original hbs equation : @xmath300 from this the binding energy is calculated as @xmath301 these values are in good agreement with the values @xmath302 and @xmath303 derived from eq . ( [ eq : positronium_mass ] ) . this shows that our numerical method works well to obtain the mass of the ground state . in this appendix we show the explicit forms of the bispinor bases for the scalar , pseudoscalar , vector , and axial - vector bound states . here , we use the notation @xmath304 with @xmath190 being the mass of the bound states , and @xmath305 \equiv a[b , c ] + b[c , a ] + c[a , b]$ ] . bispinor base for the scalar bound state ( @xmath306 ) is given by @xmath307 , \end{aligned}\ ] ] and that for the pseudoscalar bound state ( @xmath308 ) is given by @xmath309\                      \gamma_5 \ .\end{aligned}\ ] ] furthermore , for the vector bound state ( @xmath310 ) we use @xmath311                              ( p \cdot \hat q ) , \ \ \gamma_{(v)}^3 = \frac{1}{2 } [ \fsl{\epsilon},\fsl{\hat q } ] , \ \        \gamma_{(v)}^4 = \frac{1}{3!}[\fsl{\epsilon},\fsl{p},\fsl{\hat       q } ] ,        \\    & &   \gamma_{(v)}^5 = ( \epsilon \cdot p ) , \ \        \gamma_{(v)}^6 = \fsl{p } ( \epsilon \cdot p ) , \ \        \gamma_{(v)}^7 = \fsl{\hat q}(p \cdot \hat q ) ( \epsilon \cdot p ) , \ \ \gamma_{(v)}^8 = \frac{1}{2 } [ \fsl{p},\fsl{\hat q}](\epsilon \cdot p ) , \nonumber \label{v bases}\end{aligned}\ ] ] and for the axial - vector bound state ( @xmath312 ) @xmath313                              \gamma_5   , \ \ \        \gamma_{(a)}^3 = \frac{1}{2 } [ \fsl{\epsilon},\fsl{\hat q}]\ ( p \cdot \hat q )                  \ \gamma_5   , \nonumber\\       \gamma_{(a)}^4 & = & \frac{1}{3!}[\fsl{\epsilon},\fsl{p},\fsl{\hat q } ]       \ \gamma_5 , \ \ \        \gamma_{(a)}^5 = ( \epsilon \cdot p)\ ( p \cdot \hat q)\ \gamma_5   , \ \ \        \gamma_{(a)}^6 = \fsl{p } ( \epsilon \cdot p)\ \gamma_5   ,       \nonumber\\       \gamma_{(a)}^7 & = & \fsl{\hat q}\ ( \epsilon \cdot p)\ ( p \cdot \hat q )                  \ \gamma_5   , \ \ \        \gamma_{(a)}^8 = \frac{1}{2 } [ \fsl{p},\fsl{\hat q}](\epsilon \cdot p)\       ( p \cdot \hat q)\ \gamma_5   .\end{aligned}\ ] ] in this subsection we calculate coupling constants @xmath12 , @xmath13 , and @xmath14 of the vector , axial - vector , and scalar bound states to the vector current , axial - vector current , and scalar density . they are defined by @xmath314 where @xmath174 is the flavor matrix normalized as @xmath315 = 2\delta_{a b}$ ] .    by using the bs amplitudes for the vector bound state , @xmath12 is expressed as @xmath316 \ . \label{eq : f_v}\ ] ] in the above expression , the normalization of the bs amplitudes @xmath182 are determined by the following normalization condition  @xcite : @xmath317 \ . \label{eq : normalization}\end{aligned}\ ] ] here , we notice again that @xmath318 has no singularity although the fermion propagator @xmath151 has a branch cut in the time - like region . so the integral in eq . ( [ eq : normalization ] ) is well - defined . once we have obtained @xmath249 and the corresponding bs amplitudes by solving the hbs equation , we can calculate @xmath12 from eq . ( [ eq : f_v ] ) . we can also calculate @xmath13 and @xmath14 in a similar way . in fig . [ fig : decay_const ] ( a ) we show the values of @xmath12 , @xmath13 , and @xmath14 together with @xmath97 for several values of @xmath54 .    , @xmath13 , @xmath14 , and @xmath97 and  ( b )  @xmath319 , @xmath320 , and @xmath321 for several values of @xmath54 . , title="fig:",height=188 ] , @xmath13 , @xmath14 , and @xmath97 and  ( b )  @xmath319 , @xmath320 , and @xmath321 for several values of @xmath54 . , title="fig:",height=188 ]   +  ( a )   ( b )    to see the scaling properties , we plot the ratio of @xmath12 , @xmath13 , and @xmath14 to @xmath97 in fig .  [ fig : decay_const ] ( b ) . this figure shows that @xmath12 , @xmath13 , and @xmath14 have the same scaling properties as that of @xmath97 . to solve the hbs equation for the bound states numerically , we introduced the uv and ir cutoffs and converted the hbs equation into a linear eigenvalue equation by discretizing the integral . as we have discussed in subsection [ ssec : vb ] , we adjust the values of the uv and ir cutoffs in such a way that the dominant supports of the integrands of the decay constant in eq . ( [ eq : f_v ] ) and the normalization condition in eq . ( [ eq : normalization ] ) lie in the energy region between the uv and ir cutoffs . in figs . [ fig : sft0885 ] and [ fig : sft1000 ] we show those integrands for @xmath322 and @xmath323 in the case of the vector bound state .    ) and ( b ) the normalization condition in eq . ( [ eq : normalization ] ) for @xmath324 . , title="fig:",width=226 ]  ) and ( b ) the normalization condition in eq . ( [ eq : normalization ] ) for @xmath324 . , title="fig:",width=226 ] +   +  ( a )   ( b )     ) and  ( b )   the normalization condition in eq . ( [ eq : normalization ] ) for @xmath325 . , title="fig:",width=226 ]  ) and  ( b )   the normalization condition in eq . ( [ eq : normalization ] ) for @xmath325 . , title="fig:",width=226 ] +   +  ( a )   ( b )     these figures show that the dominant supports lie in the lower energy region for smaller value of @xmath54 , and that the present choices in eqs . ( [ eq : ir cutoffs ] ) and ( [ eq : uv cutoffs ] ) covers the supports . for other values of @xmath54 used in the present analysis we have checked that the dominant supports always lie within the energy region between the uv and ir cutoffs chosen as in eqs . ( [ eq : ir cutoffs ] ) and ( [ eq : uv cutoffs ] ) .    as for the numbers of the discretization , due to the limitation of the computer resources we used @xmath326 and @xmath327 as shown in eq . ( [ nbs num ] ) . here we study the dependences of the mass and the decay constant of the vector bound state on the size of discretization . we show the typical values of the mass in fig . [ fig : nbs_dependence ] ( a ) and those of the decay constant in fig .  [ fig : nbs_dependence ] ( b ) for five choices of the size of the discretization , @xmath328 , @xmath329 , @xmath330 , @xmath331 , and @xmath332 .     and  ( b )  @xmath333 for five choices of the size of discretization , @xmath328 , @xmath329 , @xmath330 , @xmath331 , and @xmath332 . , title="fig:",height=188 ]   and  ( b )  @xmath333 for five choices of the size of discretization , @xmath328 , @xmath329 , @xmath330 , @xmath331 , and @xmath332 . , title="fig:",height=188 ]   +  ( a )   ( b )    figure  [ fig : nbs_dependence ] ( a ) clearly shows that the choice @xmath334 is large enough to obtain the mass of the vector bound state . on the other hand , fig . [ fig : nbs_dependence ] ( b ) shows that there are still uncertainties in the value of the decay constant which come from the size of the discretization . apparently , this uncertainty from the size of discretization is the dominant part of the numerical uncertainties in the present analysis . t.  banks and a.  zaks , nucl . b * 196 * , 189 ( 1982 ) . t.  appelquist , j.  terning and l.  c.  wijewardhana , phys . lett .   * 77 * , 1214 ( 1996 ) [ arxiv : hep - ph/9602385 ] . v.  a.  miransky and k.  yamawaki , phys . d * 55 * , 5051 ( 1997 ) [ erratum - ibid . d * 56 * , 3768 ( 1997 ) ] [ arxiv : hep - th/9611142 ] . t.  appelquist , a.  ratnaweera , j.  terning and l.  c.  wijewardhana , phys . d * 58 * , 105017 ( 1998 ) [ arxiv : hep - ph/9806472 ] . j.  b.  kogut and d.  k.  sinclair , nucl . b * 295 * , 465 ( 1988 ) ; f.  r.  brown , h.  chen , n.  h.  christ , z.  dong , r.  d.  mawhinney , w.  schaffer and a.  vaccarino , phys . d * 46 * , 5655 ( 1992 ) ; y.  iwasaki , k.  kanaya , s.  sakai and t.  yoshie , phys . lett .   * 69 * , 21 ( 1992 ) ; y.  iwasaki , k.  kanaya , s.  kaya , s.  sakai and t.  yoshie , nucl . suppl .   * 53 * , 449 ( 1997 ) ; prog . theor .  phys . suppl .   * 131 * , 415 ( 1998 ) . r.  oehme and w.  zimmermann , phys . d * 21 * , 471 ( 1980 ) . m.  velkovsky and e.  shuryak , phys . b * 437 * , 398 ( 1998 ) . t.  hatsuda and t.  kunihiro , phys . rept .   * 247 * , 221 ( 1994 ) ; r.  d.  pisarski , hep - ph/9503330 ; g.e . brown and m.  rho , phys . * 269 * , 333 ( 1996 ) ; f.  wilczek , hep - ph/0003183 ; g.  e.  brown and m.  rho , phys . * 363 * , 85 ( 2002 ) [ arxiv : hep - ph/0103102 ] . m.  harada and k.  yamawaki , phys . lett .   * 83 * , 3374 ( 1999 ) [ arxiv : hep - ph/9906445 ] . m.  harada and k.  yamawaki , phys . lett .   * 86 * , 757 ( 2001 ) [ arxiv : hep - ph/0010207 ] . m.  harada and k.  yamawaki , to appear in phys . [ arxiv : hep - ph/0302103 ] . m.  bando , t.  kugo , s.  uehara , k.  yamawaki and t.  yanagida , phys . lett .   * 54 * , 1215 ( 1985 ) ;   m.  bando , t.  kugo and k.  yamawaki , phys . rept .   * 164 * , 217 ( 1988 ) . r.  s.  chivukula , phys . d * 55 * , 5238 ( 1997 ) [ arxiv : hep - ph/9612267 ] . n.  nakanishi , prog . * 43 * , 1 ( 1969 ) . t.  kugo , in _ proc . of 1991 nagoya spring school on dynamical symmetry breaking , nakatsugawa , japan , 1991 _ , ed . k.  yamawaki ( world scientific , singapore , 1992 ) . t.  maskawa and h.  nakajima , prog .  theor . phys .   * 52 * , 1326 ( 1974 ) ; t.  maskawa and h.  nakajima , prog .  theor . phys .   * 54 * , 860 ( 1975 ) ; m.  bando , m.  harada and t.  kugo , prog .  theor .  phys . * 91 * , 927 ( 1994 ) [ arxiv : hep - ph/9312343 ] . t.  kugo and m.  g.  mitchard , phys . b * 282 * , 162 ( 1992 ) . t.  kugo and m.  g.  mitchard , phys . b * 286 * , 355 ( 1992 ) . m.  bando , m.  harada and t.  kugo , prog . phys .   * 91 * , 927 ( 1994 ) [ arxiv : hep - ph/9312343 ] . k.  higashijima , phys . d * 29 * , 1228 ( 1984 ) . v.  a.  miransky , sov . j.  nucl . * 38 * , 280 ( 1983 ) [ yad . fiz .   * 38 * , 468 ( 1983 ) ] . k.  i.  aoki , m.  bando , t.  kugo , m.  g.  mitchard and h.  nakatani , prog . phys .   * 84 * , 683 ( 1990 ) . k.  i.  aoki , m.  bando , t.  kugo and m.  g.  mitchard , prog . * 85 * , 355 ( 1991 ) .    c.  d.  roberts and s.  m.  schmidt , prog . * 45 * , s1 ( 2000 ) [ arxiv : nucl - th/0005064 ] . t.  kugo , m.  g.  mitchard and y.  yoshida , prog . * 91 * , 521 ( 1994 ) [ arxiv : hep - ph/9312267 ] . p.  maris , c.  d.  roberts and p.  c.  tandy , phys . b * 420 * , 267 ( 1998 ) [ arxiv : nucl - th/9707003 ] . k.  i.  aoki , t.  kugo and m.  g.  mitchard , phys . b * 266 * , 467 ( 1991 ) . k.  naito , k.  yoshida , y.  nemoto , m.  oka and m.  takizawa , phys . c * 59 * , 1095 ( 1999 ) [ arxiv : hep - ph/9805243 ] ; k.  naito and m.  oka , phys . c * 59 * , 542 ( 1999 ) [ arxiv : hep - ph/9805258 ] . r.  alkofer , p.  watson and h.  weigel , phys . d * 65 * , 094026 ( 2002 ) [ arxiv : hep - ph/0202053 ] . e.  gardi and m.  karliner , nucl . b * 529 * , 383 ( 1998 ) [ arxiv : hep - ph/9802218 ] ; e.  gardi , g.  grunberg and m.  karliner , jhep * 9807 * , 007 ( 1998 ) [ arxiv : hep - ph/9806462 ] . c.  n.  leung , s.  t.  love and w.  a.  bardeen , nucl . b * 273 * , 649 ( 1986 ) ; nucl . b * 323 * , 493 ( 1989 ) . t.  appelquist , k.  d.  lane and u.  mahanta , phys . lett .   * 61 * , 1553 ( 1988 ) ; h.  georgi , e.  h.  simmons and a.  g.  cohen , phys . b * 236 * , 183 ( 1990 ) . r.  m.  corless , g.  h.  gonnet , d.  g.  e.  hare , d.  j.  jeffrey and d.  e.  knuth , adv . math .   * 5 * , 329 ( 1996 ) . h.  pagels and s.  stokar , phys . d * 20 * , 2947 ( 1979 ) . b. holdom , phys . b * 150 * , 301 ( 1985 ) ; k. yamawaki , m. bando and k. matumoto , phys . lett .   * 56 * , 1335 ( 1986 ) ; t. akiba and t. yanagida , phys . b * 169 * , 432 ( 1986 ) ; t.w . appelquist , d. karabali and l.c.r . wijewardhana , phys . lett .   * 57 * , 957 ( 1986 ) . t.  kugo and y.  yoshida , soryushiron kenkyu * 91 * , b26 ( 1995 ) ; y.  yoshida , ph.d thesis , kyoto university ( 1995 )    m.  harada and y.  yoshida , phys .  rev . d * 53 * , 1482 ( 1996 ) [ arxiv : hep - ph/9505206 ] . f.  j.  gilman and h.  harari , phys . rev .   * 165 * , 1803 ( 1968 ) . s.  weinberg , phys . rev .   * 177 * , 2604 ( 1969 ) . s.  weinberg , phys . lett .   * 18 * , 507 ( 1967 ) . kondo , s.  shuto and k.  yamawaki , mod . a * 6 * , 3385 ( 1991 ) . k.  i.  aoki , k.  morikawa , j.  i.  sumi , h.  terao and m.  tomoyose , prog . * 102 * , 1151 ( 1999 ) [ arxiv : hep - th/9908042 ] . see , for example , c.  itzykson and j.  b.  zuber , _ quantum field theory _ ( mcgraw - hill , new york , 1980 ) .
Abstract:

Answer: we solve the homogeneous bethe - salpeter ( hbs ) equation for the scalar , pseudoscalar , vector , and axial - vector bound states of quark and anti - quark in large @xmath0 qcd with the improved ladder approximation in the landau gauge . the quark mass function in the hbs equation is obtained from the schwinger - dyson ( sd ) equation in the same approximation for consistency with the chiral symmetry . amazingly , due to the fact that the two - loop running coupling of large @xmath0 qcd is explicitly written in terms of an analytic function , large @xmath0 qcd turns out to be the first example in which the sd equation can be solved in the complex plane and hence the hbs equation directly in the time - like region . we find that approaching the chiral phase transition point from the broken phase , the scalar , vector , and axial - vector meson masses vanish to zero with the same scaling behavior , all degenerate with the massless pseudoscalar meson . this may suggest a new type of manifestation of the chiral symmetry restoration in large @xmath0 qcd .