the standard model is greatly successful but it still has many free parameters which must be small to describe nature . while its supersymmetric extensions , e.g. ,  the minimal supersymmetric standard model , are attractive scenarios , small couplings are also required to explain observed facts such as the fermion mass hierarchy and mixing angles .    in recent years , extra dimensions have cast a new perspective on physics beyond the standard model .
one of the important aspects of extra dimensional models is that bulk fields can be localized with finite - width wave - function profiles .
this fact provides us with a geometrical explanation for small numbers .
that is , with a configuration where some fields are separated from each other in the extra dimensional space , the couplings among them are generally suppressed .
then how and where fields are localized is an issue to be considered . from this viewpoint , extra dimensional models with a curved background
are interesting because fields could be localized depending on the shape of the background geometry .
one of the most famous examples of curved geometries is the randall - sundrum ( rs ) model with the ads@xmath0 warped metric  @xcite .
field theories of vectors , spinors , and scalars have been studied on this background  @xcite-@xcite .
the localization behavior of zero - mode wave functions has interesting applications to phenomenology such as the suppression of unwanted operators .
for example , hierarchical forms of yukawa couplings and proton decay were studied in  @xcite .
the localization of kaluza - klein ( kk ) excited modes also leads to interesting phenomena .
for instance , the localization of higher kk gauge bosons could realize a composite scalar ( higgs ) condensation , which induces dynamical ( electroweak ) symmetry breaking on the brane where the kk gauge bosons localize  @xcite .
in addition , models on more complicated backgrounds where a warp factor oscillates generate bulk fields which localize at some points in extra dimensions  @xcite . this type of localization might be useful in explaining the observed phenomena .
however extra dimensional theories are generally nonrenormalizable and the calculations depend on the regularization scheme that one adopts . furthermore ,
extra dimensional theories are constrained by symmetries of higher dimensions .
for example , in the supersymmetric case , bulk theories are constrained by @xmath5 supersymmetry in five dimensions .
motivated by these facts , recently a four - dimensional ( 4d ) description of extra dimensional models was proposed  @xcite . with this method
, the phenomena of higher dimensional models are reproduced in terms of 4d theories , and several interesting models have been proposed along this line  @xcite .    in this paper
, we present 4d gauge theories that describe physics on 5d curved geometries .
as will be discussed below , taking generic values of gauge couplings and gauge - symmetry - breaking vacuum expectation values ( vevs ) , the models provide vector , spinor , and scalar fields on curved extra dimensions . as a good and simple illustration
, we compare our 4d model with the rs one .
we particularly focus on the `` localization '' behaviors of mass eigenstates in `` index spaces '' of gauge groups .
it will be shown that the localization profiles and the exponentially suppressed massive spectrum are certainly reproduced .
in addition , our formulation gives a localization mechanism even for massless vector fields . as a phenomenological application ,
hierarchical yukawa matrices are derived in our approach ; that is a hierarchy without symmetries in four dimensions .
the localization behavior depends on the required conditions for gauge - symmetry - breaking vevs and gauge and other couplings .
if these values are determined in the underlying theories , it may be said that the physics on warped backgrounds is dynamically generated within a four - dimensional framework .
we consider several possibilities to realize the conditions by utilizing , for example , strongly coupled gauge theories .
thus this could provide a purely 4d dynamical approach for small numbers .
we will proceed with the argument as follows . in sec .
[ sec : dwd ] , we describe our 4d gauge theories , which have generic ( nonuniversal ) values of gauge - symmetry - breaking vevs and couplings .
the models provide vector , spinor , and scalar fields in warped extra dimensions .
it is also shown that supersymmetry multiplets in flat 4d models generate supersymmetry multiplets on warped backgrounds . in sec .
[ sec : ne ] , we then numerically determine with a finite number of gauge groups that the formulation given in sec .
2 certainly reproduces various properties of bulk fields on the rs background .
in addition , a phenomenological application to quark mass matrices is also given . finally , we discuss possibilities of dynamically realizing the conditions required for curved geometries in sec .
[ sec : ddwd ] .
we conclude the discussion in sec .
[ sec : conclusion ] .
the appendix is devoted to a brief review of 5d bulk fields on a rs background .
following refs .
@xcite , we introduce @xmath6 gauge theories with gauge couplings @xmath7 ( @xmath8 ) , and scalar fields @xmath9 [ @xmath10 which are in bifundamental representations of @xmath11 .
the system is schematized by the segment diagram in fig .
[ fig : links ] .
the gauge invariant kinetic term of the scalars @xmath9 is written by @xmath12 where the covariant derivative is given by @xmath13 .
we assume that the scalar fields @xmath9 develop vevs proportional to the unit matrix , @xmath14 , which break the gauge symmetry to a diagonal @xmath15 . from the kinetic term ( [ eq : qkin ] ) , the mass terms for the vector fields @xmath16 are obtained : @xmath17 where the @xmath18 matrix @xmath19 is defined as @xmath20    the consequence of these mass terms is that we have a massless gauge boson corresponding to the unbroken gauge symmetry , which is given by the following linear combination : @xmath21 where @xmath22 and @xmath23 is the gauge coupling of the low - energy gauge theory @xmath15 .
the profile of @xmath24 is independent of the values of @xmath25 .
it is found from this that the massless vector field is `` localized '' at the points with smaller gauge couplings .
if the gauge couplings take a universal value @xmath26 , the massless mode @xmath27 has a constant `` wave function '' along the `` index space '' of gauge groups . as seen below , this direction labeled by @xmath28
becomes the fifth spatial dimension in the continuum limit ( @xmath29 ) .
the localization behavior can easily be understood from the fact that , for smaller gauge coupling @xmath7 , the symmetry - breaking scale @xmath30 of @xmath6 becomes lower , and hence the corresponding vector field @xmath16 becomes the more dominant component in the low - energy degree of freedom @xmath27 .
it is interesting to note that this vector localization mechanism ensures charge universality .
suppose that there is a field in a nontrivial representation of @xmath6 only .
that is , it couples only to @xmath31 with strength @xmath7 .
this corresponds to a four - dimensional field confined on a brane .
if there are several such fields , they generally have different values of gauge couplings .
however , note that these fields couple to the massless modes @xmath24 with a _
universal _ gauge coupling @xmath23 defined above .
this is because , in the presented mechanism , the vector fields are localized depending on the values of the gauge couplings .    as for massless eigenstates ,
the mass eigenvalues and wave functions are obtained by diagonalizing the mass matrix ( [ eq : diffop ] ) .
the simplest case is the universal couplings @xmath32 in this case , one obtains the mass eigenvalues of @xmath19 as @xmath33 in the limit @xmath29 with @xmath34 fixed ( the limit to continuum 5d theory ) , the eigenvalues become @xmath35 these are the same spectrum as that of the bulk gauge boson in the @xmath36 extra dimension with radius @xmath37 .    with generic values of vevs @xmath25 and gauge couplings
@xmath7 , the situation is rather complicated . in this case , the mass term ( [ eq : gau - mass ] ) of the vector fields becomes @xmath38(a_\mu^{i+1})^2    \nonumber \\[1 mm ]    & & + \frac{1}{2}v_n^2g_n(g_{n+1}-g_n)(a_\mu^n)^2    -\frac{1}{2 } v_1 ^ 2g_1(g_2-g_1)(a_\mu^1)^2 .
\label{eq : gbbm}\end{aligned}\ ] ] the first term becomes the kinetic energy transverse to the four dimensions in the continuum limit . on the other hand ,
the second and third terms are bulk and brane mass terms , respectively .
it should be noted that these mass terms vanish in the case of universal gauge coupling , which corresponds to a flat massless vector field in 5d theory as discussed above . in other words ,
nonuniversal gauge couplings generate bulk / brane mass terms and cause a localization of the wave function .
first we consider the series of vevs @xmath25 and couplings @xmath7 that generates a vector field on the rs warped background , namely , the ads@xmath0 background .
this model can be obtained by choosing a universal @xmath7 and by varying @xmath25 as @xmath39 substituting this and taking the continuum limit , eq .
( [ gkt ] ) becomes @xmath40 ^ 2 , \label{eq : contlimrs}\end{aligned}\ ] ] where @xmath41 represents the coordinate of the extra dimension : @xmath42 ( @xmath8 ) and @xmath43 , etc .
it is found that eq .
( [ eq : contlimrs ] ) successfully induces the kinetic energy term along the extra dimension and mass terms for the vector field on the warped background metric @xmath44 where @xmath45 with @xmath46 .
we here conclude that we can obtain the vector field on a rs warped background by varying only the vevs @xmath25 . in the following
we will see that nonuniversal gauge couplings @xmath7 induce other interesting results beyond the effects from the background metric .
now let us compare the 4d model with generic couplings ( [ eq : gbbm ] ) to extra dimensional ones .
we define the dimensionless parameters @xmath47 and @xmath48 as @xmath49 first we restrict ourselves to the case that the gauge group is @xmath50 , namely , abelian theory with _ no vector self - couplings_. similarly substituting eq .
( [ eq : deffh ] ) and taking the continuum limit , eq .
( [ gkt ] ) becomes @xmath51\bigr]^2 .
\label{eq : contlim}\end{aligned}\ ] ] equation ( [ eq : contlim ] ) induces the kinetic energy term along the extra dimension and mass terms for the vector field on the warped background metric : @xmath52 ^ 2 \eta_{\mu\nu}dx^\mu dx^\nu - dy^2 .
\label{eq : generalbkg}\ ] ] the bulk and boundary mass terms are @xmath41 dependent and proportional to the derivatives of @xmath53 .
this is also seen from the 4d model [ the second and third terms in eq .
( [ eq : gbbm ] ) ] .
the above is a generic correspondence between our 4d case and continuum 5d theory . as an example , consider the following forms of the parameters : @xmath54 where @xmath55 is a positive constant with mass dimension @xmath56 .
equation ( [ eq : contlim ] ) leads to @xmath57    -\frac{1}{2}\biggl[\zeta ke^{-2(\zeta+\eta)ky }    ( a_\mu)^2\biggr]_0^{l/2}.   \label{eq : wgmbbm}\end{aligned}\ ] ] the first term on the right hand side is the kinetic term of the gauge boson along the extra dimension with the warped background @xmath58 the second and third terms correspond to the bulk and boundary masses announced before .
as easily seen , the above equation includes the expression for vector fields on the rs background . in the 5d rs model
, the lagrangian for vector fields is written as ( see the appendix ) @xmath59 where the @xmath60 gauge fixing condition is chosen . by comparing eq .
( [ eq : wgmbbm ] ) with eq .
( [ rsgauge ] ) , we find that the case with @xmath61 realizes vector fields on the rs background .
also , a special limit , @xmath62 , produces the flat zero - mode solution .
that corresponds to the form of the parameters ( [ eq : rsparaf ] ) in the previous special argument .
the other solutions which satisfy @xmath63 correspond to nonflat wave functions of the zero - mode vector field on the rs background .
it is clearly understood in our formulation that such nonflat wave functions are caused by introducing bulk and/or boundary mass terms in the rs model .
for example , in the case of @xmath64 , the vector field has bulk and boundary mass terms , and is localized with a peak at the @xmath65 point .
it should be noticed that with these bare mass terms the zero mode is still massless .
this is understood from our formulation where the gauge symmetry @xmath66 is left unbroken in the low - energy effective theory .      in the above abelian case we discussed interpretation of the nonuniversal @xmath47 as @xmath41-dependent bulk or boundary masses in the warped extra dimension .
next we treat the non - abelian theory with vector self - couplings .
since in this case we also have @xmath41-dependent vector self - couplings in addition to the @xmath41-dependent bulk or boundary masses , it may be convenient and instructive to see @xmath47 as a @xmath41-dependent coefficient of the vector kinetic term .
to this end , we define the four - dimensional field @xmath67 , @xmath68 where @xmath69 is the lattice spacing , which goes to zero in the continuum limit . rescaling the gauge fields @xmath70 , the kinetic term @xmath71 and eq .
( [ eq : qkin ] ) become @xmath72   \nonumber \\ & & \qquad \qquad \qquad \qquad   \qquad \qquad               + i \hat{g } a_5^{i+1/2 } ( a^{i+1}_\mu - a^i_\mu )              + { \cal o}(a^{1/2 } ) \bigg|^2 ,   \label{eq : gqkin2}\end{aligned}\ ] ] where @xmath73 and @xmath74 . in the continuum limit @xmath75 with @xmath76 and @xmath77 fixed , eq .
( [ eq : gqkin2 ] ) results in @xmath78 ^ 2 \eta^{\mu \nu } f_{\mu 5 } f_{\nu 5 }         \big\},\end{aligned}\ ] ] where @xmath79 $ ] .
this completely reproduces the 5d yang - mills kinetic term with a @xmath41-dependent coefficient @xmath80 on the warped - background metric ( [ eq : generalbkg ] ) , provided that @xmath81 .
this is the generic correspondence between the present 4d model and continuum 5d theory . from eq .
( [ eq:5dym ] ) , we thus find the @xmath41-dependent factor @xmath82 in front of the canonical yang - mills term , which corresponds to a 5d dilaton vev .
the factor does carry the origin of the massless vector localization shown in eq .
( [ zero - a ] ) . with the constant gauge coupling @xmath83 ( @xmath47 = 1 )
, one obtains a bulk vector field with a constant zero mode on the warped metric ( [ eq : generalbkg ] ) .
a field redefinition @xmath84 in eq .
( [ eq:5dym ] ) gives the previous bulk and boundary mass terms but one then has @xmath41-dependent vector self - couplings in non - abelian cases .
we next consider spinor fields by arranging fermions of fundamental or antifundamental representation in each gauge theory @xmath6 .
we introduce two weyl ( one dirac ) spinors to construct a 5d bulk fermion .
the orbifold compactification in continuum theory requires that one spinor obeys the neumann boundary condition and the other the dirichlet one . in the present 4d model ,
this can be achieved by having asymmetrical numbers of fundamental and antifundamental spinors , resulting in chiral fermions in the low - energy gauge theory . here
we consider the fundamental weyl spinors @xmath85 ( @xmath8 ) in the @xmath6 theory and the antifundamental @xmath86 ( @xmath87 ) .
as seen below , @xmath88 corresponds to the bulk fermion with the neumann boundary condition and @xmath89 to that with the dirichlet one .
the generic gauge - invariant mass and the mixing terms of @xmath85 and @xmath86 are written as @xmath90 where @xmath91 and @xmath92 are dimensionless coupling constants .
we assume that @xmath9 develop vevs @xmath93 . the mass matrix for spinors
is then given by @xmath94 the spinor mass eigenvalues and eigenvectors ( wave functions ) are read from this matrix .
one easily sees that the massless mode is contained in @xmath88 and given by the following linear combination : @xmath95 therefore the localization profile of zero mode depends on the ratio of dimensionless couplings @xmath91 and @xmath92 .
a simple case is @xmath96 for all @xmath28 . in this case
, @xmath97 corresponds to a chiral zero mode obtained from a 5d bulk fermion on the flat background . if @xmath98 , the system describes a fermion with a curved wave - function profile .
for example , if @xmath99 ( @xmath100 ) , @xmath97 has a monotonically increasing ( decreasing ) wave - function profile . as another interesting example
, taking @xmath101 ( @xmath102 , @xmath103 are constants and @xmath104 ) , @xmath97 has a gaussian profile with a peak at @xmath105 .
other profiles of massless chiral fermions could also be realized in our approach .
let us discuss the 5d continuum limit .
the relevant choice of couplings @xmath91 and @xmath92 is @xmath106 the parameters @xmath107 give rise to a bulk bare mass in the continuum limit as will be seen below .
the only difference between the vector and spinor cases is the existence of possible bulk mass parameters [ see eqs .
( [ eq : diffop ] ) and ( [ eq : mass - spi ] ) ] .
the mass and mixing terms  ( [ eq : fermionmoose ] ) then become @xmath108 + \textrm{h.c.},\end{aligned}\ ] ] where @xmath109 and @xmath110 are the same as defined in the case of vector fields  ( [ eq : deffh ] ) .
similar to the vector case , this form is compared with the bulk spinor lagrangian in the rs space - time ( see the appendix ) @xmath111 .
\label{eq : rsspinor}\end{aligned}\ ] ] here the kinetic terms have been canonically normalized in order to compare them to the 4d model . in eq .
( [ eq : rsspinor ] ) , @xmath102 is a possible 5d dirac mass , and the `` 1/2 '' contribution in the mass terms comes from the spin connection with the rs metric .
it is interesting that the 5d spinor lagrangian ( [ eq : rsspinor ] ) is reproduced by taking the exact same limit of parameters as that in the vector case , defined by eq .
( [ eq : rsparaf ] ) .
furthermore , the relation between the mass parameters @xmath102 should be taken as @xmath112 that is , the @xmath107 s take a universal value .
now the localization behavior of the spinors is easily understood . in the present 4d model ,
the spinor mass matrix ( [ eq : mass - spi ] ) becomes with eq .
( [ eq : spinorcond ] ) @xmath113 a vanishing bulk mass parameter @xmath114 corresponds to @xmath115 , that is , @xmath116 in our model .
then the mass matrix @xmath117 is exactly the same as @xmath19 , and their mass eigenvalues and eigenstates are the same . in particular
, the massless mode @xmath97 has a flat wave function with universal gauge couplings as considered here .
this is consistent with the expression ( [ zero - eta ] ) , where the ratio @xmath118 determines the wave - function profile . on the other hand , in the case of @xmath119 ( @xmath120 )
, the rs zero - mode spinor is localized at @xmath121 ( @xmath65 )  @xcite , which in turn corresponds to @xmath122 ( @xmath123 ) in our model .
one can see from the spinor mass matrix ( [ eq : mass - rsspi ] ) that the zero - mode wave function is monotonically increasing ( decreasing ) with respect to the index @xmath28 .    in this way
, we have a 4d localization mechanism for the spinor fields .
nonuniversal gauge couplings or nonuniversal masses give rise to a nonflat wave function for a chiral massless fermion .
the latter option is not realized for vector fields .
notice that the charge universality still holds in the low - energy effective theory .
that is , with any complicated wave - function profiles , zero modes interact with a universal value of the gauge coupling .
this is because curved profiles of vector fields depend only on the gauge couplings .
finally we consider scalar fields .
we may introduce two types of scalar field @xmath124 and @xmath125 in the fundamental and antifundamental representations of @xmath6 , respectively .
in addition , for each type of scalar , there are two choices of the @xmath126 parity assignment in the continuum limit .
this orbifolding procedure is incorporated by removing @xmath127 or @xmath128 .
the gauge invariant mass and mixing terms for @xmath129 and @xmath130 are written as @xmath131 where @xmath132 , @xmath133 , and @xmath134 are the dimensionless coupling constants .
it is implicitly assumed that nonintroduced fields are appropriately removed in the sums .
we have included the mixing mass terms up to the nearest - neighbor interactions .
other invariant terms such as @xmath135 or terms containing @xmath9 correspond to nonlocal interactions in 5d theories , and we do not consider these in this paper .
notice , however , that for a supersymmetric case , these terms may be suppressed due to renormalizability and holomorphy of the superpotential .
the zero - mode eigenstates are given in the same form as that of the spinor shown in the previous section , replacing @xmath132 and @xmath133 by @xmath136 and @xmath137 ( @xmath138 and @xmath139 ) .
therefore the ratio @xmath140 ( @xmath141 ) determines the zero - mode wave function .
let us consider the continuum 5d limit .
in what follows , we remove @xmath128 , which corresponds to the @xmath126 assignment @xmath142 and @xmath143 .
the 5d limit can be achieved by taking the following choices of couplings : @xmath144 & & \bar\alpha'_i \,=\ , g_{i+1 } , \qquad \bar\beta'_i \,=\ , g_{i+1}(1-\bar c'_{i+1}),\end{aligned}\ ] ] where @xmath145 and @xmath146 correspond to the bulk mass parameters , as in the spinor case
. then the mass terms ( [ eq : nsp ] ) and ( [ eq : dsp ] ) for the scalars take the following forms with the parametrization ( [ eq : deffh ] ) : @xmath147\{f(y)\phi(x , y)\}\bigr|^2 + [ gv\gamma(y)]^2|f(y)\phi(x , y)|^2 \biggr],\ ] ] @xmath148\{h(y)\varphi(x , y)\ } \bigr|^2 + [ gv\bar\gamma(y)]^2|h(y)\varphi(x , y)|^2 \biggr].\ ] ]    as a special case , we compare @xmath129 and @xmath130 with the scalar fields in the rs space - time .
the scalar lagrangian on the rs background is ( see also the appendix ) @xmath149\end{aligned}\ ] ] where the 4d kinetic term is canonically normalized , and @xmath150 and @xmath151 are the bulk and boundary mass parameters , respectively , defined in the appendix [ eq .
( [ eq : smass ] ) ] . by substituting the rs limit in our model given by eq .
( [ eq : rsparaf ] ) , we find the relations between the mass parameters in 4d and 5d : @xmath152      in this subsection , we discuss 5d supersymmetry on warped backgrounds .
generally a supersymmetric theory may be obtained by relevant choices of couplings from a nonsupersymmetric theory .
we here examine whether it is possible to construct supersymmetric 4d models which describe 5d supersymmetric ones on warped backgrounds .
this is a nontrivial check for the ability of our formalism to properly describe 5d nature . in ref .
@xcite , the 5d theory on the ads@xmath0 rs background was studied . there
, supersymmetry on ads@xmath0 geometry was identified and then the conditions on the mass parameters imposed by this type of supersymmetry were derived ( also given in the appendix here ) . as seen below
, these relations among mass parameters for ads@xmath0 supersymmetry are indeed satisfied in our _
4d formalism_. this fact seems remarkable in the sense that the present analyses do not include gravity .
first consider vector supermultiplets in 5d .
the scalar fields @xmath9 and the gauge bosons @xmath31 are extended to chiral and vector superfields in 4d , respectively . notice that the vevs that were discussed above , @xmath153 are in the ( baryonic ) @xmath154-flat direction .
we start with the following 4d supersymmetric lagrangian @xmath155 the bilinear terms of the component fields are written in the unitary gauge ( we follow the conventions of @xcite ) : @xmath156 where we have rescaled @xmath157 for canonical normalization of the kinetic terms . the mass matrix @xmath158 is defined in eq .
( [ gkt ] ) .
the first term is nothing but eq .
( [ gkt ] ) , that is , the mass terms for vector fields . by
also canonically normalizing @xmath159 and @xmath160 and integrating out the auxiliary fields , we find the mass terms for the adjoint spinors and scalar @xmath161 these masses have the same forms as that of the vector field because we started from a supersymmetric theory .
we thus have a model with @xmath162 for the spinors and @xmath163 for the scalar .
( note that @xmath164 , which originate from @xmath9 , have @xmath126 odd parity . )
it is a nontrivial check to see whether the above mass terms satisfy the conditions for 5d ads@xmath0 supersymmetry .
we find from the relations ( [ eq : ci - spi ] ) and ( [ eq : massd ] ) that the mass terms for @xmath165 and @xmath166 imply @xmath167 indeed , these relations are just those required by ads@xmath0 supersymmetry  @xcite . in this way , _
5d vector supermultiplets on a rs warped background are automatically derived from a 4d supersymmetric model on a flat background_.    we also construct a 5d hypermultiplet in the warped extra dimension starting from a 4d supersymmetric theory . in order to have a hypermultiplet we introduce the chiral superfields @xmath124 and @xmath125 in the fundamental and antifundamental representations of the @xmath6 gauge theory . in the following , @xmath128 is removed to implement @xmath126 orbifolding which leaves a chiral zero mode of the fundamental representation .
the fermionic components of @xmath124 and @xmath125 , then correspond to @xmath85 and @xmath168 , respectively , in sec .
[ sec : spinor ] .
the generic renormalizable superpotential is written as @xmath169 this superpotential just leads to a spinor mass term of the form ( [ eq : fermionmoose ] ) .
in addition , the mass and mixing terms of the scalars @xmath129 and @xmath130 also have the same form as those of the spinors : @xmath170 supersymmetry induces equivalence between the boson and fermion mass matrices . in turn , this implies in our formulation given in the previous sections that the mass parameters are equal , @xmath171 and also @xmath172 .
thus , there is only one parameter @xmath102 left .
it is found from eqs .
( [ eq : ci - spi ] ) , ( [ eq : massn ] ) , and ( [ eq : massd ] ) that if one take the continuum limit the relations @xmath173 are generated .
these mass relations are exactly those imposed by supersymmetry on the ads@xmath0 geometry  @xcite .
thus hypermultiplets on the rs background are properly incorporated in our formalism with a flat background .
it may be interesting that the mass relation for vector multiplets is the one for chiral multiplets with dirichlet boundary conditions ( [ eq : abchypd ] ) with @xmath114 .
this value of @xmath102 is the limit of vanishing bulk mass parameters .
it should be noticed that our analyses have been performed for generic warped backgrounds , including the rs case as a special limit .
we thus found that even in generic warped backgrounds the conditions on the bulk mass parameters required for 5d warped supersymmetry should be the same as for the rs case .
here we perform a numerical study to confirm our formulation of the curved extra dimension discussed in the previous sections . we will also give a phenomenological application to the hierarchy among yukawa couplings .
in the following , we consider the case that corresponds to the rs model in the continuum limit , as a good and simple application . the gauge couplings and vevs are specified as given in sec .
[ sec : dwd ] ; @xmath174 the universal gauge coupling @xmath26 implies that vector zero modes have flat wave functions as shown in eq .
( [ zero - a ] ) .
the following is a summary of the mass terms for various spin fields , which were derived in the previous sections : @xmath175 { \cal l}_{fm } & = & -\psi d_c \eta + \textrm{h.c . } , \\[2 mm ] { \cal l}_{sm}^\phi & = & -\bigl|d_{3/2-b } \phi\bigr|^2 -|m \phi|^2 , \\[1 mm ] { \cal l}_{sm}^\varphi & = & -\bigl|d_{-(3/2-b)}^\dagger \varphi\bigr|^2 -|m^\dagger \varphi|^2.\end{aligned}\ ] ] the parameters @xmath151 and @xmath102 represent the bulk mass parameters for scalars and spinors , respectively .
the @xmath18 mass matrices @xmath176 and @xmath177 are defined as follows : @xmath178    m & = & \sqrt{a+4b - b^2}\,\frac{k}{v }    \pmatrix { v_1 & & \cr & \ddots & \cr & & v_{n-1 } }    \pmatrix { 1 &     0    &          &    \cr      & \ddots & \ddots   &    \cr      &         & 1 & 0 } , \end{aligned}\ ] ] where @xmath179 for supersymmetric cases , the mass matrices @xmath176 for bosons and fermions take the same form and , moreover , @xmath180 , as discussed previously .
we define the matrices @xmath181 that diagonalize the mass matrices for gauge , fermion , and scalar fields , respectively . for example
, @xmath182 satisfies @xmath183 where @xmath184 are the mass eigenvalues which should correspond to the kk spectrum of vector fields . in the following , we use the notation @xmath185 that is , the coefficients of @xmath31 in the @xmath186th massive eigenstates @xmath187 . in the continuum limit
, this corresponds to the value of the wave function at @xmath188 for the @xmath186th kk excited vector field .
similar definitions are made for spinors and scalars .    for vector fields ,
we illustrate the resultant eigenvalues @xmath189 and eigenvectors @xmath190 in figs .
[ fig : bulkgaugesys](a)[fig : bulkgaugesys](f ) . for comparison
, we also show in the figures the wave functions and kk mass eigenvalues of vector fields on the rs background .
it is found from the figures that our 4d model completely reproduces the mode function profiles [ figs .
[ fig : bulkgaugesys](b ) , [ fig : bulkgaugesys](d ) , and [ fig : bulkgaugesys](f ) ] . localization becomes sharp as @xmath191 increases ; this situation is similar to the continuum case .
the warp - suppressed spectra of kk excited modes are also realized [ figs .
[ fig : bulkgaugesys](a ) , [ fig : bulkgaugesys](c ) , and [ fig : bulkgaugesys](e ) ] . for a larger @xmath192 ( the number of gauge groups )
, the model leads to a spectrum more in agreement with the continuous rs case .
note , however , that the localization profiles of wave functions can be seen even with a rather small @xmath192 .
it is interesting that even with a finite number of gauge groups the massive modes have warp - suppressed spectra and localization profiles in the index space of gauge theory .    for fermion fields
, there is another interesting issue to be examined .
it is the localization behavior via dependence on the mass parameters @xmath102 , which was discussed in sec.[sec : spinor ] .
we show the @xmath102 dependence of the zero - mode wave function @xmath193 in fig .
[ fig : cdepwf ] .
the figure indicates that the zero - mode wave functions surely give the expected localization nature of the continuum rs limit [ eq .  ( [ eq : hzero ] ) in the appendix ] .
we find that the values of the wave functions are exponentially suppressed at the tail of localization profile even with a finite number of gauge groups .
the profiles of massive modes can also be reproduced .
now we apply our formulation to phenomenological problems in four dimensions .
let us use the localization behavior , which has been shown above , to obtain the yukawa hierarchy .
this issue has been studied in the 5d rs framework  @xcite .
we consider a model corresponding to the ( supersymmetric ) standard model in the bulk .
the yukawa couplings for quarks are given by @xmath194 where @xmath195 , @xmath196 , and @xmath197 denote the left - handed quarks and the right - handed up and down quarks , respectively , and @xmath198 are the family indices .
for simplicity , we study a supersymmetric case and introduce two types of higgs scalars @xmath199 and @xmath200 .
then the mass parameters of the higgs scalars satisfy eq .
( [ eq : abchypn ] ) and they are denoted by @xmath201 in the following .
similarly , the quark behaviors are described by their mass parameters @xmath202 .
we assume @xmath203 . generally , in supersymmetric 5d models , yukawa couplings such as eq .
( [ eq : yukawaorg ] ) are prohibited by 5d supersymmetry . however , since the present model is 4d , one may apply 5d - like results to yukawa couplings without respecting 5d consistency .
this is one of the benefits of our scheme .
we are now interested in the zero - mode part of eq .
( [ eq : yukawaorg ] ) , which generates the following mass terms @xmath204 where the fields with tildes @xmath205 stand for the @xmath186th mass eigenstate given by @xmath206 ( similarly for @xmath196 , @xmath197 , and @xmath207 ) .
the effective yukawa couplings are @xmath208 and similarly for @xmath209 . a typical behavior of @xmath210 is shown in fig .  [ fig : cdepwf ] for several values of the bulk mass parameter @xmath102 . in fig .
[ fig : cdepy ] , we show the behaviors of the zero - mode yukawa couplings @xmath211 against the quark mass parameters .
two limiting cases with @xmath212 and 1 are shown .
the former corresponds to a bulk higgs scalar localized at @xmath65 and the latter to one at @xmath121 in the continuum rs limit . from the figures
, we see that if there is a @xmath213 difference of mass ratio among the generations , it generates a large hierarchy between yukawa couplings .
combined with the mechanisms that control mass parameters discussed in the next section , one obtains a hierarchy without symmetries within the four - dimensional framework .
\(a ) @xmath212    \(b ) @xmath214    \(c ) @xmath212
\(d ) @xmath214    in the case of @xmath214 , the yukawa coupling depends exponentially on the quark bulk mass parameters @xmath202 when @xmath215 .. for @xmath216 , the higgs scalars have a peak at @xmath217 ( @xmath121 ) .
this situation is different from the one discussed in ref .
@xcite where the higgs field is localized at @xmath218 ( @xmath65 ) . ]
this implies that if @xmath202 exist in this region one obtains the following form of the yukawa matrices : @xmath219 where their exponents satisfy @xmath220 this form is similar to that obtained by the froggatt - nielsen mechanism  @xcite with a @xmath50 symmetry . as an illustration ,
let us take the following mass parameters @xmath221 and @xmath222 and @xmath223 , which generates the low - energy yukawa matrices @xmath224 where @xmath225 .
this pattern of quark mass textures leads to realistic quark masses and mixing angles  @xcite with a large value of the ratio @xmath226 . if the above analysis were extended to su(5 )
grand unified theory , realistic lepton masses and mixing may be derived .
other forms of yukawa matrices that may be realized by the froggatt - nielsen mechanism are easily incorporated in our formulation .    for more complicated patterns of mass parameters , we could realize yukawa matrices that are different from those derived from the froggatt - nielsen mechanism . in general ,
off - diagonal entries tend to be rather suppressed , that is , we have @xmath227 for the yukawa matrix ( [ yij ] )
. such a form may lead to realistic fermion masses and mixing angles .
for example , one could derive the yukawa matrix @xmath228 if initial values of @xmath229 are sufficiently suppressed . in this case
, the @xmath230 submatrix for the second and third generations does not satisfy eq .
( [ fntype ] ) . the yukawa matrix ( [ eq : compli ] )
may be relevant to the down - quark sector , indeed studied in ref .
we do not pursue further systematic studies on these types of yukawa matrices in this paper .
we have shown that 4d models with nonuniversal vevs and gauge and other couplings can describe 5d physics on curved backgrounds , including the rs model with an exponential warp factor . in the continuum 5d
theory , this factor is derived as a solution of the equation of motion for gravity . on the other hand , in the 4d viewpoint ,
warped geometries are generated by taking the couplings and vevs as appropriate forms . in the previous sections ,
we have just assumed their typical forms and examined its consequences .
if one could identify how to control these couplings by the underlying _ dynamics _ , the resultant 4d theories turn out to provide attractive schemes to discuss low - energy physics such as tiny coupling constants .
first we consider the scalar vevs @xmath231 . a simple way to dynamically control them is to introduce additional strongly coupled gauge theories  @xcite .
consider the following set of asymptotically free gauge theories : @xmath232 where @xmath233 and @xmath234 denote the dynamical scales .
we have , for simplicity , assumed common values of @xmath235 and @xmath233 for all @xmath6 .
in addition , two types of fermions are introduced : @xmath236 where their representations under @xmath237 gauge groups are shown in parentheses . if @xmath238 , the @xmath239 theories are confined at a higher scale than @xmath6 , and the fermion bilinear composite scalars @xmath240 appear .
their vevs @xmath25 are given by the dynamical scales @xmath234 of the @xmath239 gauge theories through a dimensional transmutation as @xmath241 where @xmath242 is a universal one - loop gauge beta function for @xmath239 ( @xmath243 ) .
the gauge couplings @xmath244 generally take different values and thus lead to different values of @xmath25 .
for example , a linear dependence of @xmath245 on the index @xmath28 is amplified to an exponential behavior of @xmath25 .
that is , @xmath246 which reproduces the bulk fields on the rs background as shown before .
the index dependences of the gauge couplings are actually generic situations , and may also be controlled , for example , by some mechanism fixing dilatons or the radiatively induced kinetic terms discussed below .
a supersymmetric extension of the above scenario is achieved with quantum - deformed moduli spaces  @xcite .
another mechanism that dynamically induces nonuniversal vevs is obtained in supersymmetric cases .
consider the gauge group @xmath247 and the chiral superfields @xmath9 with charges @xmath248 under @xmath249 .
it is assumed that the scalar components @xmath250 of @xmath9 develop their vevs @xmath251 .
the @xmath154 term of each @xmath252 is given by @xmath253 where @xmath254 is the coefficient of the fayet - iliopoulos ( fi ) term , and the ellipsis denotes contributions from other fields charged under @xmath252 , which are assumed not to have vevs .
given nonvanishing fi terms , @xmath255 , the @xmath154-flatness conditions mean @xmath256 and nonuniversal vevs @xmath25 are indeed realized . in this case
, the dynamical origin of nonuniversal vevs is the nonvanishing fi terms .
these may be generated at the loop level .
furthermore , if the matter content is different for each gauge theory , the @xmath254 themselves have complicated forms .
above , we supposed that the charges of @xmath9 are @xmath248 under @xmath257 .
alternatively , if @xmath9 have charges @xmath258 under @xmath257 and other matter fields have integer charges , the gauge symmetry @xmath247 is broken to the product of a diagonal @xmath50 gauge symmetry and the discrete gauge symmetry @xmath259
. such discrete gauge symmetry would be useful for phenomenology  @xcite .
models with nonuniversal gauge couplings @xmath7 are also interesting in the sense that they can describe the localization of massless vector fields .
a nonuniversality of gauge couplings is generated , e.g. ,  in the case that the @xmath6 gauge theories have different matter content from each other .
then radiative corrections to gauge couplings and their renormalization - group running become nonuniversal , even if their initial values are equal .
this fact is also applicable to the above - mentioned mechanism for nonuniversal @xmath25 .
suppose that the @xmath239 theory contains ( @xmath260 ) vectorlike quarks which decouple at @xmath261 .
the gauge couplings @xmath262 are then determined by @xmath263 where we have assumed that the @xmath239 theories are strongly coupled at a high - energy scale @xmath264  ( @xmath265 ) .
tuning of the relevant matter content thus generates the desired linear dependence of @xmath245 . with these radiatively induced couplings ( [ eq : nonuni - g ] ) at hand , the vevs are determined from eq .
( [ eq : dtvev ] ) : @xmath266
we have formulated 4d models that provide 5d field theories on generic warped backgrounds .
the warped geometries are achieved with generic values of symmetry - breaking vevs , gauge couplings , and other couplings in the models .
we focused on field localization behaviors along the index space of gauge theory ( the fifth dimension in the continuum limit ) , which is realized by taking relevant choices of the mass parameters .
as a good and simple application , we constructed 4d models corresponding to bulk field theories on the ads@xmath0 randall - sundrum background .
the localized wave functions of massless modes are completely reproduced with a finite number of gauge groups .
in addition , the exponentially suppressed spectrum of the kk modes is also generated .
these results imply that most properties of brane world models can be obtained within 4d gauge theories .
supersymmetric extensions were also investigated . in 5d warped models ,
the bulk and boundary mass terms of spinors and scalars satisfy complicated forms imposed by supersymmetry on the rs background .
however , we show in our formalism that these forms of the mass terms are derived from a 4d _ global _ supersymmetric model on a _
flat _ background .    as an application of our 4d formulation
, we derived hierarchical forms of yukawa couplings .
the zero modes of scalars and spinors with different masses have different wave - function profiles as in the 5d rs cases .
therefore by varying the @xmath213 mass parameters for each generation , one can obtain realistic yukawa matrices with a large hierarchy from the overlaps of the wave functions in a purely 4d framework .
other phenomenological issues such as proton stability , grand unified theory ( gut ) symmetry breaking , and supersymmetry breaking can also be discussed .
the conditions on the model parameters should be explained by some dynamical mechanisms if one considers the models from a fully 4d viewpoint .
one interesting way is to include additional strongly coupled gauge theories . in this case , a small @xmath213 difference between gauge couplings is converted to exponential profiles of symmetry - breaking vevs via dimensional transmutation , and indeed generates a warp factor of the rs model . a difference of gauge couplings
is achieved by , for example , the dynamics controlling dilatons , or radiative corrections to gauge couplings .
supersymmetrizing models provide a mechanism for dynamically realizing nonuniversal vevs with @xmath154-flatness conditions .
our formulation makes sense not only from the 4d points of view but also as a lattice - regularized 5d theory . in this sense , effects such as the ads / conformal field theory ( cft ) correspondence might be clearly seen with our formalism . as another application , it can be applied to construct various types of curved backgrounds and bulk or boundary masses .
for example , we discussed massless vector localization by varying the gauge couplings @xmath7 .
furthermore , one might consider models in which some fields are charged under only some of the gauge groups .
these seem not like bulk or brane fields , but `` quasi - bulk '' fields .
applications including these phenomena will be studied elsewhere .
this work is supported in part by the japan society for the promotion of science under the postdoctoral research program ( grants no .  @xmath267 and no .
@xmath268 ) and a grant - in - aid for scientific research from the ministry of education , science , sports and culture of japan ( no .  @xmath269 ) .
here we briefly review the field theory on a rs background , following ref .
one of the original motivations for introducing a warped extra dimension by randall and sundrum is to provide the weak planck mass hierarchy via the exponential factor in the space - time metric .
this factor is called `` warp factor , '' and the bulk space a `` warped extra dimension '' .
such a nonfactorizable geometry with a warp factor distinguishes the rs brane world from others .
consider the fifth dimension @xmath41 compactified on an orbifold @xmath270 with radius @xmath271 and two three - branes at the orbifold fixed points @xmath121 and @xmath272 .
the einstein equation for this five - dimensional setup leads to the solution  @xcite @xmath273 where @xmath55 is a constant with mass dimension @xmath56 .
let us study a vector field @xmath274 , a dirac fermion @xmath275 , and a complex scalar @xmath129 in the bulk specified by the background metric ( [ metric ] ) .
the 5d action is given by @xmath276 , \label{kin}\ ] ] where @xmath277 and the covariant derivative is @xmath278 where @xmath279 is the spin connection given by @xmath280 and @xmath281 . from the transformation properties under @xmath282 parity , the mass parameters of scalar and fermion fields are parametrized as is taken as @xmath283 .
here we adopt @xmath284 , and then the boundary mass parameter @xmath151 in eq .
( [ eq : smass ] ) is different from that in ref .
@xcite by the factor @xmath285 . ] @xmath286 where @xmath150 , @xmath151 , and @xmath102 are dimensionless parameters .
referring to @xcite , the vector , scalar and spinor fields are cited together using the single notation @xmath287 .
the kk mode expansion is performed as @xmath288 by solving the equations of motion , the eigenfunction @xmath289 is given by @xmath290 ,    \label{eq : modefunction}\end{aligned}\ ] ] where @xmath291 , @xmath292 , and @xmath293 for each component in @xmath294 .
@xmath295 is the normalization factor and @xmath296 and @xmath297 are the bessel functions .
the corresponding kk spectrum @xmath298 is obtained by solving @xmath299    a supersymmetric extension of this scenario was discussed in  @xcite .
the on - shell field content of a vector supermultiplet is @xmath300 where @xmath274 , @xmath301 , and @xmath302 are the vector , two majorana spinors , and a real scalar in the adjoint representation , respectively .
also a hypermultiplet consists of @xmath303 , where @xmath304 @xmath305 are two complex scalars and @xmath275 is a dirac fermion . by requiring the action ( [ kin ] ) to be invariant under supersymmetric transformation on the warped background ,
one finds that the five - dimensional masses of the scalar and spinor fields have to satisfy @xmath306 where @xmath102 remains as an arbitrary dimensionless parameter .
that is , @xmath307 , @xmath308 , and @xmath114 for vector multiplets and @xmath309 and @xmath310 for hypermultiplets .
there is no freedom to choose the bulk masses for vector supermultiplets and only one freedom parametrized by @xmath102 for the bulk hypermultiplets .
it should be noted that in warped 5d models fields contained in the same supermultiplet have different bulk and boundary masses .
that is in contrast with the flat case .
the @xmath126 even components in supermultiplets have massless modes with the following wave functions : @xmath311 \frac{e^{(1/2- c)\sigma}}{n_0\sqrt{2\pi r } } & & \textrm { for } \;h^{1\ , ( 0 ) } \ , \textrm { and } \,\psi^{(0)}_l .
\label{eq : hzero}\end{aligned}\ ] ] the subscript @xmath76 means the left - handed ( @xmath126 even ) component .
the massless vector multiplet has a flat wave function in the extra dimension .
on the other hand , the wave function for massless chiral multiplets involves a @xmath41-dependent contribution from the space - time metric , which induces a localization of the zero modes . the zero modes with masses @xmath119 and @xmath120 localize at @xmath121 and @xmath312 , respectively .
the case with @xmath114 corresponds to the conformal limit where the kinetic terms of the zero modes are independent of @xmath41 .
l.  randall and r.  sundrum , phys .
lett . * 83 * , 3370 ( 1999 ) [ hep - ph/9905221 ] ; phys .
* 83 * , 4690 ( 1999 ) [ hep - th/9906064 ] .
w.  d.  goldberger and m.  b.  wise , phys .
d * 60 * , 107505 ( 1999 ) [ hep - ph/9907218 ] ; phys .
lett . * 83 * , 4922 ( 1999 ) [ hep - ph/9907447 ] ; h.  davoudiasl , j.  l.  hewett and t.  g.  rizzo , phys .
b * 473 * , 43 ( 2000 ) [ hep - ph/9911262 ] ; a.  pomarol , phys .
b * 486 * , 153 ( 2000 ) [ hep - ph/9911294 ] ; s.  chang , j.  hisano , h.  nakano , n.  okada and m.  yamaguchi , phys .  rev .
d * 62 * , 084025 ( 2000 ) [ hep - ph/9912498 ] ; b.  bajc and g.  gabadadze , phys .
b * 474 * , 282 ( 2000 ) [ hep - th/9912232 ] ; e.  shuster , nucl .
b * 554 * , 198 ( 1999 ) [ hep - th/9902129 ] .
t.  gherghetta and a.  pomarol , nucl .
b * 586 * , 141 ( 2000 ) [ hep - ph/0003129 ] .
t.  gherghetta and a.  pomarol , nucl .
b * 602 * , 3 ( 2001 ) [ hep - ph/0012378 ] .
s.  j.  huber and q.  shafi , phys .
b * 498 * , 256 ( 2001 ) [ hep - ph/0010195 ] .
h.  abe and t.  inagaki , phys .
d * 66 * , 085001 ( 2002 ) [ hep - ph/0206282 ] .
h.  abe , k.  fukazawa and t.  inagaki , prog .
* 107 * , 1047 ( 2002 ) [ hep - ph/0107125 ] ; h.  abe , t.  inagaki and t.  muta , in _ fluctuating paths and fields _ , edited by w.  janke , a.  pelster , h .- j .
schmidt , and m.  bachmann ( world scientific , singapore , 2001 ) [ hep - ph/0104002 ] .
i.  i.  kogan , s.  mouslopoulos , a.  papazoglou and g.  g.  ross , nucl .
b * 615 * , 191 ( 2001 ) [ hep - ph/0107307 ] .
s.  dimopoulos , s.  kachru , n.  kaloper , a.  e.  lawrence and e.  silverstein , phys .
d * 64 * , 121702 ( 2001 ) [ hep - th/0104239 ] .
n.  arkani - hamed , a.  g.  cohen and h.  georgi , phys .
* 86 * , 4757 ( 2001 ) [ hep - th/0104005 ] .
c.  t.  hill , s.  pokorski and j.  wang , phys .
d * 64 * , 105005 ( 2001 ) [ hep - th/0104035 ] .
h.  c.  cheng , c.  t.  hill , s.  pokorski and j.  wang , phys .
d * 64 * , 065007 ( 2001 ) [ hep - th/0104179 ] ; a.  sugamoto , prog .
theor .  phys . *
107 * , 793 ( 2002 ) [ hep - th/0104241 ] ; h.  c.  cheng , c.  t.  hill and j.  wang , phys .
d * 64 * , 095003 ( 2001 ) [ hep - ph/0105323 ] ; n.  arkani - hamed , a.  g.  cohen and h.  georgi , phys .
b * 513 * , 232 ( 2001 ) [ hep - ph/0105239 ] ; h.  c.  cheng , d.  e.  kaplan , m.  schmaltz and w.  skiba , phys .
b * 515 * , 395 ( 2001 ) [ hep - ph/0106098 ] ; m.  bander , phys .  rev .
d * 64 * , 105021 ( 2001 ) [ hep - th/0107130 ] ; c.  csaki , g.  d.  kribs and j.  terning , phys .  rev .
d * 65 * , 015004 ( 2002 ) [ hep - ph/0107266 ] ; h.  c.  cheng , k.  t.  matchev and j.  wang , phys .
b * 521 * , 308 ( 2001 ) [ hep - ph/0107268 ] ; n.  arkani - hamed , a.  g.  cohen and h.  georgi , hep - th/0108089 ; c.  t.  hill , phys .
rev .  lett . * 88 * , 041601 ( 2002 ) [ hep - th/0109068 ] ; n.  arkani - hamed , a.  g.  cohen and h.  georgi , jhep * 0207 * , 020 ( 2002 ) [ hep - th/0109082 ] ; i.  rothstein and w.  skiba , phys .  rev .
d * 65 * , 065002 ( 2002 ) [ hep - th/0109175 ] ; t.  kobayashi , n.  maru and k.  yoshioka , hep - ph/0110117 ; n.  arkani - hamed , a.  g.  cohen , d.  b.  kaplan , a.  karch and l.  motl , hep - th/0110146 ; c.  csaki , j.  erlich , v.  v.  khoze , e.  poppitz , y.  shadmi and y.  shirman , phys .
d * 65 * , 085033 ( 2002 ) [ arxiv : hep - th/0110188 ] ; g.  c.  cho , e.  izumi and a.  sugamoto , hep - ph/0112336 ; e.  witten , hep - ph/0201018 ; w.  skiba and d.  smith , phys .
d * 65 * , 095002 ( 2002 ) [ hep - ph/0201056 ] ; r.  s.  chivukula and h.  j.
he , phys .
b * 532 * , 121 ( 2002 ) [ hep - ph/0201164 ] ; n.  arkani - hamed , a.  g.  cohen , t.  gregoire and j.  g.  wacker , jhep * 0208 * , 020 ( 2002 ) [ hep - ph/0202089 ] ; k.  lane , phys .
d * 65 * , 115001 ( 2002 ) [ hep - ph/0202093 ] ; z.  berezhiani , a.  gorsky and i.  i.  kogan , jetp  lett .
* 75 * , 530 ( 2002 ) ; pisma  zh .
.  fiz . * 75 * , 646 ( 2002 ) [ hep - th/0203016 ] ; a.  falkowski , c.  grojean and s.  pokorski , phys .
b * 535 * , 258 ( 2002 ) [ hep - ph/0203033 ] ; c.  t.  hill and a.  k.  leibovich , phys .
d * 66 * , 016006 ( 2002 ) [ hep - ph/0205057 ] . c.  csaki , j.  erlich , c.  grojean and g.  d.  kribs , phys .
d * 65 * , 015003 ( 2002 ) [ hep - ph/0106044 ] .
k.  sfetsos , nucl .
b * 612 * , 191 ( 2001 ) [ hep - th/0106126 ] .
a.  kehagias and k.  tamvakis , phys .
b * 504 * , 38 ( 2001 ) [ hep - th/0010112 ] . c.  d.  froggatt and h.  b.  nielsen , nucl .
b * 147 * , 277 ( 1979 ) .
g.  altarelli and f.  feruglio , phys .
b * 451 * , 388 ( 1999 ) [ hep - ph/9812475 ] .
p.  ramond , r.  g.  roberts and g.  g.  ross , nucl .
b * 406 * , 19 ( 1993 ) [ hep - ph/9303320 ] .
l.  e.  ibanez and g.  g.  ross , nucl .
b * 368 * , 3 ( 1992 ) ; phys .  lett .
b * 260 * , 291 ( 1991 ) .
j.  wess and j.  bagger , _ supersymmetry and supergravity _
( princeton university press , princeton , nj , 1992 )


the dynamical density  matrix renormalization group ( dmrg )  @xcite and the closely related correction  vector dmrg  @xcite have been widely used in the last decade to compute the dynamical correlation functions and spectral functions of low - dimensional strongly correlated quantum systems .
@xcite although more powerful dmrg approaches have been developed recently ,  @xcite dynamical dmrg ( ddmrg ) often remains the method of choice because it offers two practical advantages over the other approaches : it is simpler and it can be easily parallelized .
for instance , it has been recently shown that ddmrg allows us to investigate features with small spectral weights such as power - law pseudo - gaps in luttinger liquids .
@xcite the main drawback of ddmrg is that it always yields the convolution of the desired spectrum with a lorentzian distribution of finite width .
therefore , the true spectrum can only be obtained through a deconvolution of the ddmrg spectrum .
( in principle , there are some methods to get around this problem  @xcite but they are rarely used in practice . )    deconvolution is a typical ill - conditioned inverse problem , however .
@xcite a direct solution of the deconvolution equation usually yields a very noisy and thus useless spectrum . nevertheless , various regularization methods have been successfully used to deconvolve ddmrg spectra for one - dimensional systems and quantum impurity problems .
@xcite astonishingly , some of these deconvolution methods even allow us to bypass the finite - size scaling analysis and to obtain the piecewise smooth spectrum of an infinite systems directly from a broadened finite - system ddmrg spectrum .
unfortunately , regularization also smooths out the sharp features of the true spectrum .
this is a serious issue as the spectra of one - dimensional systems and quantum impurities often exhibit very interesting ( power - law ) singularities .    in this paper
we present a method , which allows us to determine sharp spectral features in the thermodynamic limit starting from a broadened finite - system ddmrg spectrum . for this purpose
we consider the extrapolation to the thermodynamic limit and the deconvolution for the lorentzian kernel to be a single blind deconvolution ,  @xcite i.e.  an inverse problem with an unknown kernel including both the lorentzian broadening and the finite - size effects . the key idea to preserve sharp spectral features in a piecewise smooth spectrum is to impose a minimal distance @xmath0 between extrema of the deconvolved spectrum . to illustrate our method we investigate the single - particle density of states ( dos ) of one - dimensional paramagnetic mott insulators represented by the half - filled hubbard model .
@xcite we confirm that this dos has the step - like onset predicted by field - theoretical studies  @xcite at least at weak to intermediate coupling up to @xmath1 .
the hubbard model  @xcite with on - site interaction @xmath2 and nearest - neighbor hopping @xmath3 is a basic lattice model for the physics of strongly interacting electrons , in particular the mott metal - insulator transition .
@xcite at half filling ( i.e. , the number of electrons equals the number of sites @xmath4 ) the ground state is a mott insulator for strong interaction @xmath5 , while it is a fermi gas in the non - interacting limit @xmath6 .
the hamiltonian of the hubbard model is defined by @xmath7 where the operator @xmath8 ( @xmath9 ) creates ( annihilates ) an electron with spin @xmath10 on the site @xmath11 , @xmath12 , and @xmath13 .
the first sum runs over all pairs @xmath14 of nearest - neighbor sites while the other two sums run over all sites @xmath15 .
here we will only consider half - filled systems and thus set the chemical potential @xmath16 to have electron - hole - symmetric spectra and a fermi energy @xmath17 .
the bulk single - particle dos @xmath18 can be measured experimentally using photoemission spectroscopy or scanning tunneling spectroscopy .
theoretically , it can be defined as the average of the local dos @xmath19 where the sum runs over both spins and all sites @xmath15 in the lattice , while @xmath20 is the local single - particle dos at site @xmath15 for spin @xmath10 and can be calculated using @xmath21 for @xmath22 and @xmath23 for @xmath24 . here
@xmath25 denotes the eigenstates of the hamiltonian @xmath26 and @xmath27 their eigenenergies in the fock space .
the ground state for the chosen number of particles corresponds to @xmath28 . the total spectral weight is @xmath29 we will consider only lattice geometry for which the hamiltonian  ( [ eq : hamiltonian ] ) is invariant under the electron - hole transformation @xmath30 .
therefore , for half filling the density of states is symmetric , @xmath31 . if the system is translation invariant , the bulk dos and the local dos are identical . for dmrg simulations , however , open boundary conditions are preferred to periodic boundary conditions . in this case
, the bulk dos can be identified with the local dos on one of the two equivalent middle sites of the system , i.e.  as far as possible from the system boundaries .
inverse problems such as ( blind ) deconvolutions  @xcite occur in many scientific fields and are among the most challenging numerical computations .
experimental measurements and computer simulations often yield approximations of the true quantities which are measured or computed , respectively .
it is often assumed that the deviations from exact results can be modelled by a convolution with a smoothing function and an additive noise due to the finite accuracy and resolution of the measurement or simulation process .
a typical example of a blind deconvolution is the reconstruction of an original signal from a degraded copy using incomplete information about the degradation process .
@xcite here we want to compute sharp spectral features in the piecewise smooth spectrum of an infinite system from a broadened finite - system spectrum calculated with ddmrg . in this section
we first show that this task can be formulated as a blind deconvolution problem , then present an algorithm for solving it .
let @xmath32 be a spectrum of a finite lattice model with @xmath4 sites .
this spectrum is a dirac - comb ( a finite sum of dirac - peaks ) @xmath33 where the sum runs over all hamiltonian eigenstates @xmath25 which contributes to the spectrum , i.e.  with a nonzero spectral weight @xmath34 . here
@xmath35 denotes the corresponding excitation energies .
this spectrum can be broadened with a lorentzian distribution of width @xmath36 @xmath37 to obtain a smooth spectral function @xmath38 with the ddmrg method we can calculate this spectrum for a discrete set of excitation energies @xmath39 . as numerical calculations are always affected by errors , ddmrg actually yields values @xmath40 which are related to the true spectral function by @xmath41 for @xmath42 , where @xmath43 represents the unknown errors .
( it should be noted that ddmrg errors @xmath43 include significant systematic contributions , for instance due to the variational nature of the procedure .
@xcite ) in principle , one could determine the true spectrum , i.e. , the excitation energies @xmath35 and the corresponding weights @xmath44 , through this system of equations . in practice
, however , this is an ill - conditioned problem except for simple discrete spectra .
moreover , we are not interested in resolving the discrete peaks of small systems but in calculating the piecewise smooth spectra of macroscopic systems .
the spectrum in the thermodynamic limit is given by @xmath45 note that , generally , the order of the two limits can not be exchanged .
typically , the spectral function @xmath46 is piecewise smooth , i.e. , it exhibits one or more continua as well as isolated sharp features such as steps , power - law singularities or cusps . in principle
, one should carry out several ddmrg simulations with varying system size @xmath4 and broadening @xmath36 and then extrapolate the numerical data to obtain @xmath46 . in most cases ,
a simultaneous extrapolation for @xmath47 and @xmath48 is possible  @xcite using a constant value of @xmath49 .
nevertheless , the computational cost of ddmrg simulations increases very rapidly with smaller @xmath36 and the overall cost of this approach is prohibitive for a full spectrum .
indeed , this approach has been mostly used to study isolated spectral features in the thermodynamic limit such as power - law singularities and steps .
@xcite    as all operations used to define @xmath46 from @xmath50 are linear , the broadened spectrum of the finite system can also be written explicitly as a function of the infinite system spectrum @xmath51 the kernel @xmath52 includes both the finite - size effects and the lorentzian smoothing .
its form is not known but it is clear that we must recover a pure lorentzian smoothing in the thermodynamic limit @xmath53 combining eqs .
( [ eq : deconv ] ) and  ( [ eq : spec - inf ] ) we obtain a system of equations @xmath54 for @xmath55 , relating the ddmrg data set @xmath56 to the infinite system spectrum @xmath46 .
determining @xmath46 from these equations is a so - called inverse problem .
@xcite this kind of problem is also called blind deconvolution since our knowledge of the kernel is incomplete .
[ strictly speaking , it is not a deconvolution because eq .
( [ eq : inverse - problem ] ) is not a convolution .
however , as the kernel approaches the form @xmath57 in the thermodynamic limit , we will use the terminology of deconvolution problems . ]
it should be obvious that this is an ill - posed problem .
first , the errors @xmath43 and the kernel @xmath52 are not known .
second , the problem is sorely underdetermined as we try to reconstruct the function of a continuous variable from a finite number @xmath58 of data points .
finally , a convolution with a lorentzian is a smoothing operation and thus the corresponding deconvolution is an extremely ill - conditioned inverse problem : the solution will be extremely sensitive to small changes or errors in the input .
various deconvolution methods have been used successfully to deduce piecewise smooth spectra from the broadened finite - system spectra calculated with ddmrg .
they include , direct inversion at low resolution ,  @xcite linear regularization methods ,  @xcite fourier transform with low - pass filtering ,  @xcite nonlinear regularization methods such as the maximum entropy method ,  @xcite parametrization with piecewise polynomial functions ,  @xcite and a deconvolution ansatz for the self energy .
@xcite however , this task has not been viewed as a blind deconvolution so far .
instead , it has been considered as the deconvolution of a perfectly known kernel .
the need for regularization or filtering techniques has been viewed as the consequence ill - conditioning and under - determination of the problem  ( [ eq : inverse - problem ] ) with a lorentzian kernel .
all of these methods offer some advantages for particular spectral forms . however , their common drawback is that they are ill - suited for sharp spectral features , such as steps or power - law singularities , within or at the edge of a continuum .
either the regularization procedure smooths out true sharp features excessively or it allows the occurrence of deconvolution artifacts ( artificial sharp structures , rapid oscillations or negative spectral weight ) , especially in the vicinity of the true spectrum singularities .
naturally , better results can be obtained if we can use _ a priori _ knowledge about the properties of the spectrum  @xcite but , in practice , this is a rare occurrence .
therefore , we need a better method for solving the inverse problem  ( [ eq : inverse - problem ] ) which allows us to determine isolated sharp spectral features accurately while preserving the positivity and the piecewise smoothness of @xmath46
.    let the ddmrg data ( [ eq : dmrg - data ] ) be evenly distributed in the energy interval @xmath59 $ ] .
the difference between two consecutive energies is @xmath60 .
additionally , consider a set of equidistant energies @xmath61 in the interval @xmath62 \subset [ \epsilon_a,\epsilon_b]$ ] .
the distance between these energies is @xmath63 . as we will always use @xmath64 , we have @xmath65 .
( typical values are @xmath66 and @xmath67 . ) as in a least - square approach we define a cost function @xmath68 as the sum of the squares of the differences between the ddmrg data and an approximate representation parametrized by a discrete set of variables @xmath69 @xmath70 the absolute minimum of @xmath68 is zero and the corresponding parameters @xmath71 are determined by a linear system of @xmath58 equations @xmath72 using this equation system to determine the parameters @xmath71 would be an unconstrained least - square fit .    in the limit @xmath73 ( followed by @xmath74 and @xmath75 )
this equation system becomes equivalent to the inverse problem  ( [ eq : inverse - problem ] ) with vanishing errors @xmath76 .
thus the absolute minimum of @xmath68 yields the spectrum @xmath46 through @xmath77    if we substitute a lorentz kernel @xmath78 for the unknown kernel , @xmath79 in  ( [ eq : linear - system ] ) , we recover the finite - system deconvolution problem defined by equations  ( [ eq : spec - finite ] ) and  ( [ eq : deconv ] ) for vanishing errors @xmath43 .
thus the absolute minimum of the cost function corresponds to the discrete finite - system spectrum @xmath32 through @xmath80 if @xmath81 and @xmath82 otherwise .
physically , the solution of the deconvolution problem is unique for vanishing errors @xmath43 and thus the cost function should have a unique absolute minimum . from a mathematical point of view , however , the equation system  ( [ eq : linear - system ] ) could have no solution or infinitely many solutions . then any small error @xmath43 can generate wildly different ( and mostly unphysical ) solutions .    therefore , as  ( [ eq : kernel - limit ] ) holds in the thermodynamic limit , it is possible and preferable to obtain a reasonable approximation of the infinite - system spectrum @xmath46 from the minimization of the cost function  ( [ eq : cost - function ] ) with a lorentz kernel under the constraint that the spectral function @xmath46 is physically allowed .
for instance , @xmath46 should be positive semidefinite and piecewise smooth .
generally , this solution does not correspond to the absolute minimum or even a local minimum of @xmath68 .
indeed , the solution of the inverse problem  ( [ eq : inverse - problem ] ) corresponds to the value @xmath83 if we assume that the relation  ( [ eq : solution ] ) holds .
of course , it could be possible to lower the cost function with other configurations @xmath84 but in that case the relation  ( [ eq : solution ] ) would no longer hold .
note that , this idea has been implicitly assumed in all previous deconvolution schemes of ddmrg data aiming at piecewise smooth spectra so far .
however , in these approaches the agreement between solution @xmath84 and numerical data @xmath85 in eq .
( [ eq : linear - system ] ) with a lorentzian kernel is considered essential while the regularization of the solution and the errors @xmath43 are seen as perturbations which should deteriorate the agreement as little as possible .    yet a blind deconvolution requires equal balancing of the agreement between solution and numerical data and of the smoothness and stability of the solution .
@xcite hence we must take a different point of view : the lorentzian kernel  ( [ eq : kernel - limit ] ) is only an approximation of the true kernel @xmath52 and the physical constraints on the deconvolved spectrum @xmath46 are essential in the minimization of the cost function . thus we accept significant deviations of @xmath84 from the conditions  ( [ eq : linear - system ] ) yielding the absolute minimum of the cost function , or , equivalently , we assume that the errors @xmath43 can be substantial .      therefore , the blind deconvolution problem can be formulated as a least - square optimization under non - linear constraints .
we want to minimize the cost function  ( [ eq : cost - function ] ) with the lorentz kernel under the constraints that the spectrum @xmath46 has the following properties :    1 .
finite band width , i.e @xmath86 for all @xmath87 and all @xmath88 for some finite @xmath89 , 2 .
positive semi - definite , @xmath90 , and 3 .
piecewise smooth .
the first two conditions are easily expressed for the parameters @xmath91 .
the somewhat fuzzy concept of a piecewise smooth spectrum must now be formulated more precisely . in principle
, we wish that @xmath46 is piecewise continuous and that the distance between discontinuities is larger than a minimal energy difference @xmath0 .
as we must work with a finite number @xmath92 of points @xmath93 , we only have a discrete representation of @xmath46 and we have to formulate a `` continuity '' condition for the discrete set of variables as well .
therefore , we require that the distance between two significant extrema is larger than a parameter @xmath94
. two neighboring extrema at energy @xmath95 and @xmath96 are significant if there relative height difference is larger than a parameter @xmath97 , @xmath98 this condition can easily be formulated for the parameters @xmath91 .
the minimal extremum distance @xmath0 must be chosen carefully .
it should be smaller than the distance between actual singularities in the spectrum @xmath46 but a too small value allows many artificial peaks in a deconvolved spectrum . in practice , we have found that we can obtain reasonable solutions to the blind deconvolution problem which look piecewise smooth using @xmath99 . in all examples discussed in this paper
every local extremum is considered to be significant ( i.e. , we have used the precision of floating - point arithmetic @xmath100 ) .
the cost function is minimized iteratively .
iterations are repeated until the procedure converges .
each iteration consists in two steps . in the first step
the cost function @xmath101 is minimized with respect to each variable @xmath102 successively .
this minimization under constraint does not present any difficulty as @xmath101 is a second - order polynomial in each variable @xmath103 . in the second step , we first find the positions @xmath104 of all significant extrema pairs in @xmath105 which are separated by less than a distance @xmath0
. then we interpolate the data @xmath105 linearly from @xmath106 to @xmath107 to smooth out the spectrum around the extrema . in doing
so we take care to preserve the total spectral weight @xmath108 the search for extrema and their smoothing is repeated until there is no more close significant extrema in @xmath109 .
then we start the next iteration . by design the first step results in a decrease of @xmath101 .
the second step nearly always results in an increase of @xmath101 . without the second step , however , we would perform an under - determined ( @xmath110 ) deconvolution devoid of any regularization mechanism and thus obtain a completely useless result .
typically , we observe a rapid and monotonic decrease of @xmath101 in the initial iterations followed by a saturation or oscillations in further iterations .
therefore , we monitor the changes in the parameters @xmath103 and the normalized cost function @xmath111 to determine converged configurations @xmath71 .
convergence requires typically @xmath112 to @xmath113 iterations depending on the quality of the ddmrg data and the complexity of the spectrum .
finally , the solution @xmath91 can be smoothened using a narrow lorentzian distribution to obtain a continuous function @xmath114 with @xmath115 .
alternatively , we can use a gaussian distribution  @xcite of width @xmath116 .
the second approach yields sharper ( real or artificial ) features because the tail of a gaussian distribution decreases faster than that of the lorentz distribution . as our minimization problem possesses many local minima , the final results @xmath71 depend somewhat on the criteria for convergence . however , if the final smoothening function is broad enough , the differences are canceled out .
if the ddmrg data  ( [ eq : dmrg - data ] ) are not evenly distributed in the interval @xmath59 $ ] or if they already exhibit numerous close extrema , it is useful to regularize them before starting the deconvolution iterations using an interpolation and the smoothening procedure described above .
the computational effort required by this procedure is negligible compared to the computational cost of the ddmrg simulations yielding the original data .
( our code in the programming language c contains less than 400 lines of instructions and the deconvolution of one spectrum takes less than 30 minutes on a single cpu . )
therefore , we have not bothered to optimize the algorithm .
nevertheless , it should be implemented in such a way that it only requires @xmath117 operations rather than the @xmath118 operations of a straightforward implementation .
the method described here can be generalized in several ways .
for instance , it is possible to use a variable spacing @xmath119 of the ddmrg data points , such as a finer mesh close to sharp spectrum features .
however , this does not seem to improve the results in practice because the broadening parameter @xmath36 , not @xmath119 , is the limiting scale .
a generalization to variable @xmath36 and @xmath119 , as proposed in ref .   for quantum impurity problems , should also be possible but we have not tested it yet . to introduce information about the variation of the spectrum with @xmath36 one could combine ddmrg data obtained for different values of @xmath36 by defining an overall cost function as the sum of the cost functions for each @xmath36 .
these generalizations will be tested in future works .
as an illustration of our deconvolution procedure we discuss its application to the dos of one - dimensional paramagnetic mott insulators .
the nature of mott insulators is a long - standing open problem in the theory of strongly correlated quantum systems .
@xcite in a paramagnetic mott insulator quantum fluctuations or frustration of the antiferromagnetic spin exchange coupling prevents the formation of a long - range magnetic order .
experimentally , non - magnetic mott insulators have been found in layered organic insulators  @xcite as well as in quasi - one - dimensional cuprate chains  @xcite and ladders  @xcite . despite decades of extensive research the properties of mott insulators , in particular their single - particle dos , are still poorly understood and thus actively investigated .
the half - filled hubbard model  @xcite with repulsive on - site interaction @xmath2 is a basic lattice model for describing mott insulators and the mott metal - insulator transition .
here we consider the case of the one - dimensional hubbard model , which is exactly solvable by bethe ansatz .
@xcite at half- filling it describes a paramagnetic mott insulator with a charge gap ( mott - hubbard gap ) @xmath120 for @xmath121 .
however , the dos can not be calculated directly from the bethe ansatz .
all ddmrg spectra used here have been calculated with a variable number of density - matrix eigenstates kept ( up to 512 ) to reach a discarded weight lower than @xmath122 and to check dmrg truncation errors .
typically , convergence was reached after three sweeps for each frequency interval of size @xmath36 .
the ddmrg method is presented in detail in ref .  .
for @xmath6 the exact dos of the tight - binding chain in the thermodynamic limit is @xmath123 for @xmath124 while @xmath125 vanishes for larger @xmath126 . for finite coupling @xmath127
the spectrum consists in two symmetric hubbard bands separated by the gap @xmath128 .
low - order strong - coupling perturbation theory  @xcite predicts a square - root divergence at the dos threshold for @xmath129 , namely @xmath130 for @xmath131 and @xmath132 otherwise with @xmath133 . because of this strong - coupling result and the result of the hartree - fock ( hf ) approximation it has often been assumed that the dos of one - dimensional mott insulators exhibits a square - root divergence at the spectrum onset @xmath134 like in a one - dimensional band insulator .
indeed , in the unrestricted hf approximation the one - dimensional half - filled hubbard model is an antiferromagnetic mott insulator for @xmath121 .
its dos is given by @xmath135 for @xmath136 and vanishes otherwise . here
@xmath137 is the hf gap . for @xmath6 , @xmath138 and
this dos reduces to the dos of the tight - biding chain  ( [ eq : dos - tb ] ) . for @xmath139 , @xmath140 and the hf dos shows a square - root divergence at the onset of the spectrum .
however , in the weak - coupling limit @xmath141 a field - theoretical analysis  @xcite predicts that the dos of one - dimensional mott insulators is constant above the threshold energy @xmath142 , @xmath143 at least for @xmath144 .
thus there is a discrepancy between the field - theoretical and strong - coupling predictions for the behavior of @xmath18 just above the threshold energy @xmath142 .
( color online ) dos of a tight - binding chain : ( a ) ddmrg spectrum for a @xmath145-site chain with a lorentzian broadening @xmath146 ( red dashed line ) and result of our deconvolution method with a gaussian broadening @xmath147 ( black line ) .
( b ) enlarged view close to the singularity at @xmath148 on a double logarithmic scale : exact results ( black long - dashed line ) , ddmrg data ( red short - dashed line ) , and the deconvolved spectra @xmath109 for @xmath149 ( black solid line ) and @xmath150 ( blue dash - dot line).,title="fig:",scaledwidth=48.0% ]   ( color online ) dos of a tight - binding chain : ( a ) ddmrg spectrum for a @xmath145-site chain with a lorentzian broadening @xmath146 ( red dashed line ) and result of our deconvolution method with a gaussian broadening @xmath147 ( black line ) .
( b ) enlarged view close to the singularity at @xmath148 on a double logarithmic scale : exact results ( black long - dashed line ) , ddmrg data ( red short - dashed line ) , and the deconvolved spectra @xmath109 for @xmath149 ( black solid line ) and @xmath150 ( blue dash - dot line).,title="fig:",scaledwidth=48.0% ]    figure  [ fig : tb](a ) shows the dos of a tight - binding chain calculated with ddmrg and the result of our deconvolution procedure .
the ddmrg spectrum has been calculated in the middle of an open chain with @xmath151 sites using a broadening @xmath146 , which is just broad enough to hide its discreteness .
we see that the square - root divergences at @xmath152 have been smoothed into two broad peaks and that there is substantial spectral weight at energies @xmath153 .
the deconvolved dos has been determined from these same ddmrg data using a minimal extremum distance @xmath154 and a final gaussian broadening with @xmath155 .
we see now that the singularities at @xmath152 are clearly visible as sharp peaks and that there is not any spectral weight at @xmath153 .
overall the deconvolved dos is in excellent agreement with the exact spectrum in the thermodynamic limit  ( [ eq : dos - tb ] ) . in particular
, we do not observe any unphysical artefact such as negative spectral weight .
however , in fig .  [ fig : tb](a ) we observe two shoulders in the deconvolved dos at energies @xmath156 , which are not present in the exact solution  ( [ eq : dos - tb ] ) .
an enlarged view close to the singularity at @xmath157 is shown in fig .  [
fig : tb](b ) on a double logarithmic scale .
we see that the ddmrg data agree with the exact result only at some distance from the singularity . in this figure
we also show deconvolved spectra @xmath109 for two different values of the normalized cost function @xmath158 . clearly , they reproduce the square - root divergence at @xmath157 much better than the original ddmrg data .
the overall divergent behavior is visible on a broader energy scale for the smaller value of @xmath158 but we see that the reduction of the cost function is also accompanied by stronger oscillations around the exact result .
these oscillations correspond to the shoulder seen in fig .
[ fig : tb](a ) .
the occurrence of artificial shoulder - like structures is the main drawback of our deconvolution procedure .
any deconvolution method magnifies the noise ( numerical errors ) which is present in the original data .
this the main issue that existing methods try to solve in different ways .
@xcite we have systematically tested our deconvolution procedure using exact results for non - interacting systems and purposely adding random numerical errors .
we have found that by preventing the formation of local maxima in the deconvolved spectrum our procedure allows us to control the noise magnification only partially .
unfortunately , it can not handle extrema ( @xmath159 oscillations ) in the spectrum derivative .
thus the magnified noise shows up as shoulder - like structures ( but not as local maxima , discontinuities , or sharp angles ) on an energy scale @xmath0 and gives a rough appearance to some deconvoled spectra presented here . in principle , we should be able to correct this deficiency with a higher - order interpolation procedure in the smoothening step or with a smoothening of the derivative of @xmath46 ( i.e. , the finite differences between the parameters @xmath103 ) .
however , we have not yet succeeded in developing a practical algorithm based on these ideas .
( color online ) ( a ) dos of the one - dimensional half - filled hubbard model at @xmath160 calculated with ddmrg in a 128-site chain using a broadening @xmath146 ( red dashed line ) and the result of our deconvolution procedure ( black line ) with a minimal extremum distance @xmath161 and a gaussian broadening @xmath162 .
the vertical dashed lines show the exact position of the dos threshold calculated from the bethe ansatz solution .
( b ) enlarged view of the same data around the dos threshold for @xmath163 .
additionally , the result of our procedure with @xmath164 ( blue dot - dash line ) and of a deconvolution with linear regularization  @xcite ( green dots ) are also shown.,title="fig:",scaledwidth=48.0% ]   ( color online ) ( a ) dos of the one - dimensional half - filled hubbard model at @xmath160 calculated with ddmrg in a 128-site chain using a broadening @xmath146 ( red dashed line ) and the result of our deconvolution procedure ( black line ) with a minimal extremum distance @xmath161 and a gaussian broadening @xmath162 . the vertical dashed lines show the exact position of the dos threshold calculated from the bethe ansatz solution .
( b ) enlarged view of the same data around the dos threshold for @xmath163 .
additionally , the result of our procedure with @xmath164 ( blue dot - dash line ) and of a deconvolution with linear regularization  @xcite ( green dots ) are also shown.,title="fig:",scaledwidth=48.0% ]    figure  [ fig : onedim4](a ) shows the single - particle dos of the half - filled one - dimensional hubbard model at @xmath160 calculated with ddmrg and the result of our deconvolution procedure .
the ddmrg spectrum has been computed  @xcite in the middle of an open chain with @xmath151 sites and a broadening @xmath146 .
the deconvolved spectrum has been obtained from these ddmrg data using a minimal extremum distance @xmath165 and a final gaussian broadening @xmath166 .
the effects of the broadening are clearly visible in the ddmrg data .
for instance , although one can recognize the opening of the mott - hubbard gap @xmath128 , spectral weight is clearly visible inside the gap .
a point - wise analysis  @xcite of the scaling for @xmath47 with @xmath167 is required to confirm that the spectral weight jumps from @xmath168 to a finite value at the onset @xmath169 and that the gap width agrees with the exact results @xmath170 calculated with the bethe ansatz .
@xcite however , the behavior for @xmath171 remains uncertain because of the relatively large broadening used in the ddmrg calculation . on the contrary
the deconvolved dos clearly shows a gap with the step - like onset  ( [ eq : dos - ft ] ) predicted by field theory  @xcite at the position @xmath172 given by the bethe ansatz solution .
we obtain similarly unambiguous results for @xmath127 up to @xmath173 ( see below ) .
therefore , our numerical investigation confirms the field - theoretical prediction  @xcite for the onset of the dos in one - dimensional mott insulators .
the superiority of the deconvolved spectrum over the original ddmrg data is even more obvious in fig .
[ fig : onedim4](b ) which shows an enlarged view of the dos around @xmath172 . in this figure
we also show the result of a deconvolution of the ddmrg data with a standard linear regularization method  @xcite .
( note that this method yields negative spectral weights for some energies @xmath174 but we show the positive parts only . )
we see that the result of the deconvolution procedure proposed in this work is much superior to that of the standard one , which is too blurred to allow us to determine the true form of the dos at the onset @xmath172 .
in addition , fig .
[ fig : onedim4](b ) shows the result of our deconvolution procedure for a minimal extremum distance @xmath175 which is deliberately too small . in that case ,
artificial oscillations on energy scales @xmath176 are clearly visible in the deconvolved spectrum .
( color online ) upper hubbard band in the dos of the one - dimensional half - filled hubbard model at @xmath177 calculated with ddmrg in a 64-site chain using a broadening @xmath178 ( red dashed line ) and the result of our deconvolution procedure ( black line ) with a gaussian broadening @xmath179 .
, scaledwidth=48.0% ]    in the strong - coupling limit @xmath180 our results are less conclusive .
for instance , we show the ddmrg and deconvolved upper hubbard and for a very strong coupling @xmath177 in fig .
[ fig : onedim40 ] .
the ddmrg spectrum has been calculated in the middle of an open chain with @xmath181 sites using a broadening @xmath178 .
the deconvolved spectrum has been obtained from these ddmrg data using a minimal extremum distance @xmath182 and a final a gaussian broadening @xmath179 .
two broad peaks are clearly visible at energies @xmath183 as predicted for the strong - coupling limit .
however , the widths and heights of these peaks after deconvolution are not compatible with the square - root divergences  ( [ eq : dos - scpt ] ) predicted by the low - order strong - coupling expansion .
actually , higher - order corrections indicate  @xcite that some spectral weight is present below the peak at @xmath184 on a scale set by the effective spin exchange coupling @xmath185 .
this width is compatible with our deconvolved spectra for @xmath186 .
thus our numerical results agree at least qualitatively with strong - coupling perturbation theory and suggest that the square - root divergences in the dos  ( [ eq : dos - scpt ] ) is an artifact of a truncated strong - coupling expansion .
nevertheless , for strong coupling such as @xmath177 we are not able to determine the shape of the dos at the onset @xmath172 . in particular
, it is not clear if the field - theoretical prediction  ( [ eq : dos - ft ] ) is still valid . indeed ,
if the distance @xmath187 between onset at @xmath169 and peak at @xmath188 becomes smaller than the minimal extremum distance @xmath0 , we can no longer distinguish both structures in the deconvolved spectrum . in practice ,
the distance @xmath0 must be comparable to the broadening @xmath36 of the ddmrg data to obtain piecewise smooth spectra .
therefore , although we can reduce the broadening of isolated spectral structures by several orders of magnitude , we can not resolve distinct spectral features on a scale lower than the original broadening of the ddmrg data .
this is a limitation of our deconvolution method that one has to keep in mind .
( color online ) deconvolved dos of the half - filled hubbard model for energies @xmath189 for @xmath6 ( dashed line ) , and from left to right @xmath190 ( red line ) , @xmath191 ( blue line ) and @xmath192 ( green line ) .
the result for @xmath177 is shifted to the left by @xmath193 .
, scaledwidth=48.0% ]    finally , fig .
[ fig : onedim ] recapitulates the evolution of the dos as a function of the interaction strength @xmath5 . as the spectrum is symmetric @xmath194 ,
we show only the deconvolved spectra for positive energies ( i.e. , the upper hubbard band ) . at @xmath6
the spectrum consists in a single band with two clearly visible square - root singularities at the band edges @xmath195 . for weak coupling @xmath127
the band splits into two symmetric hubbard bands separated by a gap @xmath128 , which agrees perfectly with the charge gap calculated from the bethe ansatz solution . at the dos onsets
@xmath196 the spectrum exhibits the step - like behavior  ( [ eq : dos - ft ] ) predicted by field theory .
@xcite there is an apparent plateau between the onset and a strong first peak , which evolves from the square - root singularities at @xmath195 for @xmath6 .
additionally , we observe substantial spectral weight and a small second peak at higher excitation energy .
however , for weak enough @xmath127 most of the spectral weight lies between the spectrum onset and the first peak .    as @xmath127 increases ( compare the spectra for @xmath160 and @xmath1 in fig .
[ fig : onedim ] ) , the spectrum and all its features shift to higher excitation energy and the spectral weight becomes more concentrated between the visible peaks . in addition , we note that the separation between onset energy @xmath142 and the strong first peak becomes systematically smaller until it is no longer resolvable with our method , the peak separation increases monotonically from about @xmath197 for @xmath198 to approximately @xmath199 for @xmath200 , and the strength of both peaks become more equal .    comparing the dos with the momentum - resolved spectral function and the bethe ansatz dispersion ( see figs . 4 and 5 in ref .  )
we note that the strong first peak corresponds to the edge of the spinon branch at momentum @xmath201 , the weak second peak corresponds to the edge of the holon branch at @xmath202 , and the upper edge of the dos spectrum coincide with the edge of the single spinon - holon continuum .
finally , we do not observe any spectral weight outside the first lower and upper hubbard bands and these two bands account for the full spectral weight  ( [ eq : sum ] ) .
thus we conclude that higher - energy hubbard bands do not carry any spectral weight in the bulk single - particle dos .
we have presented a blind deconvolution procedure which allows us to obtain piecewise smooth spectral functions for infinite - size systems from the ddmrg spectra of finite systems .
it involves a trade - off between the agreement of the deconvolved spectrum to the original ddmrg data and the piecewise smoothness and positivity of spectral functions . in practice
, the method reduces to a least - square optimization under non - linear constraints which enforce the positivity and piecewise smoothness .
we have tested this deconvolution method on many spectra which are known exactly in the thermodynamic limit , such as the single - particle density of states and the optical conductivity of correlated one - dimensional insulators .
@xcite we have found that our method works well for several kinds of singularities ( e.g.  power - law band edges , steps , excitonic peaks ) in piecewise smooth spectra . in particular , it allows us to reduce the broadening by orders of magnitude and even to substitute the lorentzian broadening by a gaussian one .
its main drawback is the frequent appearance of artificial shoulder - like structures on energy scales @xmath203 .
we have demonstrated the deconvolution procedure on the single - particle dos in the one - dimensional hubbard model at half filling .
our results show that the dos has a step - like shape but no square - root singularity at the spectrum onset in agreement with a field - theoretical prediction for one - dimensional paramagnetic mott insulators .
@xcite in addition , the deconvolution procedure has allowed us to detail the evolution of the dos from the non - interacting limit @xmath6 to the strong - coupling limit @xmath180 .
we thank karlo penc and fabian eler for helpful discussions .
the gotoblas library developed by kazushige goto was used to perform the ddmrg calculations .
some of these calculations were carried out on the rrzn cluster system of the leibniz universitt hannover .
e. jeckelmann , f. gebhard , and f.h.l .
essler , phys .
lett . * 85 * , 3910 ( 2000 ) .
e. jeckelmann , phys .
b * 66 * , 045114 ( 2002 ) . t.d .
khner and s.r .
white , phys . rev .
b * 60 * , 335 ( 1999 ) .
pati , s. ramasesha , z. shuai , and j.l .
brdas , phys .
b * 59 * , 14827 ( 1999 )
. e. jeckelmann , prog .
. suppl . * 176 * , 143 ( 2008 ) .
e. jeckelmann and h. benthien , in _ computational many particle physics _
( lecture notes in physics * 739 * ) , edited by h. fehske , r. schneider , and a. weie ( springer - verlag , berlin , heidelberg , 2008 ) , p. 621 .
white and i. affleck , phys .
b * 77 * , 134437 ( 2008 ) .
a. weichselbaum , f. verstraete , u. schollwck , j.i .
cirac , and j. von delft , phys .
b * 80 * , 165117 ( 2009 ) .
a. holzner , a. weichselbaum , i.p .
mcculloch , u.  schollwck , and j. von delft , phys .
b * 83 * , 195115 ( 2011 )
dargel , a. honecker , r. peters , r.m .
noack , and t.  pruschke , phys .
b * 83 * , 161104 ( 2011 ) .
e. jeckelmann , j. phys . : condens .
matter * 25 * , 014002 ( 2013 ) .
r. peters , phys .
b * 84 * , 075139 ( 2011 ) . w.h . press , s.a .
teukolsky , w.t .
vetterling , and b.p .
flannery , _ numerical recipes in c++  the art of scientific computing _
, second edition ( cambridge university press , cambridge , 2002 ) , chap .
oleary , _ scientific computing with case studies _
( siam , philadelphia , 2009 ) , chap . 6 and 14 .
f. gebhard , e. jeckelmann , s. mahlert , s. nishimoto , and r.m .
noack , eur .
j. b * 36 * , 491 ( 2003 ) .
s. nishimoto and e. jeckelmann , j. phys . : condens .
matter * 16 * , 613 ( 2004 ) .
s. nishimoto , f. gebhard , and e. jeckelmann , j. phys . : condens .
matter * 16 * , 7063 ( 2004 ) . c. raas , g.s .
uhrig , and f.b .
anders , phys .
b * 69 * , 041102(r ) ( 2004 ) . c. raas and g.s .
uhrig , eur .
j. b * 45 * , 293 ( 2005 ) .
e. jeckelmann and h. fehske , rivista del nuovo cimento * 30 * , 259 ( 2007 ) .
p. schmitteckert , journal of physics : conferences series * 220 * , 012022 ( 2010 ) . t. ulbricht and p. schmitteckert , euro .
* 89 * , 47001 ( 2010 ) .
p. campisi and k. egiazarian , _ blind image deconvolution  theory and applications _
( crc press , boca raton , 2007 ) . j. hubbard , proc .
london a * 276 * , 238 ( 1963 ) .
f.h.l . essler and a.m. tsvelik , phys .
b * 65 * , 115117 ( 2002 ) .
mott , _ metal - insulator transitions _ ( taylor and francis , london , 1990 ) .
f. gebhard , _ the mott metal - insulator transition _
( springer , berlin , 1997 ) .
y. kurosaki , y. shimizu , k. miyagawa , k. kanoda , and g. saito , phys .
lett * 95 * , 177001 ( 2005 ) .
kim , j.p .
hill , h. benthien , f.h.l .
essler , e. jeckelmann , h.s .
choi , t.w .
noh , n. motoyama , k.m .
kojima , s. uchida , d. casa , and t. gog , phys .
lett . * 92 * , 137402 ( 2004 ) .
m. azuma , z. hiroi , m. takano , k. ishida , and y. kitaoka , phys . rev
73 * , 3463 ( 1994 ) .
e. h. lieb and f. y. wu , phys .
20 * , 1445 ( 1968 ) .
essler , h. frahm , f. ghmann , a. klmper , and v. korepin , _ the one - dimensional hubbard model _
( cambridge university press , cambridge , 2005 ) .
a. parola and s. sorella , phys .
b * 45 * , 13  156 ( 1992 ) .
k. penc , k. hallberg , f. mila , and h.  shiba , phys .
b * 55 * , 15  475 ( 1997 ) . k. penc ( private communication ) .
essler , f. gebhard , and e. jeckelmann , phys .
b * 64 * , 125119 ( 2001 ) .