Local increase of symmetry on cooling in KNi2Se2

James R. Neilson,1,2∗ Natalia Drichko,2 Anna Llobet,3 Mali Balasubramanian,4 Matthew R. Suchomel,4 Tyrel M. McQueen1,2† 1Department of Chemistry, Johns Hopkins University, Baltimore, MD 2Department of Physics and Astronomy and the Institute for Quantum Matter, Johns Hopkins University, Baltimore, MD 3Lujan Neutron Scattering Center, Los Alamos National Laboratory, Los Alamos, NM 4Advanced Photon Source, Argonne National Laboratory, Argonne, IL ∗[email protected] and †[email protected]

KNi2Se2 exhibits an increase of symmetry on cooling below T ≤ 50 K, as observed by Raman spectroscopy and synchrotron x-ray diffraction. X-ray absorption spectroscopy confirms that the symmetry increase is due to changes in nickel-nickel interactions and suppression of charge density wave fluctuations. Density functional theory calculations reveal a zone- boundary lattice instability that provides a model of the room-temperature x-ray pair distribution function data, but fails to describe the higher local symmetry observed for T ≤ 50K. Together, these results support many- body correlation effects as drivers for the unusual heavy fermion electronic ground state in KNi2Se2.

Strong electron-electron and electron-phonon interactions give rise to broken symmetry states. A prototypical exam- ple is the charge density wave (CDW) state, which lowers the symmetry and introduces an additional periodic modula- tion of electron density.[1–3] CDWs belong to a larger fam- ily of spontaneous charge separation phenomena, such the charge ordering observed in the manganites [4, 5] and organic conductors.[6] In these materials, the formation of charge- separated states is almost always observed at low tempera- ture, when thermal fluctuations are no longer able to disrupt the electronic order. In rare situations, many-body interactions can produce a re-entrant transition in which the charge period- icity is destroyed and higher symmetry is restored on further cooling.[6, 7] In fact, a decrease of charge-order fluctuations on cooling has been observed in metallic organic conductors with 1/4-filled bands that are close to charge order; the ex- FIG. 1. (Color online) Temperature dependence of Raman spectra perimentally observed response is rather weak since the effect of KNi2Se2 illustrate an additional phonon mode then that expected results from the electronic repulsion within a single band.[8] from the I4/mmm spacegroup (arrow) which softens into the neigh- boring B (Ni) mode on cooling. Spectra are shifted along y axis for We have recently discovered that KNi2S2 [9] and KNi2Se2 1g [10, 11] have signatures of CDW fluctuations at room tem- clarity. Right: In-plane illustration of the increase of symmetry on cooling. Top inset: the average I4/mmm unit cell of KNi Se . perature that disappear below T ∼ 50 K. This surprising ob- 2 2 servation has stimulated extensive experimental [12–16] and theoretical [7, 17–20] studies aimed at understanding the un- usual nature of charge coupling with spin, orbital, and lattice degrees of freedom. the increase of local point group symmetry and loss of CDW Most strikingly, KNi2Se2 shows an increase of carrier mo- fluctuations in KNi2Se2 on cooling. An analysis of the struc- bility (µ) below T ∼ 50 K that occurs at the same temperature ture derived from high-resolution powder synchrotron x-ray at which CDW fluctuations disappear. The compound enters diffraction (SXRD) shows signatures of the CDW fluctua- a heavy fermion state below T = 20 K with an effective elec- tions, due to out-of-plane displacements of Ni atoms at room ∗ tronic mass of m ∼ 6 to 18mb, and superconductivity below temperature, that then disappear on cooling. X-ray absorption arXiv:1310.6828v1 [cond-mat.str-el] 25 Oct 2013 Tc=0.8 K.[10] Our studies showed that the magnetic response spectroscopy (XAS) and the extended x-ray absorption fine is rather weak and temperature independent.[10] An explana- structure (EXAFS) confirm that the CDW fluctuations arise tion of a heavy fermion state was proposed in terms the fluc- from displacements of the positions of nickel atoms. Den- tuation of a CDW involving electron-electron correlations.[7]. sity functional theory (DFT) calculations predict an intrinsic The theory suggests that the loss of CDW fluctuations at low lattice instability driven by electron-phonon coupling. A cal- temperature results from both many-body physics and the hy- culated imaginary frequency phonon mode at the zone bound- bridization between the localized, 1/4-filled band and a dis- ary qualitatively reproduces the room-temperature distortion persive conduction band. observed by synchrotron x-ray total scattering, but fails to In this letter, we determine that there is an increase in sym- describe the low temperature ground state of KNi2Se2. In metry on cooling in KNi2Se2, likely driven by electronic cor- agreement with the theoretical work on the canonical CDW relations. Raman spectroscopy provides direct evidence for materials,[2], the nature of the CDW fluctuations in this ma- 2 terial implies that the fluctuations are not driven by Fermi sur- face nesting, but are instead driven by local bonding interac- tions. Raman spectra were obtained from the (ab) plane of single micro-crystals in a parallel-sided fused silica capillary (evac- uated to pHe = 10 torr at RT) in a T64000 Jobin-Yvone spec- trometer equipped with an Olympus microscope. The laser power was below 0.3 mW to avoid sample heating, estimated by comparing Stokes and anti-Stokes spectra at T = 200 K. From a symmetry analysis of the low-temperature crystal structure presented in Ref. [10] (space group I4/mmm), we expect four Raman-active optical phonons.[21] In particular, Ni atoms in (4d Wyckoff positions, Ref. [11]) contribute B1g and Eg modes, of which only B1g is observed in the (ab) plane. Se atoms (4e Wyckoff positions, Ref. [11]) show A1g and Eg Raman active modes, of which only A1g is observed in the (ab) plane. Indeed, in the Raman spectra at the lowest tem- FIG. 2. (a) Unit cell volume of KNi2Se2 from high-resolution syn- peratures we observe two strong bands, located at 171 and chrotron powder x-ray diffraction (SPXRD, 11-BM, APS) with grad- ual change of slope near T ∼ 50 K (Dashed vertical line: T = 50 K. 219 cm−1 (21.2 and 27.2 meV; Fig. 1), which we assign to the Dotted line: aT 3 + bT + c). (b) Negative thermal expansion below B1g vibration of Ni atoms and the A1g vibration of Se atoms, T ∼ 18 K, measured on () cooling, (4) warming, ( ) re-cooling, respectively, based on our calculations and Ref. [12]. The low and (◦) re-warming, coincides with the increase in the carrier mobil- temperature spectrum of KNi2Se2 has a shape very similar to ity (µ, reproduced from Ref. [10]). (c) The displacement of the Ni the non-stoichiometric compound K0.95Ni1.86Se2, for which atom (z||c) illustrates the apparent increase in the Ni site symmetry no distortion is reported,[12] even though the frequencies in on cooling near T ∼ 50 K. The shaded area contours a 5% deviation the spectra of KNi2Se2 are higher. from the best Rietveld g.o.f. at each T. Solid lines guide the eye. On increasing the temperature above T > 50 K, the lower- Error bars fall within the symbols. (d) Fourier difference maps from frequency Raman band shifts to higher frequencies, and an- data collected at T = 300 K, viewed on the (110) plane, for nickel on the 4d (bottom) vs. 8g (top) Wyckoff positions. other phonon becomes well resolved at 186 cm−1. This shape of the spectra persists to room temperature. The appearance of an additional phonon mode for T > 50 K is a direct evi- dence for a decrease in local point symmetry of KNi2Se2 on transformed EXAFS signal (in real space) for the Ni K edge as warming. a function of temperature [Fig.3(a)], indicative of a change in At T = 50 K Rietveld analysis of the high-resolution SXRD the Ni coordination environment. In contrast, there is no shift data [9, 22–24] reveals a change in the average structure that in the Se K edge EXAFS [Fig.3(b)]. The EXAFS confirms gives the best description of the data [Fig 2(c)]. For T ≤ 50 K, that it is the nickel atoms that are responsible for the apparent the nickel atom is found at the 4d Wyckoff position, (0, 0.5, local increase of symmetry on cooling. 0.25). However, for T > 50 K, the goodness-of-fit (g.o.f.) is The EXAFS data were quantitatively modeled using the av- improved by displacing the Ni atom from this ideal position to erage I4/mmm crystal structure to simulate the path lengths the 8g Wyckoff position (0,0.5,z) with a 50% fractional occu- out to a cluster size of ∼ 4.5 A.˚ The quality of fit of the model, pancy, as illustrated by the T = 300 K Fourier difference maps as judged by the R-factor (R = [∑(data−fit)2]/[∑(data)2]), is in Fig. 2(d). This demonstrates an apparent decrease in the lo- temperature independent from 8.5 ≤ T ≤ 75K (∼ 0.7±0.1%); cal symmetry of Ni positions on heating above T > 50 K. In above that temperature the R-factor increases nearly linearly fact, if we describe the Ni atoms using the 8g Wyckoff posi- to 2.2% by T = 300 K. This reduced fit quality is also visi- tions in a I4/mmm unit cell, we would expect two Ni Raman- ble in Fig.3(c) and (d) and indicates that the average crystal active modes (A1g and B1g) in the (ab) plane, in agreement structure does not accurately describe the local Ni and Se en- with our experimental Raman spectra at T > 50 K. vironments for T > 75 K. Attempts to fit a distorted model to XAS data (beamline 20-BM-B [25, 26]) within the mea- the higher temperature EXAFS were unreliable. The fits of sured temperature range from T=300 to 8.5 K show no signif- the ideal undistorted structure provide physically meaningful 1.5+ 2− icant changes of the formal valence (Ni , Se ), as neither temperature dependences of the path length changes (r − r0) Ni nor Se K edges move in energy within the uncertainty of and mean-squared displacements (σ2). the XAS experiment (∆E ∼ 0.05 eV). This confirms our pre- The Ni-Se (Ni edge, Ni to Se path; 4× degenerate) and Se- vious speculation that the charge fluctuations are not due to Ni (Se edge, Se to Ni path; 4× degenerate) paths show very single-ion effects, but instead result from nickel-nickel inter- little temperature dependence [Fig.3(e)] and the σ2 exhibit actions. However, several different methods indicate slight the expected gradual increase on warming. The Se-Se path changes in the unit cell and details of crystal structure as (Se edge, Se to Se nearest neighbor) also does not show any a function of temperature, showing that the local symmetry temperature dependent inflections. In contrast, the Ni-Ni path change is due to the Ni atoms. (the 4× degenerate nearest neighbor nickel separation) has a There is a significant redistribution of the Fourier- significant temperature dependence, showing that the distance 3

FIG. 3. Temperature dependence of the EXAFS (20-BM-B, APS) of KNi2Se2 measured at (a) the Ni K edge and (b) the Se K edge show a shift in spectral weight with temperature at the Ni K edge FIG. 4. (Color online) (a) DFT calculated phonon dispersion curves predominantly around ∼2.6 A.˚ Selected fits of (c) the Ni K edge with a negative energy mode [~q = ( 1 0 0)]; the illustrated mode cor- and (d) the Se K edge EXAFS illustrate the better quality of fit for 2 responds to out-of-plane Ni displacements [Ni: blue (dark) circles; T ≤ 75 K. (e) The temperature dependence of the shifts of different Se: red (light) circles; amplitude exaggerated for clarity]. (b) Calcu- scattering path lengths illustrate that most of changes reside on the Ni lated vibrational densities of states (VDOS). (c) X-ray pair distribu- sublattice [paths: Ni-Se/Se-Ni r = 2.40 A;˚ Ni-Ni r = 2.76 A;˚ Se- 0 0 tion function (PDF; 11-ID-B, APS) of KNi Se (T = 300 K; filled Se and Ni- Ni(2) r = 3.91 A].˚ (f) The mean squared displacements 2 2 0 circles; reproduced and offset for clarity). Lines through the data of the path lengths indicate the same trend. Error bars are from the are the calculated PDF from the I4/mmm crystal structure (top), the statistical error of the fit. The vertical bars at r ∼ 2.1 A˚ indicate the imaginary frequency eigenvector from the DFT phonon calculation Ni-Se or Se-Ni path length, Ni-Ni paths at r ∼ 2.6 A,˚ and Se-Se paths (middle), and its least-squares relaxation within space group symme- at r ∼ 3.6 A.˚ try C2/m (bottom); differences are below. (d) Illustration of Ni-Se and Ni-Ni pairwise correlations. Insets illustrate the corresponding [Ni2Se2] plane (bonds denote nearest neighbors). between Ni atoms decreases on cooling from T = 300 K, sat- urating around T ∼ 50 K, which is consistent with the Ni atom position from Rietveld analysis. The next-nearest Ni-Ni sep- calculations lacks the ability to adequately describe the struc- aration [Ni-Ni(2); 4× degenerate] shows a sharp inflection at ture. However, this predicted lattice instability does resem- T = 75 K, consistent with the out-of-plane displacement of ble the distortion of KNi2Se2 observed above T ∼ 50 K. The the Ni position determined by X-ray diffraction. While the primary polarization components of the imaginary frequency significant error bars and increase in σ2 prevent an accurate eigenvector are out-of-plane Ni displacements, as illustrated assignment of the exact disposition of that path length at high in Fig. 4(a) [±z, with a small component along x, such that temperatures, it is clear that the microscopic Ni environment z ∼ 3x]. changes at T ∼ 75 K. To determine the validity of the distortion predicted by DFT calculations were carried out to gain insight to the na- DFT, X-ray total scattering experiments were performed to −1 ture of the observed structural distortions. A linear-response identify the distortion (APS, 11-ID-B, 90 ke, Qmax= 30 A˚ DFT calculation of the zone center phonons [27, 28] for the [31–34]) The DFT-predicted distortion (after scaling the am- low-temperature I4/mmm crystal structure of KNi2Se2 gener- plitude by a least-squares fit) does not directly improve the fit ally agrees with the results of the Raman scattering measure- to the x-ray PDF at room temperature where we see the ev- ments [29], as well as with the mode energies and assignments idence for CDW [Fig.4 (c) and (d); [35]]. However, using reported in Ref. [12]. Calculation of the in-plane phonon dis- ISODISTORT [36], the imaginary-frequency zone-boundary persion reveals an imaginary frequency eigenvector at the X phonon polarizations were decomposed into the subgroup 1 point in the Brillouin zone [~q = ( 2 0 0), E = 5.2i meV, Fig. 4]. symmetry, C2/m [37], which allowed for an excellent de- Imaginary frequency modes from DFT can be indicators of scription to the PDF out to r = 10 A˚ after a least-squares lattice instabilities, as exemplified by the localized, aperiodic refinement [Fig. 4(c) and (d)]. The resulting configuration distortions in Bi 2Ti2O7.[30] In the present case, this predic- contains short and long Ni-Ni bonds, akin to the long-range tion for a distorted ground state contrasts with the experimen- ordered CDW modulations in the structurally related com- tally observed undistorted ground state at the lowest tempera- pounds, SrRh2As2[38] and KCu2Se2.[39] While the bond dis- tures measured (see Figures 2 and 3, and Ref. [10]). This con- tances from the room temperature PDF can be described by tradiction suggests that the level of theory used for the DFT a C2/m distortion, the overall crystal symmetry determined 4 by high-resolution SXRD remains I4/mmm. While the space lowering of local symmetry at high temperatures. Generally, group setting that we used with C2/m should have more a raising of symmetry on cooling is unusual, as it naively modes visible in the Raman spectrum, the true atomistic con- suggests an increase in configurational entropy, when con- figuration may possess pseudo-symmetry, such that the local sidering only the degeneracy of lattice modes. Hybridization point group symmetries do not provide additional modes than of localized orbitals with the conduction band, as previously beyond what is observed. The Ni atoms exhibit the most sig- proposed,[7] may provide a state in which the total entropy is nificant distortions. Therefore, we expect that only the Ni decreased by such a raising of lattice symmetry and explain point group symmetry is significantly changed with temper- why our DFT calculations fail to correctly predict the ground ature, and the distortions are oriented out of the plane. These state of this material. Such interactions between localized and results are in agreement with the Rietveld analysis of SXRD delocalized electronic systems have been proposed as an ori- data [Fig.2(c)]and Raman scattering. gin of nematicity, coherence/incoherence crossover and heavy electrons in the iron-based superconductors,[6, 16, 17, 41– An additional piece of evidence for the importance of elec- 44] . The work presented here unambiguously demonstrates tronic correlations is an increase in unit cell volume on cool- the presence and importance of charge-lattice coupling in ing below T ∼ 20 K, yielding a negative coefficient of thermal ThCr Si -type materials. expansion [Fig. 2(a,b)]. This result is highly reproducible on 2 2 We thank James Murray and Zlatko Tesanovic for inspira- several heating and cooling cycles, akin to KNi S .[9] A neg- 2 2 tion and motivation, as well as Igor Mazin for helpful discus- ative coefficient of thermal expansion can be a signature of sions and Kevin Bayer and Karena Chapman for data collec- the formation of a heavy fermion state (confer CeAl [40]). 3 tion at 11-ID-B. This research and the Raman spectrometer The observed temperature dependence of the unit cell coin- are principally supported by the US Department of Energy cides with changes in the carrier mobility and suggests that (DoE), Office of Science, Office of Basic Energy Sciences the structure and electronic properties are highly correlated. (BES), Division of Materials Sciences and Engineering under In conclusion, this letter shows a direct evidence of a de- Award DE-FG02-08ER46544, as well as start-up funds from crease of local symmetry in KNi2Se2 with increasing temper- the Johns Hopkins University and funding from the David and ature. Out-of plane Ni atom displacements are a signature of Lucile Packard Foundation. This research has benefited from CDW fluctuations driven by local Ni-Ni interactions, rather the use of beamlines 11-BM-B, 11-ID-B, and 20-BM-B at than Fermi surface nesting, that do not exhibit any long range the Advanced Photon Source at Argonne National Laboratory, order, and do not reduce the average I4/mmm symmetry of supported by the U.S. DoE, Office of Science, Office of Basic the unit cell. Raman spectroscopy directly demonstrates the Energy Sciences, under Contract No. DE-AC02-06CH11357.

[1] G. Gruner, Density waves in solids (Addison-Wesley Publishing [13] H. Lei, K. Wang, H. Ryu, D. Graf, J. B. Warren, and C. Petro- Company, Advanced Book Program, 1994). vic, arXiv , 1211.1371 (2012). [2] M. D. Johannes and I. I. Mazin, Phys. Rev. B 77, 165135 [14] H. Wang, C. Dong, Q. Mao, R. Khan, X. Zhou, Z. Li, B. Chen, (2008). J. Yang, Q. Su, and M. Fang, arXiv (2013). [3] H.-M. Eiter, M. Lavagnini, R. Hackl, E. A. Nowadnick, A. F. [15] H. Wang, C. Dong, Q. Mao, R. Khan, X. Zhou, C. Li, B. Chen, Kemper, T. P. Devereaux, J.-H. Chu, J. G. Analytis, I. R. Fisher, J. Yang, and M. Fang, arXiv , 1305.1033 (2013). and L. Degiorgi, Proc. Natl. Acad. Sci. USA 110, 64 (2013). [16] F. Hardy, A. E. Boehmer, D. Aoki, P. Burger, T. Wolf, [4] J. Mitchell, D. Argyriou, A. Berger, K. Gray, R. Osborn, and P. Schweiss, R. Heid, P. Adelmann, Y. X. Yao, G. Kotliar, U. Welp, J. Phys. Chem. B 105, 10731 (2001). J. Schmalian, and C. Meingast, arXiv , 1302.1696 (2013). [5] E. Bozin, M. Schmidt, A. DeConinck, G. Paglia, J. Mitchell, [17] V. Stanev and P. B. Littlewood, Phys. Rev. B 87, 161122 (2013). T. Chatterji, P. Radaelli, T. Proffen, and S. J. L. Billinge, Phys. [18] V. Bannikov and A. Ivanovskii, Physica B: Cond. Matt. 418, 76 Rev. Lett 98, 137203 (2007). (2013). [6] J. Merino and R. H. McKenzie, Phys. Rev. Lett. 87, 237002 [19] F. Lu, J. Z. Zhao, and W.-H. Wang, J. Phys.: Condens. Matter (2001). 24, 495501 (2012). [7] J. M. Murray and Z. Tesanoviˇ c,´ Phys. Rev. B 87, 081103 [20] V. Bannikov and A. Ivanovskii, Physica C: Superconductivity (2013). (2013). [8] J. Merino, A. Greco, N. Drichko, and M. Dressel, Phys. Rev. [21] D. Orobengoa, C. Capillas, M. I. Aroyo, and J. M. Perez-Mato, Lett. 96, 216402 (2006). J. Appl. Crystallogr. 42, 820 (2009). [9] J. R. Neilson, T. M. McQueen, A. Llobet, J. Wen, and M. R. [22] J. Wang, B. H. Toby, P. L. Lee, L. Ribaud, S. M. Antao, Suchomel, Phys. Rev. B 87, 045124 (2013). C. Kurtz, M. Ramanathan, R. B. V. Dreele, and M. A. Beno, [10] J. R. Neilson, A. Llobet, A. V. Stier, L. Wu, J. Wen, J. Tao, Rev. Sci. Instrum. 79, 085105 (2008). Y. Zhu, Z. B. Tesanovic, N. P. Armitage, and T. M. McQueen, [23] A. C. Larson and R. B. V. Dreele, Los Alamos National Labo- Phys. Rev. B 86, 054512 (2012). ratory Report LAUR (2000). [11] J. R. Neilson and T. M. McQueen, J. Am. Chem. Soc 134, 7750 [24] B. H. Toby, J. Appl. Crystallogr. 34, 210 (2001). (2012). [25] B. Ravel and M. Newville, J. Synch. Rad. 12, 537 (2005). [12] N. Lazarevic,´ M. Radonjic,´ M. Sˇcepanovi´ c,´ H. Lei, [26] M. Newville, Journal of Synchrotron Radiation 8, 322 (2001). D. Tanaskovic,´ C. Petrovic, and Z. V. Popovic,´ Phys. Rev. B [27] “ELK, http://elk.sourceforge.net/,” . 87, 144305 (2013). [28] J. P. Perdew, A. Ruzsinszky, G. I. Csonka, O. A. Vydrov, G. E. 5

Scuseria, L. A. Constantin, X. Zhou, and K. Burke, Phys. Rev. [37] See EPAPS, Table II. Lett. 100, 136406 (2008). [38] V. Zinth, V. Petricek, M. Dusek, and D. Johrendt, Phys. Rev. B [29] See EPAPS for calculation details and mode energies. 85, 014109 (2012). [30] D. P. Shoemaker, R. Seshadri, A. L. Hector, A. Llobet, T. Prof- [39] O. Tiedje, E. E. Krasovskii, W. Schattke, P. Stoll, C. Nather,¨ fen, and C. J. Fennie, Phys. Rev. B 81, 144113 (2010). and W. Bensch, Phys. Rev. B 67, 134105 (2003). [31] P. J. Chupas, X. Qiu, J. Hanson, P. Lee, C. Grey, and S. J. L. [40] K. Andres, J. E. Graebner, and H. R. Ott, Phys. Rev. Lett. 35, Billinge, J. Appl. Crystallogr. 36, 1342 (2003). 1779 (1975). [32] A. Hammersley, S. Svensson, M. Hanfland, A. Fitch, and [41] R. H. McKenzie, J. Merino, J. B. Marston, and O. P. Sushkov, D. Hausermann, High Pressure Research 14, 235 (1996). Phys. Rev. B 64, 085109 (2001). [33] X. Qiu, J. W. Thompson, and S. J. L. Billinge, J. Appl. Crys- [42] Z. Yin, A. Kutepov, and G. Kotliar, Phys. Rev. X 3, 021011 tallogr. 37, 678 (2004). (2013). [34] C. L. Farrow, P. Juhas, J. W. Liu, D. Bryndin, E. S. Bozin, [43] J. M. Caron, J. R. Neilson, D. C. Miller, A. Llobet, and T. M. J. Bloch, T. Proffen, and S. J. L. Billinge, J. Phys.: Condens. McQueen, Phys. Rev. B 84, 180409(R) (2011). Matter 19, 335219 (2007). [44] J. M. Caron, J. R. Neilson, D. C. Miller, K. Arpino, A. Llobet, [35] See EPAPS, Table II. and T. M. McQueen, Phys. Rev. B 85, 180405(R) (2012). [36] B. J. Campbell, H. T. Stokes, D. E. Tanner, and D. M. Hatch, J. Appl. Crystallogr. 39, 607 (2006).