arXiv:0906.5492v2 [hep-th] 19 Nov 2009 n lcrn nahnyobltiegvrstephase the governs lattice honeycomb a interact- on attractively electrons of the instance, ing strength For as interaction serves transition. inverse phase which the super- the of of a term parameter e.g., mass control parameter, the order condensate, bosonic fluid governed composite is correlated a breaking strongly by symmetry very here, in is systems: This transitions fermionic term. phase mass superpotential bosonic to the a is reminiscent of of role systems parameter the supersymmetric plays control which other a many by and governed this in break- ing Dynamical spacetime. Euclidean with models phenomena. critical and transitions phase in of breaking light supersymmetry the explores statistical to paper analogies this from parti- systems, profited and often theory has symmetry field the of cle quantum in understanding in flow the phenomena RG breaking As the a by regime. near determined fixed-point system ad- are the In transition of (RG). to phase properties group related quantitative often the the dition, are of symmetry states points realized ground fixed differently between a transitions of phase and breaking fluctuation- the dynamics. of induced understanding nonperturbative a re- requir- often ing is phenomena, condensation breaking collective symmetry to low-energy lated Generically, the of formation model. the standard in role decisive a play [1] can breaking break- supersymmmetry of and dynamical variety gauge sectors, hidden a particle further to involving often addition mechanisms In ing standard- the data. attractive by precision these imposed model na- violating constraints other without in way and broken features specific be a must in supersymmetry ture hand, On Poincar´e other group. the unification combina- the with the intricate symmetries internal an supports of facilitates tion scales, and couplings, of gauge hierarchy of problem large the model parametrizing a helps standard of it the as physics, beyond particle models of of construction the sasml xml,w ocnrt nWess-Zumino on concentrate we example, simple a As symmetry physics, statistical of viewpoint the From for principle guiding important an is Supersymmetry hoeic-hsklshsIsiu,Friedrich-Schill Institut, Theoretisch-Physikalisches ftefil cln eainwihi o nw nstandard in known not 05.10.Cc,12.60.Jv,11.30.Qc is numbers: which PACS relation scaling p a control – the field with the associated of exponent Supersymmetry critical the systems. direct p between fermion quantum relevant a correlated common as in understood a transitions be is thus can which breaking superpotential supersymmet Supersymmetry of the breaking of us dynamical eter the supercharge relate one We (RG). with group models Wess-Zumino supersymmetric N eepoesprymtybekn ntelgto ihfixed rich a of light the in breaking supersymmetry explore We uecagsi w iesosin dimensions two in supercharges 1 = uesmer raiga unu hs transition phase quantum a as breaking Supersymmetry ogrGe,FazsaSntck n nra Wipf Andreas and Synatschke Franziska Gies, Holger rUiesta ea a-inPaz1 -74 ea Ger Jena, D-07743 1, Max-Wien-Platz er-Universit¨at Jena, h eiaiewt epc to respect with derivative the ftebr ueptnili odd. is superpotential power bare highest the the of if breaking supersymmetry dynamical the for resentation fields superpotential with r noe diinlsrcue hc obyn those physics. beyond statistical go in which known structures additional invokes supersymme- correlated try similarities, strongly these Beyond in similari- transition systems. will fermion phase many quantum we has a symmetry, to breaking ties this supersymmetry from supersym- that preserve Apart show on to freedom manifestly. mandatory of even metry degrees of is these degree footing of equal high treatment a A with symmetry. systems bosonic-fermionic bined proved recently has theory [6]. foot- Hertz-Millis successful equal the on not freedom beyond is of ing fermionic strategy degrees fundamental bosonic this of composite As and treatment a dealt approximation. sufficient, is local always which a description in bosonic with purely a leaving integrated intro- are out, freedom are of Sub- degrees fields fermionic transformation. the bosonic sequently, Hubbard-Stratonovich composite a Hertz-Millis by where duced the 5] is [4, systems theory fermion correlated transition. strongly phase quantum marks a parameter of control point the critical quantum of temper- a value zero [2] at critical atoms occurs the transition ultracold this ature, of As context [3]. superfluid the graphene or in the discussed and as semimetallic phase the between transition l a edfie yteLgag est nteoff-shell the in density Lagrange the formulation by defined be can els euetefntoa Gt aclt h effective the calculate to RG functional the use We ulda two-dimensional Euclidean com- are definition by systems supersymmetric Now, in transitions phase quantum to approach standard A L L kin int ψ ie iet e uesaigrelation superscaling new a to rise gives rmtradteaoaosdimension anomalous the and arameter and = = o falfie onso h system. the of points fixed all of ion yt nR eeatcnrlparam- control relevant RG an to ry 2 2 1 1 pnsystems. spin pitsrcueo two-dimensional of structure -point aetasto nlgu osimilar to analogous transition hase n h ucinlrenormalization functional the ing W ∂ L ψ µ = φ∂ ′′ n uiir field auxiliary and ( φ L µ ) kin φ ψγ γ + W + ∗ arcs hsmdlalw for allows model This matrices. ψ 4 i ( L ψ φ − int ∂ψ ,bsncfield bosonic ), / W with − ′ N φ ( φ n eueacia rep- chiral a use we and , 4 i ) esZmn mod- Wess-Zumino 1 = F, ∂ F µ h rm denotes prime The . ψγ µ ψ − φ 2 1 fermionic , F many 2 , (1) (2) 2 action Γ of the theory. The functional RG can be formu- Let us start with the lowest-order derivative expansion, lated as a flow equation for the effective average action Γk corresponding to a local-potential approximation for the [7]. This is a scale-dependent action functional which in- superpotential. In this truncation, the effective action is terpolates between the classical action S = ( + ) approximated by Lkin Lint and the full quantum effective action Γ. TheR interpola- 1 ′′ ′ tion scale k denotes an infrared (IR) regulator scale, such Γk = dτ kin + Wk (φ)ψγ∗ψ Wk(φ)F (4) that no fluctuations are included for k Λ (with Λ be- Z L 2 −  → ing the UV cutoff), implying Γ → S. For k 0, all k Λ → → with a k-dependent superpotential Wk. Projecting fluctuations are taken into account such that Γ → Γ k 0 → the flow equation (3) onto this truncation leads to is the full quantum solution. The effective average action the flow equation for the superpotential. For a sim- can be determined from the Wetterich equation [7] plified regulator choice with r1 = 0 and r2 = 2 2 −1 ( k/p 1) θ 1 (p /k ) , the superpotential flow re- 1 (2) k | |− − ∂tΓk = STr Γk + Rk ∂tRk , t = ln , (3) duces to  2 h i  Λ 2 ′′ k Wk (φ) which defines an RG flow trajectory in theory space ∂tWk(φ)= 2 ′′ 2 . (5) (2) − 4π k + Wk (φ) with S serving as initial condition. Here, Γk = −→ ←−  ab This regulator choice is convenient for analytical studies. δ Γ δ where the indices a,b summarize field com- δΨa k δΨb We have verified explicitly by using a variety of differ- ponents, internal and Lorentz indices, as well as space- ent regulators that all conclusions drawn below from this T ¯ time or momentum coordinates, i.e., Ψ = (φ, F, ψ, ψ). simple choice also hold for other regulators. The regulator function Rk is introduced in the form of The corresponding flow equation for the dimensionless a quadratic contribution to the classical action, ∆Sk = potential w = W /k admits a variety of nontrivial fixed- T k k Ψ RkΨ. Rk can thus be viewed as a momentum- point solutions: already a polynomial expansion for Z2 dependentR mass term, implementing the IR suppression antisymmetric superpotentials to order n, of modes below k. Different functional forms of Rk cor- respond to different RG schemes. Physical quantities do ′ ∂wk 2 2 4 6 2n w = λ(φ a )+b4φ +b6φ + +b2nφ , (6) not depend on the regulator. k ≡ ∂φ − ··· For a manifestly supersymmetric flow equation, we de- reveal n independent solutions to the fixed-point equa- rive the regulator from a D term of a quadratic superfield tion ∂ w = 0, in addition to the Gaußian fixed point[20]. operator, similarly to our construction for supersymmet- t k The latter has infinitely many relevant directions in ric quantum mechanics [8]. The most general supersym- agreement with perturbative power-counting in two di- metric cutoff action in the off-shell formulation in mo- mensions. The number of relevant RG directions de- mentum space is given by creases for the other fixed points, culminating in a maxi- 1 2 2 2 2 mally IR-stable fixed point with only one relevant direc- ∆Sk = d p(r1ψγ∗ψ 2r1φF +p r2φ +r2ψpψ/ r2F ). 2Z − − tion. This relevant direction exactly corresponds to the parameter a2 in the superpotential, and it is common 2 The regulator shape functions r1,2 = r1,2(p ) determine to all other fixed points as well. As a2 is related to the the precise form of the Wilsonian momentum shell inte- bosonic mass term, we introduce a νW gration. The form of ∆Sk guarantees that the whole RG for this direction, such that a2 scales as a2 k−1/νW at trajectory remains in the hypersurface of supersymmet- a fixed point. This critical exponent thus corresponds∼ to ric action functionals even within approximative calcula- the negative inverse of the associated eigenvalue of the tions; see [9, 10] for further supersymmetric RG studies. stability matrix. The exponent νW plays a similar role In general, the full quantum effective action consists for the superpotential, as the exponent ν for the effective of all possible operators which are compatible with the potential in Ising-type systems. At leading order in the symmetries of the theory. In this paper, we use a super- derivative expansion, we obtain from Eq. (5): νW LO = 1. covariant derivative expansion to systematically classify Fixed-point potentials and critical exponents at| next-to- these operators. Truncating the expansion at a given or- leading order will be discussed below. der constitutes a consistent approximation scheme. In a In order to clarify the role of this relevant direction variety of bosonic and fermionic systems, such expansions and its relation to supersymmetry breaking, we study have been proven to capture the physics of phase transi- the phase diagram of a Gaußian Wess-Zumino model at tion and critical phenomena quantitatively, see, e.g. [11]. infinite volume. This is defined in terms of a quadratic For our quantitative results presented below, we observe perturbation of the Gaußian fixed point at the UV cut- a satisfactory convergence: leading-order results generi- ′ ¯ 2 2 off Λ: WΛ = λΛ(φ a¯Λ). The dimensionless control cally receive (1 10%) corrections at next-to-leading parameter δ and coupling− γ can be associated with order, increasingO − up to 50% for observables in the ∼ −1 ¯ 2 −1 ¯ infinite-coupling limit. δ := Λ λΛa¯Λ, γ := 2Λ λΛ, (7) 3

2 of higher orders of Eq. (6), we observe a rapid con- φ Λ γφ2 δ2 t= 0.0 W'Λ( )= (0.5 - ) vergence of the phase boundary; see Fig. 1 (bottom). δ γ t= - 0.6 ] =0.03, =0.2

2 t= - 3.0 1.5 There exists a critical value for δcr characterizing the Λ t= - 10 [ in the strong coupling limit γ . -3 → ∞ 1 Our best estimate from a numerical higher-order solu- )*10 φ

( tion is δcr 0.263; (the lowest-order result from Eq. (9) 2 k ≃

W' 0.5 is π/96 0.181). We conclude that supersymmetry ≃ canp never be broken dynamically above this critical value. 0 This agrees qualitatively with earlier results in the litera- -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 φ ture [12, 13, 14, 15]. Of course, the numerical value for δcr is not universal but regulator-scheme dependent, simply 0.26 because the quadratic initial conditions near the Gaußian 0.24 fixed point are nonuniversal. A quantitative comparison unbroken SuSy 0.22 with other estimates of this critical value, e.g., taken from 0.2 the lattice, is thus only meaningful if this scheme depen- 0.18 0.16 dence is accounted for. δ 0.14 The mass of the lowest bosonic excitation is given by 0.12 10 the curvature of the full renormalized effective potential 0.1 broken SuSy φ φ8 at the minimum. In the supersymmetric phase, fermionic 0.08 φ4 and bosonic masses are identical and nonzero. In the bro- 0.06 φ2 0.04 ken phase, the lowest fermionic excitation is a Goldstino 0 2 4 6 8 10 γ/2 with vanishing mass. The latter holds also at finite k in the supersymmetry-broken regime. As the flow in the FIG. 1: Top: Typical flow of the superpotential in the IR is attracted by the maximally IR-stable fixed point, supersymmetry broken phase. Bottom: Phase diagram in the bosonic mass in the fixed-point regime is governed by the coupling–control-parameter plane for the Gaußian Wess- the only relevant direction a2. As the minimum of the Zumino model. The lowest set of data points corresponds to effective potential in the broken regime is necessarily at the lowest-order result (9). Higher orders converge rapidly. a vanishing field, the bosonic mass yields

1 2− m2 =2k2λ2 a2 k νW , (10) where the scale is set by the UV cutoff Λ. In the follow- k | |∼ ing, we determine the region for the parameters δ, γ for where the latter proportionality holds in the vicinity of which supersymmetry is dynamically broken in the IR. the IR fixed point where λ const. and a2 is governed A criterion for supersymmetry breaking is provided by by the critical exponent ν .→ For ν > 1/2, the bosonic a nonvanishing ground state energy, given by the min- W W ′ mass scales to zero upon attraction by the maximally IR- imal value of V (φ) = 1 (W (φ))2. As a true phase 2 k→0 stable fixed point. We conclude that the broken phase transition only occurs in the infinite-volume limit [12], remains massless in both bosonic and fermionic degrees the ground state energy can only be identified at k 0. of freedom. Four snapshots of a typical flow of the superpotential→ in On the other hand, the limit k 0 is an idealized IR the broken phase are depicted in Fig. 1 (top) for δ =0.03 → and γ =0.2. limit. Any real experiment as well as any lattice simu- For a simple estimate of the phase diagram, we trun- lation will involve an IR cutoff scale km characterizing cate the superpotential at order φ2 in Eq. (6). From the measurement, e.g., the scale of momentum trans- Eq. (5), we obtain the flow of the dimensionless couplings fer, detector size or lattice volume. Any measurement is therefore not sensitive to k 0 but to k k > 0. → → m 1 6λ2a2 6λ3 We conclude that a measurement of the bosonic mass in ∂ a2 = , ∂ λ = λ + , (8) t 2π − π t − π the broken phase will give a nonzero answer proportional to the measurement scale, whereas the goldstino will be which can be solved analytically. The phase transition truly massless. The bosonic mass in both phases as a curve in the (γ,δ) plane is given by the initial values for function of the control parameter δ is depicted in Fig. 2 which a2 → changes sign. This condition yields |k 0 (top) for a coupling of γ =0.2. At next-to-leading order in the derivative expansion, arcsin(α) 2 π δ = , α =1 . (9) the fields acquire a wave function renormalization Z , √24π α − 3γ k which is a multiplicative factor of the kinetic term in 2 Beyond this simple estimate, the phase transition curve Eq. (1), kin Zk kin. The flow of Zk gives rise to an L → L 2 can, of course, be identified by a full numerical inte- anomalous dimension, η = ∂t ln Zk . Convergence of the gration of the flow equation. Already upon inclusion derivative expansion requires− η . (1). O 4

# nodes η νW 6 0 0.4386 1.2809 t= - 9 5 t= -10 broken SuSy unbroken SuSy 1 0.20 1.11 t= -11 t= -12 2 0.12 1.06 4 ] Λ ] t= -13 Λ

[ 3 k /25 [ k m

TABLE I: Critical exponents of the first fixed points. m 2

1 From the polynomially expanded superpotential flow 0 at next-to-leading order we can derive the following flow 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 of the renormalized parameter a2: δ

1 η η a2 5 ∂ a2 = 1 1 a2 t ∂ λ . (11) η=0.4386 t t t t t η 2π − 4 − − 2 − λt 4 =0.20     η=0.12 At any fixed point with ∂tλ 0, we can read off that 3 2 → 2 a again denotes a relevant direction with scaling a ) χ −1/νW ∼ 2 k , where u( 1 2 ν = . (12) W 2 η 0 − This remarkable superscaling relation connects the su- -1 0 0.5 1 χ 1.5 2 2.5 perpotential exponent νW with the anomalous dimension η. Recall that in Ising-like systems the thermodynamic main exponents (i.e., α, β, γ and δ) are related among FIG. 2: Top: Scale dependence of the renormalized bosonic mass. Bottom: Fixed-point superpotentials at next-to- each other by scaling relations, and can be deduced from ′′ 2 the correlation exponents ν and η by hyperscaling rela- leading order derivative expansion, χ = Zkφ, u = wk /Zk . tions. Beyond that, there is no general relation between ν and η. The superscaling relation (12) thus represents a special feature of the present supersymmetric model and of our numerical results. From our regulator studies, we is an exact relation to next-to-leading order in the super- found a maximum variation of 10% on the positive criti- covariant derivative expansion. Beyond next-to-leading cal exponents. Subleading negative exponents can show order, Eq. (12) might receive corrections, since higher- larger variations and thus are difficult to estimate reli- derivative operators involving the auxiliary field can still ably, whereas the leading exponents vary even less, if at take influence on the flow of the superpotential. Never- all [19]. theless, even if Eq. (12) receives higher-order corrections, Each fixed point constitutes a universality class and the most important result that a potentially more in- defines a distinct asymptotically safe UV completion of volved scaling relation among νW and η exists still holds. the Wess-Zumino model, i.e., RG trajectories emanating Finally, we determine the possible nontrivial fixed- from different fixed points correspond to different theo- point superpotentials. Solving the coupled fixed-point ries. Our findings are reminiscent of those in 2D scalar equations for the dimensionless superpotential as well as field theories [16, 17], where the fixed points of the effec- the anomalous dimension self-consistently, we find a dis- tive potential can be related to conformal field theories crete set of fixed-point superpotentials with an increas- [18]. The connection between the present supersymmet- ing number of nodes, see Fig. 2 (bottom). The number ric models at their fixed points and conformal field theo- of RG relevant directions at the fixed points in addition ries remains an interesting question. to the a2 direction is equal to the number of nodes. The To summarize, we have constructed a manifestly su- superpotential with no nodes is the next-to-leading or- persymmetric flow equation for 2D = 1 Wess-Zumino N der analogue of the maximally IR-stable fixed point dis- models. These models have nontrivial fixed-point super- cussed above. Our quantitative next-to-leading order es- potentials which can be classified by their relevant direc- timates of the critical exponent νW and the anomalous tions. All fixed points share a relevant direction in the dimension η for the fixed-point superpotentials of Fig. 2 form of a bosonic mass term. At next-to-leading order, (bottom) are summarized in Tab. I. As these quantities the associated critical exponent νW and the anomalous are universal, these results should constitute regulator- dimension η satisfy a superscaling relation which has no independent predictions. Of course, the truncation of the analogue in Ising-type spin systems. effective action introduces an artificial regulator depen- The fixed point of this relevant direction corresponds dence, the size of which can be used for an error estimate to the critical point separating the supersymmetric from 5 the broken phase. The initial condition for this rele- [5] A. J. Millis, Phys. Rev. B 48, 7183 (1993). vant direction defines a control parameter, such that su- [6] P. Strack, S. Takei, and W. Metzner (2009), persymmetry breaking can be understood as a quantum arXiv:0905.3894 [cond-mat.str-el]. B301 phase transition. For the Gaußian Wess-Zumino model, [7] C. Wetterich, Phys. Lett. , 90 (1993). [8] F. Synatschke, G. Bergner, H. Gies, and A. Wipf, JHEP we observe that the control parameter stays finite even at 03, 028 (2009). arbitrarily large coupling in accord with a general argu- [9] H. Sonoda and K. Ulker (2008), arXiv:0804.1072. ment by Witten. In the broken phase, our superscaling [10] O. J. Rosten (2008), arXiv:0808.2150. relation predicts that the measured bosonic mass is pro- [11] B. Delamotte (2007), arXiv:cond-mat/0702365. portional to the momentum scale set by the detector. [12] E. Witten, Nucl. Phys. B202, 253 (1982). Helpful discussions with G. Bergner, C. Wozar, T. Fis- [13] M. Beccaria, G. F. De Angelis, M. Campostrini, and D70 chbacher and T. Kaestner are gratefully acknowledged. A. Feo, Phys. Rev. , 035011 (2004). [14] M. Beccaria, M. Campostrini, and A. Feo, Phys. Rev. This work has been supported by the Studienstiftung D69, 095010 (2004). des deutschen Volkes, and the DFG under GRK 1523, [15] J. Ranft and A. Schiller, Phys. Lett. B 138, 166 (1984). Wi 777/10-1, FOR 723 and Gi 328/5-1. [16] T. R. Morris, Phys. Lett. B345, 139 (1995). [17] R. Neves, Y. Kubyshin, and R. Potting (1998), hep- th/9811151. [18] A. Zamolodchikov, Yad. Fiz. 44, 529 (1986). [19] F. Synatschke, H. Gies, and A. Wipf (2009), [1] E. Witten, Nucl. Phys. B188, 513 (1981). arXiv:0907.4229 [hep-th]. [2] E. Zhao and A. Paramekanti, Phys. Rev. Lett. 97, 230404 [20] Beyond the polynomial expansion, we find fixed-point so- (2006). lutions corresponding to oscillating sine-Gordon-type po- [3] B. Uchoa and A. H. C. Neto, Phys. Rev. Lett. 98, 146801 tentials, bouncing potentials as well as potentials which (2007). confine the field to a compact target space; see [19] for a [4] J. A. Hertz, Phys. Rev. B 14, 1165 (1976). detailed discussion.