Bulk Topological Superconductors
DEPARTMENT OF PHYSICS, GRADUATE SCHOOL OF SCIENCE, KYOTO UNIVERSITY, JAPAN
This review introduces known candidates for bulk topological superconductors and categorizes them with time-reversal symmetry (TRS) and gap structures. Recent studies on two archetypal topological superconductors, TRS-broken Sr2RuO4 and TRS-preserved CuxBi2Se3, are described in some detail.
In recent years, research on topological materials has been explosively growing. This explosion was triggered by the prediction and discovery of quantum spin-Hall systems and topological insulators (TIs), with non-trivial topology of the valence-band wave function in the reciprocal space. Similarly, it has been recognized that some fully-gapped superconductors also possess non-trivial topological properties in their superconducting (SC) wave functions[1-6]. Such superconductors are called topological superconductors (TSCs).
One of the most intriguing properties of gapped topological materials is bulk-edge correspondence. Because it is by definition impossible to change a topologically non-trivial electronic state to a trivial state by a continuous deformation, there must be a certain kind of "singularity" at the boundaries between topologically distinct materials. Such considerations explain the closing of the gap, i.e. existence of gapless states, on the boundary between a gapped topological material and a gapped non-topological material (including vacuums).
Fully-gapped topological superconductivity is also accompanied by gapless surface states. Furthermore, because of the particle-hole symmetry in the low-energy Hamiltonian, surface states of a TSC should have the Majorana nature, i.e., particles being identical to their own antiparticles. Intriguingly, the surface state at zero energy (i.e., at the chemical potential), called the Majorana zero modes, do not obey ordinary (Fermionic or Bosonic) statistics in two dimensions. Instead, the Majorana zero modes exhibit non-Abelian statistics, which can be potentially utilized for quantum computing.
Although the concept of topology was originally considered for systems with full energy gaps, it has recently expanded to gapless (i.e. gaps with nodes) states; one can evaluate a "topological invariant of a gap node" on a manifold surrounding the node. Similarly, on a certain cut of the Brillouin zone, a k-dependent topological number can be defined. Application of this extension to superconductivity gives us new insights on unconventional superconductors.
TOPOLOGICAL SUPERCONDUCTIVITY-CANDIDATE MATERIALS
Within the conservative definition based on global topological numbers in fully gapped systems, materials exhibiting topological superconductivity are rather scarce. Nevertheless, there are several known candidates for bulk TSCs and superfluids as listed in Table I. With the recent expansion of the topology concept to nodal superconductors, a number of nodal unconventional superconductors studied over 35 years now can be added to the list, although topological phenomena in most of the new members have not been experimentally explored. Gap structures of representative topological SC and superfluid states are schematically shown in Fig. 1.
Table 1: List of known candidate materials for bulk topological superconductors and superfluids.
Fig. 1: Schematics of representative bulk topological superconducting and superfluid states. The color depicts the superconducting gap amplitude and phase as shown in the left bottom panel, with black corresponding to gap nodes, except for the BW state (corresponding to the B phase of 3He), where the color represents the d-vector direction. For all structures, cross sectional views are also shown, in which the gray tubes or spheres indicate the normal-state Fermi surface. For the parity-mixed p + s state, positive and negative helicity gaps are shown.
Prototypical superfluid 3He phases
The most classical and established cases are the superfluid phases in 3He. Both the A and B phases of the superfluid 3He have spin-triplet odd-parity Cooper pairs. The former phase, the 3D chiral p-wave state, spontaneously breaks TRS and is topologically non-trivial because of topological monopole-like point nodes (Weyl nodes). The B phase, in which TRS is preserved, is also topologically non-trivial, because of the hedgehog-like spin-orbital texture on its Fermi surface. The 3He phases have been providing prototypical playgrounds to investigate properties of topological pairings[9,10].
Classification based on TRS and gap nodes
Here, I classify bulk topological pairing states based on accompanying TRS breaking as shown in Table I. Each class is then divided into states with full gaps and those with nodal gaps. I believe that this way is convenient to grasp the current experimental situations, because TRS breaking and gap structures can be clearly detectable by experiments and indeed, have been investigated in most of the candidate materials.
TRS-broken topological superconductivity
Topological SC states with spontaneous TRS breaking due to orbital phase winding can be termed as "chiral" states. The layered perovskite oxide Sr2RuO4 is a leading candidate for quasi-two-dimensional (Q2D) chiral-p-wave superconductors[11-14]. The current research situation on this superconductor is described in §3. In the recently-discovered superconductor SrPtAs with a hexagonal honeycomb structure, spontaneous TRS breaking was reported based on μSR experiments, suggesting possible chiral-d-wave states. A fully-gapped state is indicated by the NMR study. However, the observed coherence peak in 1/T1 may favor s-wave superconductivity.
Among the nodal superconductors, two heavy-fermion superconductors, UPt3 and URu2Si2, are leading candidates for TRS-broken topological superconductivity. TRS breaking in the B-phase of UPt3 had been indicated 23 years ago by a μSR experiment but with controversy, and has been recently confirmed by the Kerr effect study. Since the odd-parity nature of superconductivity in this compound has been revealed by using the NMR technique and polarized neutron scattering, the pairing in UPt3 can be most promisingly interpreted as the chiral-f-wave kz(kx ± iky)2 or (5kz2 - 1)(kx ± iky) states (see Fig. 1), which are accompanied by topological Weyl nodes[22,23]. We note that some fundamental issues such as existence of tetra critical point have not been explained within the scenario of chiral states. The heavy fermion URu2Si2 also exhibits TRS breaking in the SC state as evidenced by the Kerr effect. Another signature for the chiral state is also obtained by the Nernst effect, which is found to be exceptionally huge just above Tc of this compound and is attributed to superconducting fluctuation with chiral nature. Together with the Pauli-limit-like Hc2 behavior, it has been believed that superconductivity of this compound is of a 3D chiral-d-wave kz(kx ± iky) (= dzx ± idyz) state (Fig. 1) with Weyl nodal points in its SC gap, although direct evidence for the spin-singlet SC nature is actually still controversial[28,29]. More recently, the half-Heusler compounds YPtBi and RPdBi (R = rare earth) with non-centro-symmetric crystal structures were found to be superconducting possibly with a line-nodal gap[30-32]. Noble "j = 3/2" pairing states with TRS breaking are theoretically proposed, although experimental evidence for the TRS breaking has not been reported yet.
TRS-preserved topological superconductivity
In 2010, Hor et al. discovered that the topological insulator Bi2Se3 exhibits superconductivity below ~ 3.8 K when Cu ions are intercalated. Theories and some experiments have indicated fully-gapped, topological, and odd-parity superconductivity in CuxBi2Se3[3,4], as described in more detail in §4. The doped topological crystalline insulator (TCI) Sn1-xInxTe was actually known to exhibit superconductivity below around 2 K, far before SnTe was recognized as a TCI. Point-contact spectroscopy in the SC state revealed the possibility of unusual surface states. A fully-opening gap and preserved TRS are indicted by specific-heat and μSR studies. These results satisfy a condition for TRS-preserved fully-gapped topological SC states.
It should be noted that two heavy-fermion superconductors, namely the spin-singlet superconductor CeCu2Si2 and the spin-triplet superconductor UBe13 (not listed in Table I), may have fully-opened gaps according to recent specific-heat studies[39,40], although they had been believed to be nodal superconductors. Although the gap functions are not clear, the possibility of fully-gapped topological superconductivity is worth investigating.
Candidates for TRS-preserved nodal TSCs are also listed in Table I. It has been established that high-Tc cuprate superconductors exhibit Q2D dx2-y2 superconductivity (Fig. 1), and the nodes of those superconductors are protected due to a topological reason. Indeed, as early as 1995, surface states originating from the SC gap had been detected as a sharp zero-bias peak in tunneling spectroscopy, due to flat-band surface states originating from the non-trivial k-dependent topology. Similarly, Q2D or Q1D d-wave or d-wave-like superconductors such as CeTIn5 (T = Co, Rh, Ir), κ-(ET)2X (X = Cu(NCS)2 etc.), and (TMTSF)2X (X = ClO4, etc.) can be regarded as topological nodal superconductors. Recently, the Cu-doped (PbSe)5(Bi2Se3)6, a "naturally-made" superlattice between a TI and an ordinary insulator, is found to become a superconductor with Tc ~ 3 K, probably with a line nodal SC gap.
Some superconductors with non-centro-symmetric crystal structures (i.e., structures without inversion symmetry) are also candidates for topological nodal superconductors. In such superconductors, mixing of spin-triplet and spin-singlet order parameters are possible because of the lack of the parity in their crystal structures, and gap nodes emerge when the spin-triplet component is dominant (see Fig. 1). The tetragonal CePt3Si and CeTSi3 (T = Rh and Ir) exhibit large upper critical fields attributed to parity-mixed superconducting states[49-51]. In the cubic Li2Pt3B, spin-triplet-dominant and nodal superconductivity is reported based on NMR studies.
Ferromagnetic superconductivity with broken TRS
In the present classification, another interesting class is the ferromagnetic superconductors UGe2, URhGe, and UCoGe, in which the TRS is already broken in the normal state due to the ferromagnetic ordering. In these superconductors, both superconductivity and ferromagnetism originate from electrons in the same (probably uranium) bands. Judging from this, as well as from the large upper critical fields, non-unitary (spin-polarized) spin-triplet states are likely to be realized. Although the gap structures are not yet clearly known, non-trivial topology is highly anticipated.
TRS-BROKEN TOPOLOGICAL SUPERCONDUCTIVITY IN Sr2RuO4
In this section, I review the current situation of research on superconductivity in Sr2RuO4, one of the leading candidates for TRS-broken TSCs[11-14]. In particular, key experiments reported in recent years are described in the last half of this section.
Evidence for the chiral-p-wave SC state
Sr2RuO4 has a layered perovskite structure consisting of conductive Ru-O layers and insulating Sr-O layers. Reflecting the crystal structure, this oxide exhibits Q2D conductivity with three cylinder-like Fermi surfaces labeled as α, β, and Γ. The normal state of this oxide can be well explained by the multi-orbital Fermi-liquid picture. This is actually rare among known unconventional superconductors, whose normal states mostly exhibit non-Fermi-liquid behavior.
The superconductivity of Sr2RuO4 was discovered in 1994. Soon afterwards, based on the similarity of its Landau parameters to those of 3He, the possibility of spin-triplet superconductivity in Sr2RuO4 was proposed . Although the actual pairing mechanism may be more complicated than that of superfluid 3He, this prescient suggestion stimulated research to study the SC nature of Sr2RuO4, along with research to improve sample qualities.
The spin-triplet nature is directly confirmed by the NMR [56-59] and polarized neutron scatterings; the spin susceptibility does not decrease in the SC state. Actually, the spin susceptibility can be even increased by ~2% in the SC state. Such an increase is attributed to weak spin polarization of spin-triplet Cooper pairs assisted by strong energy dependence of the density of states. The unusual spin-lattice relaxation rate 1/T1 using a 17O nuclear quadrupole resonance (NQR) technique at zero field suggests the existence of the collective motion of the d-vector order parameter. The recent discovery of half-quantum fluxoid states also provides strong support for the remaining spin-degree of freedom in the SC state.
The spontaneous time-reversal symmetry breaking was revealed by μSR measurements, and later by magneto-optical Kerr-effect measurements. These experimental results are interpreted as the consequence of chiral-p-wave superconductivity expressed by the d-vector order parameter d = Δz(kx ± iky). Here, the vector z means that the spins of Cooper pairs align perpendicular to the z axis, i.e. S = 1 and Sz = 0; and the orbital part kx ± iky, consisting of a complex combination of two order-parameter components kx and ky, indicates an angular-momentum polarized state (L = 1 and Lz = ±1) with TRS breaking.
SC phase diagram and pair breaking
One of the most important unresolved issues regarding superconductivity of Sr2RuO4 is the unusual SC phase diagram. It has been known that the upper critical field Hc2 exhibits strong suppression when the magnetic field is applied parallel to the conductive ab plane. In addition, physical quantities such as magnetization, specific heat, and thermal conductivity exhibit rapid recoveries to the normal-state values just below Hc2[68,69].
In 2013, a magnetocaloric-effect study of ultra-pure Sr2RuO4 samples revealed that the phase transition between the SC and normal states at Hc2 is actually the first-order transition, rather than the second-order transition for ordinary type-II superconductors (Fig. 2). In this study, a finite jump in the entropy as well as clear supercooling and/or superheating are observed. The first-order nature of the transition has been confirmed in specific-heat and magnetization studies.
Fig. 2: Superconducting phase diagram of Sr2RuO4 for magnetic fields parallel to the ab plane, revealed by the magneto-caloric effect measurements. The first-order transition was newly revealed below TFOT ~ 0.8 K. This fact indicates the existence of an unusual pair-breaking mechanism. Copyright (2013) American Physical Society.
For the ordinary orbital effect, in which superconductivity is destroyed by the kinetic energy cost of supercurrents flowing around quantum vortices, the SC transition at Hc2 is expected to be second order. Thus, the observed first-order transition strongly indicates that the superconductivity is destroyed by a mechanism other than the orbital effect. This scenario is supported by the small-angle neutron-diffraction study, in which the intrinsic anisotropy of the coherence lengths ξab / ξc deduced from the vortex lattice configuration is revealed to be 60, much larger than the observed anisotropy in Hc2; i.e., Hc2 // ab / Hc2 // c ~ 20. These ratios would be equal to each other if Hc2 were purely determined by the orbital effect. The anisotropy of 60 is also deduced from analysis of magnetic torque.
The first-order transition resembles the Pauli pair-breaking effect occurring in spin-singlet superconductors or possibly in spin-triplet superconductors with finite spin-susceptibility reduction. Indeed, there are theoretical suggestions attributing the first order transition to the Pauli effect. On the other hand, the Pauli effect is incompatible with NMR[56-59] and neutron-scattering  experiments: The Pauli effect should be absent in Sr2RuO4, since the characteristic field of the Pauli effect HP ~ 1 / (χn - χs)1/2 is infinitely large if the spin susceptibility in the SC state χs is the same as that in the normal state χn. I should emphasize that the Knight shift of Sr2RuO4 is dominated by the Pauli spin susceptibility, as evidenced by the observed negative hyper-fine coupling at the Ru site and by the theoretical calculation. Thus, the situation of Sr2RuO4 is clearly distinct from some conventional spin-singlet superconductors (such as vanadium) that do not exhibit change in the Knight shift in the SC state; the absence of the Knight-shift change in such materials is merely because the Pauli-susceptibility contribution in the Knight shift is just too small.
So far, there is no convincing scenario explaining the experimental results, including the observed first-order transition. One possibility is that the pairing glue is directly affected by the magnetic field.
It is worth commenting here that another candidate for chiral spin-triplet superconductivity UPt3 faces a similar situation for fields oriented along the c axis; an unusual limitation of Hc2(T) is seen despite the absence of the Knight-shift change in the SC state[20,21]. The origin has not been clarified, as in the case for Sr2RuO4. Interestingly, URu2Si2, which is believed to be a chiral spin-singlet superconductor, also exhibits limited Hc2(T), but an absence of clear NMR Knight-shift change in the SC state has been reported. Since this NMR measurement of URu2Si2 was performed with polycrystalline samples, it is highly necessary to perform NMR and neutron-scattering studies using high-quality single crystals.
Enhancement of superconductivity under strain
To understand pairing glue as well as superconducting symmetry, the study of strain-effects is of primary interest. It has been known that Sr2RuO4-Ru eutectic crystals, available under Ru-rich growth conditions, exhibit non-bulk superconductivity with Tc of up to around 3.5 K (often called the "3-K phase"). This enhancement has been attributed to the local strain near the interface induced by the lattice mismatch between Sr2RuO4 and Ru. Indeed, superconductivity with Tc ~ 3.5 K with a bulky nature can be induced by in-plane uniaxial pressure, and also by in-plane "bi-axial" pressure. Similar Tc increases have been also reported for out-of-plane uniaxial pressure but this increase may have been due to induced in-plane strain. It is also revealed that lattice dislocations of Sr2RuO4 hosts superconductivity of Tc of, again, around 3.5 K.
More recently, a uniaxial strain effect study using a novel strain-application instrument consisting of piezoelectric devices revealed that Tc of Sr2RuO4 increases substantially under both compressive and tensile strains along the direction. This is consistent with the two-component SC order parameter. An additional interesting observation is the substantial in-plane anisotropy; the strain effect along the direction is found to be much less substantial. This anisotropy is attributed to the existence of the Van Hove singularities along the direction and an provides important clue for the superconducting mechanism of Sr2RuO4.
Chiral edge currents and domains
The chiral superconductivity is characterized by chiral edge states emerging at the surface of the superconductor. Naively, it is anticipated that the chiral edge state carries chiral edge current circulating around the edge of the superconductor. In Sr2RuO4, although the existence of unusual edge states is confirmed by the tunneling spectroscopy with junctions attached to surfaces perpendicular to the ab plane, observation of edge current has never been reported. One direct method to detect the edge current is to observe local magnetic fields or moments generated by the circulating edge current. However, according to careful investigations of surface magnetic-field distribution by using scanning Hall and SQUID sensors, the spontaneous magnetic field generated by the edge current is negligibly small. Detection of the magnetic moment of Sr2RuO4 microscopic rings and discs induced by the edge current have also failed [64,87].
The absence of positive experiments for the chiral edge current has been one of the fundamental challenges for the chiral-p-wave scenario of Sr2RuO4. It is worth noting, however, that the chiral edge current in the 3He-A phase, the most established chiral condensate, has not been observed either. Also, the chiral edge current has never been reported in other chiral-superconductor candidates. There may be an overlooked mechanism to suppress the anticipated chiral edge current.
Another interesting consequence of the chiral order is the possible formation of multiple domains with positive and negative chiralities. In many junction-based experiments of Sr2RuO4 performed by various groups, unusual transport properties attributable to chiral domains and domain walls have been observed[88-91]. However, direct observation of chiral domains has not been reported in Sr2RuO4. It is also worth considering the energy cost to form chiral domain walls; based on the London theory, single-domain states are revealed to be energetically more favorable than multi-domain states. Thus, energetics and dynamics of the chiral domains and domain walls are still an interesting and unsettled issue.
Toward observation of topological phenomena
After the topological aspect of the superconductivity in Sr2RuO4 was widely recognized, a number of topological phenomena were predicted. Quantization of the thermal Hall conductivity of Sr2RuO4 in zero magnetic field would provide conclusive evidence for the chiral SC state, if experimentally observed. Electrical Hall conductivity, which is not quantized, can also take a finite value at zero field. However, capacitive techniques are probably required to detect such a zero-field Hall effect.
In 2D chiral-p-wave superconductors, it has been known that a Majorana zero mode, exhibiting non-Abelian statistics, exists in the core of a half-quantum vortex (HQV) . In integer-quantum vortices, there always exist an even number of Majorana zero modes corresponding to the two-fold spin degeneracy, and these two modes interfere to lose the non-Abelian nature. Thus, half vortices, in which an odd number of Majorana zero modes exist, are required to avoid such loss of the interesting statistics. When there are several HQVs with Majorana zero modes, adiabatic braiding of HQVs (thus zero modes) results in a change in the quantum state. This change can be utilized for quantum calculations. Signatures of HQVs (strictly speaking, half-quantum fluxoids) have been observed in ultra-sensitive cantilever torque magnetometory of Sr2RuO4 micro rings. In this experiment, half-height magnetization steps were found to emerge on the border of neighboring integer fluxoid states, when an additional magnetic field along the ring was applied. As a next step toward understanding the physics of HQVs and Majorana zero modes, it is desirable to detect HQVs with other techniques. Recently, magnetoresistance studies of Sr2RuO4 micro rings have been reported, but evidence for HQV has not been observed. Nevertheless, such efforts are important toward experimental demonstration of non-Abelian statistics by using, e.g., the Aharonov-Casher interference as suggested in Ref..
Another interesting direction is to fabricate new artificial structures toward observation of fascinating phenomena. As a promising example, a hybrid system between Sr2RuO4 and the ferromagnetic metal SrRuO3 has been recently reported (Fig. 3). In a ferromagnetic metal, ordinary spin-singlet Cooper pairs are expected to be unstable due to the strong magnetic exchange field. However, spin-triplet Cooper pairs can survive deep into the ferromagnet depending on the relative direction between the Cooper pairs and the magnetization, the latter of which is controllable by external magnetic fields. Thus, utilization of spin degree of Cooper pairs can be achieved using this new hybrid system. Demonstration of such novel functionality would in turn provide firm evidence for the superconducting order parameter in Sr2RuO4.
Fig. 3: New hybrid system of SrRuO3 (ferromagnetic metal) and Sr2RuO4 . A thin film of SrRuO3 was epitaxially grown on the ab surface of a Sr2RuO4 single crystal. The system exhibits ferromagnetism below ~160 K and the electrical contact between the two oxides is found to be really good. Because of these conditions, this system is suitable for study of the interference between ferromagnetism and spin-triplet superconductivity. Copyright (2015) The Japan Society of Applied Physics.
TRS-PRESERVED TOPOLOGICAL SUPERCONDUCTIVITY IN CuxBi2Se3
In this section, I explain superconductivity in CuxBi2Se3, one of the leading candidates for TRS-preserved TSCs, a counterpart of the TRS-broken Sr2RuO4.
The mother compound, Bi2Se3, is known as a typical topological insulator. The crystal structure consists of triangular-lattice layers of Bi and Se. The five layers of the stacking order Se-Bi-Se-Bi-Se, the quintuple layer, is a unit of the crystal structure. Cu can be intercalated in the van der Waals gap between quintuple layers to make the compound a superconductor (Tc ~ 3.8 K).
The normal-state electronic structure has been revealed with ARPES and quantum oscillation experiments . According to these experiments, superconducting CuxBi2Se3 has electron carriers of 1020 cm-3 with a cylindrical Q2D Fermi surface, rather than an ellipsoidal 3D surface.
Superconductivity in CuxBi2Se3
Superconductivity in CuxBi2Se3 has been reported in 2010. Point-contact spectroscopy reveals a conductance peak at zero bias voltage in the tunneling spectra, indicating the existence of unusual surface states. This observation is consistent with topological superconductivity. On the other hand, the scanning tunneling microscopy (STM) on the ab surface reported fully-gapped s-wave-like spectra, which in fact may be inconsistent with the conventional pairing. The temperature dependence of the specific heat suggests full-gapped superconducting states. The Hc2 curve determined by resistivity reaches 5-6 T at low temperatures for H // ab without any sign of the Pauli effect.
The possibility of odd-parity topological SC states was proposed just after the discovery of superconductivity of CuxBi2Se3. Owing to the multi-orbital nature as well as the strong spin-orbit coupling, topological superconductivity can be realized even with simple pairing channels. Among the proposed superconducting states named Δ1, Δ2, Δ3, Δ4x, and Δ4y, all except for Δ1 are odd-parity superconductivity and topologically non-trivial [106,105,3,4]. Various SC properties are calculated for these predicted states. The observed zero-bias peak in the tunneling spectra is revealed to be consistent with the Δ3, Δ4x, and Δ4y states. The SC states Δ4x and Δ4y, possessing two gap nodes or minima along the conduction plane, are of particular interest since the gap amplitudes break the trigonal rotational symmetry of the lattice. In analogy to the nematic liquid-crystal phases spontaneously breaking rotational symmetry, the Δ4x and Δ4y states can be called "nematic superconductivity".
Since 2015, several new reports on bulk experiments indicating anomalous SC features have been posted to the arXiv server. In particular, an NMR study using single crystalline samples revealed spin-triplet superconductivity and the spin-rotational-symmetry breaking in the SC state. Also, other dopants (such as Sr, Nb, and Tl) driving Bi2Se3 to superconductors have been reported[111-114]. These new experiments will certainly trigger further research to reveal the true nature of superconductivity in this compound.
Currently known candidate materials for bulk topological superconductivity are reviewed in this article, with some more in-depth explanations regarding the two leading candidates Sr2RuO4 and CuxBi2Se3. The recent extension of the concept of topology to nodal superconductivity allows expansion of the candidate list toward various long-studied unconventional superconductors, including cuprates, heavy fermions, and organics. Experimental efforts toward observation of topological phenomena in both gapped and gapless topological superconductors are important and interesting tasks for the near future.
Acknowledgements: I would like to thank A. Furusaki, Y. Yanase, S. Kittaka, S. Kasahara, K. Ishida, Y. Ando, Y. Nagai, and Y. Maeno for their insights and for our productive discussions. This work, as well as some studies explained in the text, have been supported by the Grants-in-Aids for Scientific Research on Innovative Areas on Topological Quantum Phenomena (KAKENHI 22103002) and Topological Materials Science (KAKENHI 15H05852), and the Grants-in-Aids for Scientific Research (KAKENHI 21740253, 23540407, 23110715, and 26287078) from the Japan Society for the Promotion of Science (JSPS).
Note added in proof: a theory describing the first-order transition of Sr2RuO4 with the so-called inter-orbital pair-breaking effect is proposed very recently (A. Ramires and M. Sigrist, arXiv:1605.03827).
 F. Wilczek, Nat. Phys. 5, 614 (2009).
 Y. Tanaka et al., J. Phys. Soc. Jpn. 81, 011013 (2012).
 Y. Ando and L. Fu, Annu. Rev. Condens. Matter Phys. 6, 361 (2015).
 S. Sasaki and T. Mizushima, Physica C 514, 206 (2015).
 A. P. Schnyder and P. M. R. Brydon, J. Phys.: Condens. Matter 27, 243201 (2015).
 M. Sato and S. Fujimoto, arXiv:1601.02726.
 A. Ivanov, Phys. Rev. Lett. 86, 268 (2001).
 A. J. Leggett, Rev. Mod. Phys. 47, 331 (1975).
 G. E. Volovik, The Universe in a Helium Droplet (Oxford University Press, 2003).
 T. Mizushima et al., J. Phys. Soc. Jpn. 85, 022001 (2016).
 A. P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 (2003).
 Y. Maeno et al., J. Phys. Soc. Jpn. 81, 011009 (2012).
 C. Kallin, Rep. Prog. Phys. 75, 042501 (2012).
 Y. Liu and Z.-Q. Mao, Physica C 514, 339 (2015).
 P. K. Biswas et al., Phys. Rev. B 87, 180503 (2013).
 K. Matano et al., Phys. Rev. B 89, 140504(R) (2014).
 G. M. Luke et al., Phys. Rev. Lett. 71, 1466 (1993).
 P. D. de Réotier et al., Phys. Lett. A 205, 239 (1995).
 E. R. Schemm et al., Science 345, 190 (2014).
 H. Tou et al., Phys. Rev. Lett. 80, 3129 (1998).
 W. J. Gannon et al., Phys. Rev. B 86, 104510 (2012).
 K. Izawa et al., J. Phys. Soc. Jpn. 83, 061013 (2014).
 P. Goswami and A. H. Nevidomskyy, Phys. Rev. B 92, 214504 (2015).
 E. R. Schemm et al., Phys. Rev. B 91, 140506(R) (2015).
 T. Yamashita et al., Nature Phys. 11, 17 (2015).
 H. Sumiyoshi and S. Fujimoto, Phys. Rev. B 90, 184518 (2014).
 Y. Kasahara et al., Phys. Rev. Lett. 99, 116402 (2007).
 Y. Kohori et al., J. Phys. Soc. Jpn. 65, 1083 (1996).
 E. Knetsch et al., Physica C 186-188, 300 (1993).
 N. P. Butch et al., Phys. Rev. B 84, 220504 (2011).
 Y. Nakajima et al., Sci. Adv. 1, 1500242 (2015).
 H. Kim et al., arXiv:1603.03375.
 P. M. R. Brydon et al., Phys. Rev. Lett. 116, 177001 (2016).
 Y. S. Hor et al., Phys. Rev. Lett. 104, 057001 (2010).
 G. S. Bushmarina et al., Sov. Phys. Solid state 28, 612 (1986).
 S. Sasaki et al., Phys. Rev. Lett. 109, 217004 (2012).
 M. Novak et al., Phys. Rev. B 88, 140502(R) (2013).
 M. Saghir et al., Phys. Rev. B 90, 064508 (2014).
 S. Kittaka et al., Phys. Rev. Lett. 112, 067002 (2014).
 Y. Shimizu et al., Phys. Rev. Lett. 114, 147002 (2015).
 C. C. Tsuei and J. R. Kirtley, Rev. Mod. Phys. 72, 969 (2000).
 S. Kashiwaya et al., Phys. Rev. B 51, 1350 (1995).
 M. Sato, Phys. Rev. B 73, 214502 (2006).
 J. L. Sarrao and J. D. Thompson, J. Phys. Soc. Jpn. 76, 051013 (2007).
 K. Kuroki, J. Phys. Soc. Jpn. 75, 051013 (2006).
 D. Jérome and S. Yonezawa, C. R. Phys. 17, 357 (2016).
 S. Sasaki et al., Phys. Rev. B 90, 220504(R) (2014).
 M. Sigrist et al., J. Mag. Mag. Mat. 310, 536 (2007).
 E. Bauer et al., Phys. Rev. Lett. 92, 027003 (2004).
 N. Kimura et al., Phys. Rev. Lett. 95, 247004 (2005).
 R. Settai et al., J. Phys. Soc. Jpn. 77, 073705 (2008).
 M.Nishiyama et al., Phys. Rev. Lett. 98, 047002 (2007).
 D. Aoki and J. Flouquet, J. Phys. Soc. Jpn. 83, 061011 (2014).
 Y. Maeno et al., Nature 372, 532 (1994).
 T. M. Rice and M. Sigrist, J. Phys.: Condens. Matter 7, L643 (1995).
 K. Ishida et al., Nature 396, 658 (1998).
 K. Ishida et al., Phys. Rev. B 63, 060507(R) (2001).
 H. Murakawa et al., Phys. Rev. Lett. 93, 167004 (2004).
 K. Ishida et al., Phys. Rev. B 92, 100502 (2015).
 J. A. Duffy et al., Phys. Rev. Lett. 85, 5412 (2000).
 K. Miyake, J. Phys. Soc. Jpn. 83, 053701 (2014).
 H. Mukuda et al., Phys. Rev. B 65, 132507 (2002).
 K. Miyake, J. Phys. Soc. Jpn. 79, 024714 (2010).
 J. Jang et al., Science 331, 186 (2011).
 G. M. Luke et al., Nature 394, 558 (1998).
 J. Xia et al., Phys. Rev. Lett. 97, 167002 (2006).
 S. Kittaka et al., Phys. Rev. B 80, 174514 (2009).
 K. Deguchi et al., J. Phys. Soc. Jpn. 71, 2839 (2002).
 K. Tenya et al., J. Phys. Soc. Jpn. 75, 023702 (2006).
 S. Yonezawa et al., Phys. Rev. Lett. 110, 077003 (2013).
 S. Yonezawa et al., J. Phys. Soc. Jpn. 83, 083706 (2014).
 S. Kittaka et al., Phys. Rev. B 90, 220502(R) (2014).
 C. Rastovski et al., Phys. Rev. Lett. 111, 087003 (2013).
 N. Nakai and K. Machida, Phys. Rev. B 92, 054505 (2015).
 M. Udagawa and Y. Yanase, J. Phys. Soc. Jpn. 79, 123601 (2010).
 A. M. Clogston et al., Rev. Mod. Phys. 36, 170 (1964).
 H. Dijk et al., J. Low Temp. Phys. 93, 101 (1993).
 Y. Maeno et al., J. Phys. Soc. Jpn. 66, 1405 (1997).
 H. Taniguchi et al., J. Phys. Soc. Jpn. 84, 014707 (2015).
 Y. Takanuki et al., JPS Conf. Proc. 3, 015018 (2014).
 S. Kittaka et al., Phys. Rev. B 81, 180510(R) (2010).
 H. Yaguchi et al., presented at JPS Annual Meeting (2016).
 Y. A. Ying et al., Nature Comm. 4, 2596 (2013).
 C. W. Hicks et al., Science 344, 283 (2014).
 S. Kashiwaya et al., Phys. Rev. Lett. 107, 077003 (2011).
 C. W. Hicks et al., Phys. Rev. B 81, 214501 (2010).
 P. J. Curran et al., Phys. Rev. B 89, 144504 (2014).
 F. Kidwingira et al., Science 314, 1267 (2006).
 H. Kambara et al., Phys. Rev. Lett. 101, 267003 (2008).
 M. S. Anwar et al., Sci. Rep. 3, 2480 (2013).
 K. Saitoh et al., Phys. Rev. B 92, 100504(R) (2015).
 V. Vakaryuk, Phys. Rev. B 84, 214524 (2011).
 K. Nomura et al., Phys. Rev. Lett. 108, 026802 (2012).
 A. Furusaki et al., Phys. Rev. B 64, 054514 (2001).
 X. Cai et al., Phys. Rev. B 87, 081104(R) (2013).
 E. Grosfeld et al., Phys. Rev. Lett. 83, 104513 (2011).
 M. S. Anwar et al., Appl. Phys. Express 8, 015502 (2015).
 P. Gentile et al., Phys. Rev. Lett. 111, 097003 (2013).
 E. Lahoud et al., Phys. Rev. B 88, 195107 (2013).
 S. Sasaki et al., Phys. Rev. Lett. 107, 217001 (2011).
 N. Levy et al., Phys. Rev. Lett. 110, 117001 (2013).
 T. Mizushima et al., Phys. Rev. B 90, 184516 (2014).
 M. Kriener et al., Phys. Rev. Lett. 106, 127004 (2011).
 T. Bay et al., Phys. Rev. Lett. 108, 057001 (2012).
 L. Fu and E. Berg, Phys. Rev. Lett. 105, 097001 (2010).
 M. Sato, Phys. Rev. B 81, 220504(R) (2010).
 S. Takami et al., J. Phys. Soc. Jpn. 83, 064705 (2014).
 Y. Nagai et al., Phys. Rev. B 86, 094507 (2012).
 L. Fu, Phys. Rev. B 90, 100509(R) (2014).
 K. Matano et al. arXiv:1512.07086 (2015).
 Shruti et al., Phys. Rev. B 92, 020506(R) (2015).
 Z. Liu et al., J. Am. Chem. Soc. 137, 10512 (2015).
 Y. Qiu et al., arXiv:1512.03519 (2015).
 Z. Wang et al., Chem. Mater. 28, 779 (2016).
Shingo Yonezawa is an assistant professor at the Graduate School of Science, Kyoto University. He received a D.Sc. from Kyoto University for his study on properties of quasi-one-dimensional organic superconductors. Since 2008, he has been working at the Graduate School of Science in Kyoto University. His main research interest is to discover interesting phenomena related to unconventional and topological superconductivity, mainly by means of thermodynamic, magnetic, and transport techniques. He also works on various unusual electronic phenomena.