Intrinsic Optical Dichroism in the 2d Model of Chiral Superconducting State

Recently, we have found a new intrinsic mechanism of dichroism in chiral superconductors. It relies on the existence of many orbital character of the superconducting state as in Sr2RuO4. Using three orbital but a two-dimensional model relevant for superconducting strontium ruthenate, we calculate temperature T and frequency ω dependent ac Hall conductivity σxy(ω,T) for a number of parameters. In particular, we study the changes of σxy due to changes in the light frequency ω, spin–orbit coupling λ and effective interaction parameters Uij between electrons occupying in-plane d2g orbitals. Our calculations qualitatively agree with the measured Kerr rotation angle in Sr2RuO4 and have a potential to describe other superconductors.


Introduction
The understanding of the many puzzling properties of Sr 2 RuO 4 superconductor has been a challenge for a long time [1]. In the normal state, this material seems to fulfill the criteria of a weakly correlated Fermi liquid [2][3][4]. Its superconducting state is fragile to the smallest number of impurities. Many of the thermodynamic properties including specific heat show power law temperature dependence requiring nodes in the order parameter. Other measurements indicate realization of the broken time reversed state [5]. In the bct crystal structure, reconciling consistently all of these properties requires a many band description of the system. The gaped but chiral order parameter of k x ± ik y variety exists on some of the bands and the order parameter with line nodes (k x ± ik y ) cos(k z c/2) on the other bands. Many of these puzzling properties of the system have been accounted for by a three-dimensional, three-band model [6,7]. Recently, we have shown that the existence of interorbital [8] components of the order parameter leads to natural explanation of the polar Kerr effect [9,10] in strontium ruthenate and possibly other superconductors [11].
In this work, we analyze the dependence of the dichroic signal resulted from interorbital or interband coupling using simple three band model of 2d RuO plane. Such a model with "hidden quasi-one-dimensional superconductivity" has been recently introduced [12] and argued to lead to spin triplet ground state. Here, we adopt our previous tight binding energy spectrum and in accordance with [12] assume that the dominant interactions are between electrons occupying ruthenium d xz and d yz orbitals in the spin triplet channel. The weak interactions between carriers on d xy orbitals lead to exponentially small gaps on the γ band. On the other hand, the interactions between carriers occupying d xz and d yz orbitals lead to both intraorbital and interorbital order parameters responsible for the Kerr effect.

Model and Approach
The Fermi surface of strontium ruthenate is known to appreciable details [13]. We shall use here our previous tight binding parameterization limited to the RuO plane. We neglect very weak hoppings leading to small corrugations of In both cases, the hybridization parameter t ab ≈ 0.1. Note that the alpha Fermi surface sheet has only two-fold symmetry, because of the shape of the Brillouin zone boundary the Fermi surface cylinders along the c crystallographic direction. However, we shall use all three partially occupied orbitals of Ru ions.
The normal state spectrum in the γ band derived from d xy orbitals (called c orbitals in the following) reads ε cc (k) = ε 0c − 2t (cos k x + cos k y ) − 4t cos k x cos k y , (1) while the α and β bands are derived from hybridize d xz and d yz orbitals. By symmetry, these orbitals do not hybridized with d xy . Thus, the full normal state Hamiltonian in the orbital representation and neglecting spin-orbit coupling readŝ where, in units in which the in-plane lattice constant a = 1, ε aa (k) = ε ab − 2(t ax cos k x + t ay cos k y ), ε bb (k) = ε ab − 2(t bx cos k x + t by cos k y ) and ε ab (k) = −4t ab sin k x sin k y . For the actual calculations, we use t as our energy unit. The other parameters fitted to the known experimental Fermi surface read: Our methodology is to solve the Bogoliubov-de Gennes equations self-consistently for eigenvectors (u N (k), v N (k)) T and eigenenergies E N at each k point of the two-dimensional Brillouin zone. The order parameter matrix is in general 6 × 6 matrix in the spin and orbital space. This is also true for H 0 (k) in Eq. (2) if spin indices are taken into account. This extension is necessary if spin-orbit coupling is taken into account [6]. To obtain Hall conductivity, we use Fermi golden rule to calculate the polarization dependent absorption of electromagnetic radiation [14,15], which is directly related to Im σ xy (ω, T ) [16].
In the following calculations, we shall only consider the superconducting states with chiral symmetry as only these states are expected to lead to nonzero dichroic signal observed in strontium ruthenate. For the spin triplet superconductor, the orbital character of the order parameter is odd and we expect p-wave component to dominate. The order parameters have the following form The complex coefficients Δ x ij (T ) and Δ y ij (T ) depend on temperature. a is the in-plane lattice constant.

Results
Let start the discussion of the results with presenting silent features of the underlying electron spectrum which Fermi surface consists of three sheets. They are shown in Fig. 1 in two different representations.
The Fermi surfaces were calculated assuming t ay = 0 (left panel-plusses) and for t ay = 0.06t ax (left panelcrosses). On the scale of the figure only small differences can be seen. Obviously, the spectrum of those bands calculated for t ay = 0 shows more one-dimensional character. On the right panel, the same Fermi surface sheets are shown in the expanded Brillouin zone with X point at the center and Γ point in its corner to highlight the differences between α and β bands.
In Fig. 2, the temperature dependence of the various parameters Δ ν ij (T ) and the Hall conductivity Im σ xy (ω 0 , T ) are shown. The interaction parameters are assumed to be:  Other order parameters are either exactly equal to those shown by symmetry or vanish. It is interesting to note the general symmetries between the Δ x ij and Δ y ij values in our model. They read Δ x aa = ±iΔ y bb and Δ x ab = ±iΔ y ab . The symmetry with respect to exchange of k x and k y is also obvious in the Hamiltonian. The symmetries become more complicated if spin-orbit interaction is taken into account. The diagonal order parameters in Fig. 2 have different phases φ aa = φ bb with the difference φ aa − φ bb = π . The same figure (right panel) shows the temperature dependences of the Hall conductivities normalized to their low T values calculated for two different frequencies of scattered light ω 0 = 0.003t and ω 0 = 0.0001t. Close to T c , Im σ xy (T , ω) calculated for lower frequency is strongly increased. We attribute this increase to coherence factors. It is typically observed for microwave frequencies of the order of the superconducting gaps or smaller. This increase is similar to the well-known Hebel-Slichter peak [17] and we shall use this name in the following. The zero temperature values of Fig. 4 The temperature dependence of the Hall conductivities normalized to their low T values and calculated for the same set of parameters as in Fig. 3 the ω 0 Im σ xy shown in the figure read 2.1557 10 −4 e 2 / for ω 0 = 0.003t and 5.1712 10 −6 e 2 / for ω 0 = 0.0001t. The peaks in the temperature dependence of Im σ xy (ω 0 , T ) are again shown in Fig. 4 for other set of parameters.
In Fig. 3 (left panel), we show the x components of the order parameters calculated for the couplings U ij = 0.6t for i, j = a, b, and U cc = 0.1t, and λ = 0. For these couplings the order parameter Δ cc is exponentially small. Realistic value of the spin-orbit coupling λ = −0.1t induces the cc component of the order parameter. For nonzero λ both order parameters and T c (not shown) increase. The corresponding T dependences of the Hall conductivities calculated with λ = 0, ω 0 = 0.0003t and λ = −0.1t, ω 0 = 0.0003t and ω 0 = 0.00003t are shown in the right panel of the figure.
The most important influence of spin-orbit interaction on the superconductor at hand is the coupling of all three bands. Due to λ the superconducting order parameter is induced in the γ band and persists up to the common transition temperature T c . Spin-orbit interaction also changes the symmetries between off-diagonal components of Δ ν,↑↓ ij (not shown). To understand such behavior one has to note that λ enters various matrix elements of the Hamiltonian with different signs, e.g., ε ↑↑ ab (k) is supplemented by the term −iλ, while ε ↑↑ ba (k) by +iλ. To the ac and bc components of the full spin dependent Hamiltonian matrix the terms ±λ are added. All this changes the effective Fermi level and the partial densities of states in all bands, and thus differently influences various components of the order parameters.
In the many orbital approaches, there are various contributions to the Hall conductivity. At zero temperature and for two orbital model Im σ xy (ω) is proportional 1 to the sum over wave vectors (suppressed in the formula below) of the expression containing inter alia the following terms For a nonzero Hall conductivity, it is enough to have two order parameters with nonvanishing imaginary part of their product. It may be two diagonal or one diagonal and one offdiagonal in orbital space. The asymmetry between elements v bb and v aa of velocities and nonvanishing v ab is in this approach a necessary condition for nonzero dichroic signal. In the model at hand, we have checked that for U aa = 0 the phase of Δ bb and Δ ab is the same and the Hall conductivity vanishes.
For some parameter sets, the Hall conductivity takes on large values. Its magnitude is obviously related to the frequency of light, inter-orbital couplings and other parameters. There seem to be no general rule for prediction if Hall conductivity will be large or small. The nonmonotonic dependence on temperature of Im σ xy (ω, T ) of more complicated character than the single peak is also observed (see Fig. 3). Typically, the peak appears in the dichroic spectrum for microwave frequencies comparable or smaller than the 1 Similar ideas of interband contribution to the Kerr effect as in [8] have been developed by Taylor and Kallin [9] gap function. This is in accordance with our previous 3d results and shows that to some extent the 2d model captures the physics of the real system. Of course, to have order parameters of symmetry compatible (in the group theoretical sense) with the underlying crystal structure with horizontal nodes one needs three dimensional spectrum.
In a general, three orbital case and at nonzero temperature each of the terms similar to those shown in Eq. (5) is additionally multiplied by the T dependent combination of Fermi functions and depends on T through the temperature dependence of the order parameters Δ ij (T ). The resulting expression depends on the frequency ω stemming from the conservation of energy and the eigenvalues of Bogoliubovde Gennes equation. Being a sum of many terms, each of which depends on k it may also change sign as a function of frequency. This aspect makes it difficult to assign a particular orbital as the cause of dichroism and will be discussed elsewhere.
As stated earlier [8], the Hall conductivity changes sign with chirality (i.e., the sign in front of i in k x ± ik y ) of the state. The mere existence of the spin-orbit coupling is not enough to have nonzero Hall conductivity. However, both spin-orbit coupling and the magnetic field B breaking time reversal symmetry will induce Kerr signal even in one band model in analogy to magnetooptic effects [18]. The signal depends on the magnitude of the field.
In summary, we have studied in some details the novel mechanism of dichroism operating in multiorbital superconductors breaking time reversal symmetry. The main features agree with previous numerically more involved 3d calculations. In particular, the role of spin-orbit interaction has been elucidated and found to lead to qualitative and quantitative changes of the temperature dependence of the Hall conductivity. We also established the existence of the Hebel-Slichter like coherence peak in a 2d model.