## Hope for SUSY?

Supersymmetry (SUSY), the symmetry between fermions, particles that form matter, and bosons, which mediate the interactions between them, has been proposed as one of the more likely extensions of the Standard Model of particle physics; however, it has so far received little experimental support. Condensed matter systems, such as the superfluid helium-3, may save the concept. In preparation for experimentation, **Grover et al.** (p. 280, published online 3 April) develop a theoretical approach that suggests SUSY describes the quantum phase transition on the boundary of a topological superconductor between a magnetic phase characterized by a bosonic order parameter and a neighboring phase hosting Majorana fermions.

## Abstract

In contrast to ordinary symmetries, supersymmetry (SUSY) interchanges bosons and fermions. Originally proposed as a symmetry of our universe, it still awaits experimental verification. Here, we theoretically show that SUSY emerges naturally in condensed matter systems known as topological superconductors. We argue that the quantum phase transitions at the boundary of topological superconductors in both two and three dimensions display SUSY when probed at long distances and times. Experimental consequences include exact relations between quantities measured in disparate experiments and, in some cases, exact knowledge of the universal critical exponents. The topological surface states themselves may be interpreted as arising from spontaneously broken SUSY, indicating a deep relation between topological phases and SUSY.

Since the 1970s, space-time “supersymmetry” (SUSY) has been actively pursued by particle physicists to attack the long-standing hierarchy problem of fundamental forces (*1*–*4*). Unlike any other symmetry, SUSY interchanges bosons and fermions and, when applied twice, generates translations of space and time, which ultimately leads to the conservation of momentum and energy (*3*). But despite sustained effort, SUSY has yet to be experimentally established in nature.

Here, we theoretically show that certain condensed matter systems display phenomena of emergent SUSY; that is, space-time SUSY naturally emerges as an accurate description of these systems at low energy and at long distances, although the microscopic ingredients are not supersymmetric. The physical systems we mainly consider are topological superconductors (TSCs) (*5*, *6*), in which pairs of fermions, which may be electrons or fermionic atoms such as helium-3, pair together in a special way. The resulting state has an energy gap to fermions in the bulk but gapless excitations at the surface. We consider the quantum phase transition at which the surface modes acquire a gap and establish emergent supersymmetry D = 1 + 1 and D = 2 + 1 dimensional surfaces (Fig. 1) by using a combination of numerical and analytical techniques (D denotes the space-time dimensionality).

The study of SUSY at a phase transition point was initiated in (*7*), where it was shown that the 1 + 1 dimensional tricritical Ising model, which can be accessed by tuning two parameters, is supersymmetric. A few other proposals that realize SUSY by fine-tuning two or more parameters have been made as well (*8*–*11*). Here, we require that SUSY be achievable by tuning only a single parameter, akin to a conventional quantum critical point (*12*). This is crucial for our results to be experimentally realizable. We also require that our theory has full space-time SUSY rather than only a limited “quantum-mechanical” SUSY (*13*). This automatically ensures translation invariance in space and time and will lead to experimental consequences, as we discuss below. Perhaps most interesting, in contrast to the strategy adopted in (*7*), our approach is not restricted to 1 + 1 dimensional theories.

There has been an explosion of activity in the field of topological phases since the discovery of topological insulators (TIs) (*14*–*17*). We will focus on a set of closely related phases, the time-reversal invariant TSCs (*5*, *6*), which include the well-known B-phase of superfluid helium-3 (*18*). These phases exist in both two and three spatial dimensions (*19*, *20*) and support Majorana modes at their boundary; the modes are protected by time-reversal symmetry from acquiring an energy gap. Spontaneous breaking of this symmetry provides a natural mechanism to gap them out. For example, electron-electron interactions at the surface could lead to magnetic order, which breaks time reversal symmetry. A natural question is how the surface modes evolve as the magnetic order sets in. We will see that space-time SUSY naturally emerges at the onset of magnetic order.

The D = 2 + 1 dimensional TSC, protected by the time-reversal symmetry, provides the simplest setting to address this question. Whereas the bulk of the superconductor is gapped, the boundary, a D_{edge} = 1 + 1 dimensional system, contains a pair of Majorana modes that propagate in opposite directions. The aforementioned instability of the edge may be described by introducing an Ising field that changes sign under time reversal. The action is given by
(1)with , , and we have used the conventional Dirac gamma matrices for the relativistic fermion . *v*_{ϕ} and *u* are, respectively, the velocity and self-interaction of ϕ; *g* is the coupling between the fermions and ϕ, whereas *r* is the tuning parameter for the transition. The symmetry-broken phase is characterized by 0, which leads to a mass gap for the fermions.

The mode count in the action is favorable for SUSY in D = 1 + 1 (*21*), with the bosons and Majorana fermions as superpartners. We now show that this is indeed the case by using a numerical simulation of a D = 1 + 1 lattice model that reproduces the action in Eq. 1 at low energies. The model is given by where(2)Here, is a single Majorana fermion at site , whereas the Ising spins sit on bond centers. When the transverse magnetic field , and lattice translation symmetry ensures that the Majorana fermions are gapless. As decreases, at some point orders antiferromagnetically, leading to a mass gap for the fermions through the coupling , reproducing the field theory in Eq. 1 at and near the critical point. tunes the relative bare velocities between the boson and the fermion modes, similar to in Eq. 1.

We numerically simulate a spin version of Eq. 2 by using the density matrix renormalization group (DMRG) method (*22*).(Fig. 2). At larger , a gapless phase is obtained that is separated by a critical line from a gapped ordered phase at small . To characterize the critical theory, consider crossing the phase boundary along fixed , say, , while monitoring the central charge , which quantifies the amount of entanglement of a 1 + 1 – D critical system (Fig. 2B). At small , *c* is almost equal to zero, which indicates a gapped phase. At , the central charge saturates at , indicating that the symmetry is restored and a gapless Majorana mode is present. At the transition, , we find , which precisely corresponds to the SUSY tricritical Ising theory (*7*). Note that we access this transition by tuning a single variable. Even when the parameter , which controls relative bare velocities of boson and fermion, is varied between , and , we obtain a second-order transition in the same universality class (*23*). We now turn to 3 + 1 dimensional TSCs, whose two-dimensional (2D) surface supports gapless Majorana fermions, protected by time-reversal symmetry. As advertised earlier, we study spontaneous breaking of time-reversal symmetry that gaps out the surface states. In the absence of numerics, we resort to an analytical calculation within a expansion to unravel the nature of quantum criticality and argue that 2 + 1 dimensional SUSY arises here. The relevant action is again given by Eq. 1 with . At , this action has SUSY (*21*) in 2 + 1D only if one sets and . We ask whether SUSY emerges in the low energy limit as we tune to the quantum critical point, without enforcing these conditions microscopically. We find (*23*) that this is indeed the case: The terms that break SUSY are irrelevant at the critical point leading to a stable superconformal fixed point in the infrared. Additional support has recently emerged from a rather different, conformal bootstrap approach (*24*) that obtains the scaling dimension , which matches well with our expansion result, .

To explore an experimental realization, consider a film of superfluid He_{3}-B, with surface perpendicular to the direction and a magnetic field parallel to the surface. Despite breaking time-reversal symmetry, this does not introduce an energy gap (*24*, *25*) because of a symmetry , that combines time reversal with 180° rotation about the axis. However, the system can lower its energy by spontaneously breaking rotation symmetry and gapping the surface modes. Indeed, the order parameter of He_{3}-B includes a rotation matrix that depends on an angle and a direction . In the absence of a magnetic field, microscopic dipolar interactions pin the vector . However, if the vector develops an in-plane component , it can break the protecting symmetry . The energetic incentive for an in-plane component increases with the applied magnetic field, because the surface gap *m* is proportional to the field: . Thus, for fields we have a gapless symmetric state, whereas for spontaneous symmetry breaking occurs, leading to a gap for surface Majorana modes. In (*24*), the critical field was predicted to be G, which is readily accessible. However, we note the challenges to realizing SUSY in this physical system (*23*). The low temperatures of the superfluid transition restricts the temperature window for experiments, and the bare boson and fermion velocities are expected to be rather different, with a ratio controlled by .

We now discuss the possibility of emergent SUSY in 3D TIs. The surface states of these materials (*14*–*17*) consist of Dirac fermions at a chemical potential . Let us fine-tune to zero and consider the instability of these surface modes to an s-wave superconductor. This multicritical point can potentially be driven by intrinsic interactions and may also be realized by patterning the surface with a Josephson junction array. We restrict ourselves to particle-hole symmetric superconducting fluctuations. The effective near the transition is given by
(3)where is the superconducting order parameter, the fermion field , and is the 2D antisymmetric tensor. *g*, *r*, and *u* have meaning analogous to Eq. 1, whereas *c* is the velocity of . The above action is known to flow to a superconformal fixed point (*26*, *27*).

Because SUSY rotates fermions into bosons, the universal critical exponents satisfy nontrivial relations between them, thus relating results from disparate experiments. In particular, the scaling dimensions of satisfy . In an electronic system, may be determined via tunneling experiments, whereas may be obtained from neutron scattering. As an example, consider the boundary of a 2D TSC, for which we showed that the transition lies in the tricritical Ising universality class. The exact critical exponents in this case are . Similarly, at the boundary of 3D TSC, we find at order in the renormalization group. Lastly, in the case of 3D TI, because of the larger () SUSY, one can determine the scaling dimensions of the boson and Fermi fields exactly, with (*28*) . Another consequence of SUSY is that the velocity for the fermion and the boson at low energies become equal, despite the fact that their bare velocities will be generically different.

SUSY has nontrivial consequences for the proximate phases as well. At the TSC boundary in the time-reversal broken phase, SUSY implies that the masses of the fermion and the boson are equal to each other as one moves away from the critical point. Mathematically, *m*_{ϕ} and *m*_{χ} are equal to α(*r* – *r _{c}*)

*, where is the critical value of the tuning parameter, α is a constant, and, owing to SUSY, .*

^{v}Emergent SUSY also provides insight into the topological excitations on either side of the transition. First, consider the time-reversal invariant boundary of a TSC in either two or three spatial dimensions, which supports a gapless Majorana fermion. This is conventionally viewed as resulting from the topological band structure. However, a completely different viewpoint is obtained by realizing that the flow from SUSY critical point () to a phase where is gapped () corresponds to spontaneous breaking of SUSY. Because SUSY is a fermionic symmetry, this leads to the generation of a massless Majorana fermion, the “goldstino” mode (*23*, *29*, *30*), which precisely corresponds to the topologically protected Majorana fermion. Thus, one can model the boundary modes with a continuum model in the same dimension if SUSY is assumed, which serves to protect their gaplessness. We provide details of this phenomenon in (*23*).

The dichotomy between SUSY and topology runs deeper: It has robust, physical manifestation in the time-reversal broken phase as well. As discussed in (*23*) and hinted above, the SUSY is unbroken as the boson condenses, resulting in equal mass for and . It is well known that, in the symmetry broken boundary phase, topological defects—kinks in 1D and line domain walls in 2D—carry fermion modes on them (*31*). We interpret these also as Goldstino modes of a type of SUSY breaking. First consider kinks, which spontaneously break translational symmetry in 1D. Because of the deep connection between translation symmetry and SUSY, this ultimately leads to a spontaneous breaking of half of the SUSY (*23*, *32*, *33*). Thus, in addition to the coordinate of the domain wall, , which can be viewed as a gapless 0 + 1 – D Goldstone boson of translation breaking, there is also a single Majorana Goldstino mode localized on the kink, which is the superpartner of . The low-energy theory of the kink is , which one might have anticipated independently on topological grounds (*31*). This argument also applies to TSCs coupled via Josephson junction discussed in (*34*).

The generalization to line domain walls at the time-reversal breaking boundary of 3D TSC is straightforward. The partial SUSY breaking leads to 1 + 1 – D chiral Majorana fermion propagating along the domain line; is a superpartner of the Goldstone mode associated with the translational symmetry breaking, . SUSY not only predicts the anticipated chiral Majorana fermions on the domain walls, it also predicts equal velocity for the superpartners and an exact “Fermi-Bose” degeneracy associated with the whole excitation spectra of the domain wall or kink. These predictions hold in the ordered phase close to the criticality. In general, we expect that the boundary modes in a topological phase can be thought of as Goldstinos unless the tuning parameter breaks SUSY explicitly; an example of the latter possibility is provided by the SUSY critical point in 3D TIs mentioned above.

## Supplementary Materials

## References and Notes

- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
Quantum mechanical SUSY corresponds to SUSY only in the time-direction; see, e.g., (
*35*) and (*10*). - ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
Supplementary materials are available on
*Science*Online. - ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
A slightly different theory, which might also exhibit emergent SUSY, arises in the context of time-reversal breaking 1 + 1 – D topological superconductors coupled via a Josephson junction (
*36*). - ↵
- ↵
- ↵
- ↵
**Acknowledgments:**We thank M. Fisher, T. Mizushima. D. Huse, R. Mong, and S. Trivedi for discussions. This work is supported by Army Research Office Multidisciplinary University Research Initiative grant W911-NF-12-0461 and NSF Division of Materials Research grant 0645691 (A.V.) and by U.S. Department of Energy Office of Basic Energy Sciences under grant DE-FG02-06ER46305 (D.N.S.).