preprint(June 28, 2021)

Time-reversal symmetry-breaking charge order in a correlated kagome superconductor

C. Mielke III,1, 2, ∗ D. Das,1, ∗ J.-X. Yin,3, ∗ H. Liu,4, 5, ∗ R. Gupta,1 C.N. Wang,1 Y.-X. Jiang,3 M. Medarde,6 X. Wu,7 H.C. Lei,8 J.J. Chang,2 P. Dai,9 Q. Si,9 H. Miao,10 R. Thomale,11, 12 T. Neupert,2 Y. Shi,4, 5, 13 R. Khasanov,1 M.Z. Hasan,3, 14, 15, 16, † H. Luetkens,1 and Z. Guguchia1, ‡ 1Laboratory for Muon Spin Spectroscopy, Paul Scherrer Institute, CH-5232 Villigen PSI, Switzerland 2Physik-Institut, Universit¨atZ¨urich,Winterthurerstrasse 190, CH-8057 Z¨urich,Switzerland 3Laboratory for Topological Quantum Matter and Advanced Spectroscopy (B7), Department of Physics, Princeton University, Princeton, New Jersey 08544, USA 4Beijing National Laboratory for Condensed Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100190, China 5School of Physical Sciences, University of Chinese Academy of Sciences, Beijing 100190, China 6Laboratory for Multiscale Materials Experiments, Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland 7Max-Planck-Institut fur Festkorperforschung, Heisenbergstrasse 1, D-70569 Stuttgart, Germany 8Department of Physics and Beijing Key Laboratory of Opto-electronic Functional Materials and Micro-nano Devices, Renmin University of China, Beijing 100872, China 9Department of Physics and Astronomy, Rice University, Houston, Texas 77005, USA 10Material Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA 11Institut fur Theoretische Physik und Astrophysik, Universitat Wurzburg, 97074 Wurzburg, Germany 12Department of Physics, Indian Institute of Technology Madras, Chennai 600036, India 13Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, China 14Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08540, USA 15Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA 16Quantum Science Center, Oak Ridge, Tennessee 37831, USA

The kagome lattice, which is composed of a net- der intertwining with unconventional supercon- work of corner-sharing triangles, is a structural ductivity in the correlated kagome lattice. motif in quantum physics first recognized more 1 The observation of orbital currents is a long stand- than seventy years ago . It has been gradually ing motive in both topological and correlated quan- realized that materials which host such special tum matter. They have been suggested to produce lattice structures can exhibit quantum diversity, the quantum anomalous Hall effect when interacting ranging from spin-liquid phases, topological mat- 25 2–13 with Dirac fermions in a honeycomb lattice (Fig. 1a) ter to intertwined orders . Recently, charge and as the hidden phase of high-temperature cuprate sensitive probes have suggested that the kagome 26,27 14–16 superconductors (Fig. 1b). In both cases, orbital superconductors AV3Sb5 (A = K, Rb, Cs) 17–24 currents run through the lattice and break time-reversal exhibit unconventional chiral charge order , symmetry. Recently, the tantalizing visualization of which is analogous to the long-sought-after quan- such exotic order has been reported17–19 in the kagome tum order in the Haldane model25 or Varma 26 superconductor AV3Sb5 (A = K, Rb, Cs) (Fig. 1c). model . However, direct evidence for the time- Scanning tunneling microscopy observes a chiral 2×2 reversal symmetry-breaking of the charge order charge order (Fig. 1d) with an unusual magnetic field re- remains elusive. Here we utilize state-of-the-art sponse. Theoretical analysis17–24 also suggests that this muon spin relaxation to probe the kagome charge chiral charge order can not only lead to a giant anoma- order and superconductivity in KV3Sb5. We ob- lous Hall effect28 but also be a precursor of unconven- serve a striking enhancement of the internal field tional superconductivity23,29. However, the broken time- width sensed by the muon ensemble, which takes reversal symmetry nature of the charge order and its in- place just below the charge ordering temperature terplay with superconductivity has not been explicitly arXiv:2106.13443v1 [cond-mat.mtrl-sci] 25 Jun 2021 and persists into the superconducting state. Re- demonstrated by experiments. markably, the muon spin relaxation rate below the charge ordering temperature is substantially To explore unconventional aspects of supercon- enhanced by applying an external magnetic field. ductivity and the possible time-reversal symmetry We further show the multigap nature of super- breaking nature of charge order and superconductivity −2 in KV Sb , it is critical to measure the superconducting conductivity in KV3Sb5 and that the Tc/λab ratio 3 5 is comparable to those of unconventional high- order parameter and weak internal fields of KV3Sb5 temperature superconductors. Our results point on the microscopic level. Thus, we concentrate on to time-reversal symmetry breaking charge or- muon spin relaxation/rotation (µSR) experiments30,31 of the normal state depolarization rate and the magnetic 2

a Haldane Model b Varma Model c Orbital currents d in graphene in cuprates in kagome lattice

Cu O

1nm Low High

e Electronic clock f Positron detectors

3

�+ 1 Muon spin 2

Sample Muon detector 4

g h i

Figure 1: Indication of time-reversal symmetry-breaking of the charge order in KV3Sb5. a, Orbital currents (red arrows) proposed in a honeycomb lattice. b, Orbital currents (red arrows) proposed in the CuO2 lattice of cuprates. c, Schematic of the orbital currents (red arrows) in the kagome lattice. d, Scanning tunneling microscopy of the Sb surface showing 2×2 charge order as illustrated by black lines. The inset is the Fourier transform of this image, shown lattice Bragg peaks marked by blue circles and 2×2 vector peaks marked by red circles. The three pairs of 2×2 vector peaks feature different intensities, denoting a chirality of the charge order. e, A schematic overview of the experimental setup. Spin polarized muons ◦ with spin Sµ, forming 60 with respect to the c-axis of the crystal, are implanted in the sample. The sample was surrounded by four detectors: Forward (1), Backward (2), Up (3), and Down (4). An electronic clock is started at the time the muon passes the muon detector and is stopped as soon as the decay positron is detected in the positron detectors. f, The ZF µSR time spectra for KV3Sb5, obtained at different temperatures, all above the superconducting transition temperature Tc. The solid black curves in panel a represent fits to the recorded time spectra, using the Eq. 1. g, The temperature dependences of the relaxation rates ∆ and Γ, obtained in a wide temperature range. h, The temperature dependence of Γ from two sets of detectors across the charge ordering temperature T ∗ ' 80 K. i, Temperature dependences of the muon spin relaxation rates ∆ and Γ, which can be related to the nuclear and electronic system respectively, in the temperature range across Tc.

penetration depth λ in KV3Sb5. Importantly, zero-field netic fields as small as 0.1 G without applying external µSR (ZF-µSR) has the ability to detect internal mag- magnetic fields, making it a highly valuable tool for 3 probing spontaneous magnetic fields due to time-reversal We note that a similar value of internal magnetic field symmetry breaking32–36 within the superconducting strength is reported in several time-reversal symmetry- and charge ordered states. Additionally, λ is one of breaking superconductors32–34,36. Increase of the inter- the fundamental parameters of a superconductor, since nal field width visible from the electronic contribution to 2 it is related to the superfluid density ns via 1/λ = the ZF-µSR relaxation rate corresponds to an anomaly 2 ∗ ∗ µ0e ns/m (where m is the effective mass). seen also in the nuclear contribution to the relaxation rate ∆12; namely, a peak coinciding with the onset of While long-range magnetism has not been reported in the charge order, which decreases to a broad minimum 37 KV3Sb5 , zero-field µSR experiments have been carried before increasing again towards lower temperatures. out above and below T to search for any weak magnetism c In order to substantiate the above zero-field µSR re- (static or slowly fluctuating). A schematic overview of sults, systematic high field µSR39 experiments were car- the experimental setup with the muon spin forming 60◦ ried out (Fig. 2a). A non-monotonous behavior of the with respect to the c-axis of the crystal is shown in the relaxation rate is clearly seen in the µSR data, measured inset of Figure 1e. The sample was surrounded by four in magnetic field of 1 T, applied parallel to the c-axis, as detectors: Forward (1), Backward (2), Up (3), and Down (4). Figure 1f displays the zero-field µSR spectra from shown in Fig. 2b. We note that in high magnetic fields detectors 3 & 4 collected over a wide temperature range. we cannot discriminate between the nuclear and the elec- We see that the muon spin relaxation shows a reason- tronic contribution to the relaxation rate. Since the data able temperature dependence. Since the zero-field relax- at 1 T looks similar to the temperature dependence of ation is decoupled by the application of a small external the zero-field nuclear rate ∆12, it seems to be dominated magnetic field applied longitudinal to the muon spin po- by the nuclear response. However, at higher fields such as 3 T, 5 T and 8 T, the rate not only shows a broad larization, BLF = 50 G (see Fig. 1f), the relaxation is ∗ therefore due to spontaneous fields which are static on bump around T , but also shows a clear and stronger the microsecond timescale. The zero-field µSR spectra increase towards low temperatures within the charge or- were fitted using the gaussian Kubo-Toyabe depolariza- dered state, similar to the behavior observed for the re- tion function38, which reflects the field distribution at laxation rates Γ12 and Γ34 in zero-field. As the nuclear the muon site created by the nuclear moments of the contribution to the relaxation cannot be enhanced by an sample, multiplied by an additional exponential exp(-Γt) external field, this indicates that the low-temperature re- term (see also the Supplementary information): laxation rate in magnetic fields higher than 1 T is dom- inated by the electronic contribution. Remarkably, we find a factor of six enhancement of the relaxation rate 1 2 h ∆2t2 i ∗ P GKT (t) = + (1 − ∆2t2) exp − exp(−Γt) with the onset of T in 8 T compared to zero-field, show- ZF 3 3 2 ing a strong field-induced enhancement of the electronic (1) response. Moreover, we find that the magnitude of the where ∆/γµ is the width of the local field distribution Knight shift, defined as Kexp = (Bint − Bext)/Bext (Bint due to the nuclear moments and γµ/2π = 135.5 MHz/T and Bext are the internally measured and externally ap- is the muon gyromagnetic ratio. A gaussian Kubo Toy- plied magnetic fields) and obtained in 3 T, 5 T and 8 T abe shape is expected due to the presence of the dense shows a sharp increase just below T ∗, as shown in Fig- system of nuclear moments with large values of nuclear ure 3c. This indicates the enhancement of local magnetic spins (I = 3/2 for 39K, I = 7/2 for 51V, and I = 5/2 susceptibility concurrent with the charge order. We also 121 for Sb) in KV3Sb5 and a high natural abundance. We note that the increase of the relaxation rate arises from therefore conclude that the additional exponential relax- nearly the entire sample volume (see the supplementary ation rate ΓZF is electronic in origin. In Fig. 1g, we information), indicating the bulk nature of the transition see the temperature dependence of both the muon spin below T ∗. Such observation is consistent with charge sen- relaxation ∆12 and Γ12,34 due to nuclear and electronic sitive probe results that the magnetic field switching of moments, respectively, over a broad temperature range the chiral charge order is observed in the impurity-free from the base temperature to 300 K. There is a notewor- region17. Both the µSR and charge probe results attest thy increase immediately visible in the relaxation rates to the intrinsic nature of the time-reversal breaking in the Γ12 and Γ34 upon lowering the temperature below the kagome lattice. One plausible phenomenological scenario charge ordering temperature T ∗, which is better visible in is that the charge order has a complex chiral order param- Fig. 1h. This observation indicates the enhanced spread eter, which exhibits a phase difference between three sub- of internal fields sensed by the muon ensemble concurrent lattices of the kagome plane. The existence of phase dif- with the onset of charge ordering. The enhanced mag- ference, if not π, breaks time-reversal symmetry. Recent netic response that sets in with the charge order persists theoretical modeling of the charge ordering in the kagome all the way down to the base temperature, and remains lattice at van Hove filling and with extended Coulomb in- constant across the superconducting transition, as seen teractions (that is close to the condition of the AV3Sb5) in Fig. 1i. The increase in the exponential relaxation also suggests that time-reversal symmetry-broken charge below T ∗ is estimated to be ' 0.025 µs−1, which corre- order with orbital currents is energetically favorable21–24. sponds to a characteristic field strength Γ12/γµ ' 0.3 G. The suggested orbital current was reported to be homo- 4

Positron detectors T* a b 0.45 c -0.011 T* KV3Sb5 )

1 0.40 -0.012 - 1 T s ) µ 3 T (

% -0.013 ( e

0.35 5 T Sample t t f a i

B r

8 T h -0.014 s n

t o 0.30 i t h a g -0.015 i x n a l

0.25 K Muon spin e -0.016 3 T R 5 T 0.20 -0.017 8 T 0 50 100 150 200 250 300 0 50 100 150 200 250 300 T (K) T (K) Figure 2: Enhanced magnetic response of the charge order with applying external magnetic fields. a, Schematic overview of the high field µSR experimental setup for the muon spin forming 90◦ with respect to the c-axis of the crystal. The sample was surrounded by 2 times 8 positron detectors, arranged in rings. The specimen was mounted in a He gas-flow cryostat with the largest face perpendicular to the muon beam direction, along which the external field was applied. Behind the sample lies a veto counter (in orange) which rejects the muons that do not hit the sample. The temperature dependence of b, the muon spin relaxation rate and c, the Knight shift Kexp (local susceptibility) for KV3Sb5, measured under the c-axis magnetic fields of µ0H = 1 T, 3 T, 5 T, and 8 T. The vertical grey lines mark the charge ordering temperature, determined from magnetization measurements (see the supplementary information). Kexp also shows a peak at around 15 K, which may point to some transition/crossover deep in the charge order state.

p 2 2 geneous on the lattice, however alternating in its flow, the total gaussian relaxation rate σtotal = σSC + σnm which would produce inhomogeneous fields at the muon for KV3Sb5 for the 5 mT and 10 mT fields applied both site. Within this framework, muons may couple to the within and out-of the kagome plane. In order to extract ∗ closed current orbits below T , leading to an enhanced the σSC contribution due only to superconductivity, the internal field width sensed by the muon ensemble concur- average value of the normal state depolarization rate σnm rent with the charge order. Thus, we conclude that the estimated from six temperature points just above the on- current results provide key evidence for a time-reversal set of the superconducting transition has been quadrati- symmetry broken charge order in KV3Sb5. These re- cally subtracted, since above Tc there is only the normal sults also indicate that the magnetic and charge chan- state contribution to σtotal. Figure 3c shows the tem- nels of KV3Sb5 appear to be strongly intertwined, which perature dependences of the superconducting relaxation can give rise to complex and collective phenomena. The rates σSC,ab and σSC,ac, determined from the data with time-reversal symmetry breaking charge order can open the field applied along the c-axis and within the kagome a topological gap on the Dirac nodal lines at the Fermi plane, respectively. The c-axis relaxation rate can be ex- 2 40 level, introducing large anomalous Hall effect. It can also tracted as σSC,c = σSC,ac/σSC,ab , which is shown as a be a strong precursor of unconventional superconductiv- function of temperature in Fig. 3c. ity as we study below. We note that the magnetic penetration depth λ(T ) The time-reversal symmetry-breaking charge order can (right axis of Fig. 3c) is related to the relaxation rate arise from extended Coulomb interactions of the kagome σSC(T ) in the presence of a triangular (or hexagonal) lattice with van Hove singularities, where the same in- vortex lattice by the equation41: teractions and instabilities can lead to correlated super- conductivity. Thus, we next focus on the low transverse- σSC (T ) Φ0 = 0.06091 2 , (2) field µSR measurements performed in the superconduct- γµ λ (T ) ing state. With a superconducting transition tempera- ture Tc of ' 1.1 K, the TF-µSR spectra above (1.25 K) where γµ is the gyromagnetic ratio of the muon and Φ0 and below (0.25 K) the superconducting transition tem- is the magnetic-flux quantum. Since the applied field is perature Tc are shown in Fig. 3a. Above Tc, the oscilla- a factor of 20-60 times smaller than the critical magnetic tions show a damping essentially due to the random local field in KV3Sb5, the Eq. 2 is valid to estimate both the fields from the nuclear magnetic moments. The damp- λab and λc. The value of the in-plane penetration depth ing rate is shown to be nearly constant between 10 K λab at 0.3 K, determined from σSC,ab (superconducting and 1.25 K. Below Tc the damping rate increases with screening currents flowing parallel to the kagome plane), decreasing temperature due to the presence of a nonuni- is found to be λab ' 877(20) nm. The value of the out- form local magnetic field distribution as a result of the of-plane penetration depth, determined from σSC,c (su- formation of a flux-line lattice in the superconducting perconducting screening currents flowing perpendicular state. Figure 3b depicts the temperature evolution of to the kagome plane), is found to be λc ' 730(20) nm. 5

a c ) 1 T - c

s 0.20 2.0

KV Sb 0.27 K µ 3 5 ( s

1.0 SC,ab c ) 1.2 K , 2 C s -

S SC,ac m s 0.15 1.5

µ , s

0.5 (

) SC,c

1 - 2 - s c ) t

µ µ H = 5 mT l

( 0 (

, 4 c 0.0 0.10 1.0 ) 3 a 2 , - P C S m s µ

( ,

-0.5 ) 2 1 -

- 0.05 0.5 b s a µ l 5 mT - H || kagome plane (

-1.0 b a

Transverse field ,

C 0.00 0.0 S

0 2 4 6 8 s 0.0 0.5 1.0 1.5 2.0 t (µs) T (K) b 0.32 d 5 mT - H || c-axis ) T -MoTe s

1 d 2 e

- t 5 mT - H || kagome plane s

s a 0.30 2H-NbSe r e 2 t µ 10 mT - H || kagome plane p

( a T u c 10 r l LaRu Si 3 2 c p

a t 0.28 d u o

t c KV3Sb5 e p s d

) o e , Nb 0.26 d p K - e

( o t

e Pb l d a c o - r n

Ta T H 0.24 o

n r T t

o c c i

t e 0.22 l a E x Tc a l 0.20 e

R 1 KV3Sb5 0.18 0 2 4 6 8 10 0.1 1 10 100 1000 T (K) -2 -2 lab (T=0.3K) (µm )

Figure 3: Correlated kagome superconductivity. a, The transverse field µSR spectra are obtained above and below Tc 6 (after field cooling the sample from above Tc). Error bars are the standard error of the mean (s.e.m.) in about 10 events. The error of each bin count n is given by the standard deviation (s.d.) of n. The errors of each bin in A(t) are then calculated by s.e. propagation. The solid lines in panel a represent fits to the data by means of Eq. 3. The dashed lines are a guide to the eye. b, The temperature dependence of the total muon spin relaxation rate σtotal measured measured in the magnetic fields of 5 mT and 10 mT applied both parallel to the c-axis and parallel to the kagome plane. The dashed lines mark the average value of σtotal estimated from few data points above Tc. c, The superconducting muon depolarization rates σSC,ab, σSC,ac, σSC,c as −2 −2 well as the inverse squared magnetic penetration depth λab and λc as a function of temperature, measured in 5 mT, applied −2 parallel and perpendicular to the kagome plane. d, Plot of Tc versus the λab (0) obtained from our µSR experiments in KV3Sb5. The dashed red line represents the relation obtained for kagome superconductor LaRu3Si2 as well as for the layered transition 31 46 metal dichalcogenide superconductors Td-MoTe2 and 2H-NbSe2 . The relation observed for underdoped cuprates is also shown (solid line for hole doping45 and the dashed black line for electron doping47,48). The points for various conventional Bardeen-Cooper-Schrieffer superconductors are also shown.

17 The λ(T ) in the applied field of 5 mT shows a well pro- by dxz/dyz orbitals of V . Different observed supercon- nounced two step behavior. This suggests that at least ducting gaps are likely associated with Sb pz and V 3d two bands are involved in superconductivity and that the Fermi surfaces. Quantitative analysis of λ(T ) and the interband coupling is nearly zero which would result in discussions on the superconducting gap symmetry31,44 42,43 different values of Tc for different bands . The multi- are shown in the Supplementary information. gap superconductivity is consistent with the presence of To place the system KV3Sb5 in the context of other multiple Fermi surfaces revealed by electronic structure superconductors, in Fig. 3d we plot the critical tempera- calculations and tunneling measurements19. The low en- −2 ture Tc against the superfluid density λ . Most uncon- ergy band structure of KV Sb consists of three nearly ab 3 5 ventional superconductors have T /λ−2 values of about independent features: a quasi-two-dimensional electron c ab 0.1−20, whereas all of the conventional Bardeen-Cooper- pocket around the Γ point formed by p orbitals from z Schrieffer (BCS) superconductors lie on the far right in Sb; a Fermi surface exhibiting a van Hove singularity the plot, with much smaller ratios45. In the other words, at the M point, formed by the d orbitals of V; and a xy unconventional superconductors are characterized by a pair of Dirac-cone like bands near the M point, formed dilute superfluid (low density of Cooper pairs) while con- 6 ventional BCS superconductors exhibit a dense super- and high-field HAL-9500 instruments at the Swiss Muon fluid. Moreover, a linear relationship between Tc and Source (SµS) at the Paul Scherrer Institut, in Villigen, −2 λab is expected only on the Bose Einstein Condensate Switzerland. Zero field is dynamically obtained (compen- (BEC)-like side of the phase diagram and is considered a sation better than 30 mG) by a newly installed automatic hallmark feature of unconventional superconductivity45, compensation device49. When performing measurements where (on-site or extended) Coulomb interactions plays in zero-field the earth magnetic field or any stray fields −2 a role. For KV3Sb5, the ratio is estimated to be Tc/λab are tabulated and automatically compensated by the au- ' 0.7, which is far away from conventional BCS super- tomatic compensation device. conductors and approximately a factor of two greater Because the KV3Sb5 samples were rather thin (0.13 than that of charge density wave superconductors 2H- - 0.35 mm along c-axis), a mosaic of several crystals NbSe2 and 4H-NbSe2 as well as Weyl-superconductor stacked on top of each other was used for these mea- 31 Td-MoTe2 and kagome superconductor nLaRu3Si2, as surements. The individual crystals were attached to a shown in Fig. 3d. The point for KV3Sb5 is close to the 25 µm thick Cu foil mounted on the Cu sample fork, electron-doped cuprates, which are well-known correlated and the entire ensemble was held together by small superconductors with poorly screened Coulomb interac- droplets of GE varnish. The crystals were aligned with tions. the same in-plane orientation, which was achievable Our work points to a time-reversal symmetry-breaking because the thin sheet-like crystals grow with a clear charge order, intertwined with correlated superconduc- hexagonal shape (see Supplementary Figure 1c). These tivity in the kagome superconductor KV3Sb5. While low- multilayer crystal mosaics were then wrapped in a temperature time-reversal symmetry-breaking supercon- single layer of 60 µm PE polyester tape. The magnetic ductivity has been discussed for many systems, the high- field was applied both in-plane (along the ab-plane) temperature time-reversal symmetry-breaking charge or- and out-of-plane (along the crystallographic c-axis). der is extremely rare, and finds a direct comparison with A schematic overview of the experimental setup for the fundamental Haldane and Varma models. The com- zero-field and low transverse field measurements in GPS plex intertwining of such a charge ordered state with cor- is shown in Figure 2a. The muon spin is forming 60◦ related superconductivity highlights the quantum diver- with respect to the c-axis of the crystal. The sample was sity of the correlated kagome lattice and hints for other surrounded by four detectors: Forward (1), Backward hitherto unknown hybrid phenomenon resulting from the (2), Up (3), and Down (4). A schematic overview of kagome quantum diversity. the experimental setup for high-field µSR instrument is shown in the inset of Figure 3a. The crystals were mounted on 10 mm circular silver sample holder by I. METHODS small droplets of GE varnish. The muon spin forms 90◦ with respect to the c-axis of the crystal. The sample was surrounded by 2 × 8 positron detectors, arranged Sample preparation: Single crystals of KV3Sb5 in rings. The specimen was mounted in a He gas-flow were grown from KSb2 alloy as a flux. K, V, Sb elements cryostat with the largest face perpendicular to the and KSb2 precursor were sealed in a Ta crucible in a molar ratio of 1:3:14:10, which was finally sealed in a muon beam direction, along which the external field was highly evacuated quartz tube. The tube was heated up applied. to 1273 K, dwelt for 20 hours and then slowly cooled down to 773 K. Single crystals were separated from µSR experiment: In a µSR experiment nearly 100% the flux by centrifuging. Crystals obtained from flux spin-polarized muons µ+ are implanted into the sample are thin, metallic platelets with high luster and the one at a time. The positively charged µ+ thermalize largest size of approximately 5 mm × 5 mm. The at interstitial lattice sites, where they act as magnetic obtained crystals with natural hexagonal facets were microprobes. In a magnetic material the muon spin easily exfoliated. The XRD pattern of a single crystal of precesses in the local field Bint at the muon site with the KV3Sb5 was collected using a Bruker D2 Phaser X-ray Larmor frequency νint = γµ/(2π)Bint (muon gyromag- −1 diffractometer with Cu Kα radiation (λ = 0.15418 nm) netic ratio γµ/(2π) = 135.5 MHz T ). Using the µSR at room temperature. The single-crystal diffraction was technique, important length scales of superconductors implemented on a Bruker D8 Venture system equipped can be measured, namely the magnetic penetration with a Mo Kα (λ = 0.71073 A).˚ The crystal structure depth λ and the coherence length ξ. If a type II su- was solved and refined using the Bruker SHELXTL perconductor is cooled below Tc in an applied magnetic Software Package. Critical magnetic fields for KV3Sb5 field ranging between the lower (Hc1) and the upper for the field applied along the c-axis and within the (Hc2) critical fields, a vortex lattice is formed which kagome plane are 100 mT and 300 mT, respectively. in general is incommensurate with the crystal lattice, with vortex cores separated by much larger distances Experimental details: Zero field (ZF) and trans- than those of the crystallographic unit cell. Because the verse field (TF) µSR experiments on the single crystalline implanted muons stop at given crystallographic sites, samples of KV3Sb5 were performed on the GPS, Dolly, they will randomly probe the field distribution of the 7 vortex lattice. Such measurements need to be performed II. ACKNOWLEDGMENTS in a field applied perpendicular to the initial muon spin polarization (so-called TF configuration). The µSR experiments were carried out at the Analysis of TF-µSR data: The TF µSR data were Swiss Muon Source (SµS) Paul Scherrer Insitute, Vil- analyzed by using the following functional form44: ligen, Switzerland. The magnetization measurements were carried out on the MPMS device of the Labora- h (σ2 + σ2 )t2 i tory for Multiscale Materials Experiments, Paul Scher- P (t) = A exp − SC nm cos(γ B t + ϕ) TF s 2 µ int,s rer Institute, Villigen, Switzerland (SNSF grant no. (3) 206021 139082). Z.G. acknowledges the useful discus- Here As denotes the initial asymmetry of the sample. sions with Robert Scheuermann, Elvezio Morenzoni and ϕ is the initial phase of the muon-spin ensemble and Alex Amato. M.Z.H. acknowledges visiting scientist sup- Bint represents the internal magnetic field at the muon port from IQIM at the California Institute of Technology. site. The relaxation rates σSC and σnm characterize the This work is also supported by Beijing Natural Science damping due to the formation of the flux-line lattice in Foundation (Grants No. Z180008), the National Key Re- the superconducting state and of the nuclear magnetic search and Development Program of China (Grants No. dipolar contribution, respectively. As indicated by the 2017YFA0302900), the National Natural Science Foun- solid lines in Figs. 2b,c the µSR data are well described dation of China (Grants No. U2032204). The work of by Eq. (1). R.G. was supported by the Swiss National Science Foun- dation (SNF Grant No. 200021 175935).

∗ These authors contributed equally to the paper. 15 Ortiz, B., Sarte, P., Kenney, E., Graf, M., Teicher, S., † Electronic address: [email protected] Seshadri, R., & Wilson, S. Superconductivity in the Z2 ‡ Electronic address: [email protected] kagome metal KV3Sb5. arXiv:2012.09097v1 (2020). 1 Syˆozi,I. Statistics of Kagome Lattice. Prog. Theor. Phys. 16 Yin, Q., Tu, Z., Gong, C., Fu, Y., Yan, S., & Lei, H. 6, 306 (1951). Superconductivity and normal-state properties of kagome 2 Barz, H. Ternary transition metal phosphides: High- metal RbV3Sb5 single crystals. arXiv:2101.10193 (2021). temperature superconductors. Mater. Res. Bull. PNAS 15, 17 Jiang, Y.-X. et al. Discovery of topological charge order 1489 (1980). in kagome superconductor KV3Sb5. arXiv:2012.15709 (Na- 3 Zhou, Y., Kanoda, K., & Ng, T.-K. Quantum spin liquid ture Materials, on-line). states. Rev. Mod. Phys. 89, 025003 (2017). 18 Shumiya, N. et al. Tunable chiral charge order in kagome 4 Guguchia, Z. et al. Tunable anomalous Hall conductivity superconductor RbV3Sb5. arXiv:2105.00550 (2021). through volume-wise magnetic competition in a topological 19 Wang, Z. et al. Anomalous transport and chiral charge or- kagome magnet. Nature Communications 11, 559 (2020). der in kagome superconductor CsV3Sb5. arXiv:2105.04542 5 Yin, J.-X., Pan, S.H., & Hasan, M.Z. Probing topological (2021). quantum matter with scanning tunneling microscopy. Na- 20 Feng, X., Jiang, K., Wang, Z., & Hu, J. Chiral flux phase ture Review Physics 3, 249-263 (2021). in the Kagome superconductor AV3Sb5. arXiv:2103.07097 6 Guguchia, Z. et al. Multiple quantum phase transitions (2021). of different nature in the topological kagome magnet 21 Denner, M., Thomale, R., & Neupert, T. Analysis of Co3Sn2−xInxS2. npj Quantum Materials 6, 50 (2021). charge order in the kagome metal AV3Sb5 (A = K, Rb, 7 Mielke III, C. et al. Nodeless kagome superconductivity in Cs). arXiv:2103.14045v1 (2021). 22 LaRu3Si2. Phys. Rev. Mat. 5, 034803 (2021). Lin, Y.-P., & Nandkishore, R. Complex charge den- 8 Yin, J-X. et al. Giant and anisotropic spin-orbit tunability sity waves at Van Hove singularity on hexagonal lattices: in a strongly correlated kagome magnet. Nature 562, 91-95 Haldane-model phase diagram and potential realization in (2018). kagome metals AV3Sb5. arXiv:2104.02725 (2021). 9 Yin, J.-X. et al. Negative flatband magnetism in a spin-orbit 23 Wu, X., et al. Nature of unconventional pairing in coupled kagome magnet. Nature Physics 15, 443 (2019). the kagome superconductors AV3Sb5. arXiv:2104.05671v1 10 Nakatsuji, S., Kiyohara, N. & Higo, T. Large anomalous (2021). Hall effect in a non-collinear antiferromagnet at room tem- 24 Chandan Setty, C., Hu, H., Chen, L., Si, Q. Electron cor- perature. Nature 527, 212-215 (2015). relations and T -breaking density wave order in a Z2 kagome 11 Ye, L. et al. Massive Dirac fermions in a ferromagnetic metal, https://arxiv.org/abs/2105.15204 (2021). kagome metal. Nature 555, 638-642 (2018). 25 Haldane, F.D.M. Model for a quantum Hall effect without 12 Ghimire, N. & Mazin, I. Topology and correlations on the Landau levels: condensed-matter realization of the parity kagome lattice. Nature Materials 19, 137-138(2020). anomaly. Phys. Rev. Lett. 61, 2015-2018 (1988). 13 J-X. Yin, et. al., Quantum-limit Chern topological mag- 26 Varma, C.M. Non-Fermi-liquid states and pairing instabil- netism in TbMn6Sn6. Nature 583, 533-536 (2020). ity of a general model of copper oxide metals. Phys. Rev. B 14 Ortiz, B. et al. CsV3Sb5:AZ2 Topological Kagome Metal 55, 14554-14580 (1997). with a Superconducting Ground State. Phys. Rev. Lett. 27 Chakravarty, S., Laughlin, R., Morr, D., & Nayak, C. Hid- 125, 247002 (2020). den order in the cuprates. Phys. Rev. B 63, 094503 (2001). 8

28 Yang, S., & Wang, Y. et al. Giant, unconventional anoma- Supplementary information lous Hall effect in the metallic frustrated magnet candidate, KV3Sb5. Sci. Adv. 6, 1-7s (2020). 29 Kiesel, M.L. & Thomale, R. Sublattice interference in the kagome Hubbard model. Phys. Rev. B 86, 121105 (2012). 30 Sonier, J. E., Brewer, J.H., & Kiefl, R.F. µSR studies of the vortex state in type-II superconductors. Rev. Mod. Phys. 72, 769 (2000). 31 Guguchia, Z. et al. Signatures of the topological s+− super- conducting order parameter in the type-II Weyl semimetal Td-MoTe2. Nature Communications 8, 1082 (2017). 32 Luke, G.M. et al. Time-reversal symmetry breaking super- A. Crystal structure of KV3Sb5 conductivity in Sr2RuO4. Nature 394, 559 (1998). 33 Hillier, A.D., Jorge, Q. & Cywinski, R. Evidence for Time- Reversal Symmetry Breaking in the Noncentrosymmet- ric Superconductor LaNiC2. Phys. Rev. Lett. 102, 117007 Additional characterization information is provided (2009). here on the kagome superconductor KV3Sb5 which crys- 34 Biswas, P.K. et al. Evidence for superconductivity with tallizes in the novel AV3Sb5-type structure (space group broken time-reversal symmetry in locally noncentrosym- P 6/mmm, where A = K, Rb, Cs). The crystallographic metric SrPtAs. Phys. Rev. B 87, 180503(R) (2013). structure of prototype compound KV Sb shown in panel 35 Grinenko, V. et al. Split superconducting and time-reversal 3 5 (a) of Figure 4 illustrates how the V atoms form a kagome symmetry-breaking transitions in Sr2RuO4 under stress, Nature Physics (2021), DOI: 10.1038/s41567-021-01182-7. lattice (medium beige circles) intertwined with a hexag- 36 Grinenko, V. et al. Unsplit superconducting and time re- onal lattice of Sb atoms (small red circles). The K atoms versal symmetry breaking transitions in Sr2RuO4 under (large purple circles) occupy the interstitial sites between hydrostatic pressure and disorder, Nature Communications the two parallel kagome planes. In panel (b) the vana- 12, 3920 (2021) dium kagome net has been emphasized, with the inter- 37 Kenney, E., Ortiz, B., Wang, C., Wilson, S., & Graf, M. penetrating antimony lattice included to highlight the Absence of local moments in the kagome metal KV3Sb5 as unit cell (see dashed lines). Figures 4c and d show an op- determined by muon spin spectroscopy. J. Phys.: Condens. tical microscope image of a 3 x 2 x 0.2 mm single crystal Matter 33, 235801 (2021). of KV Sb on millimeter paper and Scanning Transmis- 38 Kubo, R. & Toyabe, T. Magnetic Resonance and Relax- 3 5 ation North Holland, Amsterdam, (1967). sion Microscope (STM) image of the V kagome lattice 39 Sedlak, K., Scheuermann, R., Stoykov, A., Amato, A. from a cryogenically cleaved sample, respectively. GEANT4 simulation and optimisation of the high-field mSR spectrometer. Physica B 404 970-973 (2009). 40 Gupta, R. et al. Self-consistent two-gap approach in study- ing multi-band superconductivity of NdFeAsO0.65F0.35. Frontiers in Physics 8, 2 (2020). 41 Brandt, E.H. Flux distribution and penetration depth measured by muon spin rotation in high-Tc superconduc- tors. Phys. Rev. B 37, 2349 (1988). 42 Kogan, V. G., Martin, C., & Prozorov, R. Phys. Rev. B 80, 014507 (2009). 43 Khasanov, R. et al. Evolution of Two-Gap Behavior of B. Laboratory X-ray diffraction experiments on a the Superconductor FeSe1−x. Phys. Rev. Lett. 104, 087004 single crystal of KV3Sb5 (2010). 44 Suter, A., & Wojek, B.M. Musrfit: A free platform- independent framework for µSR data analysis. Physics Pro- cedia 30, 69 (2012). A single crystalline sample has been selected and 45 Uemura, Y.J. et al. Phys. Rev. Lett. 62, 2317 (1989). X-ray diffraction was performed on it. The crystal was 46 von Rohr, F. O. et al. Unconventional Scaling of the Super- oriented such that the incident X-rays scattered off of the fluid Density with the Critical Temperature in Transition ab-plane. The resultant diffraction pattern shows clear Metal Dichalcogenides. Science Advances 5(11), eaav8465 diffraction peaks, which have been indexed (see Figure (2019). 5a) and fitted using the SHELX-2018/3 program. The 47 Luetkens, H. et al. Nature Materials 8, 305 (2009). 48 obtained crystallographic information is summarized in Shengelaya, A. et al. Phys. Rev. Lett. 94, 127001 (2005). Table 2. The Laue X-ray diffraction image (see Figure 49 Amato, A. et al. The new versatile general purpose surface- muon instrument (GPS) based on silicon photomultipliers 5b) demonstrates the single crystallinity of the samples for µSR measurements on a continuous-wave beam. Review used for µSR experiments. The magnetic susceptibility of Scientific Instruments 88, 093301 (2017). measurements show the abrupt drop in susceptibility at T ∗ ' 80 K (see Figure 5c) which comes from the charge ordering below this temperature. 9

C. Crystallographic Structural Information D. Zero-field µSR spectrum for KV3Sb5

T able 1: summary of the atomic positions obtained In Figure 6 is displayed the zero-field µSR spectrum from fitting to the diffraction data. recorded at 5 K. The red solid curve represents the fit to the recorded time spectra, using only the Gaussian Kubo 38 Atomic Positions Toyabe (GKT) function : 2 Atom Wyckhoff Positions x y z Ueq/A˚ Occ Sb 1a 1 1 0 0.16 1 1 2 h ∆2t2 i 01 P GKT (t) = + (1 − ∆2t2) exp − (4) Sb02 4h 2/3 1/3 0.25381 0.019 1 ZF 3 3 2 V01 3f 1/2 1/2 0 0.015 1 K01 1b 1 0 1/2 0.045 0.977 where ∆/γµ is the width of the local field distribution due to the nuclear moments and γµ/2π = 135.5 MHz/T is the T able 2: summary of the crystallographic, refinement, muon gyromagnetic ratio. The Gaussian Kubo-Toyabe and data collection information from the refinements to depolarization function (Eq. 4) reflects the field distri- the single crystal X-ray diffraction data. bution at the muon site created by the nuclear moments of the sample. A GKT shape is expected in KV3Sb5 due Crystallographic Data to the presence of the dense system of nuclear moments 39 Chemical Formula K0.98V3Sb5 with large values of nuclear spins (I = 3/2 for K, I = g 7/2 for 51V, and I = 5/2 for 121Sb) in KV Sb and a Formula Weight 799.81 mol 3 5 Crystal System hexagonal high natural abundance. It is clear that the GKT func- Space Group P 6/mmm tion alone is not sufficient to fully describe the zero-field a 5.4831(3) A˚ µSR spectrum. Multiplying the GKT function by the c 8.9544(9) A˚ additional exponential exp(-Γt) term (Eq. 5), which is α 90◦ of electronic in origin, is essential to fully describe the γ 120◦ spectrum as shown in Figure 5 and in the inset. Thus, Z 1 A˚ ZF-µSR consists of the muon spin relaxations ∆ and Γ Cell Volume 233.14(3) A˚3 due to nuclear and electronic moments, respectively. g Density (calculated) 5.697 cm3 Temperature 273(2) K 1 2 h ∆2t2 i P GKT (t) = + (1 − ∆2t2) exp − exp(−Γt) Wavelength Mo Kα (λ = 0.71073 A)˚ ZF ◦ ◦ 3 3 2 θ Range 2.27 to 36.15 (5) Data Collection Diffractometer Bruker D8 Venture Index Ranges -9 ≤ h ≤ 9 E. High-field µSR spectrum for the single crystal -9 ≤ k ≤ 9 of KV3Sb5 obtained on the instrument HAL-9500 -14 ≤ l ≤ 14 Reflections Collected 7207 During the high-field experiments, a mosaic of several Independent Reflections 278 [Rint = 0.0380] crystals stacked on top of each other was used. The in- Absorption Coefficient 17.469 mm−1 dividual crystals were attached to 10 mm circular silver Structure Solution Tech. Direct methods sample holder and the entire ensemble was held together Structure Solution Prog. SHELX-2014/5 (Sheldrick,2014) by small droplets of GE varnish. Figure 7 shows the Refinement probability field distribution, measured at 5 K in the c- 2 axis magnetic field of 8 T. In the whole investigated tem- Refinement Method Full-matrix least-squares on F perature range, two-component signals were observed: a Refinement Program SHELX-2018/3 (Sheldrick,2018) −1 2 2 2 signal with fast relaxation 0.42 µs (broad signal) and Function Minimized Σ w(F o +F c ) another one with a slow relaxation 0.08 µs−1 (narrow Data 278 reflections (271 final) one). The narrow signal arises mostly from the muons Parameters 13 parameters stopping in the silver sample holder and its position is a Restraints 0 restraints 2 precise measure of the value of the applied magnetic field. Goodness-of-fit on F 1.233 The width and the position of the narrow signal is found Final R Indices I ≤ σ(I) R1=0.0162, wR2=0.039 to be temperature independent (see the inset of Fig. 7) all data; R1=0.0171, wR2=0.039 1 as expected and thus they were kept constant in the anal- Weighting Scheme w = 2 2 2 σ (F0 )+(0.0130P ) +0.4869P ysis. The relative fraction of the muons stopping in the 2 2 where P = (F o +2F c )/3 sample was fixed to the value obtained at the base-T Extinction Coefficient 0.0050(10) and kept temperature independent. The signal with the −3 −3 ∆ρMax and ∆ρsMin 1.259 eA˚ and -0.481 eA˚ fast relaxation, which is shifted towards the lower field R.M.S. Dev. from Mean 0.164 eA˚3 from the applied one, arises from the muons stopping in 10 the sample and it takes a major fraction (60-70 %) of lar lattice, where ϕ is the polar angle around the Fermi the µSR signal. This points to the fact that the sample surface. We note that the d-wave form ∆0,i = ∆i cos 2ϕ response arises from the bulk of the sample. Based on makes no real sense on a lattice with 3-fold rotational this two component signal we can determine Knight shift symmetry. The analysis certainly rules out a simple 1- −2 which is defined as Kexp = (Bint − Bext)/Bext, where Bint gap s-wave model as an adequate description of λ (T ) and Bext are the internal and externally applied magnetic for KV3Sb5. The 2-gap s+s-wave or 2-gap s+f-wave fields, respectively. models line up with the experimental data very well. Due to the lack of data points below 0.25 K, it is difficult to distinguish between s+s-wave and s+f-wave models and F. Superconducting gap symmetry to give an estimate for the zero-temperature value of λ. We also observe the diamagnetic shift of about 1 G in From theory side, several scenarios for electronically the SC phase, as shown in Fig. 8c, which indicates a T c mediated, unconventional superconductivity have been of ' 1.1 K, and furthermore supports the bulk nature of 23 discussed theoretically . In the bands of the AV3Sb5 superconductivity in KV3Sb5. materials, van Hove singularities are found close to the The value of the in-plane penetration depth λab at Fermi energy - an electronic structural motif shared with 0.3 K, determined from σSC,ab (superconducting screen- other systems such as the cuprate superconductors or ing currents flowing parallel to the kagome plane), is Sr2RuO4. As a particular feature of the kagome lattice, found to be λab ' 877(20) nm. The value of the out- 29 however, there is a sublattice interference mechanism of-plane penetration depth, determined from σSC,c (su- by which the Bloch states near each van Hove point are perconducting screening currents flowing perpendicular supported on a distinct sublattice. This promotes the to the kagome plane), is found to be λc ' 730(20) nm. relevance of longer-range interactions and unconventional We note that these are the upper limits of λab and λc, pairing states. since they are determined at 0.3 K and obtained by as- Figure 8a shows the temperature dependences of the suming a 100 % superconducting volume fraction. On superconducting relaxation rates σSC,ab and σSC,ac, de- the other hand vortex lattice disorder increases the mea- termined from the data with the field applied along the sured relaxation rate leading to an underestimation of c-axis and within the kagome plane, respectively. When a the penetration depth. magnetic field is applied along the c-axis, superconduct- ing screening currents will flow in the ab plane; whereas in the case of a magnetic field applied in-plane, the super- G. Analysis of the temperature dependence of the conducting screening currents will flow along the c and penetration depth a axes. Thus, σSC,ac consists of both the in-plane and c-axis contributions. A two-step nature of the supercon- ducting state with the onset of T c ' 1.1 K is clearly The temperature dependence of λ is particularly sen- sitive to the topology of the SC gap: while in a fully visible in the 5 mT data for both orientation of the mag- −2 −2 −2 netic field with respect to the c-axis. gapped superconductor, ∆λ (T ) ≡ λ (0) − λ (T ) The temperature dependence of λ is sensitive to the vanishes exponentially at low T , in a nodal supercon- topology of the SC gap: while in a fully gapped su- ductor it vanishes as a power of T . λ(T ) was calculated perconductor, ∆λ−2 (T ) ≡ λ−2 (0) − λ−2 (T ) vanishes within the local (London) approximation (λ  ξ) by the exponentially at low T , in a nodal SC it vanishes as following expression: a power of T . To quantitatively analyze the tempera- −2 2π ∞ 31,44 λ (T, ∆ ) 1 Z Z ∂f EdEdϕ ture dependence of the penetration depth λ(T ) for 0,i = 1+ ( ) , −2 p 2 2 KV3Sb5 we employ the empirical α-model. The λ(T ) λ (0, ∆0,i) π 0 ∆(T,ϕ) ∂E E − ∆i(T, ϕ) in the applied field of 5 mT shows a well pronounced (6) −1 two step behavior, as shown in Fig. 8a. This suggests where f = [1 + exp(E/kBT )] is the Fermi func- that at least two bands are involved in superconductiv- tion, ϕ is the angle along the Fermi surface, and ity and that the interband coupling is nearly zero which ∆i(T, ϕ) = ∆0,iΓ(T/Tc)g(ϕ) (∆0,i is the maximum would result in different values of Tc for different bands. gap value at T = 0). The temperature depen- With the application of 10 mT, the two-step transition dence of the gap is approximated by the expression 0.51 becomes somewhat smoothed out and less pronounced. Γ(T/Tc) = tanh {1.82[1.018(Tc/T − 1)] }, while g(ϕ) The results of the analysis for the µ0H = 10 mT using describes the angular dependence of the gap and it is α-model are presented in Fig. 8b. We consider two dif- replaced by 1 for both an s-wave and an s+s-wave gap, ferent possibilities for the gap function: either a constant | cos(2ϕ)| for a d-wave gap, and | cos(3ϕ)| for a f-wave gap, ∆0,i = ∆i, or an angle-dependent f-wave gap of the gap. form ∆0,i = ∆i cos 3ϕ which is appropriate for a triangu- 11

(a) (b)

KV3Sb5

1x1 2x2 (c) (d)

1nm

Figure 4: Crystal structure of KV3Sb5. Three dimensional representation (a) and top view (b) of the atomic structure of KV3Sb5. In panel (c) is displayed an optical microscope image of a 3 x 2 x 0.2 mm single crystal of KV3Sb5 on millimeter paper, with the scale shown. The hexagonal symmetry is immediately apparent. (d) Scanning Transmission Microscope (STM) image of the V kagome lattice from a cryogenically cleaved sample.

∗ These authors contributed equally to the paper. ‡ Electronic address: [email protected] † Electronic address: [email protected] 12

T* (a) (b) 10

) (c)

4 KV Sb 3 5 0

0 ( µ0H = 3 T KV3Sb5 9 ) ) 1 . - g u

. 1 - a

e 8 (

O

y u t i ) m 1 s

e ) 0 7

3 n -

0 7 0 0 ( e

1 t 0 ) ( ( )

2 n

8 c I

0

0 )

0 ) ) ( 5 0

7 6 (

0 6

0 0

0 0 ( 0 ( (

5 0 10 20 30 40 50 60 70 80 90 100 0 50 100 150 200 250 300 2q (degree) T (K)

Figure 5: Single Crystal X-Ray Diffraction for KV3Sb5. (a) X-ray diffraction image for KV3Sb5 recorded at 300 K. The well-defined peaks are labeled with their crystallographic indices. No second phase has been detected. (b) Laue X-ray diffraction image of the single crystal sample KV3Sb5, oriented with the c-axis along the beam. (c) The temperature dependence ∗ of magnetic susceptibility of KV3Sb5 above 1.8 K. It shows an anomaly at T ' 80 K, coinciding with emergence of a charge order. 13

0 . 2 2 0 . 2 0 0 . 2 0 ) t (

4 0 . 1 8 3 0 . 1 5 P 0 . 1 6 ) t (

4 0 . 1 4

3 0 1 2 3

P 0 . 1 0 t ( µs )

K V 3 S b 5 0 . 0 5 5 K G K T G K T ∗ e x p ( - Γt ) 0 . 0 0 0 2 4 6 8 t ( µs )

Figure 6: Zero-field µSR experiment for KV3Sb5. The ZF µSR time spectrum for KV3Sb5, obtained at T = 5 K. The solid curves represent fits to the recorded time spectra, using only Gaussian Kubo Toyabe (GKT) function (red) and the one with an additional exponential exp(-Γt) term (blue). The inset shows the low time part of the spectrum. 14 2 5 ) 1

- S a m p l e K V 3 S b 5 0 . 4 s

µ S i l v e r (

D a t a i n 8 T ) . e

t 0 . 3

u T w o c o m p o n e n t f i t a r . 2 0 n a S a m p l e s i g n a l o ( 0 . 2 i

S i l v e r s i g n a l t a n x a o 0 . 1 l i t e R u

b 0 . 0 i 1 5

r 0 5 0 1 0 0 1 5 0 2 0 0 2 5 0 3 0 0 t s

i T ( K ) d

d l e i

f 1 0

y t i l i b a b

o 5 r P B i n t B e x t 0 7 . 9 8 7 . 9 9 8 . 0 0 8 . 0 1 8 . 0 2 8 . 0 3 F i e l d ( T )

Figure 7: High-field µSR experiment for KV3Sb5. Fourier transform for the µSR asymmetry spectra of KV3Sb5 at 5 K for the applied field of µ0H = 8 T. The black solid line represents the fit to the data using the two component signal. Red and blue solid lines show the signals arising from the sample and the silver sample holder (mostly), respectively. The inset shows the temperature dependences of the muon spin relaxation rates arising from the sample and the silver sample holder. 15

( a ) T 0 . 2 0 c 2 . 0 ) ) 1 H | | c - a x i s - 2 - s H | | k a g o m e p l a n e µ m (

0 . 1 5 1 . 5 µ ( c

a , µ 2 -

C H = 5 m T

0 c S a σ λ

, 0 . 1 0 1 . 0 , ) ) 1 2 - - s m µ ( µ

(

b 0 . 0 5 0 . 5 2 a , - b C a S

K V S b λ

σ 3 5 0 . 0 0 0 . 0 0 . 0 0 . 5 1 . 0 1 . 5 2 . 0 ( b ) T ( K ) H | | k a g o m e p l a n e s + f o r s + s s - w a v e ) 0 . 1 0 1 -

s µ 0 H = 1 0 m T µ (

c a , C

S 0 . 0 5 σ

K V 3 S b 5 0 . 0 0 0 . 0 0 . 5 1 . 0 1 . 5 2 . 0 T ( K ) ( c ) 0 . 0 0

- 0 . 0 2 ) T m

( - 0 . 0 4

t n i

B - 0 . 0 6 ∆

- 0 . 0 8 K V 3 S b 5 - 0 . 1 0 0 . 0 0 . 5 1 . 0 1 . 5 2 . 0 T ( K )

Figure 8: Superconducting gap symmetry in KV3Sb5. (a) The SC muon depolarization rates σSC,ab, and σSC,ac as −2 −2 well as the inverse squared magnetic penetration depth λab and λac as a function of temperature, measured in 5 mT, applied parallel and perpendicular to the kagome plane. (b) The SC muon depolarization rate σSC,ac, measured in 10 mT, applied parallel to the kagome plane. The solid line represents the indistinguishable 2-gap s-wave and s + d wave model. The error bars represent the s.d. of the fit parameters. (c) (c) Temperature dependence of the difference between the internal field µ0HSC measured in the SC state and the one measured in the normal state µ0HNS at T = 5 K for KV3Sb5.