Signatures of the Kondo effect in VSe2

VSe2 is a transition metal dichaclogenide which has a charge- density wave transition that has been well studied. We report on a low-temperature upturn in the resistivity and, at temperatures below this resistivity minimum, an unusual magnetoresistance which is negative at low fields and positive at higher fields, in single crystals of VSe2. The negative magnetoresistance has a parabolic dependence on the magnetic field and shows little angular dependence. The magnetoresistance at temperatures above the resistivity minimum is always positive. We interpret these results as signatures of the Kondo effect in VSe2. An upturn in the susceptibility indicates the presence of interlayer V ions which can provide the localized magnetic moments required for scattering the conduction electrons in the Kondo effect. The low-temperature behaviour of the heat capacity, including a high value of γ, along with a deviation from a Curie-Weiss law observed in the low-temperature magnetic susceptibility, are consistent with the presence of magnetic interactions between the paramagnetic interlayer V ions and a Kondo screening of these V moments.

transition at the lower temperature 16,17 . The magnetic susceptibility shows a Curie-Weiss like behaviour at low temperatures which is supposed to originate because of an excess of V ions sitting in the van der Waal spaces between the VSe 2 layers in the metal rich VSe 2 4, 14-16, 18 . It is suggested that each interstitial V ion produces a net paramagnetic moment of 2.5 Bohr magnetons and the atomic percentage of the interstitial V ions can be determined using this value and the value of the Curie constant 14,16,18 . According to DiSalvo and Waszczak the CDW transition temperature is suppressed at a rate of 13 ± 3 K/at % of V interstitial ions 16 . The CDW transition is also reported to be suppressed with decreasing thickness in nanoflakes of VSe 2 19 . The resistivity of VSe 2 at low temperatures has variously been reported to show both T 2 and T 4 dependence 4,14,17 . The Hall constant is n-type and shows a large increase with temperature below the CDW transition temperature 4,14,15,17 . The transverse magnetoresistance in VSe 2 , in fields either perpendicular or parallel to the layers of the sample shows a monotonic rise with field and does not saturate 17 .
At low temperatures, quantum effects sometimes become visible in the macroscopic properties such as the resistivity of a material. One such example is the low-temperature upturn in resistivity which is seen in a variety of systems and can originate from quantum corrections to resistivity due to weak localization or electron-electron interactions [20][21][22][23][24][25] . A similar upturn in resistivity at low temperatures is also seen in the Kondo effect which is due to the scattering of conduction electrons by localized magnetic impurities through the s-d exchange interaction 25,26 . In this paper, we report a low-temperature upturn in the resistivity of single crystal samples of VSe 2 and a negative magnetoresistance that is seen at low fields and temperatures below the minimum in resistivity. The VSe 2 single crystal samples were characterized using single crystal X-ray diffraction, Laue X-ray diffraction, resistivity, Hall effect, DC magnetization measurements and specific heat measurements. The low-temperature rise in resistivity and the negative magnetoresistance seen at low fields and temperature are analysed in conjunction with the magnetization and heat capacity data to determine if they are due to the Kondo effect arising from a scattering of the conduction electrons by the paramagnetic V ions present in the van der Waal spaces in between the layers of VSe 2 .

Results
Polycrystalline powder samples of VSe 2 as well as single crystals of VSe 2 were synthesized. First, the results of the polycrystalline VSe 2 are presented, followed by those for the single crystals of VSe 2 .
Polycrystalline VSe 2 . Room temperature powder X-ray diffraction measurements were performed to determine the phase purity of the polycrystalline material. The X-ray data were fitted to the P m 3 1 (No. 164) space group using the FullProf software suite 27 . The X-ray diffraction profile for the VSe 2 polycrystalline material displayed in Fig. 1(a) shows no trace of any impurity and reveals that the hexagonal structure of the 1T polytype is formed 4 . The relative intensities of some of the Bragg peaks (hkl = (001), (002), (003) and (004)) are significantly increased for polycrystalline samples, due to the preferred orientation present in the sample. The lattice parameters (a = 3.35844(3) Å and c = 6.10141(5) Å) were found to be similar to the published values for VSe 2 1, 4 . Figure 2(a) shows the susceptibility and the inverse susceptibility versus temperature for the polycrystalline VSe 2 sample. The data reveal a prominent upturn at low temperatures and a change in slope around 112 K corresponding to the CDW transition. We have fitted both the susceptibility data and its inverse from 1.8-40 K to a modified Curie-Weiss law 0 which consists of a temperature independent term in addition to the Curie-Weiss part. The results of the fit are shown in the inset of Fig. 2(a). From the fit to the inverse susceptibility data, we obtain the values of C = (0.00797 ± 0.00003) emu-K/mol, Θ = (−1.04 ± 0.02) K and χ 0 = (1.51 ± 0.01) × 10 −4 emu/mol. We discuss this in more detail later, in the section on the susceptibility versus temperature data for the single crystal sample. Figure 2(b) shows the low-temperature susceptibility for different applied magnetic fields for the polycrystalline sample. The magnetic susceptibility decreases with increasing applied field and there is no clear signature of any magnetic order down to the lowest temperatures measured.
Single crystals of VSe 2 . Single crystals of VSe 2 were studied using electrical transport, Hall, magnetoresistance, DC magnetization, and heat capacity measurements. In Fig. 1(b) we show the X-ray Laue back scattered diffraction pattern of a single crystal of VSe 2 obtained along the [001] direction. The pattern shows the good quality of the single crystal. In order to confirm the crystal structure of our VSe 2 samples, single crystal X-ray diffraction measurements were performed on the single crystals. A refinement of the single crystal X-ray diffraction data was performed using the space group P m 3 1 (No. 164) and a goodness of fit (GooF) of 1.068 was obtained. The values of the crystal parameters obtained were a = b = 3.357(1) Å, c = 6.104(5) Å and a cell volume of 59.59(6) Å 3 which are in agreement with previously reported results 1 . We also performed compositional analysis on the single crystals using energy dispersive X-ray spectroscopy (EDS) which yielded an average composition of V = (28 ± 2)% and Se = (72 ± 2)%, showing a V deficiency in the overall stoichiometry. Below, we discuss the results of our measurements on single crystals of VSe 2 .
Low-temperature upturn in resistivity. Figure 3(a) shows the variation of resistivity with temperature, ρ(T), down to 0.45 K for a single crystal of VSe 2 . It is evident that the single crystal has a metallic nature with a plateau around 110 K which marks the CDW transition in VSe 2 . However, instead of flattening off to a residual value at  low temperatures, as is expected for a simple metal, the resistivity exhibits a minimum at around 6.5 K and then increases as the temperature decreases until about 0.5 K where ρ(T) shows signs of saturation (see the inset of Fig. 3(a)). We have measured the resistivity versus temperature of crystals taken from two separate growths and they all show an upturn in resistivity at low temperature. We show some of these results in the supplementary information. Such a low-temperature minimum in ρ(T) can occur due to electron-electron interactions, weak localization, or the Kondo effect 21,22,28,29 .
For three dimensional systems, an increase in resistivity with decreasing temperature should follow either a T −p/2 dependence if it is due to weak localization ( = p 3 2 , 2, 3 depending upon the scattering mechanism) 24, 30 , a − T dependence if it is due to electron-electron interactions 24 , or a −ln T dependence if it is due to the Kondo effect 23,25 . Weak localization and electron-electron interactions can be ruled out as possible reasons for the observed upturn in resistivity as fits to a −T 3/4 (weak localization in the dirty limit of dominant electron-electron scattering) and a −T 1/2 behaviour produced poor agreement with the experimental data. Figure 3(b) shows the low-temperature upturn in resistivity versus ln T, from which it is clear that for a certain range of temperature the upturn in resistivity is linear in ln T, as expected for the Kondo effect, and only begins to flatten off at the lowest temperatures measured. The deviation from a −ln T dependence may be the result of the formation of a Nagaoka spin-compensated state or Ruderman-Kittel-Kasuya-Yosida (RKKY) mediated interactions between the magnetic impurities responsible for the Kondo effect 31,32 .
We carefully examined the fitting of the temperature dependence of resistivity within the framework of the Kondo effect. The Kondo resistivity in the region where ρ(T) is linear in ln T can be described by the Hamann expression 32, 33 , where, ρ 0 is the unitarity limit, T K is the Kondo temperature and S is the spin of the magnetic impurity. We account for any remaining temperature dependence in ρ(T) due to electron-phonon scattering at temperatures below the resistivity minimum, by first fitting the ρ(T) data in the temperature range 10 to 20 K to an expression of the form a + bT n . We obtain the numerical values for b and n equal to (6.8 ± 0.7) × 10 −11 and 3.37 ± 0.03, respectively. The low-temperature resistivity is then fit over the range 1.5 to 10 K using an expression that includes the Hamann term, along with an electron-phonon contribution bT n and a temperature independent term ρ b , While performing this fit, the coefficient b and power n of the electron-phonon contribution are held fixed at the values obtained from the higher temperature ρ(T) data. The result is shown in Fig. 3(b) as Fit 1 and the parameters obtained are listed in Table 1.
The Hamann expression is known to give a poor fit for temperatures T < T K , where the resistivity deviates from a ln T behaviour and can give low values for S as a result of inadequacies in the Nagaoka approximation on which the Hamann expression is based 32,34 . The deviation of the resistivity from the ln T behaviour at low temperature can be accounted for by considering the RKKY interactions between the paramagnetic interlayer V ions. The low-temperature resistivity can be fit to a modified Equation (3) where the temperature T in the expression for ρ sd is replaced by an effective temperature 31 . For this second fit, S and T K along with b and n were fixed to the values obtained in Fit 1. The result of Fit 2 is also shown in Fig. 3(b) and provides a satisfactory description of the data across the entire temperature range 0.45 to 10 K. The parameters obtained from Fit 2 are listed in Table 2. The value of 1.18 ± 0.02 K obtained for the temperature T W is consistent with the temperature at which the resistivity starts deviating from a ln T behaviour.
Negative magnetoresistance. The magnetoresistance of the samples was measured with the magnetic field parallel to the c-axis of the single crystal samples. Figure 4(a) shows the relative magnetoresistance at different temperatures. It can be seen that the magnetoresistance below 7 K is negative at low fields but starts showing a positive trend as the field increases further. At 7 K and above, the magnetoresistance is positive for all fields. Thus the ρ 0 (Ω-cm) ρ b (Ω-cm) T K (K) S (5.9 ± 0.7) × 10 −7 (1.770 ± 0.005) × 10 −5 6.8 ± 0.9 0.48 ± 0.09 Table 1. Parameters obtained from a fit to the resistivity versus temperature using Equation (3) in the temperature range 1.5 to 10 K. Table 2. Parameters obtained from the fit of the resistivity versus temperature between 0.45 and 10 K using Equation (3) with T replaced by effective temperature T eff in the expression for ρ sd in order to take into account the effect of RKKY interaction between the interlayer V ions.
negative magnetoresistance seen at low fields occurs in the same temperature range at which a rise in resistivity is seen with falling temperature. This suggests that the negative magnetoresistance and the low-temperature upturn in resistivity are correlated and could arise due to the Kondo effect. We further analyse the negative magnetoresistance to see if it could be due to the Kondo effect. Negative magnetoresistance arising due to the Kondo effect is reported to show a parabolic dependence with the magnetic field at low fields 26 , 2 where S(T) is the coefficient of the magnetoresistance. In Fig. 4(b), we show the fit to the relative negative magnetoresistance up to a field of 20 kOe at different temperatures using Eq. (4). In the case of InSb, a similar negative magnetoresistance was seen due to the Kondo effect and the coefficient S(T) in that case increased with temperature and varied as (S(T)) −1/2 ∝ (T + T 0 ) where T 0 is a parameter 26 . However, for our samples, (S(T)) −1 shows a parabolic temperature dependence as shown in the inset of Fig. 4(b).
Hall concentration versus temperature. Figure 5(a) shows the variation with temperature of the Hall constant, R H , and the Hall carrier concentration, n, obtained using the relation R H = −1/ne. The Hall constant is negative, indicating that the sample is n-type. At temperatures immediately below the CDW transition seen in the resistivity versus temperature, the Hall constant shows a considerable increase in magnitude and the Hall concentration, a significant decrease. VSe 2 has been reported to have n-type charge carriers and a similar increase in the magnitude of the Hall constant at the CDW transition 14,15,17 . A CDW transition leads to a decrease in the density of states at the Fermi level which in turn reduces the carrier concentration 15  Angular dependence of magnetoresistance. We also measured the angular dependence of the transverse magnetoresistance at different temperatures which is shown as relative change in magnetoresistance at different angles in Figs. 5(b-d) for the temperatures 2, 6, and 20 K respectively. VSe 2 being a layered material cleaves perpendicular to the c-axis and hence in our samples, the c-axis is perpendicular to the plane of the sample. The magnetic field is rotated in a plane containing the c-axis and the direction perpendicular to the current. The angular dependence of the magnetoresistance at 2 and 6 K are shown in Figs. 5(b) and (c). From the figures, it is evident that the negative magnetoresistance seen at low fields does not vary much with the change in the angle the field makes with the c-axis of the sample. However, the positive trend seen in the magnetoresistance at higher fields decreases significantly in magnitude as the field becomes perpendicular to the c-axis. At the higher temperature of 20 K, (see Fig. 5(d)), where the negative magnetoresistance has already disappeared, the positive magnetoresistance decreases in magnitude as the field becomes perpendicular to the c-axis. This kind of angular dependence of the transverse magnetoresistance is another indication that the low-temperature upturn in resistivity is due to the Kondo effect with an isotropic s-d exchange interaction between the conduction electrons and the localized paramagnetic spins because in such a case the negative magnetoresistance does not show much dependence on the direction of the field 26,29 . Toriumi and Tanaka 17 reported that the magnetoresistance in VSe 2 at low fields was smaller for fields parallel to the c-axis than for fields perpendicular to the c-axis and vice versa at higher fields. However, we do not observe this effect in our sample. Magnetization measurements. Figure 6(a) shows the results of the susceptibility versus temperature measurement for several single crystals of VSe 2 , co-aligned with their c-axis pointing perpendicular to the magnetic field direction. Unlike the polycrystalline data, where a change in slope is clearly seen around 112 K at the CDW  Fits to a Curie-Weiss law along with a temperature independent term for both the field-cooled susceptibility and inverse susceptibility data above 20 K are also shown. (b) Low-temperature susceptibility at different applied magnetic fields for single crystals of VSe 2 which were co-aligned with their c-axis perpendicular to the magnetic field. The measurements at 100 Oe, 1 kOe, and 10 kOe in the temperature range 0.5 to 2 K were performed separately and the data were adjusted to remove the background signal from the iQuantum 3 He insert.
transition, there is no such prominent transition in the susceptibility data for the single crystals. However, the susceptibility shows an upturn as temperature decreases similar to that seen in the case of the polycrystalline powder. Such an upturn in the susceptibility in VSe 2 has previously been attributed to the presence of V ions in the van der Waal gaps which possess a net paramagnetic moment 16 . We analysed the susceptibility and the inverse susceptibility data for the single crystal samples by fitting it to the modified Curie-Weiss law given in Eq. (1). Figure 6(a) shows the results of these fits for both the field-cooled susceptibility and inverse susceptibility data for the single crystals. We fit only the data above 20 K, as the data clearly deviate from a Curie-Weiss behaviour at lower temperatures (see below). From the fit to the inverse susceptibility data, we obtain the values of C = (0.02773 ± 0.00007) emu-K/mol, Θ = (4.59 ± 0.06) K and χ 0 = (3.143 ± 0.004) × 10 −4 emu/mol. The temperature independent term accounts for a major part of the susceptibility at high temperatures whereas at lower temperatures the Curie-Weiss term plays an increasingly important role. It has been suggested that the V ions in the interlayer gaps each possess a net paramagnetic moment of 2.5 Bohr magnetons. Using this, the molar fraction of the V ions in the interlayer gaps can be calculated. Using the value of the Curie constant obtained from the fit to the susceptibility versus temperature, we arrive at the molar fraction of interlayer V ions of 0.0102 ± 0.0004 and 0.0357 ± 0.0001, respectively, for the polycrystalline and the single crystal samples. In Fig. 6(b), we show the low-temperature susceptibility at different fields for single crystal VSe 2 , with the data at 100 Oe, 1 kOe, and 10 kOe measured down up to 0.5 K. In the lower fields the susceptibility contains a clear anomaly at around 20 K that is not evident in the higher field curves. This feature at 20 K is not present in the χ(T) data for the polycrystalline samples that contain a lower concentration of interlayer V ions as is evident from Fig. 2(b). The measurements at 100 Oe, 1 kOe, and 10 kOe, and in particular the one at 10 kOe, show that the susceptibility deviates significantly from a Curie-Weiss law at low temperatures. The 100 Oe measurement shows a maximum at 0.5 K. The step at 20 K, the deviations from Curie-Weiss behaviour, as well as the low-temperature maximum in χ(T) are the result of the onset of magnetic correlations, short-range magnetic order, or clustering between the V ions, mediated by an RKKY interaction between the interlayer V ions, combined with a reduction in the effective moment of the V ions due to Kondo screening by the conduction electrons. These magnetic correlations and the reduction in the Kondo screened moments are also reflected in the deviation from ln T behaviour seen in the resistivity.
Specific heat measurements. We measured the specific heat of our samples to shed more light on their magnetic and electronic properties. The specific heat of several crystals of VSe 2 , which were co-aligned so that the magnetic field was perpendicular to the c-axis of the crystals, was measured from 0.5 to 10 K in different applied magnetic fields. In Fig. 7(a), we show the results of these measurements, plotting the specific heat divided by temperature versus temperature at different magnetic fields. At zero magnetic field, the specific heat divided by temperature shows a rise at low temperatures. Such a low-temperature rise in the specific heat can be due to a variety of phenomena such as the Schottky anomaly due to transitions between levels in a system of paramagnetic moments which are either intrinsic or present as impurities 36,37 , Kondo behaviour 38,39 , magnetic clusters 40 , and spin fluctuations 41 . Generally in Kondo systems, the rise in specific heat at low temperatures given by the γ value obtained by extrapolating C/T vs. T 2 to T = 0 depends on the number of Kondo scattering centres and is inversely proportional to the Kondo temperature T K . γ values of between a few mJ/mol-K 2 and 2500 mJ/mol-K 2 are not uncommon in Kondo systems 38,42,43 . This value of γ for our samples at zero field is 102 ± 2 mJ/mol-K 2 . The high field value of γ of ≈32 ± 2 mJ/mol-K 2 is also large, indicating the presence of Kondo scattering. The application of a magnetic field suppresses the low-temperature rise in specific heat in our samples and shifts the anomaly to higher temperatures, as can be seen in Fig. 7(a). This can be analysed better by subtracting the electronic and lattice contributions to specific heat from the total specific heat. To do this, we first plot the zero field C/T versus T 2 and fit the high temperature part to the expression C/T = γ + βT 2 , where γT is the high field electronic contribution and βT 3 the lattice contribution to the specific heat. This fit is shown in the inset of Fig. 7(a). We obtain a value of γ = 28.3 ± 0.2 mJ/mol-K 2 for our sample which is still larger than the values of between 3 and 16 mJ/ mol-K 2 reported for other TMDCs 44 . The Debye temperature θ D can be calculated from the β value using the expression β π θ = nR 12 /5 D 4 3 , where R is the universal gas constant and n the number of atoms per molecule. The β value of 1.119 ± 0.005 mJ/mol-K 4 for our sample gives a Debye temperature, θ D of 173.3 ± 0.3 K which is comparable to that of 213 K reported for VSe 2 35 . The values of the electronic and lattice contribution to specific heat obtained from the fit are then subtracted from the total specific heat measured at different fields including zero field, assuming that the γ and β values do not change with field. This gives the excess specific heat ΔC which is divided by temperature and plotted against temperature as shown in Fig. 7(b). From the figure, it is clear that the magnitude of the maximum in the excess specific heat at low temperatures goes down as the field is increased and that the peak shifts to higher temperatures as the field increases. This kind of behaviour of the excess specific heat is consistent with the evidence for RKKY mediated magnetic interactions seen in both the low-field low-temperature χ(T) data shown in Fig. 6(b) and the low-temperature ρ(T) curve presented in Fig. 3(b). The area under the ΔC/T versus T gives the entropy associated with this Schottky-like anomaly and for our sample at 30 kOe this is 92 ± 2 mJ/mol-K which amounts to 0.010 of the entropy Rln(2S + 1) associated with a mole of paramagnetic spins with spin S = 1, where R is the molar gas constant (8.314 J/mol-K). It should be noted that the spin corresponding to the magnetic moment of 2.5 Bohr magneton possessed by the interlayer V ions is S = 0.85. Bearing in mind the susceptibility results from which we found the molar fraction of interlayer V ions in the single crystals to be 0.0357, and the fact that the method used may underestimate the entropy associated with this anomaly, the heat capacity data are consistent with anomaly arising from magnetic cluster of the interlayer V ions, while an enhanced γ even at 90 kOe, is consistent with presence of Kondo scattering.

Discussion
The upturn in the low-temperature resistivity in our samples is best described by the Kondo effect and the fit to the low -temperature resistivity using the Hamann equation with an effective temperature that takes into account the RKKY interactions between the magnetic impurities gives a convincing description of the ρ(T) data. The temperature dependence of the resistivity cannot be fit to the behaviour expected in the case of electron-electron interactions or weak localization effects. Moreover, the negative magnetoresistance seen in our sample at low fields and at temperatures below the minimum in the resistivity is further evidence against electron-electron interactions being the source of the upturn in resistivity as the upturn would then be expected to be largely insensitive to magnetic fields 21,25 and only exhibit a small positive magnetoresistance 23,24 . For both the Kondo effect and weak localization, the low-temperature rise in resistivity is expected to disappear in a magnetic field 21,23 and both effects lead to a negative magnetoresistance 24,26 . However, the negative magnetoresistance in our samples also shows a quadratic dependence with the magnetic field at low fields, similar to what is seen in InSb due to the Kondo effect 26 . Furthermore, the magnetoresistance is not sensitive to the angle that the magnetic field makes with the plane of the sample, which is also expected in the case of a Kondo effect with an isotropic s-d exchange interaction. These facts taken together suggest that the negative magnetoresistance originates from a Kondo effect. The presence of paramagnetic interlayer V ions indicated by the magnetization and heat capacity measurements means that the localized magnetic moments needed for the scattering of the conduction electrons are present in our samples. Since, the presence of interlayer V ions is necessary for the Kondo effect it implies a certain degree of disorder is required for this effect to be seen. The localised moments appear to suppress the anomaly at the CDW transition seen in the χ(T) curves for the polycrystalline sample. A similar behaviour has been observed in the χ(T) data for 2H-Pd x TaSe 2 system, where the small local moments induced by the Pd intercalation lead to a pronounced Curie-Weiss tail in the magnetic susceptibility at low temperature and wash out the transition seen in χ(T) at the incommensurate CDW transition around ~112 K in the pure 2H-TaSe 2 material 45 . The enhanced γ, even at 90 kOe, provides further support for the Kondo effect. The low temperature deviation from a Curie-Weiss law in the susceptibility and a step in the susceptibility, along with the excess specific heat, provide further evidence for the presence of magnetic correlations at low temperature in our single crystal samples of VSe 2 . Deviations from a Curie-Weiss law may also be due to the presence of a Kondo screening of the magnetic moments on the magnetic impurities. It is to be noted that in another TMDC, 1T-TaS 2 , an increase in resistivity with decreasing temperature was reported below 60 K 46 . However, the magnetoresistance was always negative for the whole range of field, except for a small temperature range between 0.1 and 2.2 K where it was positive at low fields. This was interpreted to be the result of variable range hopping in the Anderson localized states 47 . In TMDCs generally a positive magnetoresistance is seen 48,49 . In our samples at low temperatures, the negative magnetoresistance from the Kondo effect and a linear positive magnetoresistance compete with each other, with the positive magnetoresistance finally dominating at higher fields. The fact that the Kondo effect disappears as the temperature rises, explains why only a positive magnetoresistance is seen at higher temperatures. The signatures of the Kondo effect seen in the low-temperature resistivity and magnetoresistance of our single crystal samples of VSe 2 are not widely reported in other TMDCs and this makes an interesting addition to the already diverse list of properties of the TMDCs.

Methods
Polycrystalline powder and single crystal growth. VSe 2 samples were first synthesized in polycrystalline form by reacting powders of the starting materials, vanadium (Alfa Aesar, 99.5%) and selenium (Alfa Aesar, 99.999%). Stoichiometric amounts of the powders were ground together and transferred into quartz tubes, which were evacuated and sealed. The polycrystalline mixture was then heat treated for several days at 750 °C, followed by a second heating at temperatures in the range 815-825 °C for several days. Single crystals of VSe 2 were grown by the chemical vapour transport method by using the polycrystalline VSe 2 powder as the starting material. These powders were sealed in a quartz tube under vacuum and then kept in a temperature gradient using the procedure described in Refs. 4, 15, 50. Iodine was used as the transporting agent.
Sample characterization. Phase composition analysis was carried out using a Bruker D5005 X-ray diffractometer using Cu Kα 1 and Kα 2 radiation (λ Kα1 = 1.5406 Å and λ Kα2 = 1.5444 Å). The powder X-ray diffraction patterns were collected at room temperature over an angular range of 8 to 110° in 2θ with a step size in the scattering angle 2θ of 0.02°, and a total scanning time of 17 hours. The analysis of the powder X-ray patterns was performed using the FullProf software suite 27 . Single crystal X-ray diffraction measurements were performed by mounting a small piece of single crystal measuring 12 μm × 43 μm × 126 μm, cleaved from the bulk, on a Gemini R CCD X-ray diffractometer (Agilent Technologies). The data were collected using Mo Kα radiation (λ = 0.71073 Å) at room temperature, and were then indexed and integrated using CrysAlisPro (Agilent Technologies). The structure was refined using Shelx 51 . 398 reflections were measured, of which 81 were independent reflections. An extinction parameter was refined and the obtained value was 0.04 (6). Attempts to fit the single crystal X-ray data using different models for the occupancies of the V and Se sites, showed that the best model for the refinements is one in which the occupancies of both V and Se sites are fixed. The refinement was made using isotropic atomic displacement for both V and Se sites. A Laue X-ray imaging system with a Photonic-Science Laue camera was used to investigate the quality of the crystals prepared. The composition of the single crystals was also determined using an energy dispersive X-ray spectrometer in a Zeiss SUPRA 55-VP scanning electron microscope.
Transport and Magnetic Measurements. Resistivity and magnetoresistance measurements down to 2 K were performed using an AC transport measurement option in a Quantum Design Physical Property Measurement System (PPMS). DC resistivity measurements down to 0.5 K were performed using a DC transport measurement option in a PPMS with a Helium-3 insert. Electrical contacts were made to the sample in a four probe configuration with silver wires and conductive silver paint for the resistivity and magnetoresistance measurements. The Hall resistivity was measured using a five wire Hall measurement configuration in the AC transport option. Currents of 1 to10 mA were used in the AC resistivity and magnetoresistance measurements, 0.1 to 0.5 mA for the low-temperature DC resistivity measurements, while in the Hall measurements currents of 0.05 mA were used. The angle dependence of the magnetoresistance was measured using the horizontal rotator option in the PPMS. The DC magnetization measurements were performed using a Quantum Design Magnetic Property Measurement System MPMS-5S superconducting quantum interference device (SQUID) magnetometer for temperatures between 2 and 300 K and in an MPMS magnetometer with an iQuantum 3 He system in the temperature range 0.5 to 2 K. Low temperature heat capacity measurements were made using a two-tau relaxation method in a Quantum Design PPMS equipped with a 3 He insert.