Skip to content

Superpotentials and a first interacting model

The symmetry lesson tested candidate interactions by their charges. We will now construct one and determine its dynamics. Start with one gauge-neutral chiral multiplet and the superpotential

W(ϕ)=m2ϕ2+g3ϕ3,m≠0,g≠0.W(\phi)=\frac m2\phi^2+\frac g3\phi^3,\qquad m\ne0,\quad g\ne0.

A superpotential is a holomorphic function: it depends on the chiral scalar ϕ\phi, not on the independent anti-chiral scalar ϕ~\widetilde\phi. It organizes several component interactions. It is not itself the scalar potential energy.

The target is a three-dimensional Wess–Zumino model, meaning chiral matter with a superpotential and no gauge sector. We work in flat Euclidean space, with all background sources and the central charge set to zero. We preserve the full complexified N=2\mathcal N=2 algebra, whose Lorentzian parent has four real supercharges. The kinetic normalization is that of the free chiral theory. All vacuum and mass statements below are classical, or equivalently tree level.

The dimensions are [ϕ]=1/2[\phi]=1/2, [ψ]=1[\psi]=1, [FΦ]=3/2[F_\Phi]=3/2, and [W]=2[W]=2. Thus [m]=1[m]=1 and [g]=1/2[g]=1/2. These are engineering dimensions, not a determination of an interacting scaling dimension. Aharony and collaborators describe this class of three-dimensional models in section 2.1, equation (2.2).

Keep tilded and untilded Euclidean fields independent. Write a second holomorphic function W~(ϕ~)\widetilde W(\widetilde\phi) for the anti-chiral partner. Its coefficients are initially independent too. Complex conjugation will be imposed only when discussing the physical scalar potential and an integration contour.

Primes mean derivatives with respect to the function’s own scalar argument. The interaction is

LW=W′FΦ−12W′′ψψ+W~′F~Φ+12W~′′ψ~ψ~.\boxed{ \mathcal L_W= W'F_\Phi-\frac12W''\psi\psi +\widetilde W'\widetilde F_\Phi +\frac12\widetilde W''\widetilde\psi\widetilde\psi . }

A term quadratic in fermions with a scalar-dependent coefficient is a Yukawa interaction. For our polynomial, the untilded terms are (mϕ+gϕ2)FΦ−12(m+2gϕ)ψψ(m\phi+g\phi^2)F_\Phi-\tfrac12(m+2g\phi)\psi\psi. The relative factors and the different tilded Yukawa sign are essential.

We can derive them without a superspace prerequisite. Use the contractions of the earlier component lessons:

uv=uTCv,C=(0−110),ψψ=−2ψ1ψ2.uv=u^TCv,\qquad C=\begin{pmatrix}0&-1\\1&0\end{pmatrix}, \qquad \psi\psi=-2\psi_1\psi_2.

For a commuting label ζ\zeta and an odd two-component spinor ψ\psi,

(ζψ)ψ=−12ζ(ψψ).(\zeta\psi)\psi=-\frac12\zeta(\psi\psi).

For example, the first component on either side is ζ1ψ1ψ2\zeta_1\psi_1\psi_2. This identity is where the factor 1/21/2 enters.

The bottom component W(ϕ)W(\phi) transforms as δζW=2ζ(W′ψ)\delta_\zeta W=\sqrt2\zeta(W'\psi), so its composite fermion is ψW=W′ψ\psi_W=W'\psi. Apply another transformation:

δζψW=W′′(δζϕ)ψ+W′δζψ=2W′′(ζψ)ψ+2W′ζFΦ=2ζ(W′FΦ−12W′′ψψ).\begin{aligned} \delta_\zeta\psi_W &=W''(\delta_\zeta\phi)\psi+W'\delta_\zeta\psi\\ &=\sqrt2W''(\zeta\psi)\psi+\sqrt2W'\zeta F_\Phi\\ &=\sqrt2\zeta\left(W'F_\Phi-\frac12W''\psi\psi\right). \end{aligned}

The quantity in parentheses is the composite auxiliary component FWF_W. In the anti-chiral sector, δζ~ϕ~=−2ζ~ψ~\delta_{\widetilde\zeta}\widetilde\phi=-\sqrt2\widetilde\zeta\widetilde\psi has an extra minus sign. The same calculation gives

ψ~W=W~′ψ~,F~W=W~′F~Φ+12W~′′ψ~ψ~.\widetilde\psi_W=\widetilde W'\widetilde\psi,\qquad \widetilde F_W =\widetilde W'\widetilde F_\Phi +\frac12\widetilde W''\widetilde\psi\widetilde\psi.

Thus LW=FW+F~W\mathcal L_W=F_W+\widetilde F_W. This is the F-term construction: use the highest component of a composite chiral multiplet and its anti-chiral partner as a Lagrangian. The signs agree with the chiral embeddings in P, (6.7) and (6.9) and the product rule (6.16). The paper’s (6.14) states the curved-background construction with its R-charge condition. Below we verify the flat interaction directly, including cases that do not have a continuous R-symmetry.

First calculation: eliminate the auxiliaries and find the vacua

Section titled “First calculation: eliminate the auxiliaries and find the vacua”

Add the interaction to the normalized free action:

L=∂μϕ~ ∂μϕ−iψ~γμ∂μψ−F~ΦFΦ+W′FΦ+W~′F~Φ−12W′′ψψ+12W~′′ψ~ψ~.\begin{aligned} \mathcal L={}& \partial_\mu\widetilde\phi\,\partial_\mu\phi -i\widetilde\psi\gamma^\mu\partial_\mu\psi -\widetilde F_\Phi F_\Phi\\ &+W'F_\Phi+\widetilde W'\widetilde F_\Phi -\frac12W''\psi\psi +\frac12\widetilde W''\widetilde\psi\widetilde\psi. \end{aligned}

The auxiliary fields have no derivatives. Vary them independently:

∂L∂FΦ=−F~Φ+W′=0,∂L∂F~Φ=−FΦ+W~′=0.\frac{\partial\mathcal L}{\partial F_\Phi} =-\widetilde F_\Phi+W'=0,\qquad \frac{\partial\mathcal L}{\partial\widetilde F_\Phi} =-F_\Phi+\widetilde W'=0.

Therefore

FΦ=W~′,F~Φ=W′.F_\Phi=\widetilde W',\qquad \widetilde F_\Phi=W'.

Substitution gives the scalar term +W′W~′+W'\widetilde W', not its negative. A useful way to see all signs at once is

−F~ΦFΦ+W′FΦ+W~′F~Φ=−(F~Φ−W′)(FΦ−W~′)+W′W~′.-\widetilde F_\Phi F_\Phi+W'F_\Phi+\widetilde W'\widetilde F_\Phi =-(\widetilde F_\Phi-W')(F_\Phi-\widetilde W') +W'\widetilde W'.

Now choose ϕ~=ϕ∗\widetilde\phi=\phi^* and W~(ϕ~)=W(ϕ~‾)‾\widetilde W(\widetilde\phi)=\overline{W(\overline{\widetilde\phi})}, so the tilded coefficients are the complex conjugates of m,gm,g. The scalar potential is

V(ϕ,ϕ∗)=∣W′(ϕ)∣2=∣mϕ+gϕ2∣2=∣ϕ∣2∣m+gϕ∣2.\boxed{ V(\phi,\phi^*)=|W'(\phi)|^2 =|m\phi+g\phi^2|^2 =|\phi|^2|m+g\phi|^2 . }

It appears with a plus sign in the Euclidean action. After continuation to Lorentzian signature it is the nonnegative potential energy with the canonical scalar kinetic normalization.

A constant bosonic configuration preserves all flat supercharges if its fermion variations vanish. With vanishing gradients and fermions, δψ=2ζFΦ\delta\psi=\sqrt2\zeta F_\Phi and δψ~=2ζ~F~Φ\delta\widetilde\psi=\sqrt2\widetilde\zeta\widetilde F_\Phi. Hence both auxiliaries must vanish. After elimination this becomes the F-term vacuum equation

W′(ϕ)=0,W'(\phi)=0,

together with its complex conjugate. For the stated nonzero couplings, factorization gives two distinct classical supersymmetric vacua:

v1=0,v2=−mg,V(v1)=V(v2)=0.v_1=0,\qquad v_2=-\frac m g,\qquad V(v_1)=V(v_2)=0.

These are actual zero-energy minima of the nonnegative potential. Solving the weaker stationarity equation ∂V=0\partial V=0 can also find nonsupersymmetric saddles; exercise 2 checks one explicitly.

Around either vacuum write ϕ=v+h\phi=v+h. Define μv=W′′(v)\mu_v=W''(v); it is a mass parameter for fluctuations, not a new coupling. Then

W′(v+h)=μvh+gh2,μv1=m,μv2=−m.W'(v+h)=\mu_v h+gh^2,\qquad \mu_{v_1}=m,\quad \mu_{v_2}=-m.

The quadratic scalar potential is ∣μv∣2∣h∣2|\mu_v|^2|h|^2. Writing h=(h1+ih2)/2h=(h_1+ih_2)/\sqrt2 gives two canonically normalized real scalar modes, each with mass squared ∣m∣2|m|^2. The quadratic Yukawa terms have coefficients μv,μv∗\mu_v,\mu_v^*:

Lfermion(2)=−iψ~γμ∂μψ−12μvψψ+12μv∗ψ~ψ~.\mathcal L_{\rm fermion}^{(2)} =-i\widetilde\psi\gamma^\mu\partial_\mu\psi -\frac12\mu_v\psi\psi +\frac12\mu_v^*\widetilde\psi\widetilde\psi.

Using left Grassmann variations gives

iγμ∂μψ=μv∗ψ~,iγμ∂μψ~=−μvψ.i\gamma^\mu\partial_\mu\psi=\mu_v^*\widetilde\psi,\qquad i\gamma^\mu\partial_\mu\widetilde\psi=-\mu_v\psi.

Apply the Dirac operator again and use (γμ∂μ)2=∂2(\gamma^\mu\partial_\mu)^2=\partial^2:

(−∂2+∣μv∣2)ψ=0.(-\partial^2+|\mu_v|^2)\psi=0.

The same equation holds for ψ~\widetilde\psi. The Euclidean inverse propagator has the mass scale ∣μv∣=∣m∣|\mu_v|=|m|, agreeing with the scalar modes. This is a direct tree-level mass check at both vacua; it is not a calculation of quantum mass renormalization.

The auxiliary elimination used equations of motion. To prove the action’s symmetry, return to the full off-shell fields. Constant commuting labels generate odd operations. Keep the graded product rule δ(XY)=(δX)Y+(−1)∣X∣XδY\delta(XY)=(\delta X)Y+(-1)^{|X|}X\delta Y. The relevant untilded transformations are

δζϕ=2ζψ,δζψ=2ζFΦ,δζFΦ=0,δζ~ϕ=0,δζ~ψ=−2iγμζ~ ∂μϕ,δζ~FΦ=−2iζ~γμ∂μψ.\begin{aligned} \delta_\zeta\phi&=\sqrt2\zeta\psi,& \delta_\zeta\psi&=\sqrt2\zeta F_\Phi,& \delta_\zeta F_\Phi&=0,\\ \delta_{\widetilde\zeta}\phi&=0,& \delta_{\widetilde\zeta}\psi&=-\sqrt2i\gamma^\mu\widetilde\zeta\,\partial_\mu\phi,& \delta_{\widetilde\zeta}F_\Phi&=-\sqrt2i\widetilde\zeta\gamma^\mu\partial_\mu\psi. \end{aligned}

First, δζFW=0\delta_\zeta F_W=0. Indeed,

δζ(ψψ)=22FΦζψ,\delta_\zeta(\psi\psi)=2\sqrt2F_\Phi\zeta\psi,

so the variation of W′FΦW'F_\Phi cancels that of the quadratic fermion term. The remaining candidate term is proportional to W′′′(ζψ)(ψψ)W'''(\zeta\psi)(\psi\psi) and vanishes: it contains three factors chosen from the same two odd components. This is a local classical Grassmann identity; no fermion equation is involved.

For the opposite family, use (γμζ~)ψ=−ζ~γμψ(\gamma^\mu\widetilde\zeta)\psi=-\widetilde\zeta\gamma^\mu\psi. The two contributions are

δζ~FW=−2i W′ζ~γμ∂μψ−2i W′′(∂μϕ)ζ~γμψ=−2i ∂μ(ζ~γμW′ψ).\begin{aligned} \delta_{\widetilde\zeta}F_W &=-\sqrt2i\,W'\widetilde\zeta\gamma^\mu\partial_\mu\psi -\sqrt2i\,W''(\partial_\mu\phi)\widetilde\zeta\gamma^\mu\psi\\ &=-\sqrt2i\,\partial_\mu \left(\widetilde\zeta\gamma^\mu W'\psi\right). \end{aligned}

The anti-chiral calculation must respect its different Yukawa sign. Its rules are

δζ~ϕ~=−2ζ~ψ~,δζ~ψ~=2ζ~F~Φ,δζ~F~Φ=0,δζϕ~=0,δζψ~=2iγμζ ∂μϕ~,δζF~Φ=−2iζγμ∂μψ~.\begin{aligned} \delta_{\widetilde\zeta}\widetilde\phi&=-\sqrt2\widetilde\zeta\widetilde\psi,& \delta_{\widetilde\zeta}\widetilde\psi&=\sqrt2\widetilde\zeta\widetilde F_\Phi,& \delta_{\widetilde\zeta}\widetilde F_\Phi&=0,\\ \delta_\zeta\widetilde\phi&=0,& \delta_\zeta\widetilde\psi&=\sqrt2i\gamma^\mu\zeta\,\partial_\mu\widetilde\phi,& \delta_\zeta\widetilde F_\Phi&=-\sqrt2i\zeta\gamma^\mu\partial_\mu\widetilde\psi . \end{aligned}

Under δζ~\delta_{\widetilde\zeta}, the variation −2W~′′F~Φζ~ψ~-\sqrt2\widetilde W''\widetilde F_\Phi \widetilde\zeta\widetilde\psi from W~′F~Φ\widetilde W'\widetilde F_\Phi cancels +2W~′′F~Φζ~ψ~+\sqrt2\widetilde W''\widetilde F_\Phi \widetilde\zeta\widetilde\psi from the Yukawa term. Again the triple-fermion term vanishes. Under δζ\delta_\zeta the two terms combine as

δζF~W=−2i W~′ζγμ∂μψ~−2i W~′′(∂μϕ~)ζγμψ~=−2i ∂μ(ζγμW~′ψ~).\begin{aligned} \delta_\zeta\widetilde F_W &=-\sqrt2i\,\widetilde W'\zeta\gamma^\mu\partial_\mu\widetilde\psi -\sqrt2i\,\widetilde W''(\partial_\mu\widetilde\phi)\zeta\gamma^\mu\widetilde\psi\\ &=-\sqrt2i\,\partial_\mu \left(\zeta\gamma^\mu\widetilde W'\widetilde\psi\right). \end{aligned}

Together,

δLW=∂μBWμ,BWμ=−2i(ζ~γμW′ψ+ζγμW~′ψ~).\delta\mathcal L_W=\partial_\mu B_W^\mu,\qquad B_W^\mu=-\sqrt2i\left( \widetilde\zeta\gamma^\mu W'\psi +\zeta\gamma^\mu\widetilde W'\widetilde\psi \right).

The kinetic action already has the boundary expressions derived in the free-theory lesson. Add BWμB_W^\mu to them. Periodic fields on a flat torus, or sufficiently decaying fields on R3\mathbb R^3, give zero integrated boundary flux. This proves off-shell supersymmetry of the complete interacting action without requiring mm or gg to transform.

After eliminating the auxiliaries, the reduced action remains supersymmetric: the extra terms from varying their field-dependent solutions are multiplied by their vanishing algebraic Euler derivatives. The reduced transformation algebra, however, closes using the interacting equations of motion. Keeping an invariant action and keeping off-shell closure are distinct statements.

The algebraic saddle is not a reason to impose F~Φ=FΦ∗\widetilde F_\Phi=F_\Phi^* in the off-shell Euclidean action; that would leave a growing Gaussian from −F~ΦFΦ-\widetilde F_\Phi F_\Phi.

At fixed scalar, set P=W′=A+iBP=W'=A+iB, with real A,BA,B and W~′=P∗\widetilde W'=P^*. On the free auxiliary contour

FΦ=u+iv,F~Φ=−u+iv,u,v∈R,F_\Phi=u+iv,\qquad \widetilde F_\Phi=-u+iv,\qquad u,v\in\mathbb R,

the auxiliary Lagrangian becomes

Laux=u2+v2+2i(Bu+Av).\mathcal L_{\rm aux}=u^2+v^2+2i(Bu+Av).

Its real part is positive. For one finite-dimensional mode with positive Gaussian coefficient aa, direct Fourier Gaussian integration gives

∫R2 ⁣du dv e−a[u2+v2+2i(Bu+Av)]=πa e−a(A2+B2)=πa e−a∣W′∣2.\int_{\mathbb R^2}\!du\,dv\, e^{-a[u^2+v^2+2i(Bu+Av)]} =\frac\pi a\,e^{-a(A^2+B^2)} =\frac\pi a\,e^{-a|W'|^2}.

This independently checks the sign of the effective potential. Equivalently, deform to the saddle-centered contour FΦ=P∗+fF_\Phi=P^*+f, F~Φ=P−f∗\widetilde F_\Phi=P-f^*. The completed expression is then ∣f∣2+∣P∣2|f|^2+|P|^2. For this finite Gaussian the integrand is entire, and shifting the contours is justified by their Gaussian decay. The interacting saddle generally does not lie on the original contour F~Φ=−FΦ∗\widetilde F_\Phi=-F_\Phi^*.

This check concerns auxiliary integration with the scalar held fixed. A full continuum path integral additionally requires a regulator and a consistent prescription for the remaining fields; we have not calculated its partition function.

Supersymmetry and R-symmetry impose different tests

Section titled “Supersymmetry and R-symmetry impose different tests”

For fixed, neutral m,gm,g, a continuous phase R-symmetry of the form taught in the preceding lesson would require

2r=2,3r=2.2r=2,\qquad 3r=2.

There is no common solution. The flat component proof nevertheless went through. The supercharge algebra can act on an action that lacks this continuous R-symmetry.

With g=0g=0, the quadratic term permits r=1r=1. With m=0m=0, the cubic permits r=2/3r=2/3. Treating m,gm,g as transforming spurions can organize formal covariance, but it does not create a symmetry at fixed nonzero couplings. The dimensions [m]=1[m]=1, [g]=1/2[g]=1/2 are consistent with each term of WW having dimension two and each component term having dimension three.

This distinction matters for curved backgrounds constructed from a conserved R-current. The generic fixed-coupling quadratic-plus-cubic model cannot simply be placed into that R-symmetric construction using one phase assignment for ϕ\phi. Our result is the explicitly stated flat theory; a curved completion would require its own symmetry and background analysis.

1. Determine the Yukawa signs instead of remembering them

Section titled “1. Determine the Yukawa signs instead of remembering them”

Replace the two composite auxiliary expressions by

FW(c)=W′FΦ+cW′′ψψ,F~W(c~)=W~′F~Φ+c~ W~′′ψ~ψ~.F_W^{(c)}=W'F_\Phi+cW''\psi\psi,\qquad \widetilde F_W^{(\widetilde c)} =\widetilde W'\widetilde F_\Phi +\widetilde c\,\widetilde W''\widetilde\psi\widetilde\psi .

Use the same-family transformations to determine c,c~c,\widetilde c. Explain why choosing the same sign for both fails in these conventions.

First hint

The three-fermion term vanishes. Compare the coefficient of W′′FΦζψW''F_\Phi\zeta\psi and its tilded counterpart.

Second hint

The variation of the untilded W′W' has a plus sign, while that of W~′\widetilde W' has a minus sign. The variation of either spinor square contains twice its auxiliary-spinor contraction.

Complete solution

The two variations are

δζFW(c)=2(1+2c)W′′FΦζψ,δζ~F~W(c~)=2(−1+2c~)W~′′F~Φζ~ψ~.\begin{aligned} \delta_\zeta F_W^{(c)} &=\sqrt2(1+2c)W''F_\Phi\zeta\psi,\\ \delta_{\widetilde\zeta}\widetilde F_W^{(\widetilde c)} &=\sqrt2(-1+2\widetilde c) \widetilde W''\widetilde F_\Phi \widetilde\zeta\widetilde\psi . \end{aligned}

They vanish for arbitrary off-shell fields only if c=−1/2c=-1/2 and c~=+1/2\widetilde c=+1/2. This remains necessary for the quadratic example, where W′′=m≠0W''=m\ne0. For constant fields the incorrect residual contains no spacetime derivative and can have nonzero integral on a periodic box. A field redefinition could change a whole convention package; changing only one Yukawa sign does not.

2. A stationary point need not be a supersymmetric vacuum

Section titled “2. A stationary point need not be a supersymmetric vacuum”

Take real positive m,gm,g and write ϕ=x+iy\phi=x+iy, with x,yx,y real. Show that ϕ∗=−m/(2g)\phi_*=-m/(2g) is a stationary point of V=∣mϕ+gϕ2∣2V=|m\phi+g\phi^2|^2. Compute its energy and the quadratic terms in real displacements about it. Is it a stable vacuum? Does it preserve supersymmetry?

First hint

At a constant configuration, ∂ϕ∗V=W′(ϕ)W′′(ϕ)‾\partial_{\phi^*}V=W'(\phi)\overline{W''(\phi)}. Stationarity can follow from W′′=0W''=0 even when W′≠0W'\ne0.

Second hint

Write ϕ=−m/(2g)+ξ+iη\phi=-m/(2g)+\xi+i\eta. The derivative becomes W′=−m2/(4g)+g(ξ+iη)2W'=-m^2/(4g)+g(\xi+i\eta)^2.

Complete solution

At this point W′′=m+2gϕ∗=0W''=m+2g\phi_*=0, whereas W′(ϕ∗)=−m2/(4g)W'(\phi_*)=-m^2/(4g). Its energy is

V(ϕ∗)=m416g2>0.V(\phi_*)=\frac{m^4}{16g^2}>0.

The full expansion is

V=m416g2−m22ξ2+m22η2+g2(ξ2+η2)2.V=\frac{m^4}{16g^2} -\frac{m^2}{2}\xi^2+\frac{m^2}{2}\eta^2 +g^2(\xi^2+\eta^2)^2.

The Hessian in the coordinates (ξ,η)(\xi,\eta) has eigenvalues −m2,+m2-m^2,+m^2. The point is a saddle, with an unstable real displacement. Its nonzero W′W' also gives nonzero auxiliary expectation values, so it fails the constant-configuration supersymmetry test. The two zero-energy minima found in the main calculation are the stable classical supersymmetric vacua; the intervening stationary point is neither another minimum nor another supersymmetric vacuum.

3. Follow the vacua when a coupling vanishes

Section titled “3. Follow the vacua when a coupling vanishes”

For real positive couplings, compare two limits: g→0g\to0 at fixed m>0m>0, and m→0m\to0 at fixed g>0g>0. Track both roots, the quadratic fluctuation masses, and the limiting potential. In each case, explain whether the two distinct roots are lost through escape or merging, and whether the limiting theory is free.

First hint

Track the root −m/g-m/g before setting either coupling to zero. Roots can escape to infinity as well as merge.

Second hint

At m=0m=0, use V=g2∣ϕ∣4V=g^2|\phi|^4. A zero quadratic term does not mean the interaction has vanished.

Complete solution

For g→0g\to0 at fixed mm, the root at zero remains and the other root tends to negative infinity. In every fixed bounded region of field space,

V⟶m2∣ϕ∣2.V\longrightarrow m^2|\phi|^2.

The remaining vacuum has scalar and fermion mass mm. The distant vacuum had the same local mass for every nonzero gg, but it leaves the bounded region. The potential limit is not uniform over the entire unbounded field space.

For m→0m\to0 at fixed gg, both roots merge at zero and

V⟶g2∣ϕ∣4,W′′(0)=0.V\longrightarrow g^2|\phi|^4,\qquad W''(0)=0.

The quadratic masses vanish at the merged classical vacuum, but a quartic scalar interaction and a scalar-dependent Yukawa interaction remain. This is not a free theory. The first limit loses a root to infinity; the second merges two roots. These are different mechanisms. Neither classical calculation determines a quantum index or a quantum scaling dimension.

Given WW, you can now construct the Yukawa and auxiliary terms, verify their boundary variation, eliminate the auxiliaries, and distinguish F-term vacua from other stationary points. For the quadratic-plus-cubic example, the same derivative W′′W'' controls the Yukawa coefficient and the bosonic fluctuation mass around each supersymmetric vacuum.

Before extending this construction to gauge theories, explain why a charge test alone would not have fixed the Yukawa signs, and why off-shell invariance survives even though the generic fixed-coupling model has no continuous phase R-symmetry. Then construct the vector multiplet and test which matter interactions can consistently be gauged.