gauge theories constitute the building blocks of our present understanding of natural phenomena .
the standard model ( sm ) of high energy particle interactions is , in fact , entirely based on a semi - simple gauge group @xmath0 .
the sm accounts for roughly four percent of the known universe .
however , several puzzles remain still unexplained .
for example , why do we observe , at least , three generations of elementary particles . besides , the remaining @xmath1 of the universe , of which @xmath2 dark matter , and @xmath3 dark energy , remains largely unknown .
it is therefore natural to expect that new gauge theories , or extensions of the present ones play a fundamental role in explaining the unknown side of the universe .
the space of four dimensional gauge theories at our disposal , without theoretical prejudice , is very large and moreover a large fraction of it is still _ terra incognita _ due to our limited methods to tackle nonperturbative dynamics .
a fascinating possibility is that different gauge theories display the same physical dynamics .
this is possible since the physical quantities are determined by the gauge singlet operators of a generic gauge theory .
therefore it is perfectly legal to imagine to have two different gauge theories leading to the same physics , at least , in certain regimes .
quantum chromodynamics ( qcd ) itself is a famous example .
in fact , at low energy the spectrum can be represented by hadrons interacting via yukawa interactions and its electric dual description is in terms of the @xmath4 gauge theory featuring quarks and gluons .
it is therefore natural to imagine that generic asymptotically free gauge theories have magnetic duals .
in fact , in the late nineties , in a series of ground breaking papers seiberg @xcite provided strong support for the existence of a consistent picture of such a duality within a supersymmetric framework .
supersymmetry is , however , quite special and the existence of such a duality does not automatically imply the existence of nonsupersymmetric duals .
one of the most relevant results put forward by seiberg has been the identification of the boundary of the conformal window for supersymmetric qcd as a function of the number of flavors and colors .
recently several analytic predictions have been provided for the conformal window of nonsupersymmetric gauge theories using different approaches @xcite .
we initiated in @xcite the exploration of the possible existence of a qcd nonsupersymmetric gauge dual providing a consistent picture of the phase diagram as a function of number of colors and flavors .
arguably the existence of a possible dual of a generic nonsupersymmetric asymptotically free gauge theory able to reproduce its infrared dynamics must match the t hooft anomaly conditions @xcite .
we have exhibited several solutions of these conditions for qcd in @xcite .
an earlier exploration already appeared in @xcite . in @xcite theories with fermions transforming according to higher dimensional representations were analyzed .    in @xcite
we have put the idea of nonsupersymmetric gauge duality on a firmer ground by showing that for certain scalarless gauge theories with a spectrum similar to the one of qcd the gauge dual passes a large number of consistency checks .
the theory studied in @xcite is qcd with @xmath5 dirac flavors and one adjoint weyl fermion .
an important feature of this theory is that it possesses the same global symmetry of super qcd despite the fact that squarks are absent .
this means that there are four extra anomaly constraints not present in the case of ordinary qcd , moreover we have shown that the potential dual can be constructed for any number of colors greater than two .
the magnetic dual is a new gauge theory featuring magnetic quarks and a weyl adjoint fermion , new gauge singlet fermions which can be identified as states composite of the electric variables , as well as scalar states needed to mediate the interactions between the magnetic quarks and the gauge singlet fermions .
the new scalars allow for a consistent flavor decoupling which was an important consistency check in the case of supersymmetry .
in @xcite it was also shown that the candidate dual allows to bound the anomalous dimension of the dirac fermion mass operator to be less than one in the conformal window , and there we also estimated the critical number of flavors below which large distance conformality is lost in the electric variables .    in these previous studies we implicitly assumed that there exist nontrivial infrared fixed points in the magnetic description matching the electric ones .
however , given that we do not start in the ultraviolet with a supersymmetric spectrum nor we impose supersymmetric coupling relations it is legitimate to ask whether these fixed point exist at all in the magnetic description , and if they do what are their universal properties .
the goal of this work is exactly to establish the existence of these fixed points assuming the most general coupling structure .
we will investigate the regime where the computations are trustable , i.e. when the magnetic theory features fixed points accessible in perturbation theory .
we discover , for the first time , a large number of calculable fixed points and classify them by analyzing their basin of attraction .
an extremely interesting result , according to us , is that among these very many fixed points we discover supersymmetric fixed points to emerge in the infrared along nonsupersymmetric flows .
our results have an impact on the ( non)supersymmetric gauge - gauge duality and even recent ideas on why the sm can be seen as the magnetic dual description of an electric strongly coupled theory featuring no scalars @xcite .
the consequences of this interpretation of the sm leads to a novel way to tackle the hierarchy problem and even to shed light on the mystery of the observed number of matter generations @xcite .
the phenomenological interest of our studies relies also on the fact that theories similar to the ones investigated here and featuring infrared fixed points have been used to construct sensible extensions of the standard model of particle interactions of technicolor type passing precision data and known as minimal walking technicolor models @xcite .
duality has also been used to infer relevant insights on the constraints on these models from precision electroweak tests @xcite .
in @xcite we investigated the possible existence of a magnetic dual for the electric theory constituted by a scalarless @xmath6 gauge theory with @xmath5 dirac fermions and @xmath7 larger than two , as in qcd , but with an extra weyl fermion transforming according to the adjoint representation of the gauge group .
the quantum global symmetry of the electric theory is therefore : @xmath8 at the classical level there is one more @xmath9 symmetry destroyed by quantum corrections due to the adler - bell - jackiw anomaly . of the three independent @xmath10 symmetries only two survive , a vector like @xmath11 and an axial - like anomaly free ( af )
one indicated with @xmath12 .
the spectrum of the theory and the global transformations are summarized in table [ electric ] .
@xmath13 & su_l(n_f ) & su_r(n_f ) & u_v(1)&u_{af}(1 ) \\ \hline \hline \lambda & { \rm adj } & 1 & 1 & ~~0 & ~~1 \\ q & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}}&{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}}}&1&~~1 & -\frac{n}{n_f } \\ \widetilde{q } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}}}&1 &   \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & -1 & -\frac{n}{n_f }    \\ g_{\mu}&{\rm adj }    & 1&1   & ~~0   & ~~0\\
\hline \end{array}\ ] ]    we then postulated the existence of a magnetic dual featuring the minimal spectrum of composite states and gauge group structure needed to :    * classify all of the t hooft anomaly conditions and show how to match them , for the first time , for _ any _ number of colors and flavors ; * allow for consistent flavor decoupling both in the electric and in the magnetic theory ; * ensure _ duality involution_. meaning that if we dualize once more the magnetic theory one recovers the gauge structure of the electric theory .    the proposed nonsupersymmetric magnetic gauge theory @xcite is summarized in table [ qcdadual ] .
@xmath14 & su_l(n_f ) & su_r(n_f ) & u_v(1 ) & u_{af}(1 ) \\ \hline   \hline   \lambda_m & { \rm adj } & 1 & 1 & 0 & 1 \\   q & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}}&\overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}}}&1&~~\frac{n_f - x}{x } & - \frac{x}{n_f }   \\ \widetilde{q } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}}}&1 &   { { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & -\frac{n_f - x}{x } & - \frac{x}{n_f }      \\   m   & 1 & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & 0 & -\frac{n_f-2x}{n_f } \\   \hline       { \widetilde}{\phi } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & 1 & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } & -\frac{n_f - x}{x } & \frac{n_f - x}{n_f}\\    \phi & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & 1 & \frac{n_f - x}{x } & \frac{n_f - x}{n_f}\\    g_\mu & \text{adj } & 1 & 1 & 0 & 0 \\     \hline \end{array}\ ] ]
we have also shown that one has to have @xmath15 .
differently from the dual of super qcd we did not impose a supersymmetric spectrum or supersymmetric coupling relations .
the spectrum above is nonsupersymmetric as it is clear from the fact that there is no complex scalar partner of @xmath16 .
moreover we showed that we could build the gauge singlet states using electric variables which do not contain squarks .
subsequently in @xcite it was argued that one could also add the complex scalar @xmath17 ( appearing now in table [ qcdadualm ] ) .
this is possible since it does not affect the anomaly conditions and can be built naturally out of the electric fermionic variable as follows @xcite : @xmath18 it is crucial to be able to construct all these states directly from the electric fermionic variables .
this demonstrates that supersymmetry is not a fundamental ingredient in order to construct these states .
@xmath17 in @xcite plays the phenomenologically relevant role of the sm - like higgs , elementary in terms of the magnetic variables .
although when adding the new complex scalar field the spectrum of the dual theory looks supersymmetric , the full theory is not since the couplings , in the ultraviolet , are not taken to respect supersymmetric relations . assuming the duality to exist we were able to make a number of relevant predictions for the nonperturbative gauge dynamics of the strongly coupled qcd gauge theory with one fermion in the adjoint representation .
the most relevant ones being :    * the anomalous dimension of the mass of the electric fermions at the lower boundary of the conformal window can not exceed unity .
* estimate the size of the conformal window , i.e. the critical number of flavors below which the theory looses large distance conformality .
* when the magnetic theory is used as possible dual of a minimal extension of the sm one can argue that this extension can be re - constructed in terms of a strongly coupled electric theory featuring only fermionic matter , and furthermore argue that the mathematical consistency of the electric dual requires , at least , three generations of ordinary matter @xcite .
what was still missing is the actual existence of the fixed points in the magnetic theory .
this is what this paper aims to accomplish using perturbation theory .
this regime is achieved by choosing the number of flavors @xmath5 , number of colors @xmath19 , and yukawa couplings where perturbation theory holds .
@xmath14 & su_l(n_f ) & su_r(n_f ) & u_v(1 ) & u_{af}(1 ) \\ \hline   \hline   \lambda_m & { \rm adj } & 1 & 1 & 0 & 1 \\   q & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}}&\overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}}}&1&~~\frac{n_f - x}{x } & - \frac{x}{n_f }   \\ \widetilde{q } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}}}&1 &   { { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & -\frac{n_f - x}{x } & - \frac{x}{n_f }      \\   m   & 1 & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & 0 & -\frac{n_f-2x}{n_f } \\   \hline      h & 1 & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & 0 & \frac{2x}{n_f}\\    { \widetilde}{\phi } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & 1 & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } & -\frac{n_f - x}{x } & \frac{n_f - x}{n_f}\\    \phi & { \raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } & \overline{{\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } } & 1 & \frac{n_f - x}{x } & \frac{n_f - x}{n_f}\\    g_\mu & \text{adj } & 1 & 1 & 0 & 0 \\     \hline \end{array}\ ] ]
we start with listing the yukawa operators:@xmath20 @xmath21-interactions do not affect our results to the perturbative order we are considering and therefore we will not include them here .
this yukawa - sector contributes to the two - loop beta function of the gauge coupling in the following manner : @xmath22 \beta_0 & = \frac{11}{3 } c_2(g ) - \frac{2}{3 } \sum_r t(r ) n_f(r)- \frac{1}{6}\sum_s t(s ) n_f(s)\\ \beta_1 & = \frac{34}{3 } c_2(g)^2- \sum_r \left [ \frac{10}{3 } c_2(g ) + 2 c_2 ( r ) \right ] t(r ) n_f(r ) -\sum_s \left [ \frac{1}{3 } c_2(g ) + 2 c_2 ( r ) \right ] t(s ) n_f(s)\\ \beta_y & =   \frac{1}{d(g ) } \sum_r\text{tr}\left [ c_2(r ) y^j { y_j}^\dagger\right]\label{betay}\end{aligned}\ ] ] where @xmath23 denotes the representation of fermions and @xmath24 denotes the representation of the _ real _ scalars .
@xmath25 is the trace normalization of the group generators , @xmath26 is the quadratic casimir of these and @xmath27 is the dimension of the gauge group .
we refer to appendix [ app : rg ] for a careful derivation of the yukawa contribution @xmath28 to the running of the gauge coupling .
the result is : @xmath29 ,   \end{aligned}\ ] ] with @xmath30 and where we used the notation @xmath31    we now consider the running of the yukawa couplings .
the one loop beta function for the yukawa couplings is given by @xcite : @xmath32 & + 2 y^k y^\dagger_j y^k +   \frac{1}{2}y^k \text{tr}\left [ y^\dagger_k y^j+y^\dagger_j y^k\right ] - 3 g^2 \ { c_2(r ) , y^j \}.\end{aligned}\ ] ] here @xmath33 is the yukawa coupling matrix defined by the particular interaction : @xmath34 where roman indices contract over the scalar gauge - flavor overall index and the greek indices @xmath35 are again gauge - flavor indices but reserved for the weyl fermions @xmath36 and @xmath37 .
@xmath38 is the group invariant : @xmath39 for each scalar contraction we multiply , row by column , the yukawa matrices over the fermion indices .
we again report the derivation of the beta function for each yukawa coupling in appendix [ app : rg ] .
to the second order in the couplings the set of beta function equations reads : @xmath40\\    \label{rgsystem2 }   \beta(\alpha_{\lambda } ) = 2 \alpha_{\lambda } & \left [ \frac{3 ( x^2 - 1)}{4x } \alpha_\lambda+ ( \alpha_\lambda + \alpha_{{\widetilde}{\lambda}})\frac { n_f } { 4 } + n_f \left ( \frac{\alpha_{{\widetilde}{m}}+\alpha_h}{2 } + \alpha_{m } \right ) - 3\alpha_g \frac{3x^2 - 1}{2x } \right ]   \nonumber\\ &    -4 n_f\frac { y_m y_{\widetilde m } y_\lambda y_{\widetilde \lambda}}{(4\pi)^2 }    \\    \label{rgsystem3 }   \beta(\alpha_{m } ) = 2 \alpha_{m } & \left [ \frac{3 n_f + x}{2 } \alpha_m+ \frac { x } { 2}\alpha_{{\widetilde}{m } } + \frac{n_f}{2}\alpha_h + \frac{x^2 - 1}{2x } \left ( \frac{\alpha_{{\widetilde}{\lambda}}}{2 } + \alpha_\lambda \right ) - 3\alpha_g   \frac{x^2 - 1}{2x}\right ]   \nonumber\\ &    -4\frac{x^2 - 1}{2x } \frac { y_m y_{\widetilde m } y_\lambda y_{\widetilde \lambda}}{(4\pi)^2 }    \\
\beta(\alpha_h ) = 2 \alpha_h & \left [ \frac{\alpha_m + \alpha_{{\widetilde}{m } } + 2\alpha_h}{2}n_f   + \left ( \frac{\alpha_{\lambda } + \alpha_{{\widetilde}{\lambda}}}{2}-6\alpha_g\right)c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) + x \alpha_h \right]\\    \beta(\alpha_{{\widetilde}{\lambda } } ) = \beta_{\alpha_{\lambda } } & \left ( y_{\lambda }   \leftrightarrow   y_{\widetilde \lambda }   ,   y_{\widetilde m }   \leftrightarrow   y_{m }   \right ) , \\   \beta(\alpha_{\widetilde{m } } ) = \beta_{\alpha_{m}}&\left (   { y}_{m }   \leftrightarrow   y_{\widetilde m } ,   y_{\widetilde \lambda }   \leftrightarrow   { y}_{\lambda }   \right ) ,   \label{rgsystem4}\end{aligned}\ ] ] where : @xmath41 we impose @xmath42 invariance and therefore set all the yukawa phases to zero .
we discuss the full set of equations , including the phases , in the appendix [ app : rg ] .    in the supersymmetric limit @xmath43 , @xmath44 and the beta function system simplifies to :
@xmath45\\   \beta(\alpha_{\lambda } )    & = - 2 \alpha_{\lambda } \left [ \beta_0 \alpha_g    \right]\\   \beta(\alpha_{m } ) & = 2 \alpha_{m } \left [ ( 2n_f + x ) \alpha_m - 4\alpha_gc_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) \right]\\    \beta(\alpha_h ) & = 2 \alpha_h \left [ ( 2n_f + x ) \alpha_h   - 4\alpha_gc_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}})\right],\label{susylimit4}\end{aligned}\ ] ] where @xmath46 is the well - known two - loop term of pure super qcd , i.e. _ without _ the gauge - singlet chiral superfield transforming as bilinear with respect to the non - abelian global symmetries and , if present , leading to the introduction of a superpotential term .
the theory with this chiral superfield included is seiberg s magnetic dual .
this result shows , as expected , that supersymmetry stays unbroken along the renormalization flow .
also note how the running of the gaugino yukawa coupling has collapsed to the one - loop result for the running of the gauge coupling .
these results are in agreement with the known result from supersymmetry .
the all - orders supersymmetric beta functions for seiberg s magnetic dual are : @xmath47}{1 - 2x\alpha_g}\nonumber\\ \beta^s(\alpha_m )   & =   \alpha_m    \left [ \gamma_m + 2\gamma_0 \right],\end{aligned}\ ] ] where the one - loop expressions for the anomalous dimensions of the chiral superfields read : @xmath48 by which it is readily seen that the one - loop expansion of @xmath49 and @xmath50 are in agreement with the expressions in equations - .
just below the critical number of flavors where the theory looses asymptotic freedom ( @xmath51 ) a banks - zaks perturbatively stable infrared fixed point ( irfp ) emerges once the yukawa interactions are set to zero . the well known expression for the value of the gauge coupling at this fixed point is : @xmath52 in this phase the theory shows large distance conformality .
we now investigate whether the perturbatively stable irfp persists when we turn on the yukawa interactions . from eq .
, the perturbative fixed point value for the gauge coupling reads : @xmath53 perturbative consistency is ensured by the smallness of the yukawa couplings at the fixed point .
there is , clearly , a symmetry between the tilded and untilded yukawa coupling constants at the beta functions level .
this symmetry persists at the fixed point effectively reducing the space of solutions .
the fixed points associated to zeros of @xmath54 are : @xmath55 a perturbative nontrivial infrared fixed point for all the couplings is achieved when the number of flavors are such that asymptotic freedom is almost lost . to express this point mathematically we introduce the small parameter @xmath56 as follows : @xmath57 we report in table  [ solutions0 ] and table  [ solutions1 ] the explicit solutions for the fixed points as a function of @xmath56 .
there are 18 physical solutions for the system - , where @xmath58 of these correspond to the class @xmath59 and the remaining @xmath60 to @xmath61 .
the solutions obtained exchanging tilded couplings with untilted ones @xmath62 are indicated in the tables  [ solutions0 ] and [ solutions1 ] by an asterisk .
@xmath63 \hline 2^ { ( * ) } &    \alpha_g^ * = \frac{2 ( 3 x - n_f ) ( x ( n_f + 3 x)-3)x } { 2 n_f^2 x ( 7 x^2 - 3 ) + 3 n_f ( 5 - 14 x^2 + x^4)-36 x^3 ( x^2 - 1 ) }   & = \frac{2 x ( 2 x^2 - 1 ) \epsilon}{5 - 16 x^2 + 11 x^4 } \\ & \alpha_{\lambda}^ * = \frac{12 ( 3 x - n_f ) ( 3x^2 - 1)x } { 2 n_f^2 x ( 7 x^2 - 3 ) + 3 n_f ( 5 - 14 x^2 + x^4)-36 x^3 ( x^2 - 1 ) }   & = \frac{4 x ( 3 x^2 - 1 ) \epsilon}{5 - 16 x^2 + 11 x^4}\\ &   \alpha_{{\widetilde}{\lambda } , m , { \widetilde}{m}}^ * = 0 & \\[2 mm ] \hline 3^ { ( * ) } &   \alpha_g^ * = \frac{2 ( 3 x - n_f ) ( 3n_f + x ) x } { 3 n_f^2 ( 13 x^2 - 5)-12 x^4 - 2 n_f x ( 3 + 11 x^2 ) }   & = \frac{20x \epsilon}{91 x^2 - 51 } \\ & \alpha_{m}^ * = \frac{6 ( 3 x - n_f ) ( x^2 - 1)x } { 3 n_f^2 ( 13 x^2 - 5)-12 x^4 - 2 n_f x ( 3 + 11
x^2 ) }   & = \frac{6 (   x^2 - 1 ) \epsilon}{x(91x^2 - 51)}\\ &   \alpha_{\lambda,{\widetilde}{\lambda } , { \widetilde}{m}}^ * = 0 & \\[2 mm ] \hline 4 & \alpha_g^ * = \frac { ( 3 x - n_f ) ( 3n_f + 2x ) x}{2 (   n_f^2 ( 9 x^2 - 3 ) -6 x^4- n_f x ( 3 + 2 x^2 ) ) }   & = \frac{11x \epsilon}{2(23 x^2 - 12 ) } \\ & \alpha_{m}^ * = \alpha_{{\widetilde}{m}}^*= \frac{3 ( 3 x - n_f ) ( x^2 - 1 ) } { 2 (   n_f^2 ( 9 x^2 - 3 ) -6 x^4- n_f x ( 3 + 2 x^2 ) ) }   & = \frac{3 (   x^2 - 1 ) \epsilon}{2x(23x^2 - 12)}\\ &   \alpha_{\lambda,{\widetilde}{\lambda}}^ * = 0 & \\[2 mm ] \hline 5 & \alpha_g^ * = \frac { ( 3 x - n_f ) ( 3x^2 + 2n_f x -3 ) x}{2 (   9 x^3 - 9 x^5 + n_f^2 x ( 7 x^2 - 3 ) + n_f ( 3 - 6 x^2 - 9 x^4 ) ) }   & = \frac{x \epsilon}{2(x^2 - 1 ) } \\ & \alpha_{\lambda}^ * = \alpha_{{\widetilde}{\lambda}}^*= \frac{3 ( 3 x - n_f ) ( 3x^2 - 1)x } { 2 (   9 x^3 - 9 x^5
+ n_f^2 x ( 7 x^2 - 3 ) + n_f ( 3 - 6 x^2 - 9 x^4 ) ) }   & = \frac{x\epsilon}{x^2 - 1 } = 2\alpha_g^ * + \mathcal{o}(\epsilon^2)\\ &   \alpha_{m,{\widetilde}{m}}^ * = 0 & \\[2 mm ] \hline 6^ { ( * ) } &    \alpha_g^ * = \frac{2 ( 3 x - n_f ) ( 3 ( n_f^2 - 1 ) x -8n_f+ 9 n_f x^2 + 3 x^3)x } { 3 n_f^3 x ( 13 x^2 - 5)-36 x^4 ( x^2 - 1 ) + 3 n_f x ( 5 + 18 x^2 - 31 x^4 ) + 2 n_f^2 ( 18 - 59 x^2 + 9 x^4 ) ) }   & = \frac{2 x ( 19 x^2 - 9 ) \epsilon}{41 - 141 x^2 + 100 x^4}\\ & \alpha_{\lambda}^ * = \frac{12 ( 3 x - n_f ) ( 3x^3 - 2n_f - x+8n_fx^2 ) } { 3 n_f^3 x ( 13 x^2 - 5)-36 x^4 ( x^2 - 1 ) + 3 n_f x ( 5 + 18 x^2 - 31 x^4 ) + 2 n_f^2 ( 18 - 59 x^2 + 9 x^4 ) ) }   & = \frac{4x ( 27 x^2 - 7 ) ) \epsilon}{41 - 141 x^2 + 100 x^4}\\ & \alpha_{{\widetilde}{m}}^ * = \frac{6 ( 3 x - n_f ) ( n_f x-2)(x^2 - 1 ) } { 3 n_f^3 x ( 13 x^2 - 5)-36 x^4 ( x^2 - 1 ) + 3 n_f x ( 5 + 18 x^2 - 31 x^4 ) + 2 n_f^2 ( 18 - 59 x^2 + 9 x^4 ) ) }   & = \frac { ( 6 x^2 - 4 ) \epsilon}{100x^3 - 41x}\\ &   \alpha_{{\widetilde}{\lambda},m}^ * = 0 & \\[2 mm ] \hline 7 & \alpha_g^ * = \frac { ( 3 x - n_f ) ( 3 ( n_f^2 - 1 ) x   -4 n_f   + 6 n_f x^2 + 3 x^3 ) } { 2(9 x^3 - 9 x^5 + n_f^3 ( 9 x^2 - 3 ) + n_f^2 ( x - 15 x^3 ) + n_f(3 + 6 x^2 - 21 x^4 ) ) }   & = \frac{x ( 16 x^2 - 5 ) \epsilon}{2(4x^4 - 5x^2 + 1)}\\ & \alpha_{\lambda}^ * = \alpha_{{\widetilde}{\lambda}}^*= \frac{3 ( 3 x - n_f ) ( 5 n_f x^2 -2 n_f - x + 3 x^3 ) } { 9 x^3 - 9 x^5 + n_f^3 ( 9 x^2 - 3 ) + n_f^2 ( x - 15 x^3 ) + n_f(3 + 6 x^2 - 21 x^4 ) }   & = \frac{x ( 18 x^2 - 7 ) ) \epsilon}{4x^4 -5x^2 + 1}\\ & \alpha_m^*=\alpha_{{\widetilde}{m}}^ * = \frac{3 ( 3 x - n_f ) ( n_f x+3x^2 - 2)(x^2 - 1 ) } { 2x(9 x^3 - 9 x^5 + n_f^3 ( 9 x^2 - 3 ) + n_f^2 ( x - 15 x^3 ) + n_f(3 + 6 x^2 - 21 x^4 ) ) }   & = \frac { ( 3 x^2 - 1 ) \epsilon}{x(4x^2 - 1)}\\[2 mm ] \hline \end{array}\ ] ] [ solutions0 ]    @xmath64 \hline 2^ { ( * ) } &    \alpha_g^ * = \frac{(3 x - n_f ) ( 6 x ( x^2 - 1 ) + n_f ( 7 x^2 + 2 x n_f-5))}{2 x n_f^2 ( 5 x^2 + 4 x n_f-13 ) -36 x^3 ( x^2 - 1 ) -3n_f ( 4 x^2 + 9 x^4 - 5 ) }   & = \frac{x ( 15 x^2 - 7 ) \epsilon}{5 - 26 x^2 + 21 x^4 } \\ & \alpha_h^ * =   \frac{3 ( 3 x - n_f ) ( x^2 - 1 ) ( 3 x^2 + 2 x n_f-5)}{x(2 x n_f^2 ( 5 x^2 + 4 x n_f-13 ) -36 x^3 ( x^2 - 1 ) -3n_f ( 4 x^2 + 9 x^4 - 5 ) ) }   & = \frac{(9x^2 - 5)\epsilon}{x(21x^2 - 5 ) } \\ & \alpha_{\lambda}^ * = \frac{12 x ( 3 x - n_f ) ( 3 x^2 + 2 x n_f-1 ) } { 2 x n_f^2 ( 5 x^2 + 4 x n_f-13 ) -36 x^3 ( x^2 - 1 ) -3n_f ( 4 x^2 + 9 x^4 - 5 ) }   & = \frac{4 x ( 9 x^2 - 1 ) \epsilon}{5- 26 x^2 + 21 x^4}\\ &   \alpha_{{\widetilde}{\lambda } , m , { \widetilde}{m}}^ * = 0 & \\[2 mm ] \hline 3^ { ( * ) } &   \alpha_g^ * = \frac{(3 x - n_f ) ( 2 x^2 + 8 x n_f + 5 n_f^2 ) } { 20n_f^3x - 12 x^4 + n_f^2(17x^2 - 15 ) - 2 n_f x ( 17x^2 + 3 ) }   & = \frac{71x \epsilon}{193x^2 - 51 } \\ & \alpha_h^ * =   \frac{3 ( 3 x - n_f ) ( 2 x + 5 n_f ) ( x^2 - 1 ) } { x(20n_f^3x - 12 x^4 + n_f^2(17x^2 - 15 ) - 2 n_f x ( 17x^2 + 3 ) ) }   & = \frac{51(x^2 - 1)\epsilon}{x(193x^2 - 51 ) } \\ & \alpha_{m}^ * = \frac{6 ( 3 x - n_f ) ( x^2 - 1 ) } { 20n_f^3x - 12 x^4 + n_f^2(17x^2 - 15 ) - 2 n_f x ( 17x^2 + 3 ) }   & = \frac{6 (   x^2 - 1 ) \epsilon}{x(193x^2 - 51)}\\ &   \alpha_{\lambda,{\widetilde}{\lambda } , { \widetilde}{m}}^ * = 0 & \\[2 mm ] \hline 4 & \alpha_g^ * = \frac { ( 3 x - n_f ) ( 2n_f + x ) ( n_f+2x)}{2 ( 4 n_f^3 x - 6 x^4 + n_f^2 ( 7 x^2 - 3 ) - n_f x ( 3 + 8 x^2 ) ) }   & = \frac{35x \epsilon}{2(47x^2 - 12 ) ) } \\ & \alpha_h^ * =   \frac{3 ( 3 x - n_f ) ( n_f+x ) ( x^2 - 1 ) } { x(4 n_f^3 x - 6 x^4 + n_f^2 ( 7 x^2 - 3 ) - n_f x ( 3 + 8 x^2 ) ) }   & = \frac{12(x^2 - 1)\epsilon}{x(47x^2 - 12 ) } \\ & \alpha_{m}^ * = \alpha_{{\widetilde}{m}}^*= \frac{3 ( 3 x - n_f ) ( x^2 - 1 ) } { 2 (   4 n_f^3 x - 6 x^4 + n_f^2 ( 7 x^2 - 3 ) - n_f x ( 3 + 8 x^2 ) ) }   & = \frac{3 (   x^2 - 1 ) \epsilon}{2x(47x^2 - 12)}\\ &   \alpha_{\lambda,{\widetilde}{\lambda}}^ * = 0 & \\[2 mm ] \hline 5 & \alpha_g^ * = \frac{(3 x - n_f ) ( 3 x ( x^2 - 1 ) + 2 n_f ( 2 x^2 + x n_f-1))}{2 ( 4 n_f x + 3 x^2 - 3 ) ( n_f^2 x - 3 x^3 - n_f ( x^2 + 1 ) ) }   & = \frac{x(11x^2 - 3 ) \epsilon}{2(5x^2 - 1)(x^2 - 1 ) } \\ & \alpha_h^ * =   \frac{3 ( 3 x - n_f ) ( n_fx-1 ) ( x^2 - 1 ) } { x ( 4 n_f x + 3 x^2 - 3 ) ( n_f^2 x - 3 x^3 - n_f ( x^2 + 1 ) }   & = \frac{(3x^2 - 1)\epsilon}{x(5x^2 - 1 ) } \\ & \alpha_{\lambda}^ * = \alpha_{{\widetilde}{\lambda}}^*= \frac{3x ( 3 x - n_f ) ( 3x^2 + 2xn_f-1 ) } { ( 4 n_f x + 3 x^2 - 3 ) ( n_f^2 x - 3 x^3 - n_f ( x^2 + 1 ) ) }   & = \frac{x(9x^2 - 1)\epsilon}{(5x^2 - 1)(x^2 - 1 ) } \\ &   \alpha_{m,{\widetilde}{m}}^ * = 0 & \\[2 mm ] \hline 6 & \alpha_g^ * = \frac { ( 3 x - n_f ) ( 2 n_f + x ) (   n_f^2 + 2 n_f x + 3 x^2 - 3 ) } { 2(4 n_f^4 x + 9 x^3 - 9 x^5 + n_f^3 ( x^2 - 3 ) - n_f^2 x ( 3 + 23 x^2 ) + 3 n_f ( 1 + 4 x^2 - 9 x^4 ) ) }   & = \frac{7x   \epsilon}{2(x^2 - 1)}\\ & \alpha_h^ * =   \frac{3 ( 3 x - n_f ) ( n_fx - n_f^2 + 1 ) ( x^2 - 1 ) } { x(4 n_f^4 x + 9 x^3 - 9 x^5 + n_f^3 ( x^2 - 3 ) - n_f^2 x ( 3 + 23 x^2 ) + 3 n_f ( 1 + 4 x^2 - 9 x^4 ) ) }   & = \frac{\epsilon}{x } \\ & \alpha_{\lambda}^ * = \alpha_{{\widetilde}{\lambda}}^ * \frac{3 ( 3 x - n_f ) ( x + 2 n_f ) ( 3 x^2 + x n_f-1 ) } { 4 n_f^4 x + 9 x^3 - 9 x^5 + n_f^3 ( x^2 - 3 ) - n_f^2 x ( 3 + 23 x^2 ) + 3 n_f ( 1 + 4 x^2 - 9 x^4 ) ) }   & = \frac{7x \epsilon}{x^2 - 1 } = 2\alpha_g^ * + \mathcal{o}(\epsilon^2)\\ & \alpha_m^*=\alpha_{{\widetilde}{m}}^ * = \frac{3 ( 3 x - n_f ) ( 3n_f x+3x^2 - 2)(x^2 - 1 ) } { 2x(4 n_f^4 x + 9 x^3 - 9 x^5 + n_f^3 ( x^2 - 3 ) - n_f^2 x ( 3 + 23 x^2 ) + 3 n_f ( 1 + 4 x^2 - 9 x^4 ) ) }   & = \frac { \epsilon}{x } = \alpha_h^ * + \mathcal{o}(\epsilon^2)\\ \hline \end{array}\ ] ] [ solutions1 ]    solution ( 5 ) in table  [ solutions0 ] corresponds to the perturbative infrared fixed point for super qcd while the last solution ( 6 ) in table [ solutions1 ] corresponds to the perturbative fixed point for seiberg s magnetic dual and it is the only one with all nonvanishing couplings . interestingly , when setting to zero @xmath59 ( solution @xmath65 in table [ solutions0 ] ) we find yet another fixed point corresponding to the infrared stable fixed point of the non - supersymmetric magnetic dual theory proposed in @xcite and shown in table  [ qcdadual ] .
this lends further support to the proposed non - supersymmetric duality .
the stability of the fixed points reported in tables  [ solutions0 ] and [ solutions1 ] is determined by first evaluating the following matrix @xmath66 at the fixed point of interest denoted in brief by @xmath67 .
if it has only real and positive eigenvalues the fixed point is said to be stable . in the tables [ exponents0 ] and [ exponents1 ] we present the eigenvalues corresponding to the fixed points of tables [ solutions0 ] and [ solutions1 ] respectively for @xmath68 and @xmath69 . for convenience ,
we have marked the cells with negative eigenvalues in gray .
@xmath70{0.8}}{\cellcolor{white}}~~1~~ & -0.049 & -0.049&-0.031&-0.016&-0.016 & { \cellcolor{white}}0.0034\\ [ 2 mm ] \hline 2 & 0.20&{\cellcolor[gray]{0.8}}-0.048 & { \cellcolor[gray]{0.8}}-0.045&{\cellcolor[gray]{0.8}}-0.015&0.0063&0.0011\\ [ 2 mm ] \hline 3 & { \cellcolor[gray]{0.8}}-0.049&{\cellcolor[gray]{0.8}}-0.044 & 0.035&{\cellcolor[gray]{0.8}}-0.029&{\cellcolor[gray]{0.8}}-0.015&0.0036\\ [ 2 mm ] \hline 4 & { \cellcolor[gray]{0.8}}-0.064 & 0.038&0.030&{\cellcolor[gray]{0.8}}-0.027&{\cellcolor[gray]{0.8}}-0.023&0.0039\\ [ 2 mm ] \hline 5 & 0.30 & 0.092 & 0.063 & { \cellcolor[gray]{0.8}}-0.060&{\cellcolor[gray]{0.8}}-0.058 & 0.0089\\ [ 2 mm ] \hline 6 & 0.20 & { \cellcolor[gray]{0.8}}-0.062&{\cellcolor[gray]{0.8}}-0.043&0.031 & 0.021 & 0.0065\\ [ 2 mm ] \hline 7 & 3.0 & 1.8 & 1.4 & 1.1 & { \cellcolor[gray]{0.8}}-0.14 & 0.033\\ [ 2 mm ] \hline \end{array}\ ] ] [ exponents0 ]    @xmath71{0.8}}-0.067&{\cellcolor[gray]{0.8}}-0.067&{\cellcolor[gray]{0.8}}-0.0073&{\cellcolor[gray]{0.8}}-0.0073&0.0058\\ [ 2 mm ] \hline 2 & 0.37&0.17 & { \cellcolor[gray]{0.8}}-0.078 & 0.026 & 0.012 & 0.00038\\ [ 2 mm ] \hline 3 & 0.12&{\cellcolor[gray]{0.8}}-0.067&{\cellcolor[gray]{0.8}}-0.064&0.015&{\cellcolor[gray]{0.8}}-0.0075&0.0058\\ [ 2 mm ] \hline 4 & 0.12&{\cellcolor[gray]{0.8}}-0.074&{\cellcolor[gray]{0.8}}-0.053&0.015&0.014&0.058\\ [ 2 mm ] \hline 5 & 0.73 & 0.26 & 0.18 & 0.16&{\cellcolor[gray]{0.8}}-0.081 & 0.018\\ [ 2 mm ] \hline 6 & 21 & 5.4 & 4.0 & 4.0 & 0.77 & 0.052\\ [ 2 mm ] \hline \end{array}\ ] ] [ exponents1 ]    there is a rich structure of fixed points classified according to the relevant directions clearly visible from the list of associated eigenvectors reported in the appendix [ tables ] .
the dimensions of the critical surfaces are dictated by the number of positive eigenvalues .
once the unstable directions , in the coupling space , are removed the remaining subset of couplings generate a critical @xmath72-dimensional surface with @xmath72 given by the number of positive eigenvalues containing the nontrivial fixed points . as an example consider the first fixed point of table  [ exponents0 ] . here
we have a critical line with an infrared stable fixed point and five unstable directions . the critical _ surface _ in this case is just a line .
another interesting example is constituted by the fixed point solution 7 of table  [ solutions0 ] . in this case
we find a five - dimensional critical surface and one unstable sixth direction parallel to the @xmath73 axis .
we plot in figure [ fp7 ] the projections of the renormalization group flow around this solution for a subset of two couplings at the time .
the other filled small circles in the figure correspond to the projections on the chosen planes of fixed points such as the seiberg s magnetic dual , indicated by the susy label on the plot , or the super qcd banks - zaks fixed point labelled by susybz .     and
@xmath69 in the planes ( from left to right and top to bottom ) : ( @xmath74 ) , ( @xmath75 , @xmath73 ) , ( @xmath76,@xmath73 ) and ( @xmath77 , @xmath78 ) . the other couplings are kept fixed at the fixed point value of solution 7 given in table  [ solutions0 ] . besides solution 7 clearly indicated in the plot we also plot the projections of the seiberg s magnetic dual indicated by the susy label on the plot , and the super qcd perturbative fixed point indicated by susybz . ,
title="fig:",scaledwidth=45.0% ]   and @xmath69 in the planes ( from left to right and top to bottom ) : ( @xmath74 ) , ( @xmath75 , @xmath73 ) , ( @xmath76,@xmath73 ) and ( @xmath77 , @xmath78 ) . the other couplings are kept fixed at the fixed point value of solution 7 given in table  [ solutions0 ] . besides solution 7 clearly indicated in the plot we also plot the projections of the seiberg s magnetic dual indicated by the susy label on the plot , and the super qcd perturbative fixed point indicated by susybz .
, title="fig:",scaledwidth=45.0% ] +   and @xmath69 in the planes ( from left to right and top to bottom ) : ( @xmath74 ) , ( @xmath75 , @xmath73 ) , ( @xmath76,@xmath73 ) and ( @xmath77 , @xmath78 ) .
the other couplings are kept fixed at the fixed point value of solution 7 given in table  [ solutions0 ] . besides solution 7 clearly indicated in the plot we also plot the projections of the seiberg s magnetic dual indicated by the susy label on the plot , and the super qcd perturbative fixed point indicated by susybz .
, title="fig:",scaledwidth=45.0% ]   and @xmath69 in the planes ( from left to right and top to bottom ) : ( @xmath74 ) , ( @xmath75 , @xmath73 ) , ( @xmath76,@xmath73 ) and ( @xmath77 , @xmath78 ) . the other couplings are kept fixed at the fixed point value of solution 7 given in table  [ solutions0 ] . besides solution 7 clearly indicated in the plot we also plot the projections of the seiberg s magnetic dual indicated by the susy label on the plot , and the super qcd perturbative fixed point indicated by susybz .
, title="fig:",scaledwidth=45.0% ]    the fixed point featuring all positive eigenvalues is the one corresponding to seiberg s magnetic dual .
we uncovered the full spectrum of perturbative fixed points associated to an @xmath79 nonsupersymmetric gauge theory featuring scalars and fermions , of the type summarized in table  [ qcdadualm ] .
although the noninteracting field theory has a supersymmetric looking spectrum we did not assume the bare couplings to respect supersymmetry .
we then analyzed the fixed points by studying the beta functions to second order in the couplings . to this order , which is renormalization scheme independent , we discovered supersymmetric and nonsupersymmetric fixed points .
we further analyzed their stability and discovered that :    * when all bare couplings are nonzero the theory flows to the seiberg s magnetic dual fixed point .
this occurs on supersymmetric and nonsupersymmetric renormalization flow directions .
* when the yukawa couplings @xmath77 , @xmath80 and @xmath73 are all set to zero the theory flows to the super qcd fixed point , as above , independently if the renormalization flow is supersymmetric .
* when all the yukawas are set to zero we achieve the perturbative nonsupersymmetric fixed point in the gauge coupling . * when @xmath77 and @xmath80 are nonzero and @xmath81 we discover a new nonsupersymmetric fixed point .
this can be identified with the magnetic dual fixed point for nonsupersymmetric gauge theories proposed in @xcite .
another relevant fixed point emerges when @xmath82 and @xmath83 . *
all the other fixed points emerge similarly by setting to zero different yukawas .
the generic conclusions we draw from our results are :    * supersymmetry naturally emerges as a fixed point theory from a nonsupersymmetric lagrangian , meaning without fine - tuning the bare couplings .
this relevant result demonstrates that supersymmetry can be viewed as an emergent phenomenon in field theory .
our results , among other things can be used to argue that one does not need to fine - tune the bare couplings when performing lattice simulations aimed to study supersymmetry on the lattice @xcite . *
duality is not a prerogative of supersymmetry given that the theory features many nonsupersymmetric fixed points , depending on the details of the interactions .
therefore these new fixed points can describe the nonperturbative physics of electric dual gauge theories like the one envisioned in @xcite .
the potential physical applications of our results range from ( supersymmetric ) extensions of the sm to providing new ideas on how to investigate nonperturbative supersymmetric dynamics via first principle lattice simulations .
we thank simon catterall , stefano di chiara and marco nardecchia for enlightening discussions and comments .
we follow the notation of machacek and vaughn@xcite , who derived the general expressions to two - loop order for the running of the gauge coupling and yukawa couplings for a general gauge theory with real scalars and majorana fermions . in ref .
@xcite is showed how to use the expressions of machacek and vaughn for complex scalars and weyl fermions .
we need only to do a slight rewriting of @xmath84 given in eq . to follow their notation : @xmath85 where @xmath86 is the gauge index reserved for the adjoint majorana fermion with @xmath27 the dimension of its representation , greek gauge indices @xmath87 are reserved for the weyl fermions and roman gauge indices @xmath88 are reserved for the complex scalars .
flavor indices have been suppressed .
the yukawa matrices are defines as follows , once we take flavor indices into account with @xmath89 denoting the indices of @xmath90 and @xmath91 indices of the @xmath92 global symmetries :    @xmath93    where @xmath94 are the generators of the gauge group .
the minus sign of @xmath95-term has been chosen to facilitate the connection with supersymmetry , i.e. for susy @xmath96 .
we will first derive the yukawa contribution to the 2-loop beta function , given in eq .
- , in particular : @xmath97 \ , \end{aligned}\ ]
] where @xmath98 is the quadratic casimir of the spinor representations @xmath23 and @xmath27 is the dimension of the gauge group .
it is instructive to consider the contribution from the yukawa sector to the two - loop beta function of the gauge coupling diagrammatically .
there is only one type of diagram that has non - vanishing contribution , which is :        where the solid line represents spinor fields , and the dashed line represents the scalar(s ) coupled via the yukawa interactions .
note that the scalar - gauge interactions do not contribute .
the fermion @xmath16 is not coupled to the gauge fields .
thus we get : @xmath99 \nonumber\\ & = ( \alpha_\lambda + \alpha_{{\widetilde}{\lambda } } ) t({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) n_f \left [ c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) + c_2(g ) \right ] + ( \alpha_m + \alpha_{\widetilde{m } } + 2 \alpha_h ) n_f^2 \frac{c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) x}{d(g ) } \nonumber \\ & =   ( \alpha_\lambda + \alpha_{{\widetilde}{\lambda } } )   \frac { n_f}{2 } \frac{3x^2 - 1}{2x } + \left ( \frac{\alpha_m + \alpha_{\widetilde{m}}}{2 } + \alpha_h \right ) n_f^2,\end{aligned}\ ] ] where we used the notation @xmath31 we now consider the running of the yukawa couplings .
the one loop beta function for the yukawa couplings is given in @xcite : @xmath32 & + 2 y^k y^\dagger_j y^k +   \frac{1}{2}y^k \text{tr}\left [ y^\dagger_k y^j+y^\dagger_j y^k\right ] - 3 g^2 \
{ c_2(r ) , y^j \}.\end{aligned}\ ] ] where @xmath38 is the group invariant : @xmath100 it is again instructive to consider each term in the beta function in terms of the diagrams they originate from .
the first two terms in the bracket corresponds to the fermion leg self - energy contribution from the yukawa sector , and the expression tells us that they each contribute with a factor of @xmath101 to the beta function :    cccc    c @xmath102 +   +   +   +   +    & &    c @xmath103 +   +   +   +   +        the next term in the beta function is the vertex correction from the yukawa sector :        the fourth term is the scalar leg self energy contribution from the yukawa sector :    ccc   &    c @xmath104 +   +   +   +        note the different directions of the arrow in the fermion loops . for dirac fermions the above contribution comes naturally with a factor of two .
finally the gauge sector has a non - vanishing contribution only to the fermion leg self energy , and each term contributes with a factor of -3 :    cccc    c @xmath105 +   +   +   +   +    & &    c @xmath106 +   +   +   +   +        we thus find the one - loop coefficient of the beta function for @xmath33 to be : @xmath107 + 2 { \widetilde}{u}^i u^\dagger_j { \widetilde}{y}^i   \nonumber\\[2 mm ] & + y^k \text{tr}\left [ y^\dagger_k y^j\right ] + y^k \text{tr}\left [ u^k u^\dagger_j\right ] - 3 g^2 \left\ { c_2(q ) y^j + y^j c_2(\lambda ) \right\}.\end{aligned}\ ] ] inserting the definitions of the yukawa - matrices , we find : @xmath108{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}})\mid y_\lambda \mid^2 + n_f \mid \tilde{y}_m \mid^2   + n_f \mid y_h \mid^2 + t({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}})n_f(\mid y_\lambda \mid^2 + \mid y_{\widetilde \lambda } \mid^2 ) \right\ } \right . \nonumber \\ & \quad \left .
+ c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}})\mid y_\lambda \mid^2 + n_f\mid y_m \mid^2 - 3 g^2 \left \ { c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}})+ c_2(g ) \right \ } \right ] y^j   -2 n_f
y_m^ * y_{\widetilde m } y_{\widetilde \lambda } { \left(t^a\right)_\alpha}^j \delta_l^{l ' } \ ] ] note that @xmath109 in the last term is the matrix part of @xmath33 , which cancels on both sides of the equation once we do the rewriting @xmath110 .
the beta function for @xmath111 is obtained via : @xmath112 @xmath113 and @xmath114 are related by a similar transformation .
we compute the former : @xmath115 + 2 { \widetilde}{y}^i y^\dagger_j { \widetilde}{u}^i   \nonumber\\[2 mm ] & + u^k \text{tr}\left [ u^\dagger_k u^j\right ] + u^k \text{tr}\left [ y^\dagger_k y^j\right ] - 3 g^2 u^j c_2({\widetilde}{q}).\end{aligned}\ ] ] using the expression for the matrices we deduce : @xmath116{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) \alpha_{{\widetilde}{\lambda } } + n_f \alpha_h \right\ }    \right.\nonumber \\ & \quad \left .
+ n_f \alpha_m + c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) \alpha_\lambda - 3 \alpha_g    c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}})\right ] u^j   - \frac{2 c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) } { ( 4\pi)^2 } y_{\widetilde \lambda } y_\lambda^ * y_{\widetilde m } \delta_\alpha^j \delta_l^{l ' } \delta_{r}^{r ' } , \ ] ] where we used @xmath117 . as before @xmath118 is the matrix part of @xmath119 , which cancels on both sides of the equation .
the renormalization of the @xmath73 coupling reads : @xmath120\nonumber \\ & + v \text{tr}\left [ v^\dagger v \right ] - 3g^2 \left \
{ c_2(q ) v + v c_2({\widetilde}{q } ) \right \}\nonumber\\[2 mm ] \beta(v)= & \left [ \frac{\alpha_m + \alpha_{{\widetilde}{m } } + 2\alpha_h}{2}n_f   + \left ( \frac{\alpha_\lambda + \alpha_{{\widetilde}{\lambda}}}{2}-6\alpha_g\right)c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) + x \alpha_h \right ] v\end{aligned}\ ] ]    we summarize now the results for beta functions , after having cancelled the matrices :    @xmath121{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) \left [ c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) + c_2(g ) \right ] ( \alpha_\lambda + \alpha_{{\widetilde}{\lambda } } ) +   \frac{c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}})x n_f^2}{d(g ) } ( \alpha_m + \alpha_{\widetilde{m } } + 2 \alpha_h )     \right]\\      & \beta(y_\lambda ) =   \mid y_\lambda \mid e^{i\theta_\lambda }    \left [    \frac{3}{2 } c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) \alpha_\lambda + \frac{t({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) n_f}{2 } ( \alpha_\lambda + \alpha_{{\widetilde}{\lambda } } )   + n_f \left ( \frac{\alpha_{{\widetilde}{m}}+\alpha_h}{2 } + \alpha_{m } \right )   - 3\alpha_g   \left \ { c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt}}}})+ c_2(\text{g } ) \right \}\right ]    \nonumber\\ & \qquad \quad    -2 n_f\frac{\mid y_m \mid \mid y_{\widetilde m}\mid \mid y_{\widetilde \lambda}\mid}{(4\pi)^2}e^{i(\theta_{{\widetilde}{\lambda}}+
\delta \theta_m    ) }    \\[2 mm ]      & \beta ( y_m ) = \mid y_m \mid e^{i\theta_m }    \left [   \frac{3}{2}n_f \left ( \alpha_m + \frac{\alpha_h}{3}\right ) + \frac{x}{2 } ( \alpha_m + \alpha_{{\widetilde}{m}})+   c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) \left ( \frac{\alpha_{{\widetilde}{\lambda}}}{2 } + \alpha_\lambda \right ) - 3\alpha_g c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } )    \right ] \nonumber\\ & \qquad \quad    - 2c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) \frac{\mid y_{\widetilde m}\mid \mid y_\lambda\mid \mid y_{\widetilde \lambda}\mid}{(4\pi)^2}e^{i(\theta_{{\widetilde}{m}}+ \delta \theta_\lambda ) }    \\[2 mm ]   & \beta(\alpha_h ) = 2 \alpha_h \left [ \frac{n_f}{2}(\alpha_m + \alpha_{{\widetilde}{m } } + 2\alpha_h )   + c_2({\raisebox{-.5pt}{{\hbox{\rule{0.4pt}{6.5pt}\hskip-0.4pt\rule{6.5pt}{0.4pt}\hskip-6.5pt\rule[6.5pt]{6.5pt}{0.4pt}}\rule[6.5pt]{0.4pt}{0.4pt}\hskip-0.4pt\rule{0.4pt}{6.5pt } } } } ) \left ( \frac{\alpha_{\lambda } + \alpha_{{\widetilde}{\lambda}}}{2}-6\alpha_g\right)+
x \alpha_h \right]\\    & \beta(\alpha_{{\widetilde}{\lambda } } ) = \beta_{\alpha_{\lambda } } \left ( y_{\lambda }   \leftrightarrow   y_{\widetilde \lambda }   ,   y_{\widetilde m }   \leftrightarrow   y_{m }   \right ) , \quad   \beta(\alpha_{\widetilde{m } } ) = \beta_{\alpha_{m}}\left (   { y}_{m }   \leftrightarrow   y_{\widetilde m } ,   y_{\widetilde \lambda }   \leftrightarrow   { y}_{\lambda }   \right),\end{aligned}\ ] ]    where we have made explicit the phases of the yukawa s couplings : @xmath122 note also that we define @xmath123 where @xmath124 is the complex conjugate of @xmath125 .
we have considered the case of real couplings in the main text .
we provide here the eigenvectors associated to the eigenvalues given in tables [ exponents0 ] and [ exponents1 ] .
shaded cells indicate unstable directions , i.e. corresponding to negative eigenvalues .
f.  sannino and k.  tuominen , phys .
d * 71 * , 051901 ( 2005 ) [ arxiv : hep - ph/0405209 ] .
d.  k.  hong , s.  d.  h.  hsu and f.  sannino , phys .
b * 597 * , 89 ( 2004 ) [ arxiv : hep - ph/0406200 ] .
h.  gies and j.  jaeckel , eur .
j.   c * 46 * , 433 ( 2006 ) [ arxiv : hep - ph/0507171 ] .
j.  braun and h.  gies , phys .
b * 645 * , 53 ( 2007 ) [ arxiv : hep - ph/0512085 ] .
j.  braun and h.  gies , jhep * 0606 * , 024 ( 2006 ) [ arxiv : hep - ph/0602226 ] .
d.  d.  dietrich and f.  sannino , phys .  rev .
d * 75 * , 085018 ( 2007 ) [ arxiv : hep - ph/0611341 ] .
t.  a.  ryttov and f.  sannino , phys
.  rev .
d * 76 * , 105004 ( 2007 ) [ arxiv:0707.3166 [ hep - th ] ] .
t.  a.  ryttov and f.  sannino , phys .
d * 78 * , 065001 ( 2008 ) [ arxiv:0711.3745 [ hep - th ] ] .
j.  braun and h.  gies , jhep * 1005 * , 060 ( 2010 ) [ arxiv:0912.4168 [ hep - ph ] ] . o.  antipin and k.  tuominen , phys .
d * 81 * , 076011 ( 2010 ) [ arxiv:0909.4879 [ hep - ph ] ] .
o.  antipin , k.  tuominen , [ arxiv:0912.0674 [ hep - ph ] ] .
j.  alanen , k.  kajantie and k.  tuominen , phys .
d * 82 * , 055024 ( 2010 ) [ arxiv:1003.5499 [ hep - ph ] ] .
h.  s.  fukano and f.  sannino , phys .
d * 82 * , 035021 ( 2010 ) [ arxiv:1005.3340 [ hep - ph ] ] . c.  pica and f.  sannino , arxiv:1011.3832 [ hep - ph ] . c.  pica and f.  sannino , phys .
d * 83 * , 035013 ( 2011 ) [ arxiv:1011.5917 [ hep - ph ] ] .
m.  t.  frandsen , t.  pickup , m.  teper , phys .
lett .   * b695 * , 231 - 237 ( 2011 ) .
[ arxiv:1007.1614 [ hep - ph ] ]            j.  braun , c.  s.  fischer and h.  gies , arxiv:1012.4279 [ hep - ph ] .
m.  jarvinen and f.  sannino , jhep * 1102 * , 081 ( 2011 ) [ arxiv:1009.5380 [ hep - ph ] ] .
f.  sannino , phys .
d * 79 * , 096007 ( 2009 ) [ arxiv:0902.3494 [ hep - ph ] ] .
t.  a.  ryttov and f.  sannino , int .
j.  mod .
a * 25 * , 4603 ( 2010 ) [ arxiv:0906.0307 [ hep - ph ] ] .
f.  sannino , phys .
d * 80 * , 065011 ( 2009 ) [ arxiv:0907.1364 [ hep - th ] ] .
g. t hooft , _ recent developments in gauge theories _ , plenum press , 1980 , 135 ; reprinted in _ unity of forces in the universe vol .
, a. zee ed . ,
world scientific 1982 , 1004 .
j.  terning , phys .
lett .   * 80 * , 2517 ( 1998 ) [ arxiv : hep - th/9706074 ] .
f.  sannino , nucl .
b * 830 * , 179 ( 2010 ) [ arxiv:0909.4584 [ hep - th ] ] .
m.  mojaza , m.  nardecchia , c.  pica and f.  sannino , phys .
d * 83 * , 065022 ( 2011 ) [ arxiv:1101.1522 [ hep - th ] ] .
f.  sannino , arxiv:1102.5100 [ hep - ph ] .
f.  sannino , phys .
d * 82 * , 081701 ( 2010 ) [ arxiv:1006.0207 [ hep - lat ] ] .
f.  sannino , phys .
lett .   * 105 * , 232002 ( 2010 ) [ arxiv:1007.0254 [ hep - ph ] ] .
s.  di chiara , c.  pica and f.  sannino , arxiv:1008.1267 [ hep - ph ] .
m.  e.  machacek and m.  t.  vaughn , nucl .
b * 222 * , 83 ( 1983 ) .
m.  e.  machacek and m.  t.  vaughn , nucl .
b * 236 * , 221 ( 1984 ) .
m.  x.  luo , h.  w.  wang and y.  xiao , phys .
d * 67 * , 065019 ( 2003 ) [ arxiv : hep - ph/0211440 ] .
s.  catterall , e.  dzienkowski , j.  giedt , a.  joseph and r.  wells , jhep * 1104 * , 074 ( 2011 ) [ arxiv:1102.1725 [ hep - th ] ] .
s.  catterall , arxiv:1005.5346 [ hep - lat ] .