key: cord-0651272-88eynwhx
authors: Lin, Shi-Zeng; Kim, Duk Y.; Bauer, Eric D.; Ronning, Filip; Thompson, J. D.; Movshovich, Roman
title: Interplay of the spin density wave and a possible Fulde-Ferrell-Larkin-Ovchinnikov state in $mathrm{CeCoIn_5}$ in rotating magnetic field
date: 2019-02-13
journal: nan
DOI: nan
sha: 819c1cb26fdc10619c991e5ccf52d9e8acb35da6
doc_id: 651272
cord_uid: 88eynwhx

The $d$-wave superconductor $mathrm{CeCoIn_5}$ has been proposed as a strong candidate for supporting the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state near the low-temperature boundary of its upper critical field. Neutron diffraction, however, finds spin-density wave (SDW) order in this part of the phase diagram for field in the $a$-$b$ plane, and evidence for the SDW disappears as the applied field is rotated toward the tetragonal $c$ axis. It is important to understand the interplay between the SDW and a possible FFLO state in $mathrm{CeCoIn_5}$, as the mere existence of an SDW does not necessary exclude an FFLO state. Here, based on a model constructed on the basis of available experiments, we show that an FFLO state competes with an SDW phase. The SDW state in $mathrm{CeCoIn_5}$ is stabilized when the field is directed close to the $a$-$b$ plane. When the field is rotated toward the $c$ axis, the FFLO state emerges, and the SDW phase disappears. In the FFLO state, the nodal planes with extra quasiparticles (where the superconducting order parameter is zero) are perpendicular to the field, and in the SDW phase, the quasiparticle density of states is reduced. We test this model prediction by measuring heat transported by normal quasiparticles in the superconducting state. As a function of field, we observe a reduction of thermal conductivity for field close to the $a$-$b$ plane and an enhancement of thermal conductivity when field is close to the $c$ axis, consistent with theoretical expectations. Our modeling and experiments, therefore, indicate the existence of the FFLO state when field is parallel to the $c$ axis.

Introduction.-A key property of a superconductor is how it responds to an external magnetic field. Generally, magnetic field suppresses superconductivity. For superconductors with singlet Cooper pairing, magnetic field destroys superconductivity in two ways. Firstly, the Lorentz force that an external magnetic field exerts on the two electrons of a Cooper pair have opposite directions. These forces tear a Cooper pair apart, thereby suppressing superconductivity via orbital limiting. Secondly, a magnetic field tends to polarize electron spins via Zeeman coupling, reducing the electrons' energy and leading to an enhanced Pauli susceptibility. Under certain circumstances, superconductivity can be destroyed by this mechanism of Pauli limiting even when Cooper pairs form with electrons of opposite spin alignment. While the maximal magnetic field that most known superconductors can sustain is defined by orbital limiting, some superconductors have been identified whose upper critical field is determined largely by Pauli limiting. These superconductors also can stabilize a spatially modulated superconducting state, known as the Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) state, before superconductivity is suppressed entirely [1, 2] . In the FFLO state, Cooper pairs acquire a nonzero momentum, and the superconducting order parameter vanishes locally in space along nodal planes. This results in an excess of quasiparticles, which can significantly modify physical properties of the superconductors. Theoretically, it is clear that the formation of an FFLO state is unavoidable when a clean Pauli limited superconductor is subjected to a strong magnetic field. However, the experimental detection of an FFLO state remains a challenge despite some encouraging experimental evidence [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] , mainly because the superconducting order parameter cannot be measured directly. Often, competing effects render the interpretation of the ex-perimental data difficult.

The discovery of the heavy-fermion superconductor CeCoIn 5 with tetragonal crystal structure and a transition temperature T c = 2.3 K has provided an exciting playground to search for the FFLO state. High quality CeCoIn 5 with a large electron mean free path l ≈ 10ξ can be achieved, where ξ is the superconducting coherence length. Various measurements have revealed that the upper critical field in CeCoIn 5 is mainly limited by the Pauli mechanism. Initial experimental measurements have shown the existence of a new phase inside the superconducting phase for field both along the crystal c and a axis, see Ref. 14 for a review. For instance, a double-peak structure has been observed in NMR spectra in the high-field and low-temperature corner of the superconducting phase diagram [4] . In the case of H a, neutron-scattering measurements have identified this new phase as a spin-density-wave (SDW) state that coexists with the superconducting state. In the SDW, magnetic moments with magnitude 0.1 µ B (µ B is the Bohr magneton) are aligned along the crystalline c axis due to crystal field effects. A similar SDW is observed in the slightly Hg-doped CeCoIn 5 [15] . This SDW is induced by magnetic field and disappears when superconductivity is destroyed by field. The SDW phase is suppressed (eventually entirely) when the field is rotated away from the a-b plane [16] . The double peak structure in NMR spectra for H a can be explained in terms of the SDW phase. However the state responsible for similar NMR features for H c requires further study.

The appearance of an SDW only inside the superconducting phase comes as a surprise, because it is generally believed that SDW competes with superconductivity for the density of states at the Fermi surface. Several scenarios have been put forward to account for the stabilization of the SDW. These scenarios highlight the importance of the vortex lattice [17] , Pauli pair breaking [18, 19] , the FFLO state [20, 21] , a superconducting pairing density wave [22] [23] [24] , and Fermi surface nesting improved by the magnetic field [25] [26] [27] . It was shown in Ref. 28 that an in-plane magnetic field enhances the transverse magnetic susceptibility in Pauli-limited d-wave superconductors, which can result in a divergence of the dressed susceptibility. As a consequence, the magnetic fluctuations condense and static SDW order sets in. This magnon condensation picture is favored by inelastic neutron-scattering measurements [29] [30] [31] .

The emergence of an SDW does not rule out the existence of an FFLO state. In the case when the SDW competes with the FFLO state, one can suppress the SDW phase by rotating the field out of the a-b plane (as is the case in CeCoIn 5 ), which in turn can favor the FFLO state. This requires understanding the interplay between the SDW and the FFLO states, with a proper model that is relevant for CeCoIn 5 . In this Letter, we study the relation between SDW and FFLO states using a theoretical model described below. We find that the SDW competes with the FFLO phase. When H a, an SDW phase emerges inside the superconducting phase due to magnon condensation triggered by the magnetic field. When H is rotated towards the c axis, the SDW phase is disfavored and the FFLO state appears. Because of the coupling between superconductivity and magnetism, the SDW order induces weak modulation in the superconducting order parameter and vice versa. Guided by the theory, we performed thermalconductivity measurements in CeCoIn 5 for a magnetic field applied within the a-c plane as a function of angle away from the c-axis. The data are consistent with expectations from the model. The synergy between modeling and experiment suggests the presence of an FFLO phase in CeCoIn 5 when H c. Model.-We construct a model Hamiltonian based on the following experimentally established facts. (1) Various experiments have shown that CeCoIn 5 is close to an SDW instability [29] [30] [31] [32] [33] [34] [35] [36] , which is consistent with the fact that a field of order of 11 T is sufficient to trigger an SDW instability. (2) The SDW is formed by gapping quasiparticles in nodes of the superconducting d x 2 −y 2 order parameter. The SDW ordering wave vector is Q = (0.44, ±0.44, 0.5). (3) The SDW moments have strong Ising anisotropy and lie along the c axis. (4) The Fermi surface relevant for superconductivity is a warped cylinder. Based on these facts, we can construct the following mean-field Hamiltonian in two dimensions:

where t i j describes electron hopping on a square lattice with the dispersion (k) = 2t 1 [cos(k x ) + cos(k y )] + 4t 2 cos(k x ) cos(k y ) + 2t 3 (NN), the second NN (along the diagonal), and the third NN (along the bond) hopping amplitudes, respectively, and their strengths have been estimated by using density functional theory [37] . The electron density for up and down spins is

is the Fermi distribution function and u i↑,l , v i↓,l , E l are the l-th eigen vector and eigen energy of the Bogoliubov-de Gennes equation associated with Eq. (1) [38] . We fix the electronic density at n = 0.72 by tuning the chemical potential µ.

The self-consistent equation for the pairing potential ∆ i j and molecular field h i are ∆ i j = V 4 l u i↑,l v * j↓,l + u j↑,l v * i↓,l tanh(E l /2T ) and h i = − j J i j (n j↑ − n j↓ ). The spin of conduction electrons is

There is a strong anisotropy in the electron g-factor [39] , and we take g c /g ab = 3. Here µ B is the Bohr magneton. To stabilize SDW order at Q, we use J i j = J 1 (NN antiferromagnetic interaction) and the third NN (competing) interaction J 3 = −J 1 /4 cos(2πQ x ). To ensure d-wave pairing symmetry, we choose the NN pairing potential ∆ i j with V = 4.5t 1 in the calculations. The d-wave order parameter is given by

wherex andŷ are the unit vectors in the x and y direction, respectively [40] . Similar models were introduced to describe the emergence of an SDW state in Pauli-limited superconductors in high field [25, 26] , as well as the phase diagram for Nd-doped CeCoIn 5 [27, 41] .

This model describes the competition between an SDW and d-wave superconductivity, which depends on values of V and J 1 [42] . When J 1 /t 1 > 3.6, SDW order develops and coexists with superconducting order [41] . To model CeCoIn 5 , the system is tuned to an SDW instability by setting J 1 /t 1 = 3.2. The phase diagram at zero temperature T = 0, when field is rotated in the a-c plane, is displayed in Fig. 1 . For H a, the model correctly describes the development of the SDW phase inside the superconducting state at high field. When field is canted towards the c axis, the SDW state is suppressed, and the FFLO state appears [43] . The transition from the uniform superconducting state to the FFLO (SDW) state is of first (second) order, while the transition from the FFLO state (SDW) to the normal state is of second (first) order. The spatial distributions of ∆ d and SDW orders are shown in Fig. 2 . In the FFLO state, the modulation of ∆ d induces modulation of the SDW; while in the SDW phase, modulation of SDW order generates weak modulation in ∆ d . Because of the competition between SDW and superconductivity, the maxima in ∆ d corresponds to the minima in the SDW order, and vice versa.

In our 2D calculations, the wavevector of the FFLO, Q FFLO , is confined to the plane. For an s-wave superconductor with a vortex lattice, Q FFLO is parallel to the vortex lines and H [44, 45] . For d-wave, one would also expect that Q FFLO H, otherwise there would be modulation of the superconducting order parameter with two different periods in the a-b plane. This is energetically disfavored both by Pauli and orbital pair breaking effects. This picture is supported by numerical calculation based on a quasiclassical theory for a single vortex in a d-wave superconductor [46] . For Q FFLO H, the thermal conductivity will be enhanced when thermal current J ⊥ H due to the excess of quasiparticles around nodal planes in the FFLO state. In contrast, some of nodal quasiparticles of the dwave superconductivity are gapped out in the SDW state and one expects a reduction of the thermal conductivity. These the crystal because the electronic mean free path is mainly limited by vortex scattering even with very large inter-vortex distances at low fields. For H close to the a-b plane, the thermal conductivity κ drops sharply with decreasing field just below H c2 . This drop originates from the development of the SDW and concomitant pair-density-wave (PDW) orders which gap out some of d-wave nodal quasiparticles [47] [48] . As field is rotated toward the c axis, the drop decreases, signaling the suppression of SDW order. The threshold angle where this sharp drop disappears is approximately 70 • , which is consistent with previous measurements [16] . For H close to the c axis, a significant increase of κ is observed as H approaches H c2 from below. We ascribe this enhancement to a signature of the FFLO state, in accord with the model calculation. Such enhancement of thermal conductivity below H c2 is not expected for a Pauli-limited superconductor with a first-order superconducting transition. Instead, thermal conductivity below H c2 should be rather field independent up to H c2 , and display a sharp step at H c2 . A contour plot of κ in the complete field and angle range is depicted in Fig. 4 . The behavior of κ near H c2 is consistent with the expectation from our model. Discussions and summary.-We stress the importance of Ising-like magnetic anisotropy in the model. The presence of this anisotropy naturally explains why the SDW phase only occurs for field close to the a-b plane within the magnon condensation picture [28] . It also leads to the competition between the SDW and FFLO states, as shown in Fig. 1 . Our approach and results are different from a model with isotropic magnetic fluctuations [20] , where it was proposed that the SDW phase is stabilized by Andreev bound states localized around FFLO nodal planes. We too observe an induced SDW oscillation in the FFLO state [see Fig. 2 ], but with a much weaker amplitude compared to that for H [100], where the full SDW develops. In addition, the picture that the SDW state for H [100] is induced by the FFLO state is not supported by recent neutron-scattering measurements, which reveal that the SDW state is induced by closing the magnon gap [29] [30] [31] . A number of measurements reveal the presence of a quantum critical point around H c2 = 5 T when H

[001] [33, 49, 50] . Strong quantum fluctuations should suppress the thermal conductivity near H c2 . There must exist another mechanism which counters this suppression and leads to the enhancement of κ, as observed experimentally. Model calculations show that κ increases with H near H c2 in orbital-limited superconductors [51] . The orbital-limited critical field H orb c2 in CeCoIn 5 is about three times larger than the Pauli limited critical field H P c2 , and therefore, the experimentally measured superconducting critical field H c2 as well. The effect of orbital limiting on thermal conductivity near H c2 , can be gleaned from calculations for the case of orbital limiting alone around a field of about 1 3 H c2 [52, 53] . The variation of thermal conductivity with field is very slow in this region, in contrast to the sharp increase in the thermal conductivity near H c2 shown in Fig. 3 . Therefore, orbital limiting cannot be causing the observed increase in κ. The vortex lattice undergoes a number of structural transitions as a function of field strength when H [001], including transition from square to rhombic to triangular between 3 T and H c2 = 5 T [54] . The vortex lattice transition is a gradual process. For instance, the apex angle of a unit cell for triangle vortex lattice is 60 • and that for a square lattice is 90 • . The evolution from triangular to the square lattices corresponds to the continuous change of the apex angle from 60 • to 90 • [55] . The smooth deformation of the vortex lattice, while the vortex lines are kept perpendicular to heat current, is therefore unlikely to give rise to a dramatic increases of κ. The Pauli pairing breaking effect, therefore, should be dominant in CeCoIn 5 at fields near H c2 .

While our modeling and experiments suggest the existence of the FFLO state when the field is orientated close to the c axis, two experimental observations deserve further attention. First, the region of enhanced κ near H c2 is wider than that identified as the FFLO region by NMR measurements [4] . However, a similar effect was observed for the field dependence of thermal conductivity in CeCoIn 5 at the SDW transition for field in the a-b plane [47] . There, the onset of the reduction in thermal conductivity at very low temperature was observed to take place at 9 T, whereas specific heat and neutron scattering measurements firmly place the SDW transition at 10 T. This may be due to fluctuations of the order parameter. Second, we are not able to resolve hysteresis in κ near H c2 when the direction of the field sweep is reversed. Such hysteresis is expected from the first order phase transition between the uniform superconducting state and the FFLO state. The lack of hysteresis may be due to a weak first order nature of the transition.

To summarize, combined modeling and thermalconductivity measurements suggest the existence of an FFLO state in CeCoIn 5 for field aligned along the c axis. The FFLO state competes with an SDW phase, and their relative stability can be tuned by rotating field in the a-c plane. Additional neutron scattering measurements with a scattering plane that includes the c axis may be promising in resolving the modulated susceptibility caused by the FFLO state. We expect the FFLO modulation wave vector to lie along the c axis when magnetic field is applied along the c axis. We note that the anomalous drop in the vortex lattice form factor, observed by the neutron scattering for field close to H c2 , was ascribed to the formation of an FFLO state [54] .

The authors would like to thank Takanori Taniguchi and Daniel Mazzone for helpful discussions. Computer resources for numerical calculations were supported by the Institutional Computing Program at LANL. The work by SZL was carried out under the auspices of the U.S. DOE Contract No. DE-AC52-06NA25396 through the LDRD program. Experimental work was supported by the LDRD program at LANL. * szl@lanl.gov † duk0@skku.edu

Superconductivity in a strong spin-exchange field

Inhomogeneous state of superconductors

Possible Fulde-Ferrell-Larkin-Ovchinnikov superconducting state in CeCoIn 5

Fulde-Ferrell-Larkin-Ovchinnikov state in a perpendicular field of quasi-twodimensional CeCoIn 5

Calorimetric evidence for a Fulde-Ferrell-Larkin-Ovchinnikov superconducting state in the layered organic superconductor κ−(BEDT−TTF) 2 Cu(NCS) 2

Zeeman-driven phase transition within the superconducting state of κ−(BEDT−TTF) 2 Cu(NCS) 2

Angle-dependent evolution of the Fulde-Ferrell-Larkin-Ovchinnikov state in an organic superconductor

Evidence of andreev bound states as a hallmark of the FFLO phase in κ−(BEDT−TTF) 2 Cu(NCS) 2

Microscopic study of the Fulde-Ferrell-Larkin-Ovchinnikov state in an all-organic superconductor

Evidence for the presence of the Fulde-Ferrell-Larkin-Ovchinnikov state in CeCu 2 Si 2 revealed using 63 Cu NMR

Thermodynamic evidence for the Fulde-Ferrell-Larkin-Ovchinnikov state in the KFe 2 As 2 superconductor

Magnetic torque studies on fflo phase in magnetic-field-induced organic superconductor λ-(BETS) 2 FeCl 4

Fulde-Ferrell-Larkin-Ovchinnikov superconductivity in the layered organic superconductor β − (BEDT − TTF) 4 [(H 3 O)Ga(C 2 O 4 ) 3 ]C 6 H 5 NO 2

Fulde-Ferrell-Larkin-Ovchinnikov State in Heavy Fermion Superconductors

From ising resonant fluctuations to static uniaxial order in antiferromagnetic and weakly superconducting CeCo(In 1−x Hg x ) 5 (x = 0.01)

Exploring the fragile antiferromagnetic superconducting phase in CeCoIn 5

Theory of an inherent spin-density-wave instability due to vortices in superconductors with strong pauli effects

Antiferromagnetic ordering induced by paramagnetic depairing in unconventional superconductors

Emergent antiferromagnetism in a d-wave superconductor with strong paramagnetic pair-breaking

Antiferromagnetic Order and π-Triplet Pairing in the Fulde-Ferrell-Larkin-Ovchinnikov State

Possible triplet superconducting order in a magnetic superconducting phase induced by paramagnetic pair breaking

Order parameter and vortices in the superconducting Q phase of CeCoIn 5

Theory for Coupled SDW and Superconducting Order in FFLO State of CeCoIn 5

Magnetic-field-induced pattern of coexisting condensates in the superconducting state of CeCoIn 5

Antiferromagnetic order in pauli-limited unconventional superconductors

Structure of magnetic order in pauli limited unconventional superconductors

Impurity-induced antiferromagnetic order in pauli-limited nodal superconductors: Application to heavyfermion CeCoIn 5

Field-induced spin-exciton condensation in the d x 2 -y 2 -wave superconductor CeCoIn 5

Spin resonance in the d-wave superconductor CeCoIn 5

Ising incommensurate spin resonance of CeCoIn 5 : A dynamical precursor of the Q phase

Magnetic field splitting of the spin resonance in CeCoIn 5

Avoided antiferromagnetic order and quantum critical point in CeCoIn 5

Field-induced quantum critical point in CeCoIn 5

Pressure study of quantum criticality in CeCoIn 5

Strong magnetic fluctuations in a superconducting state of CeCoIn 5

Anisotropic spin fluctuations in heavy-fermion superconductor CeCoIn 5 : In-NQR and Co-NMR studies

Theory of Superconductivity in CeMIn 5 (M=Co, Rh, Ir) on the Basis of the Three Dimensional Periodic Anderson Model

Bogoliubov-de Gennes Method and Its Applications

Field-tuned quantum critical point in cecoin 5 near the superconducting upper critical field

Structure of a vortex line in a d x 2 −y 2 superconductor

Impurity-induced magnetic droplet in unconventional superconductors near a magnetic instability: Application to Nd-doped CeCoIn 5

Superconductivity in narrow-band systems with local nonretarded attractive interactions

In the calculation with finite system size, there are several FFLO states with different wave-vector Q separated by the normal state at higher fields because of the commensuration of Q of the FFLO state with the system size

Fulde-Ferrell effect in type-II superconductors

Topological structure of a vortex in the Fulde-Ferrell-Larkin-Ovchinnikov state

Vortex state in a Fulde-Ferrell-Larkin-Ovchinnikov superconductor based on quasiclassical theory

Intertwined orders in heavy-fermion superconductor CeCoIn 5

The PWD is an induced order, and is of secondary importance in current discussions. The FFLO state is a dominant order when

Nonvanishing energy scales at the quantum critical point of CeCoIn 5

Quantum critical quasiparticle scattering within the superconducting state of CeCoIn 5

Quasiparticle thermal conductivity in the vortex state of high-T c cuprates

Unconventional superconductors under a rotating magnetic field. I. density of states and specific heat

Unconventional superconductors under a rotating magnetic field. II. thermal transport

Superconducting vortices in CeCoIn 5 : Toward the pauli-limiting field

Strain-induced intervortex interaction and vortex lattices in tetragonal superconductors