# Topological superfluid He-B in magnetic field and Ising variable

###### Abstract

The topological superfluid He-B provides many examples of the interplay of symmetry and topology. Here we consider the effect of magnetic field on topological properties of He-B. Magnetic field violates the time reversal symmetry. As a result, the topological invariant supported by this symmetry ceases to exist; and thus the gapless fermions on the surface of He-B are not protected any more by topology: they become fully gapped. Nevertheless, if perturbation of symmetry is small, the surface fermions remain relativistic with mass proportional to symmetry violating perturbation – magnetic field. The He-B symmetry gives rise to the Ising variable , which emerges in magnetic field and which characterizes the states of the surface of He-B. This variable also determines the sign of the mass term of surface fermions and the topological invariant describing their effective Hamiltonian. The line on the surface, which separates the surface domains with different , contains gapless fermions, which are protected by combined action of symmetry and topology.

## 1. Introduction

Recently topological insulators, semimetals, superconductors, superfluids and other topologically nontrivial gapless and gapped phases of matter have attracted a lot of attention, see e.g. [1, 2, 3, 4]. Superfluid He-B represents the fully gapped topological superfluids with time reversal symmetry [3, 4, 5]. This system contains gapless Andreev bound states on the surface of the liquid or at the interfaces between the bulk states with different topological invariants; these fermion zero modes resemble Majorana fermions in relativistic quantum field theories [6, 7, 8]. Andreev bound states on the surface of He-B were discussed theoretically (see e.g. [9, 10, 11]) and probed experimentally [12, 13, 14]. However, the Majorana signature of these states has not yet been observed. One of the possible tools for observation of the “relativistic” spectrum of the bound states is NMR, which requires an external magnetic field. The influence of magnetic field on Majorana fermions has been recently discussed in [11]. Here we consider the effect of magnetic field in more detail.

The topological invariant for He-B which gives rise to the gapless fermion zero modes on its surface is supported by the time reversal symmetry [5]. Magnetic field destroys this symmetry, as a result the Andreev bound states acquire gap (mass) proportional to the magnitude of magnetic field. The resulting massive Andreev-Majorana fermions also have nontrivial topology, described by the nonzero topological invariant in momentum space. The value of this invariant is determined by the Ising type variable , which characterizes the surface of He-B in the presence of magnetic field. The surface of He-B can be in two states, with and , and both states of the surface have been recently probed by NMR experiments with He-B confined in a single micron-size slab [15]. These experiments also detected the coexistence of the domains with different on the surface of He-B; these domains are separated by the one-dimemsional domain walls (lines) on the surface. Since the domains are characterized by different values of the topological invariant, the boundary line between the domains contains gapless fermions with “relativistic” spectrum.

## 2. Distortion of order parameter

Magnetic field destroys the rotational symmetry of He-B, leading to the anisotropic gap. The gap anisotropy is also induced by the walls. The order parameter in He-B is matrix [16]. It transforms as vector under spin rotation (the Greek index), and as vector under orbital rotation (the Latin index). For the distorted He-B it has the form

(0) |

where is the distortion parameter with in the non-deformed B-phase and in the planar state; is the gap for quasiparticles propagating in direction along the normal to the wall; is the gap for quasiparticles propagating in directions parallel to the wall; is rotation matrix, which connects spin and orbital spaces. As said, the distortion of the gap is caused both by magnetic field and by boundaries. Let us introduce the unit vector in direction of magnetic field, and the unit vector in the orbital space. The orbital vector shows the direction of the orbital angular momentum of Cooper pairs (see e.g. [17, 18]) and plays the role of the anisotropy axis for the uniaxial deformation of the gap.

## 3. Ising variable probed by NMR

Let us consider NMR experiments with superfluid He-B confined in a single micron size slab [15], where the field is normal to the plates: . Typical fields used for NMR exceed the dipole field (spin-orbit coupling). As a result the boundary conditions for -vector is that is normal to the walls, [19]. This boundary condition ensures that the anisotropy axis caused by the walls and anisotropy axis caused by magnetic field coincide.

In a given geometry with the field being normal to the plates, the projection of the -vector on the magnetic field direction

(0) |

represents the Ising variable. Thus is because only two orientations of -vector with respect to manetic field are allowed by boundary conditions: , and . The degeneracy is lifted by spin-orbit interaction. The latter is small compared with all other energies but is important for the measured NMR frequency shift from the Larmor value , which is just caused by spin-orbit interaction. The measurements of the frequency shift allows to determine experimentally the value of the Ising spin on the surface.

For the distorted He-B, the corresponding frequency shifts can be found from equations (4.5) and (4.10) of Ref. [20]. The frequency shift is positive when -vector is parallel to the field, i.e. for , and negative for -vector antiparallel to field, i.e. for . The ratio of the frequency shifts is

(0) |

Both positive and negative frequency shifts have been observed in [15]. From Fig. 1 of [15] and Eq.(3.) it follows that in the geometry of experiment, the distortion parameter at temperature is about . From this figure it also follows that in cooling one obtains three different states: the state with characterized by small positive frequency shift; the state with characterized by large negative frequency shift; and the state with coexisting domains of both orientations of the Ising variable, which is characterized by two peaks in the NMR spectrum. The latter state contains domain wall(s) – the line(s) on the surface separating surface domains with opposite orientations of .

Let us also for completeness consider the field parallel to the plates. In this case the -vector is perpendicular to the field and . The frequency shift is positive, .

We have seen that the orientation of the magnetic field with respect to the Ising spin is crucial for the NMR experiments. Let us consider the effect of orientation of the magnetic field on the spectrum of Majorana fermions.

## 4. Majorana fermions in distorted B-phase

Fermions in magnetic field are described by the following Hamiltonian:

(0) | |||

(0) |

where are Pauli matrices of Bogolyubov-Nambu spin; are Pauli matrices of He nuclear spin; is the gyromagnetic ratio of the He atom; and finally are the Pauli matrices obtained from matrices by rotation. In the absence of magnetic field there is a topological invariant for He-B [5]:

(0) |

At , when the time reversal symmetry is preserved, the integral (4.) is invariant under deformations because the matrix anti-commutes with Hamiltonian . The nonzero value of this invariant supports gapless fermions at the interface between bulk states with different values of [6] and on the surface of He-B [8]. These fermions contribute to the ground-state spin supercurrent [6].

In the presence of magnetic field, the integral (4.) does not represent the topological invariant because the Pauli term does not anti-commute with , and the gaplessness of the Andreev bound states is not guaranteed. The Pauli term can be expressed in terms of the rotated spin :

(0) |

Equation (4.) demonstrates that the Hamiltonian (4.) does not depend on the orientation of magnetic field. It depends only on the Ising variable , which describes two possible orientations of the vector on the surface. This is because in large fields used in experiments, the spin-orbit interaction can be neglected and the boundary conditions for the vector on the surface are fully determined by the orbital effects, which are invariant under rotation of the magnetic field together with the spin subsystem.

Let us, however, consider the orientation of magnetic field along the axis . This allows us to express the Hamiltonian in terms of the variable in (3.), which is regulated in NMR experiments in Ref. [15] where the orientation of the field normal to the plates is used:

(0) | |||

(0) |

The second order secular equation (10) in [8] produces the following Hamiltonian for fermion zero modes:

(0) |

which can be written as

(0) |

This the so-called helical (see e.g. [21] and references therein) Hamiltonian produces the relativistic spectrum of massive particles:

(0) |

where the ‘speed of light’ [8]

(0) |

The mass of the surface fermions

(0) |

appears due to violation of the time reversal symmetry by magnetic field. The mass (4.) was obtained in Ref. [11] only for the case of magnetic field normal to the wall. However, as distinct from Ref. [11], in our case the equation (4.) is valid for any orientation of magnetic field. This is because we consider relatively large magnetic fields used in NMR experiments, when the spin-orbit interaction is relatively small: though the spin-orbit interaction is responsible for the frequency shift, it does not influence the orientation of the vector on the surface. As a result the Pauli term (4.) expressed via the rotated spin does not depend on orientation of magnetic field.

This is contrary to the effect of magnetic field discussed in Refs. [7, 11], where the magnetic field is assumed to be so small that the boundary conditions are solely determined by spin-orbit interaction. This leads to anisotropy of magnetic susceptibility [7, 11], which however disappears in large fields. Note also that in our case the spin which enters the Hamiltonian (4.) for Majorana fermions, is obtained from the real spin by rotations provided by the order parameter matrix . This means that in large fields, when the spin-orbit interaction can be neglected, the quantization axis for the real spin of Majorana fermions is not fixed.

The effective Hamiltonian (4.) has nontrivial topology described by the topological invariant

(0) |

where is the Green’s function of fermion zero modes. The invariant (4.) was first introduced in relativistic theories [22, 23, 24]. It is responsible for quantization of Hall conductivity induced by topological structure of energy-momentum space, the so-called intrinsic quantum Hall effect. It was later independently introduced for the film of superfluid He-A in condensed matter [25, 26]. Under experimental conditions in Ref. [15], the value of this invariant is determined by the Ising variable in (3.):

(0) |

NMR experiments in Ref. [15] detected coexistence of the surface domains with different , which are thus separated by the one-dimensional domain wall on the surface. Since the domains have different values of topological invariant , then according to the index theorem (see e.g. [27]), the line separating two domains contains gapless fermion zero modes.

## 5. Discussion

The symmetry of superfluid He-B is important for its topological properties. Magnetic field violates isotropy of He-B: the gap become anisotropic with . This distortion of the gap does not change the topological invariant introduced in [5], and thus cannot destroy the massless relativistic fermions on the surface of He-B. However, the violation of the other symmetry of He-B by magnetic field – the time reversal symmetry – is crucial. The topological invariant supported by this symmetry ceases to exist and thus the gapless fermions on the surface of He-B are not protected any more by topology: they become fully gapped. Nevertheless, if perturbation of symmetry is small, the surface fermions remain relativistic with mass proportional to symmetry violating perturbation – magnetic field.

This is not the whole story. The spontaneously broken relative spin-orbit symmetry of He-B [28] – symmetry under relative rotations of spin subsystem with the respect to the orbital one [29] – gives rise to the vector of orbital anisotropy in the presence of magnetic field. This in turn leads to Ising variable which characterizes the surface of He-B. The Ising variable determines the sign of the mass term of the surface fermions and the topological invariant (4.) describing their effective Hamiltonian (4.). The line on the surface, which separates domains with different , contains gapless fermions, which are protected by combined action of symmetry and topology. So, the superfluid He-B provides many examples of the interplay of symmetry and topology.

The other examples of the interplay of symmetry and topology in topological matter and in vacua of relativistic quantum fields can be found in [27, 32, 33]. Some topological invariants exists only due to symmetry, while some values of the topological invariant are not compatible with symmetry. It is right time for the general classification of topological matter in terms of its symmetry and topology. The attempt of such classification was made in [1, 2, 3, 4]. However, only non-interacting systems were considered, and not all the symmetries were exploited. Examples are again provided by the superfluid phases of He, where the symmetry is enhanced due to relative smallness of the spin-orbit interactions. One task should be to consider symmetry classes, including crystal symmetry classes, magnetic classes, superconductivity classes, etc., and to find out what are the topological classes of Green’s function which are allowed within a given symmetry class. Then one should find out what happens when the symmetry is smoothly violated, etc. The Green’s function matrices with spin and band indices must be used for this classification, since they take into account interaction. In this way we may finally obtain the full classification of topological matter, including insulators, superconductors, magnets, liquids, and vacua of relativistic quantum fields.

The problem of the interplay of topology and symmetry exists not only for the momentum-space topology, but also for the real-space topology, which considers classification of topological defects and textures in systems with broken symmetry. For example, some values of the winding number describing the topological defect are not compatible with symmetry of this defect, which connects for example the points and . This happens with monopoles [34]; disclinations – topological defects in liquid crystals [35]; vortices in superfluids [17]; cosmic strings; etc. One of the tasks there is the classification of topological objects in terms of their symmetry. Similarly we need the symmetry classification of topological objects in momentum space, where the symmetry connects different points in momentum space, say and [1].

The special role in classification of topological systems is played by dimensional reduction. See for example Chapters 11 and 21 in [27] where the topological invariants for the 2D fully gapped systems are obtained by dimensional reduction of the topological invariants for the 3D gapless system; the dimensional reduction has been recently discussed in [36]. The dimensional reduction allows us to use for classification of the gapped systems the scheme, which was suggested by Horava for the classification of the topologically nontrivial nodes in spectrum [37]. But this also must be supplemented by the symmetry analysis, and also by analysis of the types of the emergent symmetries which necessarily appear within some topological classes.

I owe special thanks to John Saunders for discussion on the experiments [15] and am grateful to Shinsei Ryu and Liang Fu for discussion on symmetry and topology. This work is supported in part by the Academy of Finland, Centers of Excellence Program 2006–2011, and the Khalatnikov–Starobinsky leading scientific school (grant 4899.2008.2).

- 1. L. Fu, C.L. Kane and E.J. Mele, Topological insulators in three dimensions, Phys. Rev. Lett. 98, 106803 (2007).
- 2. J. E. Moore and L. Balents Topological invariants of time-reversal-invariant band structures, Phys. Rev. B 75, 121306 (2007).
- 3. A.P. Schnyder, S. Ryu, A. Furusaki and A.W.W. Ludwig, Classification of topological insulators and superconductors, AIP Conf. Proc. 1134, 10–21 (2009); arXiv:0905.2029.
- 4. A. Kitaev, Periodic table for topological insulators and superconductors, AIP Conf. Proc. 1134, pp. 22–30 (2009); arXiv:0901.2686.
- 5. G.E. Volovik, Topological invariant for superfluid He-B and quantum phase transitions, Pis’ma ZhETF 90, 639–643 (2009); JETP Lett. 90, 587–591 (2009); arXiv:0909.3084.
- 6. M.M. Salomaa and G.E. Volovik, Cosmiclike domain walls in superfluid He-B: Instantons and diabolical points in (,) space, Phys. Rev. B 37, 9298–9311 (1988).
- 7. Suk Bum Chung, Shou-Cheng Zhang, Detecting the Majorana fermion surface state of He-B through spin relaxation, Phys. Rev. Lett. 103, 235301 (2009).
- 8. G.E. Volovik, Fermion zero modes at the boundary of superfluid He-B, JETP Lett. 90, 398–401 (2009); arXiv:0907.5389.
- 9. N.B. Kopnin, P.I. Soininen and M.M. Salomaa, Parameter-free quasiclassical boundary conditions for superfluid He at rough walls. B-phase order parameter and density of states, J. Low Temp. Phys. 85, 267–282 (1991).
- 10. K. Nagai, Y. Nagato, M. Yamamoto and S. Higashitani, Surface bound states in superfluid He, J. Phys. Soc. Jap. 77, 111003 (2008).
- 11. Y. Nagato, S. Higashitani and K. Nagai, Strong anisotropy in spin suceptibility of superfluid He-3-B film caused by surface bound states, J. Phys. Soc. Japan 78, 123603 (2009).
- 12. C.A.M. Castelijns, K.F. Coates, A.M. Guénault, S.G. Mussett and G.R. Pickett, Landau critical velocity for a macroscopic object moving in superfluid He-B: evidence for gap suppression at a moving surface, Phys. Rev. Lett. 56, 69–72 (1986).
- 13. J.P. Davis, J. Pollanen, H. Choi, J.A. Sauls, W.P. Halperin and A.B. Vorontsov, Anomalous attenuation of transverse sound in He, Phys. Rev. Lett. 101, 085301 (2008).
- 14. S. Murakawa, Y. Tamura, Y. Wada, M. Wasai, M. Saitoh, Y. Aoki, R. Nomura, Y. Okuda, Y. Nagato, M. Yamamoto, S. Higashitani and K. Nagai, New anomaly in transverse acoustic impedance of superfluid He-B with a wall coated by several layers of He, Phys. Rev. Lett. 103, 155301 (2009).
- 15. L.V. Levitin, R.G. Bennett, A.J. Casey, B. Cowan, J. Parpia and J. Saunders, Superfluid He confined in a single 0.6 micron slab: A phase transition between superfluid phases with hysteresis, J. Low Temp. Phys. 158, 159–162 (2010).
- 16. D. Vollhardt and O. Wölfle, The superfluid phases of helium 3, Taylor and Francis, London (1990).
- 17. M.M. Salomaa and G.E. Volovik, Quantized vortices in superfluid He, Rev. Mod. Phys. 59, 533–613 (1987).
- 18. V.V. Dmitriev, V.B. Eltsov, M. Krusius, J.J. Ruohio and G.E. Volovik, Superflow-stabilized nonlinear NMR in rotating He-B, Phys. Rev. B 59, 165–168 (1999).
- 19. W.F. Brinkman and M.C. Cross, Spin and orbital dynamics of superfluid He, in: Progress in Low Temperature Physics, edited by D.F. Brewer (North-Holland, Amsterdam), Vol. VIIA, p. 105 (1978).
- 20. Yu.M. Bunkov, and G.E. Volovik, Homogeneously precessing domains in He-B, JETP 76, 794–801 (1993).
- 21. Jiadong Zang and Naoto Nagaosa, Monopole current and unconventional Hall response on topological insulator, arXiv:1001.1578.
- 22. K. Ishikawa and T. Matsuyama, Magnetic field induced multi component QED in three-dimensions and quantum Hall effect, Z. Phys. C 33, 41–45 (1986).
- 23. K. Ishikawa and T. Matsuyama, A microscopic theory of the quantum Hall effect, Nucl. Phys. B 280, 523–548 (1987).
- 24. T. Matsuyama, Quantization of conductivity induced by topological structure of energy-momentum space in generalized QED, Progr. Theor. Phys. 77, 711–730 (1987).
- 25. G.E. Volovik, Analog of quantum Hall effect in superfluid He film, JETP 67, 1804–1811 (1988).
- 26. G.E. Volovik and V.M. Yakovenko, Fractional charge, spin and statistics of solitons in superfluid He film, J. Phys.: Condens. Matter 1, 5263–5274 (1989).
- 27. G.E. Volovik, The Universe in a Helium Droplet, Clarendon Press, Oxford (2003).
- 28. A.J. Leggett, NMR lineshifts and spontaneously broken spin-orbit symmetry: I. General concepts, J. Phys. C 6, 3187–3204 (1973).
- 29. In relativistic theories the 4D analog of the broken relative spin-orbit symmetry is the broken relative Lorentz symmetry, see [30, 31].
- 30. G.E. Volovik, Superfluid He-B and gravity, Physica B 162, 222–230 (1990).
- 31. C. Wetterich, Gravity from spinors, Phys. Rev. D 70, 105004 (2004).
- 32. G.E. Volovik, Quantum phase transitions from topology in momentum space, in: ”Quantum Analogues: From Phase Transitions to Black Holes and Cosmology”, eds. W.G. Unruh and R. Schützhold, Springer Lecture Notes in Physics 718 (2007), pp. 31–73; cond-mat/0601372.
- 33. G.E. Volovik, Topological invariants for Standard Model: from semi-metal to topological insulator, Pis’ma ZhETF 91, 61–67 (2010); arXiv:0912.0502.
- 34. D. Wilkinson and A.S. Goldhaber, Spherically symmetric monopoles, Phys. Rev. D 16, 1221–1231 (1977).
- 35. A.A. Balinskii, G.E. Volovik and E.I. Kats, Disclination symmetry in uniaxial and biaxial nematic liquid crystals, ZhETF 87, 1305–1314 (1984); JETP 60, 748–753 (1984).
- 36. S. Ryu, A. Schnyder, A. Furusaki and A. Ludwig, Topological insulators and superconductors: ten-fold way and dimensional hierarchy, arXiv:0912.2157.
- 37. P. Hořava, Stability of Fermi surfaces and -theory, Phys. Rev. Lett. 95, 016405 (2005).