arxiv: v1 [cond-mat.supr-con] 8 Jan 2019

Transcription

1 Anisotropic shear stress σ xy effects in the basal plane of Sr 2 RuO 4 P. Contreras 1 1 Departamento de Física and Centro de Física Fundamental, arxiv: v1 [cond-mat.supr-con] 8 Jan 2019 Universidad de Los Andes, Mérida 5101, Venezuela Juan Moreno 2 2 Facultad de Ingeniería, Universidad de la Empresa, Montevideo 11300, Uruguay (Dated: January 10, 2019) In this short note we calculate the jumps for the specific heat C σxy, the elastic compliance S σxy xyxy and the thermal expansion α σxy due to a shear stress σ xy in the basal plane of Sr 2 RuO 4. Henceforth we clarify some issues regarding the elastic theoretical framework suitable to explain the sound speed experiments [10], and partially the strain experiments [12, 13] in strontium ruthenate. We continue to propose that the discontinuity in the elastic constant C xyxy of the tetragonal crystal Sr 2 RuO 4 gives an unambiguous experimental evidence that the Sr 2 RuO 4 superconducting order parameter Ψ has a two components order parameter with a broken time-reversal symmetry state, and that the γ band couples the anisotropic electronphonon interaction to the [xy] in-plane shear stress in Sr 2 RuO 4 [15, 16]. Keywords: Shear stress; time-reversal symmetry; unconventional superconductors; elastic constants, thermal expansion. PACS numbers: De; Rp, Pq I. INTRODUCTION In Sr 2 RuO 4 the electrons in the Cooper pairs are bound in spin triplets, where the spins are lying on the basal plane and the pair orbital momentum is directed along the z-direction. Henceforth, the order parameter Ψ is represented by a vector d(k) (of the type k x ±ik y ) [1, 2, 5, 6]. Based on the results of the Knight-shift experiment performed through T c [7, 8], it has been proposed that Sr 2 RuO 4 is a triplet superconductor, also it has been reported [9] that Ψ

2 2 breaks time reversal symmetry, which constitutes another key feature of unconventionality. The Sr 2 RuO 4 elastic constants C xyxy have been carefully measured as the temperature T is lowered through T c, showing the existence of small step in the transverse sound mode T[100] [10] which is related to the shear stress tensor component σ xy. This experimental result theoretically implies that Ψ has two different components with a time-reversal symmetry broken state [3]. Similar conclusions from a muon-spin relaxation (µsr) experiment were reported by Luke et. al. [9]. Recently, experiments on the effects of uniaxial strain ǫ xy, were performed [12], reporting that for Sr 2 RuO 4 the symmetry-breaking field can be controlled experimentally. Additionally, a most recently [13] experiment found that the transition temperature T c in thesuperconductor Sr 2 RuO 4 rises dramaticallyunder theapplicationofaplanaranisotropic strain, followed by a sudden drop beyond a larger strain. Furthermore, a recent theoretical work suggest that those recent experiments tuned the Fermi surface topology efficiently by applying planar anisotropic strain emphasizing again, the point of view that in plane effects (even by means of a more complicated renormalization group theory framework) also shows a clear evidence of a symmetry broken stated in Sr 2 RuO 4. Furthermore, they reported a rapid initial increase in the superconducting transition temperature T c, that the authors associated with the evolution of the Fermi surface toward a Lifshitz Fermi surface reconstruction under an increasing strain [14]. In this short preprint we aim to clarify some particular concepts and methods following an elasticity theory framework [24 29]. First: The term stress is used to express the loading in terms of force applied to a certain cross-sectional area of an object. The term strain is the response of a system to an applied stress. Second, Liu et at. [14] claim that only the γ band responses to the strain sensitively, and we emphasize that this physical phenomenon is because the γ band couples the anisotropic electron-phonon interaction to the [xy] in-plane shear stress σ xy [15, 16]. Third, we do not expand our analysis to a Lifshitz reconstruction of the Fermi surface mainly because we do not have experimental evidences that show a topological Lifshitz phase transition in Sr 2 RuO 4 even in its normal state, neither we have observed a generalized topological transition in Sr 2 RuO 4. In our opinion a 2D Fermi contour evolution under an applied external strain as the one for the γ band in Sr 2 RuO 4 needs further interpretations in terms of a topological electronic Lifshitz phase transition [18 21]. We remember that σ ik = p δ ik shows how pression and stress are in general

3 3 related, the stress σ ik becomes the delta function δ ik if a volumetric pression is applied to a sample. [22]. This set of previous experimental results and theoretical interpretations needs to be clarify moderately in order to unify several theoretical models which try to explain the changes occurring in the C xyxy elastic constant at T c, and the changes in the shear stress and shear strain components which, as far as we know, were carried out in several experiments [10, 12, 13]. Thus, the aim of our short note is to clarify that our elasticity Gizburg-Landau phenomenological framework applied to the shear stress σ xy partially demonstrates that Sr 2 RuO 4 is an unconventional superconductor with a two-component order parameter Ψ [3, 16]. We based our interpretation on a Sr 2 RuO 4 [T100] transversal response-impulse mode experimentally measured as the temperature T is lowered through T c showing only an small step change [10] (for that particular result please see bottom panel on Fig 5.7 and Fig 5.8 pp [10]). The result clearly shows a discontinuity in the T[100] mode. Here, let us mention that a different theory of Sr 2 RuO 4 elastic properties was presented [23]. However, this theory does not take into account the splitting of T c due to shear σ xy, and directly calculates the jumps at zero stress, where the derivative of T with respect to σ xy does not exist [4]. II. SHEAR STRESS σ xy ANALYSIS As we stated in [3, 16] when a shear stress σ xy is applied to the basal plane of Sr 2 RuO 4, the crystal tetragonal symmetry is broken, and a second order transition to a superconducting state occurs. Accordingly, for this case the analysis of the sound speed behavior at T c also requires a systematic study of the two successive second order phase transitions. Hence, the C xyxy discontinuity [10] at T c, can be explained in this context. Due to the absence of discontinuity in S xyxy for any of the one-dimensional Γ representations, the superconductivity in Sr 2 RuO 4 must be described by the two dimensional irreducible representation E 2u of the tetragonal point group D 4h [3]. If there is a double transition, the derivative of T c with respect to σ xy i.e. dt c /dσ xy is different for each of the two transition lines. At each of these transitions, the specific heat C σxy, the thermal expansion α σxy, and the elastic compliance S σxy xyxy show discontinuities [3],

4 4 the sum of them gives the correct expressions for the discontinuities at zero shear stress, where the Ehrenfest relations do not hold. The T c -σ xy phase diagram is showed in fig.(1) FIG. 1: Phase diagram showing the upper and lower superconducting transition temperatures, T c+ and T c, respectively, as functions of the shear stress σ xy In the case of an applied shear stress σ xy, the change for the in-plane Gibss free energy is given by G σxy = α( ψ x 2 + ψ y 2 )+σ xy d xyxy (ψ x ψy +ψxψ y )+ b 1 4 ( ψ x 2 + ψ y 2 ) 2 +b 2 ψ x 2 ψ y 2 + b 3 2 (ψ2 x ψ 2 y +ψ2 y ψ 2 x ), (1) where d xyxy couples the stress σ xy to the order parameter, the thermal expansion coefficient α = α (T T c0 ), and the minimization of G σxy is performed by substituting the general expression for Ψ as was previously calculated [16, 23]. Therefore G σxy becomes G σxy = α(η 2 x +η2 y )+2η xη y σ xy sinϕd xyxy + b 1 4 (η2 x +η2 y )2 + (b 2 b 3 )η 2 x η2 y +2b 3η 2 x η2 y sin2 ϕ. (2) In the presence of σ xy, the second order term determines the phase below T c+, which is characterized by ψ x and by ψ y = 0. As the temperature is lowered below T c, depending of the value of b 3 a second component ψ y may appear. If at T c a second component occurs, the fourth order terms in eqn. (2) will be the dominant one. Thus for very low T s, or for σ xy 0, a time-reversal symmetry-breaking superconducting state may emerge. The analysis of eqn. (2) depends on the relation between the coefficients b 2 and b 3. It also depends on

5 5 the values of the quantities η x and η y, and of the phase ϕ. If b 3 < 0, and η x and η y are both nonzero, the state with minimum energy has a phase ϕ = π/2. The transition temperature is obtained from eqn. (2), by performing the canonical transformations: η x = 1 2 (η µ + η ξ ) and η y = 1 2 (η µ η ξ ). After their substitution, eqn. (2) becomes G σxy = α + η 2 ξ +α η 2 µ (η2 ξ +η2 µ )2 +(b 2 +b 3 )(η 2 ξ η2 µ )2. (3) As it was done before, η ξ = η sinχ and η µ = η cosχ, eqn. (3) takes the form [ ] G σxy = α + η 2 sin 2 χ+α η 2 cos 2 χ+ η4 b 1 +(b 2 +b 3 )cos 2 2χ. (4) 4 G σxy is minimized if cos2χ = 1, this is, if χ = 0. Also, in order for the phase transition to be of second order, b, defined as b b 1 +b 2 +b 3, must be larger than zero. Therefore, if σ xy is non zero, the state with the lowest free energy corresponds to b 3 < 0, phase ϕ equal to π/2, and Ψ of the form: (ψ x,ψ y ) η (e iϕ 2, e iϕ 2 ). (5) In phase 1 of fig. (1), ϕ = 0, and as T is lowered below T c, phase 2, ϕ grows from 0 to approximately π/2. The two transition temperatures T c+ and T c are obtained to be: T c+ (σ xy ) = T c0 σ xy d α xyxy, (6) T c (σ xy ) = T c0 + b 2b 3 α σ xy d xyxy. (7) The derivative of T c+ with respect to σ xy, and the discontinuity in C + σ xy at T c+ are respectively: d T c+ d σ xy = d xyxy α, (8) C + σ xy = 2T c+ α 2 b. (9) After applying the Ehrenfest relations [24] the results for α σxy and S xyxy at T c+ are: α + σ xy = 2α d xyxy b, (10) S + xyxy = 2d2 xyxy b. (11)

6 6 For T c, the derivative of this transition temperature with respect to σ xy, and the discontinuities in the specific heat, thermal expansion and elastic stiffness respectively are: d T c = bd xyxy d σ xy 2b 3 α, (12) C σ xy = 4T c α 2 b 3 bb, (13) ασ xy = 2α d xyxy, b (14) S xyxy = d2 xyxy b. b b 3 (15) Since for the case of σ xy, the derivative of T c with respect to σ xy is not defined at zero stress point (see fig.(1)), theehrenfest relations do not holdat T c0. Thus, thediscontinuities occurring at T c0, in the absence of σ xy, are calculated by adding the expressions obtained for the discontinuities at T c+ and T c, Cσ 0 xy = 2T c0 α 2, (16) b S 0 xyxy = d2 xyxy, b 3 (17) ασ 0 xy = 0, (18) in this case, there is no discontinuity for the thermal expansion of Sr 2 RuO 4 at zero stress α 0 σ xy which is another physical feature encountered in our theory [3]. Further experimental studiesonthechangesofthethermalexpansioncoefficient α i belowt c inhtswerereported [32] showing that an additive lattice jump was found to appear spontaneously at T c as it is schematically shown in fig. (2).

7 7 FIG. 2: Schematic dependence of the thermal expansion on the temperature. Notice the additive jump in the thermal coefficient at T c for a superconducting transition of a compound with only one component Ψ. Since the phase diagram was determined as a function of σ xy, rather than as a function of the shear strain ǫ xy, (see fig. (1)), in this work, as in refs. [3, 16], we make use of the 6 6 elastic compliance matrix S. However, the sound speed measurements [10, 11] are best interpreted in terms of the elastic stiffness tensor C, with matrix elements C ijkl, which is the inverse of S [31]. However, usually the strain is easier to measure than the stress [30]. Therefore, it is important to be able to obtain the discontinuities in the elastic stiffness matrix in terms of the elastic compliance matrix for the shear stress case. Thus, close to the transitionline we follow[4]: C(T c +0 + ) = C(T c 0 + )+ C ands(t c +0 + ) = S(T c 0 + )+ S, where 0 + is positive and infinitesimal. By making use of the fact that C(T c +0 + ) S(T c +0 + ) = ˆ1, where ˆ1 is the unit matrix, it is shown that, to first order, the discontinuities satisfy, C C S C. In this manner, it is found that, for instance at T c+, the expressions that define the jumps for the discontinuities in elastic stiffness and compliances, due to an external stress, have either a positive or a negative value. In this way, S xyxy is negative; while, the stiffness components C xyxy are positive.

8 8 III. REMARKS The most noteworthy outcome of this short note, is that the observation of a discontinuity in the elastic constant C xyxy [10, 11] is an evidence that Ψ in Sr 2 RuO 4 has two components, and also the theoretical prediction that the difference in the in-plane thermal expansion coefficient ασ 0 xy = 0 is zero [3] at zero shear stress. Hence, the use of Sr 2 RuO 4 as a material in detailed studies of superconductivity symmetry-breaking effects has significant advantages because is described by a two-component order parameter. Nevertheless, determining from Sr 2 RuO 4 experimental measurements the magnitude of the parameters in the Ginzburg-Landau model is complicated [3]. The number of independent parameters occurring for the case of tetragonal symmetry in the compound Sr 2 RuO 4 are three linearly independent parameters: b 1, b 2, and b 3, which are required to specify the fourth order terms in G σxy. However for Sr 2 RuO 4, two independent ratios can be formed from the three independent b i parameters, and these two independent ratios could be determined, for example, by experimentally determining the ratios C σ + / Cσ in the presence of the normal stress σ xx and the shear stress σ xy [3]. In the experimental work of the sound velocity measurements [10, 11] in Sr 2 RuO 4 a discontinuity in the behavior for C xyxy below T c0, without a significant change in the sound speed slope as T goes below 1 Kelvin was understood as a signature of an unconventional transition to a superconducting phase [3, 10, 23]. Thus, this set of results and new other results [12 14], considers Sr 2 RuO 4 as a strong candidate for a detailed experimental investigation of the effects of a symmetry-breaking field by means of strain or stress experimental measurements. We also suggest that an in-plane thermal expansion measurement at zero shear stress might further clarify our theoretical interpretation [3]. Acknowledgments P. Contreras wishes to express his gratitude to Professors C. Lupien and L. Taillefer now at Universite de Sherbrooke from for stimulating discussions regarding their experimental results several years ago at University of Toronto. We also acknowledge Dr. Vladimir Grechka for clarifying certain aspects regarding the elastic framework here considered.

9 9 [1] Y. Maeno, T. M. Rice, and M. Sigrist, Physics Today 54, No. 1, 42 (2001). [2] Y. Maeno, H. Hashimoto, K. Yoshida, S. Nishizaki, T. Jujita, J. G. Bednorz, and f. Lichtenberg, Nature (London) 372, 532 (1994). [3] M. B. Walker and P. Contreras, Phys. Rev. B 66(21) (2002) [4] M. B. Walker, Phys. Rev. B 22 (1980) [5] A. Mackenzie, and Y. Maeno, Rev. of Modern Physics 75, 657 (2003). [6] T. Rice and M. Sigrist, J. Phys. Condens. Matter 7, L643 (1995). [7] K. Ishida, H. Mukuda, Y. Kitaoka, K. Asayama, Z. Q. Mao, Y. Mori, Y. Maeno, Nature (London) 396, 658 (1998). [8] J. Duffy, S. Hayden, Y. Maeno, Z. Mao, J. Kulda, and G. McIntyre, Phys. Rev. Lett. 85, 5412 (2000). [9] G. M. Luke, Y. Fudamoto, K. M. Kojima, M. I. Larkin, J. Merrin, B. Nachumi, Y. J. Uemura, Y. Maeno, Z. Q. Mao, H. Nakamura, and M. Sigrist, Nature 394, 558 (1998). [10] C. Lupien, Ph. D. Thesis. University of Toronto, [11] C. Lupien, W. A. MacFarlane, Cyril Proust, Louis Taillefer, Z. Q. Mao, and Y. Maeno, Phys. Rev. Lett. 86, 5986 (2001). [12] C. Hicks, D. Brodsky, E. Yelland, A. Gibbs, J. Bruin2, and M. Barber Science (2014). [13] A. Steppke, L. Zhao, M. Barber, T. Scaffidi, F. Jerzembeck, H. Rosner, A. Gibbs, Y. Maeno, S. Simon, A. Mackenzie, and C. Hicks arxiv (2016). [14] Y. Liu, F. Zhang, T. Rice, and Q. Wang Nature: Quantum Materials 2:12 (2017). [15] M. Walker, M. Smith, and K. Samokhin. Phys. Rev. B 65 (2002) [16] P. Contreras, Ph.D. Thesis. University of Toronto (2006) [17] P Contreras, J Flores, R Almeida Rev. Mex. Fis. 62(5) 442 (2016) [18] Lifshitz, I. Journal of Experimental and Theoretical Physics (JETP) (1960). [19] Kaganov, M. and Lifshitz, I. Sov. Phys. Usp (1979) [20] Kaganov, M. and Contreras, P., Journal of Experimental and Theoretical Physics (JETP) 79(6) 985 (1994). [21] Kaganov, M. and Nurmagarnbetov, A. Journal of Experimental and Theoretical Physics (JETP) (1982)

10 10 [22] L. D. Landau and E. M. Lifshitz, Fluid Mechanics Elsevier (2009) [23] M. Sigrist Prog. Theor. Phys. 107, 917 (2002). [24] L. D. Landau and E. M. Lifshitz, Statistical Physics Butterwoth-Heinemann (Oxford, 1980) [25] L. D. Landau and E. M. Lifshitz, Elasticity Theory Pergamon Press (London, 1970) [26] L. R. Testardi, Phys. Rev. B 3, 95 (1971). [27] A. Auld, Acoustic fields and waves in solids (Wiley ans Sons, New York, 1973). [28] Musgrave, M. Crystal Acoustics Holden-Day (1970). [29] Philip, J., Pramana-J. Phys, 29(4) 391 (1987). [30] Fiber optics yields dynamic fidelity of a fraction of a nanostrain, and tiltmeters yield comparable static fidelity, according to Dr. Grechka explanation. [31] F. Ney Physical Properties of Crystals (Oxford, London, 1979). [32] T. Asahi, H. Suzuki, M. Nakamura, H. Takano, and J. Kobayashi Phys. Rev. B 55 (1997) 9125