Establishing magneto-structural relationships in the solid solutions of the skyrmion hosting family of materials: GaV4S8−ySey

The GaV4S8−ySey (y = 0 to 8) family of materials have been synthesized in both polycrystalline and single crystal form, and their structural and magnetic properties thoroughly investigated. Each of these materials crystallizes in the F\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{\bf{4}}}$$\end{document}4ˉ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{\bf{4}}}$$\end{document}4ˉ3m space group at ambient temperature. However, in contrast to the end members GaV4S8 and GaV4Se8, that undergo a structural transition to the R3m space group at 42 and 41 K respectively, the solid solutions (y = 1 to 7) retain cubic symmetry down to 1.5 K. In zero applied field the end members of the family order ferromagnetically at 13 K (GaV4S8) and 18 K (GaV4Se8), while the intermediate compounds exhibit a spin-glass-like ground state. We demonstrate that the magnetic structure of GaV4S8 shows localization of spins on the V cations, indicating that a charge ordering mechanism drives the structural phase transition. We conclude that the observation of both structural and ferromagnetic transitions in the end members of the series in zero field is a prerequisite for the stabilization of a skyrmion phase, and discuss how the absence of these transitions in the y = 1 to 7 materials can be explained by their structural properties.

Skyrmions attract widespread attention partly as promising candidates for energy-efficient high-density data storage and in spintronic devices 37,38 . Skyrmions are magnetic vortices that have been observed in magnetic materials in two forms: Bloch skyrmions and Néel skyrmions, so-named because a cut across the diameter of the skyrmion reveals either a Bloch or Néel wall of spins. Magnetic skyrmions have been observed chiefly in chiral B20 and β-Mn ferromagnets, in the form of a lattice of Bloch-type skyrmions 39 . In contrast, lattices of Néel-type skyrmions have been found in GaV 4 S 8 and GaV 4 Se 8 . In fact, the Néel-type skyrmion lattice in the bulk has, thus far, only been observed in three materials: GaV 4 S 8 , GaV 4 Se 8 , and VOSe 2 O 5 40 . There is therefore a justifiable interest in enlarging the family of materials hosting the Néel-type skyrmion, in order to acquire a better understanding of these spin textures. Alteration of the chemical composition, e.g. substitution of a transition metal in Cu 2 OSeO 3 , has proven to be an effective way to expand families of skyrmion-hosting materials and obtain valuable insights into the physics of these systems [41][42][43] . Given the many similarities in the two end members of the series: GaV 4 S 8 and GaV 4 Se 8 (i.e. similar structural phase transitions as well as exhibiting ferromagnetic-like ground states and hosting skyrmion phases), it is interesting to explore the properties of this family further by investigating the intermediate compositions lying between the GaV 4 S 8 and GaV 4 Se 8 end members.

Results
The phase purity of polycrystalline materials used for the crystal growth was determined by powder X-ray diffraction. The diffraction profiles of the GaV 4 S 8−y Se y samples with y ranging between 0 and 8 could be indexed in a cubic F43m symmetry, as shown in Fig. 1a and Supplementary Figs. S2 and S3. The lattice parameter a increases monotonically with the increasing Se substitution levels, as expected when substituting S 2− (ionic radius of 1.70 Å) with larger Se 2− (1.98 Å) anions (Supplementary Table S1) 46 . Materials with y between 0 and 5 are phase pure, while small amounts (up to 3.3%) of VSe 2 impurity were observed in the GaV 4 SSe 7 and GaV 4 Se 8 samples. In contrast to the rest of the materials investigated, the presence of two cubic phases with slightly different lattice parameters was detected in the GaV 4 S 2 Se 6 sample. The refinement of S and Se occupancies revealed higher concentrations of sulphur in GaV 4 S 8−y Se y materials than the nominal concentrations used during the synthesis, except in GaV 4 S 3 Se 5 . Moreover, the Se 2− is found to preferentially occupy the site near the Ga cations (GaX 4 tetrahedron), while S 2− is found to prefer the site near V cations (V 4 X 4 heterocubane).
The compositions of single crystals with y between 0 and 8 were determined by refining single crystal X-ray diffraction data, and were further checked using energy dispersive X-ray measurements as shown in Table 1. The obtained lattice parameters can be seen to follow Vegard's law (Fig. 1b). Single crystal X-ray diffraction shows that GaV 4 S 8−y Se y (y = 0 to 8) adopts the same crystal structure down to 50 K, where GaX 4 and V 4 X 4 units exhibit the NaCl-like arrangement as shown in Fig. 1c. In contrast to the two end compounds, GaV 4 S 8 and GaV 4 Se 8 , substitutional disorder is observed in all the intermediate solid solutions examined. To model the substitutional disorder properly, different refinement approaches were used: either allowing the atomic positions of S and Se to refine separately or constraining them. The most reliable structural models were obtained when S2 and Se2 in the GaX 4 tetrahedra were fixed to occupy the same atomic position, and the atomic positions of S1 and Se1 in the V 4 X 4 unit were allowed to refine separately, but were constrained by the crystal symmetry. The occupancies of sulphur and selenium were allowed to be freely refined, with the total occupancy constrained to 1, and there is clear evidence that S prefers to occupy the atomic position in the V 4 X 4 cluster, while the Se prefers to be situated in the GaX 4 tetrahedron (Supplementary Table S2 and S3). These findings are slightly at odds with the hard-soft acids-bases 47 concept, where Ga 3+ is classified as a harder acid than V 3+ , and S 2− is classified as a stronger base than Se 2− , and therefore, it would be expected that Se 2− would prefer to bind to V 3+ , and S 2− to bind to Ga 3+ . However, Se 2− is considerably bigger than S 2− , hence the tetrahedral coordination in the GaX 4 tetrahedron is preferred over the trigonal coordination in the V 4 X 4 cluster. Table 1. Composition of the single crystals obtained from EDX analysis and from the refinement of single crystal X-ray diffraction (SCXRD) data.
www.nature.com/scientificreports www.nature.com/scientificreports/ (y between 1 and 7) do not show this trend (Supplementary Table S3). A close examination of the structural models reveals that V-V distances within the cluster increase much more slowly than expected as a function of y, while V-S2/Se2 distances increase slightly faster. Moreover, a higher deviation from the linearity in V-V distance can be observed in the material y = 4. This can be partially explained by the S1 and Se1 site occupancies. However, the trend of the Se1 site occupancy does not exactly follow the trend of the V-V bond distances, as shown in Fig. 2d. In contrast, the changes of the V-S2/Se2 bond distances are in very good agreement with changes of Se2 site occupancy. These findings indicate that the tendency of S 2− and Se 2− to occupy the atomic positions in the V 4 X 4 cluster and in the GaX 4 tetrahedron respectively is the driving force for the structural motifs observed in the GaV 4 S 8−y Se y solid solutions. Additionally, the volumes of the GaX 4 and V 4 tetrahedra increase with increasing Se substitution levels (as expected, the volume of V 4 tetrahedra follow the same trend as V-V bond distances), however, the volume of (S1) 4 and (Se1) 4 tetrahedra stays nearly constant throughout the series (Fig. S4).
The end compounds, GaV 4 S 8 and GaV 4 Se 8 , undergo a structural phase transition at 42 K and 41 K respectively. At this transition, one of the vanadium and one of the sulphur/selenium atoms move along the [111] crystallographic direction, as shown in Fig. 2a, reducing the F43m to R3m symmetry, where the cubic diagonal becomes the c crystallographic lattice. Lattice parameters of a = b = 6.8098(3) Å and c = 16.7100(17) Å were obtained for GaV 4 S 8 at 30 K, below the structural transition, and are in good agreement with the reported values 8 . Three of the V1-V2 bond distances (2.8927(13) Å) are elongated in comparison to the three V2-V2 distances (2.8759(10) Å) and the rhombohedral angle of 59.62(3)° deviates from 60° observed in cubic symmetry; see Fig. 2. Also, three V1-S2 (2.2920(12) Å) and three V2-S1 (2.2923(16) Å) bond distances are longer than the remaining six V2-S2 (2.2882(13) Å) bond distances as shown in Fig. 2. Additionally, the GaS 4 tetrahedron distorts along the crystallographic c direction with one Ga-S4 (2.271(2) Å) bond distance elongated compared to the remaining three Ga-S3 (2.2664(10) Å) bond distances (Fig. 2). The observed structural distortion is in very good agreement with that reported previously 8,11 . However, the changes in bond distances obtained in single crystals are less pronounced than those seen in polycrystalline samples.

Neutron diffraction.
To complement the X-ray diffraction data, powder neutron diffraction was performed on the end compound (GaV 4 S 8 ) and one intermediate solid solution (GaV 4 S 6 Se 2 ) with representative properties of the solid solutions, to obtain the nuclear structures at a range of temperatures. Powder neutron diffraction data obtained on GaV 4 S 8 measured at 50 K, above the structural transition, refine in the high temperature F43m space group (Supplementary Table S4). Due to the small coherent scattering cross section of V, the positions of the V atoms were constrained to those obtained in single crystal X-ray diffraction. The evidence of a structural phase transition from cubic to a pseudo-cubic, rhombohedral, structure can be clearly seen through the splitting of some of the cubic Bragg peaks such as the cubic (448) peak as shown in Fig. 3a. Below the structural phase transition at 41 K, the crystal structure remains the same until 1.5 K and refines in the R3m space group, as shown in Supplementary Table S4. These results are in good agreement with those reported in the literature 8 . Powder neutron diffraction data for GaV 4 S 6 Se 2 show the absence of a structural phase transition, with data refining in a F43m space group down to 1.5 K, as shown in Supplementary Fig. S5. The site location and occupancy of S 2− and www.nature.com/scientificreports www.nature.com/scientificreports/ Se 2− within the V 4 X 4 units are in good agreement with X-ray diffraction results. Similar to the X-ray diffraction data, S1 and Se1 are found to occupy different sites within the V 4 X 4 cluster with a higher fraction of S 2− located within the V 4 X 4 cluster and Se 2− preferentially occupying sites in the GaX 4 tetrahedron.
Similar to powder neutron diffraction, one single crystal of the end composition (GaV 4 Se 8 ) and one single crystal of a member of the intermediate materials (GaV 4 S 4 Se 4 ) with indicative properties of the set of intermediate solid solutions were chosen for investigation. The data obtained for GaV 4 Se 8 refines in a F43m space group at 55 K whereas the data between 30 to 1.5 K refine in a R3m space group, in very good agreement with both our X-ray data and the literature 14 (see Supplementary Table S6). As in the analysis of the powder neutron diffraction, the positions of the vanadium atoms were constrained to those obtained from single crystal X-ray diffraction. The low temperature measurements show that GaV 4 S 4 Se 4 does not undergo a structural phase transition, and remains cubic down to 1.5 K, similar to all the other intermediate compositions examined. The nominal GaV 4 S 4 Se 4 single crystal is found to be GaV 4 S 4.36 Se 3.64 ,based on the refinement of the S and Se site occupancies, which is in excellent agreement with single crystal x-ray diffraction (SCXRD) and EDX (Table 1). Moreover, the individual site occupancies show that there is a higher fraction of S atoms situated within the V 4 X 4 cluster and a lower fraction in the Ga 4 X 4 tetrahedron, also in good agreement with our powder neutron diffraction and X-ray diffraction data.

Magnetic neutron scattering.
Magnetic structures for the end compound GaV 4 S 8 and two representative members of the intermediate solid solutions (GaV 4 S 6 Se 2 and GaV 4 S 4 Se 4 ) were probed using both powder and single crystal diffraction. The magnetic scattering was isolated by taking the difference pattern between 1.5 and 15 K, above and below the magnetic transition. For GaV 4 S 8 , the magnetic peaks are located at the same positions as the nuclear Bragg peak positions and have a symmetrical peak shape indicating a k = 0 propagation vector as seen in Fig. 3c. Magnetic symmetry analysis with help of the program BASIREPS 49,50 was used to determine the allowed irreducible representations (IR) and their basis vectors (BV) for the two vanadium sites. Testing the different allowed IRs showed that the magnetic difference pattern for GaV 4 S 8 refines best in a R3m' magnetic structure, shown in Fig. 3b. The V1 atoms refine to have a magnetic moment of 0.30(6) μ B , pointing along the hexagonal c axis whereas the V2 magnetic moments have a smaller moment of 0.18(7) μ B and are canted outward from the cluster. The localization of a higher moment on the distorted V1 site is consistent with density functional theory (DFT) calculations showing a higher spin density residing on the distorted V1 atom as opposed to the V2 atoms 16,51 . It is worth noting that, due to the nature of powder neutron diffraction, it cannot be established if the orientations of the moments are in the configuration shown in Fig. 3b or pointing in the opposite direction with www.nature.com/scientificreports www.nature.com/scientificreports/ the same magnitude, i.e. moments on the V2 atoms pointing in towards the cluster with the moment on the V1 atom pointing in the negative crystallographic c axis. Attempts were made to refine the structure with only one magnetic moment per V 4 cluster to be localized either on the V1 site or in the centre of the cluster. However, a single magnetic moment per V 4 cluster does not produce the correct ratios of peak intensities, supporting the fact that both V1 and V2 carry a magnetic moment.
Powder neutron diffraction patterns of GaV 4 S 6 Se 2 show no difference between the diffractograms measured at 1.5 and 15 K, as shown in Supplementary Fig. S5, confirming the absence of any magnetic intensity. Additionally, single crystal neutron diffraction of GaV 4 S 4 Se 4 revealed no magnetic scattering intensity at 1.5 K. There is a strong indication that no long-range magnetic ordering is present in solid solutions based on our results obtained on the two representative materials.
Magnetic measurements. Bulk magnetization measurements on single crystals were carried out to identify the nature of magnetic ordering occurring within the GaV 4 S 8−y Se y family of materials. Obvious differences in the magnetic susceptibilities across the GaV 4 S 8−y Se y family are observed, as shown in Fig. 4a. The end compounds, GaV 4 S 8 and GaV 4 Se 8 , undergo long-range ferromagnetic-like ordering at 13 and 18 K, respectively, with GaV 4 S 8 having a susceptibility over twice as large as its selenium counterpart GaV 4 Se 8 8,14 . The structural phase transitions can be observed as large discontinuities in the inverse magnetic susceptibilities, as shown in Supplementary Fig. S6, as was previously reported in the powder materials 8,10,11,14,27,44 . This highlights the magneto-structural link between the cubic and rhombohedral phase, with the discontinuities likely related to the lifting of the degeneracy of V orbitals from the structural phase transition leading to a change in the location of spin density and hence magnetic interactions. The effective moment calculated from the susceptibility by using the Curie law μ χ = .
T 797 8 m eff is given in Fig. 4c. A discontinuity in the effective moment μ eff can be seen at the temperature of the structural phase transition in GaV 4 S 8 (42 K) and GaV 4 Se 8 (41 K), and a large increase in μ eff is seen with decreasing temperature. The peak in effective moment is in good agreement with the onset of long-range magnetic order. Above the structural phase transition in GaV 4 S 8 and GaV 4 Se 8 , μ eff continuously increases in the temperature region between ~50 and 300 K, as previously reported 8 . This means that the Weiss temperature cannot be accurately extracted from χ vs T measurements, although it can be seen that the interactions change from antiferromagnetic to ferromagnetic at the structural phase transition. Based on the M vs μ 0 H www.nature.com/scientificreports www.nature.com/scientificreports/ measurements shown in Fig. 4d, the saturation magnetization at applied fields of 5 T and 1.8 K was found to be 0.87(2) μ B /f.u. and 0.85(2) μ B /f.u. for GaV 4 S 8 and GaV 4 Se 8 , respectively. These values are in good agreement with those reported previously 10 . Additionally, the value of 0.87(2) μ B /f.u. for GaV 4 S 8 is also in good agreement with the value of 0.84 (7) μ B /f.u., the total magnetic moment per V 4 cluster, obtained from our neutron diffraction refinements.
A large suppression in both the magnetic transition temperature and the magnitude of the magnetization is evident in all of the intermediate compounds measured, as seen in Fig. 4a, suggesting that the absence of a structural phase transition has a dramatic effect on the magnetism. Examination of the magnetization at low temperature shows a bifurcation of the zero-field-cooled-warming (ZFCW) and field-cooled-warming (FCW) susceptibility curves, shown in Fig. 4b, indicative of either a canted antiferromagnetic ordering or a glassy transition. It is notable that similar behaviour has been seen in related compounds [10][11][12][13][14] . Our recent μ + SR measurements show no evidence for long-range magnetic ordering in the intermediate compounds and instead indicate the presence of a glass-like state 45 . In Fig. 4c, it can be seen that even though GaV 4 S 4 Se 4 does not undergo a structural phase transition, the effective moment starts to decrease rapidly at a similar temperature to where the structural phase transitions occur in GaV 4 S 8 and GaV 4 Se 8 , suggesting the onset of magnetic correlations. The effective moment behaves qualitatively very differently with temperature in the intermediate compounds, consistent with a lack of ordered magnetic response and the likely formation of a low-temperature glassy state.
In order to investigate the magnetic behaviour of one of the intermediate compositions in depth, ac susceptibility measurements on GaV 4 S 4 Se 4 single crystals were performed. Figure 5a shows behaviour again indicative of a spin glass, with the peak of the in-phase and out-of-phase components of susceptibility shifting in temperature and amplitude, depending on the applied ac frequency. The peaks observed in the ac susceptibility data are caused by a freezing of the spins in the material, and hence correspond to a freezing temperature T f . The freezing temperature for GaV 4 S 4 Se 4 can be seen to be ~2.4 K. The frequency dependence of the peaks in χ′ ac can be characterized by a shift of the freezing temperature per decade of frequency given by , no sharp lambda anomaly, indicative of the structural phase transition, can been seen for GaV 4 S 4 Se 4 , as shown in Fig. 5b. Moreover, at low temperatures a broad peak can be seen centred at ~3.1 K. The position of this peak is located at ~20% above T f , inline with known spin glass materials 52 , and indicates the slowing down and freezing of magnetic moments above T f . The phonon contributions were estimated using one Debye and three Einstein modes. Fitting was performed in the temperature region between 45 to 300 K with a Debye temperature of 201(4) K, and three Einstein frequencies corresponding to temperatures of 263(8), 465 (8), and 508(13) K. By subtracting the phonon contribution, the associated magnetic entropy can be calculated as shown in the inset of Fig. 5b. The magnetic entropy can be seen to be recovered up to approximately 42 K where the entropy is in good agreement with the value of Rln2, the entropy corresponding to a spin 1/2 system. These heat capacity results are consistent with those of magnetometry, indicating that there are short-range magnetic correlations occurring at temperatures far above the spin glass transition. The phonon subtracted heat capacity below T f shows T 3 2 like dependence deviating from the linear expected behaviour of a canonical spin glass.

Discussion
Our structural analysis of single crystals confirms that the end compounds, GaV 4 S 8 and GaV 4 Se 8 , undergo a structural phase transition from cubic (F43m) to rhombohedral (R3m) symmetry at 42 and 41 K respectively, as previously reported 8,14 . The vanadium and X atoms in the V 4 X 4 (X = S or Se) heterocubane-like entity face along all four [111] crystallographic directions in the cubic structure. At the structural phase transition one of the vanadium and the X pairs of atoms distorts, elongating along the cubic diagonal, reducing the symmetry from cubic to rhombohedral, as shown in Fig. 2. This distortion has an equal probability to occur along any of the four [111] crystallographic directions. (In addition, our powder neutron diffraction measurements indicate the presence of noticeable negative thermal expansion in GaV 4 S 8 in the rhombohedral phase (Table S4). In the temperature region between 30 and 9 K, the lattice parameter c increases by 0.04%, while the lattice parameter a contracts by 0.01%, resulting in a negative thermal expansion. The origin of this subtle negative thermal expansion is a subject of further investigation.) The structural phase transition is followed by long-range ferromagnetic-like ordering at 13 and 18 K for GaV 4 S 8 and GaV 4 Se 8 respectively. Our powder neutron diffraction measurements allowed the determination of the magnetic structure of GaV 4 S 8 revealing that, on the neutron time scale, spins are localized on the individual vanadium cations. A magnetic moment of 0.30(6) μ B on V1 points along the c axis, while a V2 cation moment of magnitude 0.18(7) μ B is canted outward from the V 4 cluster (Fig. 3b), yielding a magnetic moment of 0.84(7) μ B in the spin-polarized state. This value is in very good agreement with the magnetic moment of 0.87 μ B obtained from our magnetization measurements (M vs μ 0 H at 1.8 K) and with previously-reported values 8,10 . The spin arrangement on V cations in the V 4 S 4 cluster observed in our magnetic model is similar to the arrangement of spin density on V cations in GaV 4 S 8 obtained from DFT calculations 16,51,53 .
To date, the electronic state and the origin of long-range ferromagnetic-like ordering in GaV 4 S 8 and GaV 4 Se 8 have been explained via a recombination of d orbitals on the V cations, leading to the formation of a molecular system with S = 1/2 per V 4 cluster. In this picture the structural phase transition is driven by the Jahn-Teller effect, distorting one of the V and S atoms from an ideal cubic symmetry 8,44 . This model has been used to successfully explain a range of physical properties in this series, such as the polar relaxation dynamics observed in dielectric and THz spectroscopy 26 , whose characteristic behaviour can be described within the framework. However, it is hard to reconcile this model with the magnetic behaviour of these materials. We note that the experimentally-obtained magnetic moment (0.821 μ B ) is considerably lower than the 1 μ B anticipated for a S = 1/2 system. Additionally, the ionic radius of the V 3+ cation is 0.64 Å 46 , and so the bond distances between V cations within the cluster (~2.9 Å) are exceptionally long to permit a significant overlap of d orbitals, required for this model. Similarly, direct exchange interactions over the inter-cluster distance of ~3.9 Å are unlikely.
Here we present an alternative interpretation. It is notable that both intra-and inter-cluster V cations are bonded to bridging S anions with an angle of ~102°. This configuration permits ferromagnetic superexchange interactions between spins located on the V cations via S 2− or Se 2− 54,55 . We propose, therefore, that the low-temperature electronic state (and consequently the origin of long-range ferromagnetic ordering in GaV 4 S 8 ) follows from the structural phase transition at 42 K, where there is a freezing out of a delocalized electron onto one of the vanadium atoms, forming three V 3+ and one V 4+ cations (i.e. a charge-ordered state). Evidence for the presence of V 3+ and V 4+ in the ratio of 0.75:0.25 in both cubic and rhombohedral structures comes from bond-valence analysis performed using the SoftB software 56 . Charge ordering is also consistent with a sharp increase of the resistivity at the structural phase transition observed previously 8,27 .
Compared to the end compounds (GaV 4 S 8 and GaV 4 Se 8 ), the solid solutions (with y between 1 and 7) show quite distinct behaviour. Substitutional disorder is found to be present in all intermediate compounds (shown in Fig. 1c), where Se prefers to occupy the positions in the GaX 4 tetrahedron, while S prefers to occupy the positions in the V 4 X 4 heterocubane-like cluster. The solid solutions do not undergo a structural phase transition and do not exhibit long-range ferromagnetic magnetic ordering. However, an irreversible bifurcation between ZFCW and FCW traces in the χ vs T measurements is observed in these materials, that probably reflects the presence of a spin glass-like state. From the point of view of our interpretation, this implies that substitution of cations on either the B or A site has a profound effect on the ability of the material to undergo a structural phase transition. However, there is still an appreciable electronic response in this temperature region, demonstrated by the drastic decrease of μ eff , observed in the μ eff vs T curves depicted in Fig. 4c.
The mechanism preventing a structural phase transition in intermediate compositions (y = 1 to 7) requires further investigation. However, in general terms we can comment that the orbital degeneracy of V 3+ -based compounds would often be expected to be lifted by a Jahn-Teller (JT) distortion which, if it happens at each site via a coherent lattice distortion, implies orbital ordering 57 . Generally, three energy scales decide whether this occurs. These are the on-site Coulomb repulsion U, the electronic bandwidth W and intra-atomic Hund rule coupling J h . In an insulator,  U W and we expect a JT distortion; in a metal (where  W U), the JT distortion is suppressed. Where U and W are similar, it is also possible 57 that the Hund's rule coupling overcomes on-site repulsion and Scientific RepoRtS | (2020) 10:9813 | https://doi.org/10.1038/s41598-020-65676-9 www.nature.com/scientificreports www.nature.com/scientificreports/ leads to a charge ordering, involving double electron occupancy on some sites (V 3+ ) and single occupancy on others (V 4+ ). This charge ordering occurs via a mechanism involving a delocalization effect, where some of the electronic energy levels become band-like and consequently reduce the energy penalty for the double electron occupancy required for charge ordering 57 . In the compounds GaV 4 S 8 and GaV 4 Se 8 , we might expect that their insulating nature puts them in the JT regime. However, the possible delocalization of charges over V 4 clusters that occurs in these systems, could push the systems into the regime where U and W are similar and explain how the charge ordering mechanism leads to the distortion in these materials.
At a JT transition 58 the distortion of the structure costs elastic energy as αq 2 (where q is a displacement coordinate and α a constant) while the resulting splitting of the electronic states (through local crystal field levels and Hund's rule coupling) saves energy as −βq (where β is a constant). An equilibrium distortion (with q 0 = β/2α) leads to a net saving in energy and a JT distortion results. To avoid a JT distortion in the intermediate materials, we could imagine either (i) the distortion costs too much energy or (ii) the electronic energy saving is not large enough, or both. In the charge-ordering scenario, (i) still applies in cases where there is a distortion, but here the electronic energy saving occurs because delocalization leads to a reduction in the Coulomb energy cost of having two electrons on a site. In this case, any reduction in the ability of electrons to diffuse between sites can reduce delocalization and hence also suppress the transition.
Since the end members of the series have similar electronic behaviour, it seems unlikely that S/Se disorder radically alters the electronic structure, such that these materials lie in a different energetic regime. Indeed, the behavior of μ eff suggests that all of the compounds are susceptible to a change in the electronic structure around 40 K. (However, we note that the μ + SR results suggested that the ground states of the two end systems are quite different to each other, so a change in the electronic states with substitution cannot be ruled out.) To explain the lack of structural transition in the intermediate materials, it is more probable that either the S/Se disorder creates a large elastic energy cost for a structural distortion or, if the distortion is driven by charge ordering, that electron delocalization is suppressed by the S/Se disorder. In the first hypothesis, the uneven local distribution of S and Se in the V 4 X 4 units, where neighbouring units contain slightly different concentrations of S and Se, could prevent a structural distortion and therefore retain the material in the cubic symmetry. The second hypothesis involves the disorder shutting down electron diffusion pathways via S and Se anions, reducing electron delocalization and consequently preventing the mechanism for charge ordering of V cations, which seems to be the driving force for the structural phase transition. Finally we note that while our proposed charge-ordering scenario provides a plausible explanation for the magnetic properties discussed in this paper, whether it also provides a straightforward and complete description of other (non-magnetic) behaviour of this system remains an open question, deserving of future investigation.
All of the results of our investigation on this series, when taken together, provide a direct link between the structural phase transitions and the observation of magnetic transitions in the materials. This is the case in the two end members, GaV 4 S 8 and GaV 4 Se 8 , which exhibit both these features whilst the intermediate compositions investigated exhibit neither of these transitions. The two parent compounds still remain the only members which host a skyrmion phase, making them unique, and the observation of both of the transitions appear to be prerequisites for the stabilization of the skyrmion phases. Polycrystalline and single crystal samples of the entire series of compounds GaV 4 S 8−y Se y (y = 0 to 8) have been synthesized and single crystals of the materials obtained by the CVT technique. Both the polycrystalline and single crystal samples have been thoroughly investigated for the structural properties through the use of powder and single crystal x-ray as well as neutron diffraction. Additionally, the magnetic properties of the whole series of materials have been studied using dc and ac magnetic susceptibility, magnetization and heat capacity measurements. Our results show that while both GaV 4 S 8 and GaV 4 Se 8 both undergo structural phase transitions, followed by magnetic ordering at 13 and 18 K respectively, all the intermediate compositions retain their cubic symmetry down to the lowest temperature investigated (1.5 K). These intermediate compositions do not exhibit long-range magnetic order but are, however, found to exhibit a spin glass-like behaviour with freezing temperatures of about 2.5 K, as ascertained from ac susceptibility and heat capacity measurements. Examining the magnetic structure of GaV 4 S 8 in detail, there is a clear indication of the localization of spins on the V atoms within the V 4 X 4 cluster which appears to favour a charge ordering mechanism. The results presented form a thorough investigation of both the structural and magnetic properties of this entire family of interesting materials. In the search for additional members of this family beyond the two end compositions to host skyrmions, this study throws light on the crucial links that exist between the structure and magnetism in these materials. Given the subtle nature of the local structural distortions in the intermediate compositions as a function of temperature, one strategy that emerges from our study is the possibility that investigating materials with much smaller substitutional levels close to the two end members (y between 0.1 and 0.5, and 7.5 and 7.9) maybe promising. The study of these materials forms part of our ongoing investigations to provide additional insight into structure-property correlations in this interesting family.

Methods
Synthesis and crystal growth. Polycrystalline GaV 4 S 8−y Se y materials with y between 0 and 8 were synthesized by thoroughly grinding stoichiometric amounts of V (99.5%, metals basis, Sigma-Aldrich), S (99.99%, Acros Organics), and Se (99.999%, Alfa Aesar) powders inside an argon-filled glove box. The mixtures of powders were transferred into silica tubes along with stoichiometric amounts of Ga (99.999%, metals basis, Alfa Aesar) granules and 50 mg of I 2 (99.99%, Alfa Aesar; added to prevent the formation of V 5 S 8 ). The tubes were then evacuated and sealed. The mixtures were heated at a rate of 10 °C /h to 600 °C and at a rate of 5 °C/h to 820 °C, kept at this temperature for 200 h, followed by water quench cooling. This yielded black polycrystalline air-stable powders with metallic lustre.
Single crystals of GaV 4 S 8−y Se y were grown using the chemical vapour phase transport (CVT) technique. Polycrystalline materials (1.5 g), described above, and 3.3 mg/cc of the transporting agent I 2 or PtCl 2 , were sealed Powder and single crystal neutron diffraction. Two complementary sets of powder neutron diffraction experiments were performed at the Institute Laue Langevin (ILL) on the D2B and D20 diffractometers, to determine the nuclear and magnetic structures of two of the materials in the series, namely GaV 4 S 8 and GaV 4 S 6 Se 2 . The D2B diffractometer was used to obtain high-resolution diffractograms for nuclear structure characterization using a neutron wavelength of 1.594 Å over the temperature range 1.5−50 K. High-intensity powder neutron diffraction was performed on the D20 diffractometer using a wavelength of 2.417 Å to determine the magnetic scattering over the temperature range 1.5−50 K. The GaV 4 S 8 and GaV 4 S 6 Se 2 powder samples ranging from 3-4 g were sealed in thin-walled cylindrical vanadium containers of 8 mm diameter and placed in a standard orange cryostat. Structural and magnetic Rietveld refinements were carried out using the TOPAS academic v6.0 and FullProf software 63 .
Single crystal neutron diffraction experiments were preformed using the SXD diffractometer at the ISIS Neutron and Muon spallation source 64 . Data was collected on ~2 × 2 × 2 mm 3 crystals of GaV 4 S 4 Se 4 and GaV 4 Se 8 . The samples were mounted in an aluminum loop on top of an Al pin and were placed in a vanadium cryostat which was evacuated to ~10 −5 mbar. Data were processed using the SXD2001 software and intensities extracted using a least squares procedure 64 with 3D Gauss-ellipsoids, line fittings, and box fittings with time-of-flight asymmetry as the profile function. The nuclear structures were refined using the JANA software 65 .

Magnetic measurements. Quantum Design Magnetic Property Measurement System, MPMS-5S and
MPMS-7XL, superconducting quantum interference device (SQuID) magnetometers were used for the investigation of the magnetic properties of the polycrystalline samples and single crystals as a function of temperature and field. dc temperature-dependent magnetic susceptibility (χ) measurements were carried out in an applied field of 10 mT and 2 T in the temperature region between 1.8 and 300 K under both zero-field-cooled-warming (ZFCW) and field-cooled-warming (FCW) protocols. The magnetic moments were obtained from M vs B measurements in applied fields between 0 and 5 T at 1.8 K. ac frequency-dependent susceptibility (χ′ ac and χ′′ ac ) measurements were performed on non-oriented single crystals at zero applied magnetic field between 1.8 and 3.5 K.
Heat capacity measurements were performed in a Quantum Design Physical Property Measurement System from 1.8 ≤ T ≤ 300 K. A He 3 insert was used to obtain temperatures down to 0.5 K.