1 Introduction
The remarkable properties of the spin-triplet, p-wave phases of superfluid have stimulated research efforts to discover and identify electronic superconductors with novel broken symmetries and non-trivial ground-state topology [1–6], driven in part by predictions of novel transport properties. Chiral superfluids and superconductors are topological phases with gapless Fermionic excitations that reflect the momentum-space topology of the condensate of Cooper pairs. The A-phase of superfluid was definitively identified as a chiral p-wave superfluid by the observation of anomalous Hall transport of electrons moving through a quasiparticle fluid of chiral Fermions [7, 8]. A chiral d-wave state was proposed for doped graphene [9, 10], while a chiral p-wave state is proposed for MoS [11]. There is evidence from of broken time-reversal symmetry onsetting at the superconducting transition for the pnictide SrPtAs [12], and a chiral d-wave state has been proposed as the ground state [13]. Recent experiments also provide evidence for chiral d-wave superconductivity in the pnictide LaPt3P [14]. The perovskite, , has been studied extensively and was proposed as a promising candidate for chiral p-wave superconductivity ( pairing with ), in part based on similarities of its normal-state Fermi-liquid properties with those of liquid [15, 16]. Evidence of broken time-reversal symmetry from both and Kerr rotation measurements support an identification of as a chiral superconductor [17, 18]. However, experiments designed to detect the theoretically predicted chiral edge currents [19], or to test for the two-dimensionality of the representation that is a necessary requirement to support a chiral ground state, are so far inconclusive, or report null results [20–23]. Recent transport measurements also appear to conflict with the chiral p-wave identification based on the representation; i.e., thermal conductivity measurements at low temperatures and as a function of magnetic field, which probe the low-energy quasiparticle excitation spectrum, are consistent with the nodal line structure of a d-wave order parameter, and inconsistent with the gap structure expected based on the representation [24, 25]. The possibility that is an even parity chiral superconductor has so far not been ruled out (see, e.g., Refs. [26–30]).
The first superconductor reported to show experimental evidence of broken time-reversal symmetry was the heavy fermion superconductor, , based on linewidth measurements [31]. This experiment followed theoretical predictions of broken time-reversal symmetry in the B-phase of , i.e., the lower temperature superconducting phase [32]. Another notable signature of broken time-reversal symmetry is the onset of Kerr rotation as enters its low-temperature B-phase [33]. More recently, a neutron scattering experiment, using vortices as a probe for the superconducting state in the bulk, offers yet another piece of evidence for broken time-reversal symmetry in [6]. These results support the identification of a chiral superconducting phase of , and they also support the basic theoretical model of a multi-component order parameter belonging to a two-dimensional representation of the hexagonal point group, , in which a weak symmetry breaking field lifts the degeneracy of the two-component order stabilizing two distinct superconducting phases in zero magnetic field [32, 34]. In this theory the predicted A phase of is time-reversal symmetric with pronounced anisotropic pairing correlations in the hexagonal plane [35, 36], is preferentially selected by the symmetry breaking field, and nucleates at as the first superconducting phase. The B-phase develops as the sub-dominant partner of the two-dimensional representation nucleates at , such that the low-temperature superconducting phase spontaneously breaks both time-reversal and mirror-reflection symmetries, the latter defined by a plane containing the chiral axis which is aligned parallel (or anti-parallel) to the c-axis of .
There are four two-dimensional representations of : two even-parity representations, and , and two odd-parity representations, and , all of which allow for chiral ground states [37, 38]. The chiral ground states belonging to the and representations are defined by zeroes of the Cooper pair amplitude at points on the Fermi surface that are protected by the topology of the orbital order parameter in momentum space, i.e., , where is the azimuthal angle defining a point on the Fermi surface, and for the representations1 The bulk of the experimental evidence—thermodynamic, H-T phase diagram [37, 39, 40], thermal transport [41, 42], ultra-sound [43], Josephson tunneling [3], SANS [5] and optical spectroscopy measurements [33]—supports the identification of as an odd-parity superconductor with an order parameter belonging to the representation, and a chiral B-phase order parameter of the form, . The vector is the quantization axis along which the spin-triplet Cooper pairs have zero spin projection, i.e., an equal-spin-pairing (ESP) state [39]. A key feature of the chiral order parameter is the winding number . The Josephson interference experiment described in Ref. [3] can discriminate between and chiral ground states. Indeed the report of a phase shift in the Fraunhofer pattern for the corner-SQUID geometry, combined with the observations of broken time-reversal symmetry [31, 33], provides strong evidence in favor a chiral B-phase of . However, conclusive evidence for bulk chiral superconductivity remains elusive. A zero-field Hall transport measurement is an ideal experiment to confirm broken time-reversal and mirror symmetries in the bulk of a chiral superconductor candidate.
2 Anomalous Hall transport
The winding number of the order parameter for a chiral superconductor reflects the topology of the superconducting ground state. For a fully gapped chiral superconductor is related to the Chern number defined in terms of the Bogoliubov-Nambu Hamiltonian in 2D momentum space, or for chiral superconductors defined on a 3D Fermi surface the effective two-dimensional spectrum at fixed , , where is the Berry curvature [44]. The result for the Chern number is , which provides topological protection for a spectrum of chiral Fermions.
For 2D chiral phases there is a spectrum of massless chiral Fermions confined on the boundary (edge states) with the zero-energy state enforced by the bulk topology. However, for a chiral order parameter defined on a closed 3D Fermi surface there is also a bulk spectrum of gapless Weyl-Majorana Fermions with momenta near the nodal points , in addition to a spectrum of massless chiral Fermions confined on surfaces normal to the [1,0,0] and [0,1,0] planes [45].
2.1 Anomalous edge transport
For a fully gapped chiral p-wave ground state in two dimensions the spectrum of chiral edge Fermions is predicted to give rise to quantized heat and mass transport in chiral superfluids and superconductors [45–49]. In particular, an anomalous thermal Hall conductance is predicted to be quantized, based on the low-energy effective field theory of the chiral edge states [46, 50]. This result is also obtained from the topology of the bulk order parameter combined with linear response theory based on the Bogoliubov Hamiltonian for 2D topological superconductors [49].
For a chiral superconductor defined on a 3D Fermi surface an anamolous thermal Hall current is predicted, but is not quantized in units of a fundamental quantum of conductance. Based on the linear response theory of Qin et al. [50] Goswami and Nevidomsky obtained a result for the anomalous thermal Hall conductivity of the B-phase of for [44],
The anomalous thermal Hall conductivity reflects the number of branches of chiral Fermions confined on the [1,0,0] or [0,1,0] surface, i.e., for the chiral ground state. The non-universality of the thermal Hall conductivity is reflected by the term , which is the “distance” between the two topologically protected Weyl points at on a projected surface containing the chiral axis; e.g., for a spherical Fermi surface [44].
Thus, heat transport experiments could decisively identify the broken symmetries and topology of superconductors predicted to exhibit chiral order. The thermal conductivity depends on both the topology of the order parameter and the Fermi surface. The anomalous thermal Hall effect, in which a temperature gradient generates heat currents perpendicular to it, results from broken time-reversal and mirror symmetries—a direct signature of chiral pairing2. A zero-field thermal Hall experiment can also be used as a signature of chiral edge states. However, zero-field thermal Hall transport has remained elusive thus far.
2.2 Impurity-induced anomalous transport
Here we consider zero-field Hall transport resulting from electron-impurity interactions in the bulk of the superconductor, which we show are easily several orders of magnitude larger than the edge contribution [52]. There are earlier theoretical predictions for impurity-induced anomalous thermal Hall effects in chiral superconductors based on point-like impurities by several authors [53–55]. As we show, the point-like impurity model, which includes only s-wave quasiparticle-impurity scattering, predicts zero Hall response except for Chern number [52], i.e., only for chiral p-wave superconductors [53].
In the following we present a self-consistent theory incorporating the effects of finite-size impurities and show that such effects are essential for a quantitative description of Hall transport in chiral superconductors. Experimental observation of an impurity-induced anomalous thermal Hall effect would provide a definitive signature of chiral superconductivity. The bulk effect can easily dominate the edge state contribution to the anomalous Hall current, except in ultra-pure fully gapped chiral superconductors.
3 Transport theory
We start from the Keldysh extension [56] of the transport-like equations originally developed by Eilenberger, Larkin and Ovchinnikov for equilibrium states of superconductors [57, 58], and extended by Larkin and Ovchinnikov to describe superconductors out of equilibrium [59]. This formalism is referred to as “quasiclassical theory”. For reviews see Refs. [61–63]. The quasiclassical theory is formulated in terms of matrix propagators for Fermionic quasiparticles and Cooper pairs that describe the space-time evolution of the their non-equilibrium distribution functions, as well as the dynamical response of the low-energy spectral functions and the superconducting order parameter. Here we are interested in the response to static, or low-frequency, thermal gradients and external forces that couple to energy, mass and charge currents. We follow as much as possible the notation and conventions of theory developed for thermal transport in unconventional superconductors by Graf et al. [64].
3.1 Keldysh-Eilenberger equations
The quasiclassical transport equations are matrix equations in particle-hole (Nambu) space which describe the dynamics of quasiparticle excitations and Cooper pairs. Physical properties, such as the spectral density, currents or response functions are expressed in terms of components of the Keldysh matrix propagator,
where are the retarded (R), advanced (A) and Keldysh (K) matrix propagators.
The nonequilibrium dynamics is described by a transport equation for the Keldysh propagator,
as well as transport equations for the retarded and advanced propagators,
where
is defined in terms of the excitation energy, , the coupling to external fields, , and the self-energies, . Pairing correlations, as well as effects of scattering by impurities, phonons and quasiparticles are described by the self-energies, . The convolution product (-product) appearing in Eqs. 3, 4, in the mixed energy-time representation, is defined by,
Note that is the excitation energy and is the external time variable. The operator expansion for the convolution product is particularly useful if the external timescale, is slow compared to the typical internal dynamical timescales, and , i.e., and . In this limit we can expand Eq. 6,
The quasiclassical transport equations are supplemented by the normalization conditions [57, 58],
3.2 Quasiclassical propagators
The quasiclassical propagators are -matrices whose structure describes the internal quantum-mechanical degrees of freedom of quasiparticles and quasiholes. In addition to spin, the particle-hole degree of freedom is of fundamental importance to our understanding of superconductivity. In the normal state of a metal or Fermi liquid there is no quantum-mechanical coherence between particle and hole excitations. By contrast, the distinguishing feature of the superconducting state is the existence of quantum mechanical coherence between normal-state particle and hole excitations. Particle-hole coherence is the origin of persistent currents, Josephson effects, Andreev scattering, flux quantization, and all other nonclassical superconducting effects. The quasiclassical propagators are directly related to density matrices which describe the quantum-statistical state of the internal degrees of freedom. Nonvanishing off-diagonal elements in the particle-hole density matrix are indicative of superconductivity, indeed the onset of non-vanishing off-diagonal elements is the signature of the superconducting transition.
The Nambu matrix structure of the propagators and self energies is
The 16 matrix elements of are expressed in terms of four spin-scalars (, , , ) and four spin-vectors (, , , ). All matrix elements are functions of , , and . The spin scalars , determine spin-independent properties such as the charge, mass and heat current densities, , and , as well as the local density of states
where is the normal-state density of states at the Fermi energy. The integration is over the Fermi surface weighted by the angle-resolved normal density of states at the Fermi surface, , normalized to
The current densities are determined from Fermi-surface averages over the elementary currents, , mass, , and energy, , weighted by the scalar components of the diagonal Keldysh propagator. In particular, the charge and heat current densities are given by
The off-diagonal components, and , are the anomalous propagators that characterize the pairing correlations of the superconducting state. Spin-singlet pairing correlations are encoded in , while is the measure of spin-triplet pairing correlations. Pair correlations develop spontaneously at temperatures below the superconducting transition temperature . The anomalous propagators are not directly measurable, but the correlations they describe are observable via their coupling to the “diagonal” propagators, and , through the transport equations.
3.3 Coupling to external and internal forces
The couplings of low-energy excitations to electromagnetic fields are defined in terms of the scalar and vector potentials,
Note that encodes the charge coupling of both particle and hole excitations to the electromagnetic field. The magnetic field also couples to the quasiparticles and pairs via the Zeeman energy, , where , is the gyromagnetic ratio of the normal-state quasiparticles, and is the Nambu representation of the Fermion spin operator.
3.3.1 Mean-field self-energies
Superconductors driven out of equilibrium are also subject to internal forces on quasiparticles and Cooper pairs, originating from electron-electron, electron-phonon and electron-impurity interactions. These interactions enter the quasiclassical theory as self-energy terms, , in the transport Eqs. 3, 4, 5. We include self-energies that contribute to leading order in expansion parameters, , , that define the low-energy, long-wavelength region of validity of Landau Fermi-liquid theory, and its extension to include BCS condensation [61, 63, 65].
The leading order contributions to the self-energy from quasiparticle-quasiparticle interactions correspond the mean-field self-energies, , in the particle-hole (Landau) and particle-particle (Cooper) channels, and are represented by Eqs. 16, 17, respectively,3
Note that and represent the upper row of the Nambu matrix, . Since the mean-field self-energies are independent of , , and . The interaction vertices, , in Eq. 16 represent the quasiparticle interactions in the particle-hole channel. In the non-relativistic limit these interactions are spin-rotation invariant, in which case there are two real amplitudes: the spin-independent quasiparticle-quasiparticle interaction, , the exchange term, , describing the spin-dependent quasiparticle-quasiparticle interaction. These interactions are defined by the renormalized four-point vertex in the forward-scattering limit for quasiparticles with momenta and energies confined to the Fermi surface, i.e., and , which is a good approximation in the Fermi-liquid regime far from a quantum critical point. Thus, the propagator is integrated over the low-energy bandwidth defined by , and the corresponding self-energies depend on the direction of the quasiparticle momentum on the Fermi surface, but are independent of .
In the Cooper channel the mean-field self energy from quasiparticle interactions is given by Eq. 17. The interaction vertex separates in terms of an even-parity, spin-singlet interaction, , and an odd-parity, spin-triplet interaction, , the latter resulting from exchange symmetry in the non-relativistic limit.4 In a rotationally invariant Fermi liquid like liquid , the interactions in the Cooper channel further separates according to the irreducible representations of the rotation group in three dimensions, ,
which are labeled by the orbital angular momentum quantum number , with the basis functions given by the spherical harmonics , normalized to . The Cooper instability occurs in the pairing channel defined by the most attractive interaction, , which for is the odd-parity, spin-triplet , (p-wave) channel. For strongly correlated materials Cooper pairing is mediated by quasiparticle-quasiparticle interactions. This is necessarily the case in a single-component Fermi system like liquid , and is prevalent in strongly correlated electronic compounds such as the heavy-fermion superconductors, and , and unconventional superconductors like , all of which exhibit experimental signatures of broken time-reversal symmetry by the superconducting state. For these superconductors the pairing channel belongs to an irreducible representation of the crystal point group. Eq. 18 holds with summed over the irreducible representations of the point group, the second sum is over the set of orthogonal basis functions, , that span the irrep labeled by . For materials with hexagonal point symmetry, e.g., , we consider the four two-dimensional “E-reps”: even parity and representations and odd-parity and . All four E-reps allow for a chiral ground state with minimum Chern numbers of or .
3.3.2 Impurity self-energy
The effects of impurity disorder originate from the quasiparticle-impurity interaction, , which corresponds to the transition matrix element for elastic scattering of a quasiparticle with momentum to the point on the Fermi surface. Multiple scattering of quasiparticles and quasiholes by an impurity is described by the Bethe-Salpeter equation,
where is the t-matrix for quasiparticle-impurity scattering, and is the quasiclassical Keldysh matrix propagator for particles, holes and Cooper pairs. The leading-order contribution to the configurational-averaged self energy is then determined by scattering of quasiparticles off an uncorrelated, random distribution of statistically equivalent impurities with average density, ,
where is the t-matrix evaluated self-consistently in the forward-scattering limit. Thus, the Nambu-matrix components of the impurity Keldysh self energy, , are determined by the corresponding components of the t-matrix,
Before proceeding to non-equilibrium quasiparticle transport we need to discuss the equilibrium state, including the effects of impurity scattering, on the equilibrium states of chiral superconductors and superfluids.
4 Equilibrium
For homogeneous systems in equilibrium the transport equations for the retarded and advanced propagators reduce to
where is the mean-field order parameter and are the equilibrium self-energies resulting from quasiparticle-impurity scattering. We consider the low-temperature limit in which the thermal populations of quasiparticles and phonons are sufficiently small that we can neglect quasiparticle-quasiparticle scattering and quasiparticle-phonon scattering contributions to the self energy. Thus, we retain only the mean-field pairing self energy and impurity self energy resulting from elastic quasiparticle-impurity scattering. The propagator is also constrained by the normalization condition, which for equilibrium reduces to matrix multiplication,
A chiral superconducting ground state is defined by spontaneous breaking of time-reversal and mirror symmetries by the orbital state of the Cooper pairs. We restrict our analysis to unitary superconductors in which the Nambu matrix order parameter obeys the condition,
Unitary states preserve time-reversal symmetry with respect to the spin-correlations of the pairing state. In the clean limit is the energy gap for quasiparticles with momentum near the Fermi surface, i.e., the Bogoliubov quasiparticle excitation energy is doubly degenerate with respect to spin and given by , with and defined for on the Fermi surface. The unitarity condition is necessarily satisfied by spin-singlet pairing states, and is also the case for all known spin-triplet superconductors in which the parent state in zero external field is non-magnetic.5 An important class of unitary triplet states are the equal-spin-pairing (ESP) states defined by the spin-matrix order parameter, , in which is the direction in spin space along which the Cooper pairs have zero spin projection. Equivalently, this state corresponds to equal amplitudes for the spin projections and with . For the chiral A-phase of , the direction can be controlled by a small magnetic field, , through the nuclear Zeeman energy that orients . For chiral superconductors spin-orbit coupling and the crystalline field typically lock along a high-symmetry direction of the crystal.
We consider four classes of chiral ground states corresponding to the even-parity, spin-singlet, and , and odd-parity, spin-triplet, and representations of the hexagonal point group, . These representations all allow for chiral ground states with principle winding numbers, for the representations.6 Table 1 provides representative basis functions for these two-dimensional representations.
For even-parity, spin-singlet pairing the Nambu-matrix order parameter has the form, , where and are matrices spanning particle-hole (Nambu) space. The spin-singlet correlations are represented by the Pauli matrix , which is anti-symmetric under exchange. The orbital order parameter, , is symmetric under exchange implying . The general form of the orbital order parameter is spanned by the two-dimensional space of basis functions. For the chiral basis can be constructed from the 2D vector representation: , where the latter two forms correspond to basis functions defined on a cylindrical Fermi surface with corresponding to the azimuthal angle of . Note that chiral pairing also breaks reflection symmetry in the plane normal to the chiral axis, and has a line of nodes in the energy gap for momenta in the plane . Thus, pairing is not realized in 2D, but is defined on a 3D Fermi surface.
For chiral pairing the basis functions can be defined as , with . The latter two forms correspond to pairing defined on a cylindrical Fermi surface. Note that the chiral ground state for also breaks reflection symmetry in one or more planes containing the chiral axis, , but, in contrast to , preserves reflection symmetry in the plane normal to the chiral axis. Thus, a fully-gapped chiral ground state is possible in 2D, as well as a 3D Fermi surface that is open in the direction. For a 3D Fermi surface that is closed in the direction, the chiral ground state has topologically protected nodal points of at , and a corresponding spectrum of massless chiral Fermions in the bulk phase [44].
For odd-parity, ESP triplet states the Nambu-matrix order parameter takes the form, , where we have chosen the ESP state with .7 The ESP triplet-correlations are represented by the symmetric Pauli matrix , and the odd-parity orbital order parameter, , which is necessarily anti-symmetric under exchange, i.e., . For pairing the chiral basis is constructed from the odd-parity 2D vector representation: , the latter two forms correspond to basis functions defined on a cylindrical Fermi surface with corresponding to the azimuthal angle of . In contrast to , the chiral ground states are fully gapped in 2D, and in 3D for an open Fermi surface in the direction. For chiral pairing the basis functions are constructed from those of by multiplying by odd-parity function of . Thus, , with . These chiral states correspond to the pairing model for the B-phase of .
4.1 2D chiral superconductors
Here we consider the fully gapped and chiral ground states defined on a 2D cylindrical Fermi surface. These two cases illustrate nearly all of the key physical phenomena responsible for anomalous thermal and electrical transport mediated by non-magnetic impurity scattering in chiral superconductors. At low temperatures, thermally excited quasiparticles and phonons are dilute, therefore quasiparticles interact predominantly with quenched defects. For randomly distributed impurities, the self-energy is given by , where is the mean impurity density and is the forward scattering limit of the single-impurity -matrix in the superconducting state. This -matrix can be expressed in terms of the normal-state -matrix, and the latter can be expressed in terms of scattering phase shifts for normal-state quasiparticle-impurity scattering,
where is the normal-state density of states per spin at the Fermi surface and represents averaging over the Fermi-surface—for an isotropic 2D Fermi surface, . The superscripts denote three types of quasiclassical propagators: retarded , advanced and Keldysh . In deriving Eq. 26, the bare electron-impurity interaction is eliminated in favor of the normal-state propagator, with , and the normal-state t-matrix,
with the scattering phase shift in the cylindrical harmonic.8 Here and in the following, the directions on the Fermi surface and their corresponding azimuth angles are used interchangeably.
The mean field order parameter for unitary chiral states can be expressed in the following form, , where is the unitary matrix for singlet or triplet pairing,
in which case reduces to
for both and . Thus, in the absence of external magnetic fields, magnetic impurities or spin-dependent perturbations, the spin structure of the order parameter can be transformed away by a unitary transformation, and as previously noted the quasiparticle excitation spectrum is doubly degenerate with respect to the quasiparticle spin.
This representation of the mean-field order parameter extends to the off-diagonal components of the impurity self energy. In Eq. 29 we chose to be real. In this gauge the off-diagonal impurity self-energies reduce to
with the gauge condition, . The Nambu-matrix impurity self energy can then be expressed in terms of three functions
The term proportional to the unit Nambu matrix, , drops out of Eq. 23 for the equilibrium propagators, , and thus plays no role in determining the equilibrium properties of the superconductor. However, the unit-matrix term does contribute to the linear response of the superconductor, e.g., the a.c. conductivity [67].
The diagonal term proportional to can be combined with the excitation energy and expressed as
and similarly the impurity renormalized off-diagonal self energy is given by
Thus, for any of the chiral, unitary states described by Eq. 30, the equilibrium propagators that satisfy the transport equation and normalization condition, Eqs. 23, 24, are given by
Note that the functions and satisfy the symmetry relations: and . In equilibrium, the Keldysh propagator is determined by the spectral functions for quasiparticles and Cooper pairs, and the thermal distribution of excitations,
4.1.1 Gap equation: mean-field order parameter
The pairing interaction combined with the off-diagonal component of the Keldysh propagator determines the mean-field pairing self-energy for any of the unitary chiral states is given by the “gap equation,”
where the pairing interaction in any of the two-dimensional E-reps defined on a cylindrical Fermi surface has the form
Thus, projecting out the amplitude of the chiral mean-field order parameter we obtain the gap equation,
In practice the pairing interaction strength is eliminated in favor of the critical temperature.
The equilibrium retarded and advanced propagators are given by
where , and .9 The equilibrium spectrum is renormalized by interactions with impurities, i.e., and , where is the mean-field excitation gap from Eq. 39, and and are the diagonal and off-diagonal terms in the impurity self energy, Eq. 31.10 The self-energy is obtained from the equilibrium matrix,
Upon solving Eq. 26, we obtain
The diagonal terms and are the amplitudes for quasi-particles and quasi-holes scattering off an impurity with relative angular momentum . The off-diagonal terms, and , are the amplitudes for branch conversion scattering in which a quasi-particle (quasi-hole) scatters off an impurity and also converts to a quasi-hole (quasi-particle). The branch conversion process is accompanied by the creation (destruction) of a Cooper pair. In a chiral superconductor the Cooper pairs have angular momentum , and thus branch conversion scattering requires a corresponding change in the angular momentum of the scattered excitation, e.g., for an incident quasi-particle scattering with relative angular momentum converting to a quasi-hole and a Cooper pair of angular momentum . Thus, for branch conversion scattering to occur the quasiparticle-impurity potential must support amplitudes, , and thus non-vanishing phase shifts with , as can be seen from Eq. 42b. A direct consequence is that isotropic impurity scattering from point-like impurities does not support branch conversion scattering in chiral superconductors since the incoming and outgoing scattering states are limited to the s-wave scattering channel. As we show in what follows, the ionic radius of the impurity and branch conversion scattering are central in determining the magnitude and temperature dependence of anomalous Hall transport in chiral superconductors. Finally we note that the propagators, t-matrix and self-energies must be computed self-consistently with the gap equation, Eq. 39. In Section 4 we summarize results for thermal transport in 2D fully gapped chiral superconductors with a random distribution of finite size impurities.
4.2 3D chiral superconductors
Here we consider chiral superconductors in 3D belonging to the two-dimensional E-representations of the tetragonal and hexagonal point groups, both even- and odd-parity and representations. These groups describe the discrete point symmetries of and , respectively. See Table 1.
4.2.1 Symmetries of the order parameter
The mean-field pairing self-energy, after factoring the spin-structure using the unitary transformation in Eq. 28, has the structure,
The weak-coupling mean-field order parameter, “gap function”, is independent of energy and related to the equilibrium Keldysh pair propagator by the BCS gap equation,
where provides the pairing interaction and the energy integral is over the bandwidth of attraction, , and are the directions of the relative momentum of pairs of Fermions with zero total momentum, and is an average over the Fermi surface. The pairing interaction respects the maximal symmetry of the point group and can be expressed as a sum over bi-linear products of basis functions of the irreducible representations analogous to Eq. 18. We assume the irreducible representation, , with the most attractive pairing interaction, , dominates, in which case we can ignore the sub-dominant pairing channels,
where the interaction amplitude determines the critical temperature and denotes the basis functions for the irreducible representation, , of the relevant point group. Table 1 summarizes the basis functions, expressed in chiral basis, for the point groups and , and which are relevant for and the heavy-fermion compound , respectively.11 For a chiral superconductor the order parameter is proportional to one of the chiral basis functions, e.g., for a chiral ground state belonging to the irrep with winding number . For the analysis to follow it is sufficient to consider pairing of states near a 3D spherical Fermi surface, in which case the mean-field pairing self-energy is proportional to the spherical harmonic, i.e., , where is the orbital angular momentum corresponding to the irrep ,
where is the maximum value of the order parameter, and the normalized spherical harmonics are related to the standard spherical harmonics via .
4.2.2 Impurity self-energy
In the low temperature limit quasiparticle scattering from thermally excited quasiparticles and phonons is negligible compared to scattering off the random impurity potential. For a homogeneous uncorrelated random distribution of impurities the corresponding self-energy is a product of the mean impurity density and the forward scattering limit of the single impurity -matrix,
where the -matrices are obtained as a solution of the integral equation,
The Keldysh component of the -matrix then given by
Note that can be obtained from the symmetry relation, . In Eqs. 48, 49 we eliminated the electron-impurity matrix element, , in favor of the normal-state quasiparticle propagator, with , and the normal-state -matrix,
where is the phase shift in the relative angular momentum channel, . The normalization of the spherical harmonics is given by
An important feature of scattering theory by central force potentials, in this case the quasiparticle-impurity potential, is that the characteristic range of the potential leads to phase shifts that decay rapidly to zero for , effectively truncating the summations over and .
4.2.3 Equilibrium properties
Below we present the framework for determining the self-consistent equilibrium propagators. To highlight the effects of chiral phase winding we consider systems that are gauge-rotation invariant, i.e., invariant under a rotation around the chiral axis combined with a specific element of the . As a result the diagonal equilibrium propagator depends on only through the polar angle measured from the chiral axis, . The azimuth angle appears only in the phase factor of the pair propagator, i.e., , where is winding number of the chiral order parameter. Thus we parametrize the propagators as
where the quasiparticle and pair propagators read
with
The equilibrium Keldysh propagator is determined by and the Fermi distribution function,
Note that the retarded and advanced propagators are related by the symmetry relation, . The excitation energy and order parameter are renormalized by impurity scattering through the impurity self-energies, and , via
where is the polar-angle dependence of the normalized spherical harmonics. The functions and are defined such that
The function encodes the asymmetry in scattering rates for particles and holes, which has implications for the thermoelectric response of chiral superconductors [70].
Since the scattering potential is rotationally invariant we can expand the -matrix Equation 48 into a set of decoupled equations for each cylindrical harmonic channel. Thus, we parametrize the -matrix as follows
where the diagonal part of the -matrix is given by
and the off-diagonal part is given by
Thus, by factoring out the dependence on the azimuth angle as shown in Eq. 48, we obtain integral equations for the cylindrical harmonics of the -matrix,
where . The off-diagonal -matrix describes Andreev scattering in which an incoming particle branch converts into an outgoing hole and vice versa. This process relies on the creation and destruction of a Cooper pair. For chiral pairing, the conservation of angular momentum implies that the cylindrical harmonics of incoming and outgoing scattering states must differ by the orbital angular momentum quantum number of a Cooper pair, hence the phase winding factor in front of in Eq. 58. Since branch conversion scattering requires two distinct angular momentum channels, branch conversion scattering is absent for point-like impurities which support only s-wave scattering. Equations 61, 62 are solved for the 3D Fermi surface by expanding in the spherical harmonic basis functions, . For cylindrical Fermi surfaces the dependence on drops out and the -matrices are obtained by matrix inversion as described in Section 4.1.
5 Linear response theory
For small departures from equilibrium driven by a small temperature bias between different edges of the superconductor the heat current is proportional to the temperature gradient,
where is the thermal conductivity tensor which is constrained by the chiral symmetry of the ground state. To obtain these transport coefficients, , we determine the self-consistent, equilibrium propagators and their first-order non-equilibrium corrections to linear order in . The equilibrium propagators encode information about the bound and unbound quasiparticle pair spectrum, and are key inputs to the determination of the linear-response functions. The heat current is computed from the non-equilibrium Keldysh propagator in Eq. 14.
Here we consider the linear response functions for a static and homogeneous thermal gradient. For convenience we separate the Keldysh response into a spectral and anomalous part. The anomalous response encodes information about the non-equilibrium distribution function and is defined by,
where stands for the propagator or self-energy . We focus on the anomalous functions because the spectral response functions, , do not contribute to the thermal conductivity tensor in linear response theory to leading order in the quasiclassical expansion parameters [64]. For a uniform thermal gradient the anomalous propagator is obtained from the solution of the linearized transport equations (see Ref. [64] for the general solution),
where is the gradient of the local equilibrium distribution function. We added the subscript “eq” to denote the equilibrium propagators. We also adopt the shorthand notation,
with and defined in Eqs. 54, 57, respectively. It is also efficient to express the response functions as column vectors whose elements correspond to those of their corresponding matrices in particle-hole space,
The expression for the anomalous propagator (Eq. 65) can then be recast as
where the static thermal gradient leads to the perturbation,
The linear-response matrix is defined in terms of the equilibrium propagators,
where , and
5.1 Self-energy—vertex corrections
The r.h.s. of Eq. 69 consists of two terms. The first is the contribution that is explicitly proportional to the external field, . This term contributes only to the longitudinal thermal conductivity. Indeed the anomalous Hall conductivity arises solely from the non-equilibrium self-energy term. The self-energy corrections are the vertex corrections in the field-theoretical formulation based on Kubo response theory. These terms describe the response to perturbations by long-wavelength collective excitations of the interacting Fermi system [63]. In the context of the linear response theory developed for disordered chiral superconductors, the vertex corrections resulting from interactions of Bogoliubov quasiparticles with static impurities are obtained from the linear response corrections to the equilibrium -matrix Eqs. 48, 49 obtained from the first-order non-equilibrium corrections to the full -matrix Eqs. 21, 22. For the anomalous self-energy expressed in Nambu matrix form,
can be recast in column vector form as defined by Eq. 68,
where the impurity vertex-correction operator is given by
The retarded [advanced] -matrix elements are evaluated at , and the equilibrium -matrix elements, and , are defined such that
In general the mean-field pairing self-energy also contributes a vertex correction (i.e., ). These terms play a central role in collective mode response of the condensate, however, in the present context they contribute only to the retarded and advanced self-energies. The vertex correction contributing to anomalous heat transport arises only from the impurity-induced self-energy.
For point-like impurities, the vertex correction, and thus the anomalous Hall current, vanishes in all but chiral p-wave states. This can be shown by noting that for isotropic impurity scattering the vertex correction from Eq. 73, , is obtained from the isotropic components of the anomalous propagator. The diagonal components of the equilibrium propagators are isotropic, and thus their contribution to the linear response function involves momentum dependence only from the perturbation, . This p-wave term vanishes when averaged over the Fermi-surface, and as a result does not contribute a vertex correction. On the other hand, the off-diagonal components of the equilibrium propagators acquire the phase factor , reflecting the angular momentum of the chiral Cooper pairs. The linear response functions from these terms include a phase factor , which contributes a vertex correction only when , i.e., for chiral p-wave pairing. This is a non-universal result specific to point-like impurities. For finite-size impurities vertex corrections and thus anomalous Hall effects result for chiral superconductors with with results varying with the ionic radius of the impurity.
5.2 Cylindrical harmonic decomposition for 2D chiral superconductors
For chiral superconductors with cylindrically symmetric Fermi surfaces and pairing interactions we can parametrize the non-equilibrium corrections to the propagators and self-energies in terms of cylindrical harmonics,
The response in different cylindrical harmonic channels can decoupled such that Eq. 69 reduces to
where , , and the temperature gradient along the -axis generates the perturbation
with . The impurity self-energy correction from Eq. 74 becomes
where the vertex-correction operators are given by
with . Substituting Eq. 80 into Eq. 78 results in a linear matrix equation for the cylindrical harmonics of the anomalous response.
5.3 Spherical harmonic decomposition for 3D chiral superconductors
To exploit the axial symmetry of the Fermi surface and chiral symmetry of the order parameter, we write the anomalous propagator and self-energy as a sum of spherical harmonic components
with
The spherical harmonic components are then given by
The anomalous response in Eq. 69 can now be expressed in terms of solutions for each cylindrical harmonic component,
where the perturbation is
and the linear response matrix is given by
Similarly the vertex correction, Eq. 74, is recast as
where
Finally we use Eq. 88 to eliminate the self-energy term from Eq. 85, yielding
This equation is solved by matrix inversion.
6 Results for 2D chiral superconductors
To quantify the effects of finite-size impurities, we consider hard-disc scattering for which the scattering phase shifts are given by [71], where is the hard-disc radius and, and are Bessel functions of the first and second kind, respectively. Results presented in this section were reported in Ref. [52]. They are included here to highlight the effects of disorder on fully gapped topological chiral superconductors and to compare with new results for 3D nodal chiral superconductors. We start with the effects of impurities on the equilibrium properties and the sub-gap excitation spectrum.
6.1 Suppression of and pair-breaking
For temperatures approaching the critical temperature, , from below temperature the order parameter approaches zero continuously at the second order transition. The resulting linearized gap equation yields the transition temperature in terms of the pairing interaction, , bandwidth of attraction (“cutoff”), , and the pair-breaking effect of quasiparticle-impurity scattering. The pairing interaction and cutoff can be eliminated in favor of the clean-limit transition temperature, , with the result being a transcendental equation for the suppression of from quasiparticle scattering off the random distribution of impurities. The resulting critical temperature is given by.12
where is the digamma function, is the critical temperature and is the coherence length in the clean limit. The effects of pair-breaking by impurity scattering is determined by the pair-breaking cross section
for a chiral order parameter with winding number . For s-wave pairing , and consequently as expected from Anderson’s theorem [72]. In Figure 1, we see that is generally different from the total cross section . The two cross sections approach one another only in the point-like impurity limit . Furthermore, and both depend on the ionic radius and the winding number. A feature of the hard-disk scattering model is that for and cross at . For radii smaller (larger) than this value, pair breaking is stronger and is lower for .
6.2 Density of states
The quasiparticle spectrum, , also depends sensitively on the winding number, , as shown in Figure 2. Note the existence of multiple sub-gap impurity bound states, which are broadened into bands with increasing impurity density. These states are generated by the combination of potential scattering by impurities and multiple Andreev scattering by the chiral order parameter. As a result, the number of bound states and their sub-gap energies are determined by not only the impurity potential, e.g., the ionic radius, but also the winding number . The impurity-induced sub-gap spectrum has important implications for all quasiparticle transport processes. In the low-temperature limit, , the thermal conductivity is dominated by excitations at energies well below the clean limit gap edge. Diffusion within the lowest energy band of sub-gap states near the Fermi level determines the low temperature heat current as we discuss below.
6.3 Thermal conductivity and the anomalous Thermal Hall effect
In normal metals the thermal conductivity is limited by the transport mean free path for quasiparticles scattering off the random distribution of impurities, , where is determined by the transport cross-section.13
In the superconducting state the thermal conductivity depends on both the mean impurity density as well as the impurity cross-section via,
where we define the thermal transport lengths for the longitudinal and transverse currents by
Figure 3 shows the temperature dependence of longitudinal thermal conductivity for fully gapped chiral superconductors with . Note that the presence of impurities generally enhances the low-temperature thermal conductivity through the formation of sub-gap states, but the enhancement depends on winding number of the chiral order. For impurities with note that a band of Andreev bound states with a finite density of states at develops for a chiral order parameter with , but not for as shown in Figure 2. This is because the state with has more phase space for scattering on the Fermi surface with a nearly perfect sign change that leads to maximal pairbreaking (i.e., scattering with ) compared to the state with (scattering with ). Thus, for a gapless, diffusive, “metallic” band results in a low-temperature thermal conductivity which is linear in temperature as , i.e., . We also note that for , such behavior only occurs for sufficiently large impurity densities where the impurity bands broaden to close the gap at .
Figure 4 illustrates perhaps the most pronounced effects of finite-size impurities on transport properties. Although the longitudinal conductivity is relatively insensitive to the impurity size or the winding number, the Hall conductivity depends strongly on both and . For point-like impurities with radii smaller than the Fermi wavelength, , the thermal Hall conductivity is finite for , but is dramatically suppressed for chiral states with , as is clear in the comparison between and for shown in the lower two panels of Figure 4. This supports our previous argument that Hall currents vanish for point-like impurities, i.e., , for all chiral winding numbers except . Also note that as we increase the radius of the impurities such that , the Hall conductivity for increases dramatically and can be substantially larger than that for . Furthermore, for a fixed normal-state transport mean free path, the Hall conductivity exhibits a non-monotonic dependence on impurity size, reaching maximum at an intermediate radius. Thus, the details of the impurity potential, and thus the sub-gap spectrum, are of crucial importance for a quantitative understanding of anomalous Hall effects in chiral superconductors.
It is also instructive to compare the low-temperature limit of thermal Hall transport originating from the bulk topology in the form of chiral edge states with the bulk thermal Hall conductance from the random distribution of impurities embedded in the bulk of the superconductor. For chiral p-wave pairing the edge-state contribution to the thermal Hall conductance is universal [46, 75, 76]. By contrast the bulk impurity contribution to the low-temperature thermal Hall conductance can be expressed as [see Eq. 94], where is the effective transport length, which is non-universal and depends on the impurity density and scattering cross-section. Furthermore, is finite in a finite range of impurity density for which there is a finite density of states at , but not so disordered as to destroy superconductivity, as shown in Figure 5. At sufficiently low impurity density the spectrum is gapped at and so the edge contribution, which is linear in can dominate at sufficiently low temperatures. While above a critical impurity density superconductivity is destroyed and with it the Hall transport. However, over a significant range of impurity density both the edge and bulk impurity contributions are present for all .
To compare the edge and bulk contributions to when both are present we consider typical values of the coherence length to Fermi wavelength, , and the relative impurity size, . For example, taking , and we find at the value of that maximizes as shown in Figure 5. In general we find that the bulk contribution to the anomalous thermal Hall conductivity is generally dominant over most of the temperature range.
7 Results for chiral superconductors in 3D
We have extended the analysis for chiral states in 2D to chiral states defined on closed 3D Fermi surfaces which often include symmetry enforced line and point nodes of the excitation gap. The results reported here include anomalous thermal Hall effects in candidates for 3D chiral superconductors belonging to tetragonal and hexagonal crystalline point groups, particularly the perovskite and the heavy-fermion superconductor . To investigate the effects of ionic radius and the dependence on the ionic cross-section, we use the hard-sphere impurity potential for which the scattering phase shifts are analytically given in terms of the hard-sphere radius, , and the Fermi wavevector [77],
where and are spherical Bessel functions of the first and second kind, respectively [78].
7.1 Critical temperature
For 3D chiral superconductors we obtain a result of the same form as Eq. 91 for the suppression of by disorder, but with a pair-breaking cross-section appropriate for scattering of a 3D Fermi surface with finite-size impurities in 3D,
where the pair-breaking cross section is given by,
with are Clebsch-Gordan coefficients and is the Cooper pair angular momentum quantum number. For s-wave pairing the Clebsch-Gordan coefficient vanishes unless , and thus and , consistent with Anderson’s theorem [72]. Note also that the pair-breaking cross-section is in general different from both the total cross-section and transport cross section, which are defined by
which determine the quasiparticle scattering lifetime and transport mean-free path, respectively. For point-like impurities all of the above cross sections coincide except for pairing in the s-wave channel, in which case , but .
In the limit the total cross-section and pair-breaking cross-section both approach , i.e., four times the geometric cross section of the hard sphere impurity. However, for the pair-breaking cross section is typically smaller than the total cross section as shown in Figure 6. In the limit and approach the geometric limit, . However, in general the pair-breaking cross section is dependent on the topology of the order parameter. Chiral states with higher angular momentum are subject to stronger pair-breaking effects as we show for hard-sphere impurities. In the lower panel of Figure 6 we show the pair-breaking effects of impurity size and concentration on the critical temperature as described by Eqs. 97, 98.
7.2 Quasiparticle spectrum
Central to the interplay between chiral symmetry, topology and disorder is the impact of impurity scattering on pair-breaking and the resulting sub-gap quasiparticle spectrum. Distinct from fully gapped 2D topological states, 3D chiral ground states support symmetry protected nodes of the order parameter which leads to quasiparticle states over the entire energy range from the maximum gap on the Fermi surface down to the Fermi energy. The quasiparticle spectral function defines the angle-resolved quasiparticle density of states is determined by the retarded diagonal propagator,
The local density of states is the Fermi-surface average of the spectral function,
where is the normal-state density of states at the Fermi level. Figures 7, 8 show the effects of impurity induced scattering on the quasiparticle spectrum. The coherence peak at the maximum gap edge is broadened by impurity scattering. The spectral weight is redistributed to sub-gap energies by the formation of sub-gap resonances. The formation of sub-gap impurity bands is clearly visible in the spectral function for positions on the Fermi surface corresponding to the maximum gap as shown in the bottom panel of Figures 7, 8. These resonances correspond to Andreev bound states that hybridize with continuum states near nodal regions of the order parameter (c.f. Ref. [8]). Impurity-induced sub-gap states are formed by multiple Andreev scattering from the combined potential scattering and branch-conversion scattering by the phase-winding of the order parameter on the Fermi surface. The spectrum depends on the structure of the scattering potential as well as the topological winding number of the order parameter. These impurity-induced sub-gap states play a central role in determining the magnitude and temperature dependence of the anomalous thermal Hall conductivity because these states couple to the chiral condensate is the source of broken time-reversal and mirror symmetries.
7.3 Thermal conductivity tensor for chiral superconductors
The heat current density in Eq. 14 for chiral ground states with embedded impurity disorder reduces to
where is the spherical harmonic component of with angular momentum quantum numbers . In deriving this formula we used the symmetry relation . We also note that in linear response theory a thermal gradient does not generate a spectral response, in which case the anomalous response is equal to the Keldysh propagator, . From , we can express the longitudinal and transverse components of the thermal conductivity tensor as
where the spectral resolved transport mean free paths are defined by
with and the thermal gradient is chosen to be along the -axis, . In Eq. 104 the derivative of the Fermi distribution leads to the factor which confines the quasiparticle contribution to the heat current to excitations with . Thus, if the transport mean free paths, , vary with on a scale , then they may be approximated by , in which case the integration over the spectrum and thermal distribution yields,
In the normal state, , and the above formula yields the well known result for the normal-state thermal conductivity with given by the transport mean-free path. In particular, in the normal state the matrices that determine the anomalous response and vertex corrections are
where and given by Eqs. 99, 100. Then Eq. 90 yields the anomalous response function,
Combining Eqs. 105, 106 yields the normal-state thermal conductivity
is the transport mean free path. Furthermore, the normal state does not break time-reversal and mirror symmetries and thus vanishes.
In the chiral superconducting phase, is a low energy scale set by the width of the impurity band at the Fermi level, . When it exists a metallic-like band develops which at very low temperature in the superconducting state gives rise to diffusive heat transport that is again linear in temperature for , now for both the longitudinal and Hall conductivities. This regime is shown in Figure 9 for both components of the conductivity tensor for the pairing states and , i.e., the states with with and 2, respectively. For these states, low-energy excitations are located around the point nodes at (in the clean limit), and therefore do not contribute to low-temperature transport in the basal plane. Instead low-temperature transport relies on the impurity-induced sub-gap bands (see Figures 7, 8). The linear regime onsets for and at a threshold impurity density above which the impurity-induced resonances broaden sufficiently to generate a finite density of state at . The longitudinal conductivity increases with the impurity density as more states become available at the Fermi level even as the impurity scattering rate goes up. This behavior is due to the fact that the total cross section, which characterizes spectral broadening, is greater than the transport cross section as seen in Figure 6. Above a critical impurity density, , where is the Euler-Mascheroni constant, the number of available states no longer depends on impurity scattering, i.e., . Thus, increasing the density of impurities only decreases the thermal current by reducing the transport mean free path, , and thus . At the zero-temperature limit of is proportional to the ratio of the pair-breaking and transport cross sections, . As a result the p-wave pairing state has a universal upper bound for , i.e., independent of the structure of the impurity potential, because for (Eq. 100),
The thermal Hall conductivity also initially increases with the impurity concentration above the lower threshold density shown in Figure 9. However, peaks below as the Hall response must vanish when disorder destroys the condensate and restores time-reversal and mirror symmetries. The thermal Hall conductivity also depends strongly on the topology of the order parameter and the structure of impurity potential. The latter is highlighted by the comparison between the thermal Hall currents for the and states with decreasing impurity size. The state with is severely suppressed . This behavior results from the suppression of branch-conversion scattering which couples impurity scattering to the order parameter of the chiral condensate. In the limit of pure s-wave impurity scattering only the chiral states with can support a finite Hall conductivity (see Section 5.1).
Figure 10 shows the thermal conductivity in the zero-temperature limit for the states and , i.e., with with and 2, respectively. The presence of the line node at guarantees the availability of low-energy quasiparticles for transport in the basal plane even in the clean limit. Consequently the low-temperature limit of does not rely solely on impurity-induced sub-gap states at the Fermi level, and is finite even for as shown in the upper panels of Figure 10. Indeed approaches universal values, identical to those obtained for point-like impurities by Ref. [64], shown as the diamond symbols in Figure 10. However, the thermal Hall conductivity, does not onset at . A finite impurity density is still essential for a non-vanishing anomalous thermal Hall conductance, , at low temperatures. The reason is that the transverse component of the heat current is generated by branch-conversion scattering induced by potential scattering off the distribution of impurities. For this process to generate a finite the sub-gap Andreev resonances must be sufficiently broadened to generate a finite density of states at the Fermi level.
7.4 Comparison with the anomalous thermal Hall conductivity from Berry curvature
Anomalous Hall transport in ultra-clean topological superconductors with broken time-reversal and mirror symmetries was predicted by several authors [44, 46, 50]. In particular, anomalous Hall conductance originating from the gapless edge spectrum confined on the boundary of a topological chiral superconductor is predicted to be quantized, . This edge contribution to the anomalous thermal Hall conductance can be computed from the Berry curvature via the formula [44, 50],
where denotes the Fermi-Dirac distribution, , the quasiparticle spectrum with being the band index and the volume of the system. In the superconducting state energy eigenstates of the Bogoliubov Hamiltonian separates into two bands: above and below the Fermi level with eigenenergies where is the normal-state excitation energy measured from the Fermi level. The Berry curvature reflects the topology of the order parameter, and thus decays rapidly away from the Fermi surface (c.f. Ref. [79]). In the limit the Berry curvature confines the summation over to the Fermi surface,
where as before is the phase winding of the order parameter about the -axis. The resulting Berry phase contribution to the anomalous thermal Hall conductance for isotropic Fermi surfaces in dimensions is
in the low-temperature limit.
In Figures 9, 10 we compare the our results for the impurity-induced thermal Hall conductivity with the prediction of the edge contribution based on the Berry curvature in the low temperature limit for four different chiral ground states. The comparison is based on a typical coherence length scale . For all four chiral states the Berry phase contribution is dominant at impurity densities below the threshold for impurity-induced transverse transport in the limit . However, above this threshold the impurity-induced Hall conductivity is comparable to or much larger than the Berry phase contribution. For example, for the chiral state (bottom left panel of Figure 9), the impurity-induced Hall effect yields transverse heat currents in the zero temperature limit which are approximately an order of magnitude larger than the Berry phase contribution for typical impurity dimensions.
Figure 11 depicts the zero-field thermal Hall conductivity as a function of temperature for chiral states belonging to the spin-triplet, odd-parity and representations and the spin-singlet, even-parity and representations of the hexagonal point group, and the and, representations of . Almost all proposed chiral superconductor candidates, including the perovskite and the heavy-fermion superconductor , belong to one of these representations. The results show that the impurity-induced anomalous Hall effect (solid lines) dominates the Berry curvature contribution [44, 50] (dashed lines) over the full temperature range in all four chiral pairing states for impurities with . In this context it is worth reiterating our earlier estimate of the magnitude of the impurity-induced anomalous thermal Hall conductivity for the chiral phase of [52]. Namely, for , and , representative of [38] we estimate for for the chiral state with and impurity radius (Figure 11). Compared to the normal-state thermal conductivity at , , one needs sensitivity to transverse heat currents at the level of as shown in Figure 4.
8 Summary and outlook
We presented the theoretical framework for understanding disorder-induced anomalous Hall transport in chiral superconductors, and we reported quantitative predictions for the thermal conductivity and the anomalous thermal Hall conductivity in superconductors with phase winding for chiral superconducting ground states belonging to the 2D irreducible representations of the hexagonal and tetragonal point groups. We highlight the role of quasiparticle-impurity scattering by finite-size impurities, i.e., . Our analysis demonstrates that an anomalous thermal Hall effect is obtained for chiral superconductors with winding , provided the ionic radius of the impurities satisfies . Thus, for point-like impurities with the anomalous thermal Hall current vanishes for all but chiral p-wave ground states. We also discussed the spectrum of impurity-induced Andreev bound states, which are formed via multiple Andreev scattering. The spectrum depends sensitively on the winding number of the chiral order parameter as well as the structure of the impurity potential. Our results also show that the impurity-induced anomalous thermal Hall transport dominates the edge state contribution by an order of magnitude or more over most of the temperature range below . The impurity- and edge contributions to the thermal Hall effect both depend on broken time-reversal and mirror symmetries. Thus, they are equally good signatures of chiral superconductivity. The bulk impurity effect is likely more accessible experimentally; it produces larger Hall currents, and it is insensitive to the quality of the surfaces of a sample. In summary this work provides the theoretical framework for computing and analyzing experiments seeking to identify broken time-reversal and mirror symmetries, as well as non-trivial topology of chiral superconductors, from bulk transport measurements.
8.1 Outlook
There are a number of candidates for chiral superconductivity that have been proposed theoretically and pursued experimentally. The chiral phase of was proven to be chiral p-wave based on the observation of anomalous Hall transport of electrons embedded in superfluid [7, 8]. The heavy electron metal shows evidence of broken time-reversal symmetry based on Kerr rotation [33], Josephson interferometry [3], [31] and SANS studies of diffraction by the vortex lattice [6]. Observation of an anomalous thermal Hall effect onsetting at the A to B transition would provide a definitive bulk signature of broken time-reversal and mirror symmetries in . Analysis of the temperature- and impurity-dependences of the Hall conductivity could provide new and quantitative experimental constraints on the symmetry class of E-rep of . For a number of proposed candidates for chiral superconductivity, e.g., , doped graphene, SrPtAs, etc., observation of an anomalous thermal Hall effect would provide confirmation of broken time-reversal and mirror symmetry by the superconducting order parameter. NMR experiments revealed the existence of new superfluid phases of liquid when it is infused into low density, anisotropic, random solids - “aerogels” [80] - or confined into sub-micron cavities [81]. Analysis based on Ginzburg–Landau theory predicts that the ground state of under anisotropic confinement is a chiral phase [82]. Thus, experiments designed to measure the transverse heat current could provide a definitive test of the theory for the ground state of superfluid infused into anisotropic aerogels [83], and similarly for confined in sub-micron cavities [84].
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
VN: Writing–original draft, Writing–review and editing. JS: Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The research of VN was supported through the Center for Applied Physics and Superconducting Technologies at Northwestern University and Fermi National Accelerator Laboratory. The research of JS was supported in part by the National Science Foundation (Grant DMR-1508730), “Nonequilibrium States of Topological Quantum Fluids and Unconventional Superconductors,” and by the U.S. Department of Energy, Office of Science, National Quantum Information Science Research Centers, Superconducting Quantum Materials and Systems Center (SQMS) under contract number DE-AC02-07CH11359.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
1More complex chiral order parameters with large winding numbers are allowed by the point group symmetry, c.f. Ref. [51]; , are the loweset order harmonics consistent with the and representations, respectively.
2Note that time-reversal and mirror symmetries need not be broken simultaneously. For example, a three-band superconductor may break timereversal symmetry when the order parameter defined on each band has a different phase [60]. Mirror symmetry is however preserved and Hall effects are therefore not expected in this system.
3The sign in Eq. (17) is such that corresponds to an attractive pairing interaction.
4This separation does not apply to superconducting materials without an inversion center, i.e. non-centrosymmetric superconductors.
5The phase of superfluid , which is stabilized by an externally applied magnetic field, is a non-unitary spin-triplet state [66]. The Uranium-based ferromagnetic superconductors are also believed to be non-unitary, spin-polarized, triplet superconductors.
6, chiral ground states with any integer winding number are possible. For the discrete point group higher winding numbers with are possible for pairing basis functions exhibiting strong hexagonal anisotropy, but in general the chiral basis functions with higher winding numbers will mix with .
7Results for heat and charge transport in zero field do not depend on the choice for the direction of
8The summation over is truncated as a defect with characteristic radius leads to rapidly decaying phase shifts for .
9Hereafter the retarded (R) and advanced (A) superscripts are not shown for etc., but are implied.
10Despite its absence from spectral renormalization, encodes particle-hole asymmetry, e.g. the difference in scattering lifetimes for particles and holes which could have implications for transport properties [68], especially in thermoelectric responses [70].
11A comprehensive set of tables of basis functions for the tetragonal, hexagonal and cubic point groups is provided in Ref. [69]
12Similar results were derived for the suppression of by non-magnetic disorder in p-wave superconductors and superfluid in aerogel [73, 74].
13The transport and pair-breaking cross sections are different except for , c.f. Eq. (92).
References
1. Kycia JB, Hong JI, Graf MJ, Sauls JA, Seidman DN, Halperin WP. Suppression of superconductivity in UPt3 single crystals. Phys Rev B (1998) 58:R603–6. doi:10.1103/physrevb.58.r603
CrossRef Full Text | Google Scholar
2. Mackenzie AP, Maeno Y. The superconductivity of Sr2RuO4 and the physics of spin-triplet pairing. Rev Mod Phys (2003) 75:657–712. doi:10.1103/revmodphys.75.657
CrossRef Full Text | Google Scholar
3. Strand JD, Van Harlingen DJ, Kycia JB, Halperin WP. Evidence for complex superconducting order parameter symmetry in the low-temperature phase of UPt3 from josephson interferometry. Phys Rev Lett (2009) 103:197002. doi:10.1103/physrevlett.103.197002
PubMed Abstract | CrossRef Full Text | Google Scholar
4. Maeno Y, Kittaka S, Nomura T, Yonezawa S, Ishida K. Evaluation of spin-triplet superconductivity in Sr2RuO4. J Phys Soc Jpn (2012) 81:011009. doi:10.1143/jpsj.81.011009
CrossRef Full Text | Google Scholar
5. Gannon WJ, Halperin WP, Rastovski C, Schlesinger KJ, Hlevyack J, Eskildsen MR, et al. Nodal gap structure and order parameter symmetry of the unconventional superconductor UPt3. New J Phys (2015) 17:023041. doi:10.1088/1367-2630/17/2/023041
CrossRef Full Text | Google Scholar
6. Avers KE, Gannon WJ, Kuhn SJ, Halperin WP, Sauls JA, DeBeer-Schmitt L, et al. Broken time-reversal symmetry in the topological superconductor UPt3. Nat Phys (2020) 16:531–5. doi:10.1038/s41567-020-0822-z
CrossRef Full Text | Google Scholar
8. Shevtsov O, Sauls JA. Electron bubbles and Weyl fermions in chiral superfluid 3He-A. Phys Rev B (2016) 94:064511. doi:10.1103/physrevb.94.064511
CrossRef Full Text | Google Scholar
9. Nandkishore R, Levitov LS, Chubukov AV. Chiral superconductivity from repulsive interactions in doped graphene. Nat Phys (2012) 8:158–63. doi:10.1038/nphys2208
CrossRef Full Text | Google Scholar
10. Black-Schaffer AM, Le Hur K. Topological superconductivity in two dimensions with mixed chirality. Phys Rev B (2015) 92:140503. doi:10.1103/physrevb.92.140503
CrossRef Full Text | Google Scholar
12. Biswas PK, Luetkens H, Neupert T, Stürzer T, Baines C, Pascua G, et al. Evidence for superconductivity with broken time-reversal symmetry in locally noncentrosymmetric SrPtAs. Phys Rev B (2013) 87:180503. doi:10.1103/physrevb.87.180503
CrossRef Full Text | Google Scholar
13. Fischer MH, Neupert T, Platt C, Schnyder AP, Werner H, Goryo J, et al. Chiral d-wave superconductivity in SrPtAs. Phys Rev B (2014) 89:020509. doi:10.1103/physrevb.89.020509
CrossRef Full Text | Google Scholar
14. Biswas PK, Ghosh SK, Zhao JZ, Mayoh DA, Zhigadlo ND, Xu X, et al. Chiral singlet superconductivity in the weakly correlated metal LaPt3P. Nat Comm (2021) 12:2504. doi:10.1038/s41467-021-22807-8
CrossRef Full Text | Google Scholar
15. Rice TM, Sigrist M. Sr2RuO4: an electronic analogue of 3He? J Phys Cond Mat (1995) 7:L643–8. doi:10.1088/0953-8984/7/47/002
CrossRef Full Text | Google Scholar
17. Luke GM, Fudamoto Y, Kojima KM, Larkin MI, Merrin J, Nachumi B, et al. Time-reversal symmetry breaking superconductivity in Sr2RuO4. Nature (1998) 394:558–61. doi:10.1038/29038
CrossRef Full Text | Google Scholar
18. Xia J, Maeno Y, Peter B, Fejer M, Kapitulnik A. High resolution polar Kerr effect measurements of Sr2RuO4: evidence for broken time-reversal symmetry in the superconducting state. Phys Rev Lett (2006) 97:167002. doi:10.1103/physrevlett.97.167002
PubMed Abstract | CrossRef Full Text | Google Scholar
19. Matsumoto M, Sigrist M. Quasiparticle states near the surface and the domain wall in a px + ipy wave superconductor. J Phys Soc Jpn (1999) 68:994–1007. doi:10.1143/jpsj.68.994
CrossRef Full Text | Google Scholar
20. Kirtley JR, Kallin C, Hicks CW, Kim E-A, Liu Y, Moler KA, et al. Upper limit on spontaneous supercurrents in Sr2RuO4. Phys Rev B (2007) 76:014526. doi:10.1103/physrevb.76.014526
CrossRef Full Text | Google Scholar
21. Hicks CW, Kirtley JR, Lippman TM, Koshnick NC, Huber ME, Maeno Y, et al. Limits on superconductivity-related magnetization in Sr2RuO4 and PrOs4Sb12 from scanning SQUID microscopy. Phys Rev B (2010) 81:214501. doi:10.1103/physrevb.81.214501
CrossRef Full Text | Google Scholar
22. Curran PJ, Bending SJ, Desoky WM, Gibbs AS, Lee SL, Mackenzie AP. Search for spontaneous edge currents and vortex imaging in Sr2RuO4 mesostructures. Phys Rev B (2014) 89:144504. doi:10.1103/physrevb.89.144504
CrossRef Full Text | Google Scholar
23. Hicks CW, Brodsky DO, Yelland EA, Alexandra SG, Edkins SD, Nishimura K, et al. Strong increase of Tc of Sr2RuO4 under both tensile and compressive strain. Science (2014) 344:283–5. doi:10.1126/science.1248292
PubMed Abstract | CrossRef Full Text | Google Scholar
24. Hassinger E, Bourgeois-Hope P, Taniguchi H, René de Cotret S, Grissonnanche G, Anwar MS, et al. Vertical line nodes in the superconducting gap structure of Sr2RuO4. Phys Rev X (2017) 7:011032. doi:10.1103/physrevx.7.011032
CrossRef Full Text | Google Scholar
25. Graf MJ, Balatsky AV. Identifying the pairing symmetry in the Sr2RuO4 superconductor. Phys Rev B (2000) 62:9697–702. doi:10.1103/physrevb.62.9697
CrossRef Full Text | Google Scholar
26. Mackenzie AP, Scaffidi T, Hicks CW, Maeno Y. Even odder after twenty-three years: the superconducting order parameter puzzle of Sr2RuO4. npj Quan Mater (2017) 2:40. doi:10.1038/s41535-017-0045-4
CrossRef Full Text | Google Scholar
27. Kashiwaya S, Saitoh K, Kashiwaya H, Koyanagi M, Sato M, Yada K, et al. Time-reversal invariant superconductivity of Sr2RuO4 revealed by josephson effects. Phys Rev B (2019) 100:094530. doi:10.1103/physrevb.100.094530
CrossRef Full Text | Google Scholar
28. Pustogow A, Luo Y, Chronister A, Su YS, Sokolov DA, Jerzembeck F, et al. Constraints on the superconducting order parameter in Sr2RuO4 from oxygen-17 nuclear magnetic resonance. Nature (2019) 574:72–5. doi:10.1038/s41586-019-1596-2
PubMed Abstract | CrossRef Full Text | Google Scholar
29. Sharma R, Edkins SD, Wang Z, Kostin A, Chanchal Sow , Maeno Y, et al. Momentum-resolved superconducting energy gaps of Sr2RuO4 from quasiparticle interference imaging. Proc Nat Acad Sci (2020) 117:5222–7. doi:10.1073/pnas.1916463117
PubMed Abstract | CrossRef Full Text | Google Scholar
30. Chronister A, Pustogow A, Kikugawa N, Sokolov DA, Jerzembeck F, Hicks CW, et al. Evidence for even parity unconventional superconductivity in Sr2RuO4. Proc Nat Acad Sci (2021) 118:e2025313118. doi:10.1073/pnas.2025313118
PubMed Abstract | CrossRef Full Text | Google Scholar
32. Hess DW, Tokuyasu TA, Sauls JA. Broken symmetry in an unconventional superconductor: a model for the double transition in UPt3. J Phys Cond Matt (1989) 1:8135–45. doi:10.1088/0953-8984/1/43/014
CrossRef Full Text | Google Scholar
33. Schemm ER, Gannon WJ, Wishne CM, Halperin WP, Kapitulnik A. Observation of broken time-reversal symmetry in the heavy-fermion superconductor UPt3. Science (2014) 345:190–3. doi:10.1126/science.1248552
PubMed Abstract | CrossRef Full Text | Google Scholar
34. Machida K, Ozaki M. Splitting of superconducting transitions in UPt3. J Phys Soc Jpn (1989) 58:2244–7. doi:10.1143/jpsj.58.2244
CrossRef Full Text | Google Scholar
35. Huxley A, Rodiere P, Paul DMK, Niels van Dijk , Cubitt R, Flouquet J. Re-alignment of the flux-line lattice by a change in the symmetry of superconductivity in UPt3. Nature (2000) 406:160–4. doi:10.1038/35018020
PubMed Abstract | CrossRef Full Text | Google Scholar
36. Strand JD, Bahr DJ, Van Harlingen DJ, Davis JP, Gannon WJ, Halperin WP. The transition between real and complex superconducting order parameter phases in UPt3. Science (2010) 328:1368–9. doi:10.1126/science.1187943
PubMed Abstract | CrossRef Full Text | Google Scholar
37. Sauls JA. The order parameter for the superconducting phases of UPt3. Adv Phys (1994) 43:113–41. doi:10.1080/00018739400101475
CrossRef Full Text | Google Scholar
39. Choi CH, Sauls JA. Identification of odd-parity superconductivity in UPt3 based on paramagnetic limiting of the upper critical field. Phys Rev Lett (1991) 66:484–7. doi:10.1103/physrevlett.66.484
PubMed Abstract | CrossRef Full Text | Google Scholar
41. Graf MJ, Yip S-K, Sauls JA. Thermal conductivity of UPt3 at low temperatures. J Low Temp Phys (1996) 102:367–79. (E) 106, 727 (1997). doi:10.1007/bf00755120
CrossRef Full Text | Google Scholar
42. Graf MJ, Yip S-K, Sauls JA. Thermal conductivity of the accidental degeneracy and enlarged symmetry group models for superconducting UPt3. J Low Temp Phys (1999) 114:257–73. doi:10.1023/a:1021814206687
CrossRef Full Text | Google Scholar
43. Matthias J, Graf S-KY, Sauls JA. Identification of the orbital pairing symmetry in UPt3. Phys Rev B (2000) 62:14393–402. doi:10.1103/PhysRevB.62.14393
CrossRef Full Text | Google Scholar
44. Goswami P, Andriy H. Nevidomskyy, Topological Weyl superconductor to diffusive thermal Hall metal crossover in the B phase of UPt3. Phys Rev B (2015) 92:214504. doi:10.1103/physrevb.92.214504
CrossRef Full Text | Google Scholar
45. Tsutsumi Y, Machida K. Edge mass current and the role of Majorana fermions in A-phase superfluid 3He. Phys Rev B (2012) 85:100506. doi:10.1103/physrevb.85.100506
CrossRef Full Text | Google Scholar
46. Read N, Green D. Paired states of Fermions in two dimensions with breaking of parity and time-reversal symmetries and the fractional quantum Hall effect. Phys Rev B (2000) 61:10267–97. doi:10.1103/physrevb.61.10267
CrossRef Full Text | Google Scholar
47. Stone M, Roy R. Edge modes, edge currents, and gauge invariance in px + ipy superfluids and superconductors. Phys Rev B (2004) 69:184511. doi:10.1103/physrevb.69.184511
CrossRef Full Text | Google Scholar
48. Sauls JA. Surface states, edge currents, and the angular momentum of chiral p-wave superfluids. Phys Rev B (2011) 84:214509. doi:10.1103/physrevb.84.214509
CrossRef Full Text | Google Scholar
49. Sumiyoshi H, Fujimoto S. Quantum thermal Hall effect in a time-reversal-symmetry-broken topological superconductor in two dimensions: approach from bulk calculations. J Phys Soc Jpn (2013) 82:023602. doi:10.7566/jpsj.82.023602
CrossRef Full Text | Google Scholar
54. Li S, Andreev AV, Spivak BZ. Anomalous transport phenomena in px + ipy superconductors. Phys Rev B (2015) 92:100506. doi:10.1103/physrevb.92.100506
CrossRef Full Text | Google Scholar
55. Yip S. Low temperature thermal hall conductivity of a nodal chiral superconductor. Supercond Sci Technol (2016) 29:085006. doi:10.1088/0953-2048/29/8/085006
CrossRef Full Text | Google Scholar
56. Keldysh LV. Diagram technique for nonequilibrium processes. Zh Eskp Teor Fiz (1965) 47:1515. Sov. Phys. JETP, 20, 1018 (1965).
Google Scholar
57. Eilenberger G. Transformation of gorkov’s equation for type II superconductors into transport-like equations. Zeit F Physik (1968) 214:195–213. doi:10.1007/bf01379803
CrossRef Full Text | Google Scholar
58. Larkin AI, Ovchinnikov YN. Quasiclassical method in the theory of superconductivity. Sov Phys JETP (1969) 28:1200–5.
Google Scholar
59. Larkin AI, Ovchinnikov YN. Nonlinear conductivity of superconductors in the mixed state. Sov Phys JETP (1975) 41:960.
Google Scholar
60. Stanev V, Tešanović Z. Three-band superconductivity and the order parameter that breaks time-reversal symmetry. Phys Rev B (2010) 81:134522. doi:10.1103/physrevb.81.134522
CrossRef Full Text | Google Scholar
62. Larkin AI, Ovchinnikov YN. Nonequilibrium superconductivity. In: D Langenberg, and A Larkin, editors. Modern problems in condensed matter Physics. Amsterdam: Elsevier Science Publishers (1986).
Google Scholar
63. Rainer D, Sauls JA. Strong-coupling theory of superconductivity. In: Superconductivity: from basic Physics to new developments. Singapore: World Scientific (1994). p. 45–78. arXiv. Available at: https://arxiv.org/abs/1809.05264.
Google Scholar
64. Graf MJ, Yip S-K, Sauls JA, Rainer D. Electronic thermal conductivity and the wiedemann-franz law for unconventional superconductors. Phys Rev B (1996) 53:15147–61. doi:10.1103/physrevb.53.15147
PubMed Abstract | CrossRef Full Text | Google Scholar
65. Rainer D. Principles of ab initio calculations of superconducting transition temperatures. In: Progress in low temperature Physics, 10. Amsterdam: Elsevier Science Publishers B.V. (1986). p. 371–424.
CrossRef Full Text | Google Scholar
66. Ambegaokar V, Mermin ND. Thermal anomalies of 3He pairing in a magnetic field. Phys Rev Lett (1973) 30:81–4. doi:10.1103/physrevlett.30.81
CrossRef Full Text | Google Scholar
68. Arfi B, Bahlouli H, Pethick CJ. Transport properties of anisotropic superconductors: influence of arbitrary electron-impurity phase shifts. Phys Rev B (1989) 39:8959–83. doi:10.1103/physrevb.39.8959
PubMed Abstract | CrossRef Full Text | Google Scholar
69. Mineev VP, Samokhin K. Introduction to unconventional superconductivity. Amsterdam: Gordon and Breach Science Publishers (1999).
Google Scholar
70. Löfwander T, Fogelström M. Large thermoelectric effects in unconventional superconductors. Phys Rev B (2004) 70:024515. doi:10.1103/physrevb.70.024515
CrossRef Full Text | Google Scholar
71. Richard Lapidus I. Scattering by two-dimensional circular barrier, hard circle, and delta function ring potentials. Am J Phys (1986) 54:459–61. doi:10.1119/1.14585
CrossRef Full Text | Google Scholar
73. Larkin AI. Vector paring in superconductors in small dimensions. Sov Phys JETP Lett (1965) 2:130–2.
Google Scholar
74. Thuneberg EV, Yip S-K, Fogelström M, Sauls JA. Scattering models for superfluid 3He in aerogel. Phys Rev Lett (1998) 80:2861–4. cond-mat/9601148v2. doi:10.1103/physrevlett.80.2861
CrossRef Full Text | Google Scholar
75. Senthil T, Marston JB, Fisher MPA. Spin quantum hall effect in unconventional superconductors. Phys Rev B (1999) 60:4245–54. doi:10.1103/physrevb.60.4245
CrossRef Full Text | Google Scholar
76. Sumiyoshi H, Fujimoto S. Quantum thermal hall effect in a time-reversal-symmetry-broken topological superconductor in two dimensions: approach from bulk calculations. J Phys Soc Jpn (2013) 82:023602. doi:10.7566/JPSJ.82.023602
CrossRef Full Text | Google Scholar
77. Albert M. Quantum mechanics, I. North-Holland (1958).
Google Scholar
78. Abramowitz M, Stegun IA. Handbook of mathematical functions. New York: Dover (1970).
Google Scholar
79. Andreas PS, Brydon PMR. Topological surface states in nodal superconductors. J Phys Cond Matt (2015) 27:243201. doi:10.1088/0953-8984/27/24/243201
CrossRef Full Text | Google Scholar
80. Pollanen J, Li JIA, Collett CA, Gannon WJ, Halperin WP, Sauls JA. New chiral phases of superfluid 3He stabilized by anisotropic silica aerogel. Nat Phys (2012) 8:317–20. doi:10.1038/nphys2220
CrossRef Full Text | Google Scholar
81. Levitin LV, Bennett RG, Casey A, Cowan B, Saunders J, Drung D, et al. Phase diagram of the topological superfluid 3He confined in a nano-scale slab geometry. Science (2013) 340:841–4. doi:10.1126/science.1233621
PubMed Abstract | CrossRef Full Text | Google Scholar
82. Sauls JA. Chiral phases of superfluid 3He in an anisotropic medium. Phys Rev B (2013) 88:214503. doi:10.1103/physrevb.88.214503
CrossRef Full Text | Google Scholar
83. Sharma P, Sauls JA. Anomalous thermal Hall effect in chiral phases of 3He-aerogel. J Low Temp Phys (2022) 208:341–55. doi:10.1007/s10909-021-02657-w
CrossRef Full Text | Google Scholar
84. Sharma P, Vorontsov A, Sauls JA. Disorder induced anomalous thermal Hall effect in chiral phases of superfluid 3He. JPS Conf Proc (2023) 38:011002. Proc. of LT29, Hokkaido, Japan, 2209.04004. doi:10.7566/JPSCP.38.011002
CrossRef Full Text | Google Scholar