Two Different Mechanisms of Metal Cluster Formation in the

Apr 25, 2018 - The pyrochlore sublattice (Figure 1a) can be imagined by two types of planes with alternate stackings: one type is a two-dimensional tr...
2 downloads 6 Views 2MB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Two different mechanisms of metal cluster formation in the rhombohedral spinel structures: AlVO and CuZr S 2

4

1.86(1)

4

Mikhail V. Talanov Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.8b00151 • Publication Date (Web): 25 Apr 2018 Downloaded from http://pubs.acs.org on April 26, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Two different mechanisms of metal cluster formation in the rhombohedral spinel structures: AlV2O4 and CuZr1.86(1)S4 Mikhail V. Talanov Southern Federal University, Rostov-on-Don 344090, Russia

ABSTRACT: An unusual case of the structural mechanisms of inorganic crystals formation is considered theoretically on the example of the AlV2O4 and CuZr1.86(1)S4 isosymmetric rhombohedral spinels with the same character of the atom distribution on the Wyckoff positions. It is shown that the difference in structures is due to the critical influence of the order parameter value that describes the real (in the case of AlV2O4) and virtual (in the case of CuZr1.86(1)S4) structural phase transitions. Owing to the low value of the order parameter, the phase transition to the rhombohedral AlV2O4 structure is associated with relatively small atomic displacements and the formation of heptamer structural blocks in the V-sublattices. In contrast, the virtual phase transition to the rhombohedral CuZr1.86(1)S4 structure is associated with high value of the order parameter and giant atomic displacements. As a result, the Cu- and Zr-atoms form metal nanoclusters – a ‘‘bunch’’ of dimers. The CuZr1.86(1)S4 is a second substance (the first example is CuTi2S4) with a new type of self-organization of atoms, in which the metal clusters are formed

ACS Paragon Plus Environment

1

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 32

by the atoms from two (but not one, as in the known materials) geometrically frustrating spinel sublattices.

1. INTRODUCTION Geometrical frustration is an important feature of the electronic and structural organization of spinels giving rise to a rich variety of exotic condensed phases (such as spin liquid1,2, spiral spin liquid3, valence bond crystal4,5, orbital liquid6, orbital glass7 states and other) and some novel phenomena (such as heavy fermion behavior8,9, anomalous Hall effect induced by the chiral order10, large negative magnetoresistance11, superconductivity12,13 and so forth). Frustration is caused by competition between different interactions, and the energy cannot be simultaneously minimized on all interactions. The most important feature of the frustrated structures is the existence of highly degenerate ground states. The lifting of the degeneracy may be accompanied by structural deformation, atom-, charge- orbital-, spin-ordering. As results of such processes, unique states of material with unusual physical behavior may appear1,14,15. One of the most widespread type of crystals with a geometrically frustrated sublattices is AB2X4 spinel. The A-atoms occupy Wyckoff position 8a (site symmetry 3 m), B − Wyckoff position 16d (site symmetry 3m), X − Wyckoff position 32e (site symmetry 3m) in the cubic Fd 3 m-phase of normal spinel. The spinel structure can be viewed as interpenetrating sublattices with B atoms, forming three-dimensional network of tetrahedra (the so-called pyrochlore sublattice), while the A atoms constitute a diamond lattice (Figure 1a). The pyrochlore sublattice (Figure 1a) can be imagined by two types of planes with alternate stackings: one type is a twodimensional triangular lattice (Figure 1b,c), consisting of the apical B ions of the tetrahedra, while the other type is a two-dimensional kagome lattice, consisting of the B triangles (Figure

ACS Paragon Plus Environment

2

Page 3 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

1b)1,15,16. It has been recently demonstrated that the tetrahedral site (A-site) of the normal spinel structure, containing transition element atoms, is also strongly frustrated in relation both to spins6 and orbitals17.

Figure 1. The cubic spinel structure fragments. A- and B- cation tetrahedra in the spinel structure (a). A atoms and B atoms are shown by green and deep blue color respectively; X-anions are not shown. In the spinel structure, B-sublattices are formed by a network of corner-sharing Btetrahedra (deep blue color). The B-sublattice can be visualized as alternating kagome (b) and triangular (b, c) layers of structure along a [111] direction as designated in (b). Kagome layers are shown by deep blue, triangular layers are shown by dark yellow color.

One of the characteristic structural phenomena, particularly observed in geometrically frustrated spinel systems, is the spontaneous formation of high-order atom structures – unusual molecular metal clusters or orbital molecules18. Well known examples of such atom clusters in ordered spinel structures are Ti2 dimers in MgTi2O419,20, decagons in Na4Ir3O821,22 (such decagons are also in Na3Ir3O823 and in other ordered inorganic structures24), trimers in LiVO225, «trimerons» in Fe3O426-28, heptamers29 or trimers + tetramers30 in AlV2O4, octamers in CuIr2S431, tetrahedra in LiXY4O8 (X = Ga, Fe, In; Y = Cr, Rh), Cr4GaLiO8, Cr4InLiO8, Rh4InLiO8, Ag0.5In0.5Cr2S4 and Cu0.5In0.5Cr2S432–37 and in others. The formation of metal clusters is

ACS Paragon Plus Environment

3

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 32

accompanied by the charge and/or spin and orbital orderings. All this is a result of structural (and/or magnetostructural) transition in spinels. It is important to note that clustering in all these materials occurs only in B-sublattices of the ordered spinels. Recently, we have first found the existence of a new exotic type of high-order clusters in the CuTi2S4 rhombohedral modification with R 3 m space group. Each of these clusters is formed not only by B-atoms, but by A-atoms spinel also. These metal clusters were called a "bunch" of dimers38. One can suppose that the "bunch" of dimers should exist in CuZr1.86(1)S4 structure, since this substance is isostructural with CuTi2S439. The structure of the AlV2O4 rhombohedral phase is characterized by the same space group as the CuZr1.86(1)S4 structure40. The atom distribution on the Wyckoff positions is the same for both phases, but the structures of these phases differ noticeably. This unusual and intriguing situation needs theoretical analysis. This research has two aims. One goal is a new interpretation of the clustering in the CuZr1.86(1)S4 structure. Within this interpretation, we first prove the existence of (Cu-Zr) clusters similar to clusters in CuTi2S4 structure. Other purpose of this paper is to study the reasons for the substantial differences between isosymmetric spinel-like structures (CuZr1.86(1)S4 and AlV2O4). To realize the latter, we carry out a detailed group-theoretical and structural study of the rhombohedral phases formation mechanisms from a cubic spinel structure as well as clustering in CuZr1.86(1)S4 and AlV2O4 structures. As a result, we expect to determine two fundamentally different mechanisms of the metallic clusters formation in rhombohedral spinel modifications. In our study, we use a combination of ISOTROPY41/ISODISTORT42 and AMPLIMODES43,44 computer software to analyze the structural phase transition from the hightemperature spinel phase with the Fd 3 m space group to the low-temperature phase with the R 3 m space group. In the case of AlV2O4, this structural phase transition is real and was detected

ACS Paragon Plus Environment

4

Page 5 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

by synchrotron x-ray and electron diffraction experiments at T 20 GPa on hydrostatic46 and non-hydrostatic conditions47, and also with an increase of the Cr concentration in B-sublattice at x> 0.148. Note that the temperature-induced phase transition Fd 3 m ↔ R 3 m accompanied with charge ordering in AlV2O4 was observed by Matsuno et al. for the first time49, but in more recent works the distribution of the atoms on the Wyckoff positions of R 3 m-phase was revised. On the other hand, there are no any phase transitions in rhombohedral CuZr1.86(1)S440 and the structural phase transition from parent phase with spinel structure is virtual. By using the concepts of one critical irreducible representation (irrep)50 and the parent phase (praphase, archetype), we derive the theoretical structure (or more precisely, motif of crystal structure) of low-symmetry R 3 m-structure from high-symmetry spinel structure (parent Fd 3 m-phase) and find out the different role of critical and noncritical (improper) order parameters in the rhombohedral spinel modifications formation. 2. RESULTS AND DISCUSSION 2.1. Order parameter symmetry. The spinel structure has a face-centred cubic lattice. For this lattice, the first Brillouin zone represents a body-centred lattice and contains four points of high symmetry, namely: k11(Г), k10(X), k9(L) and k8(W)51. For these points there are stars of the following wavevectors: k11(Г) = 0; k10(X) = 1/2(b1+ b2); k9(L) = ½(b1 + b2 + b3); k8(W) = b1/4 - b2/4 + b3/2 There are ten irreps of the wavevector k11(Г) (four one-dimensional, two twodimensional and four three-dimensional); four six-dimensional irreps of the wavevector k10(X), four four- and two eight-dimensional irreps of the wavevector k9(L), and two 12-dimensional irreps of the wavevector k8(W). Thus, totally there are 22 irreps. All irreps of the wavevector

ACS Paragon Plus Environment

5

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 32

k8(W) and also two irreps of the wavevector k10(X) (τ3 and τ4) do not satisfy the Lifshitz criterion, i.e. they induce incommensurate phases50. Since the considered rhombohedral structure is commensurate phase, we exclude these irreps from our analysis. Based on the results of the group-theoretical analysis of the phase transitions occurring according to one critical irrep in the group Fd 3 m, which satisfies the Lifshitz criterion50, we find that the phase with the R 3 m space group may be generated by three-dimensional irrep k11(τ7), as well as by two six-dimensional irreps k10(τ1) and k10(τ3), two four- dimensional irreps k9(τ1) and k9(τ4). The expression kj(τi) means the star of wavevectors kj, where i is number of corresponding irrep τ for given star j (according to Kovalev51). The correspondence of notation by Kovalev51 and Miller & Love52 is as follows53. The analysis shows that only in the case of phase, generated by irrep k9(τ4), the calculated distribution of the atoms on the Wyckoff positions of R 3 m-phase is consistent with experimental data on AlV2O429 and CuZr1.86(1)S440. As a result of group-theoretical calculations, the theoretical structure formula of the R3 m

low-symmetry

rhombohedral

spinel

modification

should

be

A2c1/2(1)A2c1/2(2)B1a1/4(1)B1b1/4(2)B6h3/2(3)X2c1/2(1)X2c1/2(2)X6h3/2(3)X6h3/2(4) (the rhombohedral presentation) or A6c1/2(1)A6c1/2(2)B3a1/4(1)B3b1/4(2)B18h3/2(3)X6c1/2(1)X6c1/2(2)X18h3/2(3)X18h3/2 (the hexahedral presentation). The theoretical formula agrees with experimental data29,40. 2.2. Structural mechanisms of rhombohedral modifications formation. The physical meaning of the order parameters that are transformed according to the corresponding critical and improper irreps, is determined by their entry into the mechanical and permutation reducible representations of the high-symmetry structure of the parent Fd 3 m-phase. The critical irrep k9(τ4) enters into both the mechanical representation and permutation representation of the spinel structure on the Wyckoff positions 8a, 16d and 32e22. This means

ACS Paragon Plus Environment

6

Page 7 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

that the low-symmetry R 3 m-phase formation is the result of the tetrahedral and octahedral cations and anions displacements as well as the ordering of all atom types. Group-theoretical analysis showed that the formation of the R 3 m-rhombohedral phase is due to the following types of atom ordering: – binary tetrahedral cation ordering (the type of order is 1:1); – ternary octahedral cation ordering (the type of order is 1:1:6); – quarternary anion ordering (the type of order is 1:1:3:3). The contribution of a critical irrep in crystal structure formation defines the symmetry of the low-symmetry phase completely. However, when the low-symmetry structure is far from the phase transition temperature, then the contribution of improper representations can become essential54. Here, the interpretation of some experimental data cannot be carried out only by means of a critical irrep. A method of finding improper atomic displacements and ordering is proposed54. In case of the Fd 3 m → R 3 m phase transition, besides the critical order parameter k9(τ4), the total condensate of the order parameters includes also two improper order parameters with the wave vector k11 = 0 from the center of the Brillouin zone: k11(τ1) and k11(τ7). The improper order parameters are due to nonlinear interactions between different degrees of freedom in the crystal. They cause additional distortions of the crystal structure, which can be significant in magnitude far from the phase transition temperature. Improper irrep k11(τ7) enters into the mechanical representation on the Wyckoff positions 8a and 32e and into the permutation representation on the Wyckoff positions 16d and 32e. This means that we can expect additional displacements of tetrahedral cations and anions and ordering of octahedral cations (the type of order is 1: 3) and anions (the type of order is 1: 3).

ACS Paragon Plus Environment

7

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 32

Improper irrep k11(τ1) enters into the mechanical representation on the Wyckoff position 32e and into permutation representation on the Wyckoff positions 8a, 16d and 32e. This means that additional ordering of all types of atoms is possible, but additional displacements can be expected for anions only. 2.3. Exploring the role of critical and improper order parameters to the rhombohedral spinel modifications formation. We calculated the contribution of the critical and improper order parameters to the formation of rhombohedral phases (Table 1). The contributions of critical and improper order parameters to the low-symmetry crystal structures formation of AlV2O4 and CuZr1.86(1)S4 differ significantly.

Table 1. The role of critical and improper order parameters in the AlV2O4 and CuZr1.86(1)S4 rhombohedral phases formation. The amplitudes normalized with respect to the primitive unit cell of the parent phase.

Amplitudes of global distortions (Å) AlV2O4 CuZr1.86(1)S4

Irrep

Dimensions

Order parameter direction

kwavevectors

Isotropy subgroups

k11(τ1)

1

(a)

(0 0 0)

Fd 3 m (227)

0.06090

2.25228

k11(τ7)

3

(a a a)

(0 0 0 )

R 3 m (166)

0.07323

3.87712

k9(τ4)

4

(a 0 0 0)

(1/2 1/2 1/2)

R 3 m (166)

0.37003

11.49193

R 3 m (166)

0.38209

12.33569

Total

In the case of AlV2O4, the typical picture for second order phase transitions is observed, i.e. the contribution from the critical order parameter to the displacement amplitude is dominant, and the atom displacements are relatively small – global distortion (total displacement of all

ACS Paragon Plus Environment

8

Page 9 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

atomic positions) is 0.38209 Å. The absolutely different situation is in CuZr1.86(1)S4. Indeed, there are giant atomic displacements (global distortion is 12.33559 Å) from parent phase with a large contribution from the improper order parameters k11(τ7) and k11(τ1) at the hypothetical/virtual Fd 3 m → R 3 m phase transition. Such giant displacements practically do not occur during phase transitions in inorganic compounds. Therefore, it is no wonder that the cubic phase has never been obtained for CuZr1.86(1)S4 (the rhombohedral phase of CuTi2S4 is detected in the compound prepared by using a KCl/KI flux at low temperatures, but an irreversible phase transition to the cubic phase occurs at above 450°C)39. One should emphasize that the choice of a cubic spinel Fd 3 m phase as a parent phase for describing the transition to a CuZr1.86(1)S4 rhombohedral structure is important for comparing the mechanisms of the low-symmetry phases formation in the two title compounds under study. Therefore, the absolute values of the theoretical displacements in the formation of the CuZr1.86(1)S4 rhombohedral phase have an exclusively conditional value. Thus, in the formation of rhombohedral R 3 m-modifications of spinels, two different mechanisms compete, namely with small (AlV2O4) and giant (CuZr1.86(1)S4) atomic displacements. These two mechanisms cause the existence of two types of rhombohedral phases with different metal clusters. 2.4. Differences in the AlV2O4 and CuZr1.86(1)S4 rhombohedral structures. Let us consider the changes in the pyrochlore and diamond spinel sublattices resulting in R 3 m -spinel modifications formation (Figure 2). A kagome sublattice in the rhombohedral structure, formed by B(3)-ions, has become trimerized or breathing kagome sublattice (Figure 2a,b). A twodimensional triangular sublattice, consisting of the apical octahedral B ions in cubic phase, split into two sublattices, formed by B(1)- and B(2)-ions (Figure 2c). The diamond spinel A-

ACS Paragon Plus Environment

9

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 32

sublattices are transformed into two type A(1)- and A(2)-sublattices in the rhombohedral structure (Figure 2d,e).

Figure 2. The transformations of pyrochlore (a, b, c) and diamond (d, e) sublattices of the Fd 3 m cubic spinel structure into sublattices of R 3 m -rhombohedral spinel modification. As a result of the phase transition B-sublattice of spinels (pyrochlore sublattice) splits into layers that are mutually interconnected by B(1) and B(2)-atoms (only in AlV2O4) (a). In this case [B(1)B(3)6]and [B(2)B(3)6]-heptamer structural blocks are formed too (only in AlV2O4). Kagome sublattices in low-symmetry R 3 m-phase have two types of [B(3)2]-dimers (dark blue and red) and two types of [B(3)3]-trimers (dark blue and red) with different B(3) – B(3) interatomic distances (red-

ACS Paragon Plus Environment

10

Page 11 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

long, blue - short) (b). Triangular sublattices, formed the apical В-atoms in cubic spinel structure, are transformed into two types of triangular sublattices in rhombohedral structure (c). One type of sublattice is formed by B(2)-atoms (blue), and another type is formed by B(1)-atoms (deep blue). Triangular diamond sublattice, formed by tetrahedral А-atoms in cubic spinel structure, splits in the R 3 m-rhombohedral structure into two types of triangular sublattices too, namely one type of triangular sublattice is formed by A(1)-atoms (dark green) and another type is formed by A(2)-atoms (yellow) (d, e). The disposition of the layers formed by A(1)-and A(2)-atoms is shown in (e).

In the kagome sublattice, there are two kinds of B(3) – B(3) bonds (B-dimers: [B(3)2]) with different bond lengths: the shorter (dark blue; 2.610 Å for AlV2O4 and 3.579 Å for CuZr1.86(1)S4) and the longer (red; 3.141 Å for AlV2O4 and 3.773 Å for CuZr1.86(1)S4) (Figure 2a,b). These dimers alternate in the kagome sublattice (Figure 2a,b). It is important to underline that the B-B distances in the cubic form of AlV2O4 and CuZr1.86(1)S4 are equal to 2.896 Å and 3.670 Å, respectively29,40. The kagome-layers, between which the B(2)-atoms are located, are shifted relative to each other in CuZr1.86(1)S4 (Figure 2a). As a result, the formation of heptamers (as structural blocks, but not as orbital molecules) linking kagome layers in CuZr1.86(1)S4 structure is observed only with B(1) atoms, while in the AlV2O4 structure the formation of heptamers with both B(1)- and B(2)-binding atoms is possible (Figure 2a). In addition, the triangle layers A(1)- and A(2)-atoms are shifted relative to each other in the CuZr1.86(1)S4 structure (Figure 2d,e). Table 2 shows the absolute atom displacements at the forming of the low-symmetry phases AlV2O4 and CuZr1.86(1)S4. The formation of the rhombohedral phase of CuZr1.86(1)S4 is

ACS Paragon Plus Environment

11

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 32

associated with much more noticeable changes than in AlV2O4 compared to the parent phase structure. The largest displacements occurs with the atoms A(2), B(3), X(2) and X(4). The amplitudes of their absolute displacements in thiospinels exceed 1Å. An important feature of the phase transition Fd 3 m → R 3 m in the above compounds is the absence of displacements of the B(1)3a and B(2)3b atoms, as well as the displacement of cations and anions in the Wyckoff position 6c along only the z axis (Supporting Information).

Table 2. Atomic displacements in the transition Fd 3 m → R 3 m. In bold type are the groups of atoms with the largest displacements

Atoms

Wyckoff position

Al(1)/Cu(1) Al(2)/Cu(2) V(1)/Zr(1) V(2)/Zr(2) V(3)/Zr(3) O(1)/S(1) O(2)/S(2) O(3)/S(3) O(4)/S(4)

6c 6c 3a 3b 18h 6c 6c 18h 18h

Displacements, (Å) AlV2O4 CuZr1.86(1)S4 0.13056 0.03941 0.09214 1.54078 0.00000 0.00000 0.00000 0.00000 0.15721 6.08456 0.03755 0.09047 0.14587 3.13621 0.08197 0.04341 0.04048 3.04844

Note the giant displacements of Zr(3) atoms from neighboring kagome layers towards each other. As a result, the kagome layers shift relative to each other, and [Zr(2)Zr(3)6] – structural blocks are not formed, while [V(2)V(3)6] in AlV2O4 are formed (Figure 3a). Cu(2) and Zr(3) atoms form a «bunch» of dimers (Figure 3a). Each «bunch» consists of three [Cu(2)– Zr(3)]-dimers which are joined by common Cu(2)-atom. Each Cu(2) atom is surrounded by three Zr(3)-atoms, resulting in three Cu(2)-Zr(3) bonds with interatomic distances of 2.937 Å. The shortest Cu-Zr distance in cubic spinel is much longer and equals 4.304 Å. This distance is too

ACS Paragon Plus Environment

12

Page 13 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

large for a significant interaction, while Cu(2)-Zr(3) interatomic distance in the rhombohedral form of CuZr1.86(1)S4 is the shortest one among all the metal-metal bond lengths. Note that Zr(1)and Zr(3)-atoms also form [Zr(1)Zr(3)]-dimers, but the interatomic distance 3.787 Å (on Figure 2a they are shown blue) is much longer. Thus, the «bunches» of [Cu(2)Zr(3)]-dimers are the smallest-sized clusters. It is interesting to note, that only contracted [Zr(3)3]-trimers (with shortest interatomic distances Zr(3)-Zr(3)) belong to the «bunch» of dimers. Thus, it has been found that CuZr1.86(1)S4 is the second compound (after CuTi2S438) in which such type of metal clusters exists.

ACS Paragon Plus Environment

13

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 32

Figure 3. Cations A(1) (green) and B(2) (blue) located between the kagome layers of cations B(3) (grey) (a). Anions are not shown. Kagome layers and the B(2) coordination polyhedra [B(2)X(3)6] along the [001] direction (b). The coordination polyhedra around B(2) and B(3) octahedra: O(1)/S(1) (red), O(2)/S(2) (purple), O(3)/S(3) (orange), O(4)/S(4) (dark purple) (c).

In the CuZr1.86(1)S4 rhombohedral thiospinel structure, Zr(2) atoms occupy a position being equidistant from 12 Zr(3)-atoms with a considerable distance between them that practically

ACS Paragon Plus Environment

14

Page 15 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

excludes the possibility of a chemical bonding between them (for comparison, in CuTi2S4 distance Ti(2)-Ti(3) = 4.488 Å39) (Figure 3a,b). As a result, the octahedron Zr(2)S(3)6 occupies a position in the voids (Figure 3a) formed as a result of the kagome layers shift, and the Zr(3) atoms form a peculiar cage-skeleton, causing the rattling effect. Such local environment of B(2) atoms in the rhombohedral CuZr1.86(1)S4 and CuTi2S4 structures is assumed to result in the high values of isotropic displacement factors observed39,40. In contrast, the V(2)-atoms in the AlV2O4 rhombohedral structure, as assumed in29, are to be a part of the [V(2)V(3)6]-heptamer structural blocks, and the [Al(2)O(3)6]-octahedron is more "clamped" between the two kagome layers (Figure 3b) that leads to the limitation in their thermal vibrations, and the values of their isotropic displacement factors do not differ significantly from the remaining atoms47. The existence of voids where [B(2)S(3)6]-octahedra are located should be taken into account when designing thermoelectric materials based on the rhombohedral thiospinels. Another important difference in the rhombohedral structures of AlV2O4 and CuZr1.86(1)S4 is the mutual arrangement of coordination polyhedrons formed around the atoms B(2) and B(3) (Figure 3c). In the AlV2O4 structure, they have common edges, and in the CuZr1.86(1)S4 structure they have common vertices. This significant transformation in the local environment of transition metal atoms determines electronic interactions, and as a consequence, the formation of charge, orbital and spin degrees of freedom55. As a result, the localized ground states can differ drastically that must be taken into consideration when choosing the materials for spintronics and orbitronics. Previously, the authors40 suggested that the structural formula of CuZr1.86(1)S4 rhombohedral modification that considers atom valence states (obtained from the crystalchemical analysis of the crystal bond lengths and electronic structure calculations) should be

ACS Paragon Plus Environment

15

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 32

(Cu+)4Zr(1)4+Zr(2)3+Zr(3)3.5+6(S2-)16. This formula means that the Zr(1)4+ ion has no t2g orbitals, i.e. there are no chemical bonding between Zr(1)- and Zr(3)-atoms. Therefore, the formation of Zr-heptamers is improbable. Besides, the calculated densities of states of rhombohedral modification show that Cu(1)–Zr(3) and Zr(3)-Zr(3) interactions in «bunch» of dimers and [Zr(3)3]-trimers (with interatomic distances Zr(3)-Zr(3) 3.579 Å40) have bonding character40. We consider the formation of clusters only on the basis of the analysis of the average structure and interatomic distances obtained from the experimental diffraction data. Therefore, in this case heptamers in AlV2O4 are only the structural units, i.e. building blocks. We believe that in AlV2O4 only a fraction of the [V(2)V(3)6]-heptamers, namelly [V(3)3]-trimers with the interatomic distances V(3) – V(3) equal to 2.610 Å can be considered as orbital metal molecules. In this orbital molecules, the lengths of the V(3) – V(3) bonds are even lower than the V-V bond length in the metal vanadium (2.62 Å). According to J.B. Goodenough’s56-58, there is a critical cation-cation distance Rc in oxide, dividing the regions of localized and itinerant-electron behavior. Note that Rc for interatomic distances V–V is 2.62 Å. Therefore, we believe that the nature of the V(3)-V(3) bond in the [V(3)3]-trimers is close to the metal bond. Recently, Browne et al30,59, based on an analysis of the pair distribution functions of the AlV2O4 and its analogue GaV2O4 structures, showed that some local distortions in the above structures do not agree with that the heptamers are orbital molecules29. Thus, in the of rhombohedral spinels structure, two mechanisms of formation of metal clusters are observed. The above mechanisms are explained by various displacement amplitudes for a real or hypothetical phase transitions from a cubic structure. With small distortions in the structure of the cubic spinel, [B(2)B(3)6]-heptamer structural blocks are formed, giving rise to [V(3)3]-trimers60 or (trimers + tetramers) local molecules30 as in AlV2O4. At large amplitudes of

ACS Paragon Plus Environment

16

Page 17 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

distortions, giant displacements of B(3) atoms occur and no [B(2)B(3)6]-heptamer structural blocks are formed, whereas a new type of clusters – "bunches" of dimers [Cu(1)Zr(3)3] is formed. As mentioned above, all known clusters in the structures of the low-symmetry phase of spinels are formed by atoms only one sublattice, while «bunches» of dimers are formed by atoms of two frustrating sublattices (kagome and triangular sublattices). Therefore, the substances, in the structures of which similar clusters exist, are expected to have unique physical properties, above all magnetic and electrical ones. In particular, one should note that the interaction between the A- and B-cations will lead to an unusual spiral short-range magnetic order. This magnetic order can induce ferroelectricity as proposed in CoCr2O461.

3. CONCLUSIONS The main result of the work is to establish the general reasons (group-theoretical and structural) for the significant differences in AlV2O4 and CuZr1.86(1)S4 structures. It is proved that different values of the critical and improper order parameters explain the formation of various structures. In the case of the AlV2O4 rhombohedral modification, the atomic displacements are insignificant. The role of improper order parameters in the phase formation is also insignificant. In this case, Landau thermodynamics, which operates with small in magnitude order parameters, is quite suitable tool for describing of the structural transition and the rhombohedral phase formation mechanism. The situation changes completely in the case of the CuZr1.86(1)S4 rhombohedral modification. If a phase transition from a cubic spinel to a rhombohedrally distorted spinel occurred in this substances, then the displacements of the atoms, according to our calculations, would have to be giant. At the same time, the displacements generated by improper

ACS Paragon Plus Environment

17

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 32

(secondary) parameters of the order of k11(τ7) and k11(τ1) should contribute significantly to the calculated value of the displacements of Cu(2)-, S(2)- and S(4)-atoms. The various values of the critical and the improper order parameters describing the formation of rhombohedral structures are caused by the crystal-chemical features of the compositions and possibly by the methods of sample preparation. This is an interesting topic for further research. This study first proposed a structural mechanism of the low-symmetry rhombohedral modifications formation from high-symmetry cubic spinel phase, and showed that the calculated low-symmetry rhombohedral phase is due to displacements and the ordering of all types of atoms. An important result of the work is that a «bunch» of [Cu(1)Zr(3)]-dimers, two types of Zr(3)-trimers [Zr(3)3] exist in CuZr1.86(1)S4. The discovery of metal nanostructures, such as [Cu(1)Zr(3)]-dimers, trimers [Zr(3)3] with the shorter bond lengths and a «bunch» of [Cu(1)Zr(3)]-dimers is of particular interest. The metalic dimers and trimers in geometrically frustrated structures (for example, dimers: MgTi2O419,20,62, CuTe2O563,64, VO265, CuGeO366; trimers: LiVO225, LiVS267, NaV6O1068, BaV10O1569-72, SrV8Ga4O1973, AV13O18 (A = Ba, Sr)74, A2V13O22 (A = Ba, Sr)75, Ba4Ru3O1076) have long been known. The «bunch» of [Cu(1)-Zr(3)]dimers in CuZr1.86(1)S4 is a new type of self-organization of atoms in geometrically frustrated spinel-like structures. It is vital to stress that metal clusters in all known spinel materials are formed by the atoms with only one geometrically frustrating sublattice. CuZr1.86(1)S4

and CuTi2S4

rhombohedral phases are the only examples of the crystals in which the metal clusters are formed in two geometrically frustrating spinel .sublattices. The appearance of two sublattices means new

ACS Paragon Plus Environment

18

Page 19 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

combinatorial possibilities for the selection of cations in these sublattices. Therefore, multiple scenarios of interaction between two frustrating sublattices are possible in crystals having similar clusters. This suggests that CuZr1.86(1)S4 and CuTi2S4 are prototypes of a new class of materials with possible exotic phases as well as unexpected intriguing physical properties. Also, our study showed for the first time that the values of critical and improper (noncritical) order parameters will be considered as a new source of structural diversity of materials.

ASSOCIATED CONTENT Supporting Information Atomic representations of AlV2O4 and CuZr1.86(1)S4 structures. Calculation details (ISODISTORT: modes details for AlV2O4 and CuZr1.86(1)S4).

AUTHOR INFORMATION Corresponding Author *E-mail: [email protected] ORCID Mikhail V. Talanov: 0000-0002-5416-9579 Present Address ¶ Southern Federal University, Research institute of physics, 194 Stachki av., Rostov-on-Don, 344090. Notes The author declares no competing financial interests.

ACKNOWLEDGMENTS

ACS Paragon Plus Environment

19

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 32

The reported study was funded by RFBR, according to the research project No. 16-3260025 mol_а_dk. REFERENCES (1) Ramirez, A. P. Strongly geometrically frustrated magnets. Annu. Rev. Mater. Sci. 1994, 24, 453-480 (and references therein). (2) Krimmel, A.; Mutka, H.; Koza, M. M.; Tsurkan, V.; Loidl, A. Spin excitations in frustrated -site spinels investigated with inelastic neutron scattering. Phys. Rev. B 2009, 79, 134406. (3) Bergman, D.; Alicea, J.; Gull, E.; Trebst, S.; Balents, L. Order-by-disorder and spiral spinliquid in frustrated diamond-lattice antiferromagnets. Nature Physics 2007, 3, 487-491. (4) Anderson, P. W. Resonating valence bonds: A new kind of insulator? Mater. Res. Bull. 1973, 8, 153-160. (5) Fazekas, P.; Anderson, P. W. On the ground state properties of the anisotropic triangular antiferromagnet. Philos. Mag. 1974, 30, 423-440. (6) Fritsch, V. J.; Hemberger, J.; Büttgen, N.; Scheidt, E.-W.; Krug von Nidda, H.-A.; Loidl, A.; Tsurkan, V. Spin and orbital frustration in MnSc2S4 and FeSc2S4. Phys. Rev. Lett. 2004, 92, 116401. (7) Fichtl, R. V.; Tsurkan, V.; Lunkenheimer, P.; Hemberger, J.; Fritsch, V.; Krug von Nidda, H.-A.; Scheidt, E.-W.; Loidl, A. Orbital freezing and orbital glass state in FeCr2S4. Phys. Rev. Lett. 2005, 94, 027601. (8) Kondo, S.; Johnston, D. C.; Swenson, C. A.; Borsa, F.; Mahajan, A. V.; Miller, L. L.; Gu, T.; Goldman, A. I.; Maple, M. B.; Gajewski, D. A.; Freeman, E. J.; Dilley, N. R.; Dickey, R. P.;

ACS Paragon Plus Environment

20

Page 21 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Merrin, J.; Kojima, K.; Luke, G. M.; Uemura, Y. J.; Chmaissem, O.; Jorgensen, J. D. LiV2O4: A heavy fermion transition metal oxide. Phys. Rev. Lett. 1997, 78, 3729-3732. (9) Fujimoto, S. Hubbard chains network on corner-sharing tetrahedra: Origin of the heavyfermion state in LiV2O4. Phys. Rev. B 2002, 65, 155108. (10) Shindou, R.; Nagaosa, N. Orbital ferromagnetism and anomalous hall effect in antiferromagnets on the distorted fcc lattice. Phys. Rev. Lett. 2001, 87, 116801. (11) Furubayashi, T.; Suzuki, H.; Kobayashi, N.; Nagata, S. Large negative magnetoresistance in thiospinel CuCrZrS4. Solid State Commun. 2004, 131, 505-508. (12) Johnston, D. C., Prakash, H., Zachariasen, W. H., Viswanathan, R. Mater. Res. Bull. 1973, 8, 777-784. (13) Robbins, M., Willens, R. H., Miller, R. C. Solid State Commun. 1967, 5, 933-934. (14) Sadoc, J.-F.; Mosseri, R. Geometrical Frustration; Cambridge University Press: Cambridge, 2006. (15) Takagi, H.; Niitaka, S. Highly Frustrated Magnetism in Spinels. In: Lacroix C., Mendels P., Mila F. (eds) Introduction to Frustrated Magnetism. Springer Series in Solid-State Sciences, 164. Springer: Berlin, 2011. (16) Anderson, P. W. Ordering and antiferromagnetism in ferrites. Phys. Rev. 1956, 102, 10081013. (17) Büttgen, N.; Hemberger, J.; Fritsch, V.; Krimmel, A.; Mücksch, M.; Krug von Nidda, H.A.; Lunkenheimer, P.; Fichtl, R.; Tsurkan, V.; Loidl, A. Orbital physics in sulfur spinels: ordered, liquid and glassy ground states. New Journal of Physics 2004, 6, 191. (18) Attfield, J. P. Orbital molecules in electronic materials. APL Materials. 2015, 3, 041510.

ACS Paragon Plus Environment

21

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 32

(19) Schmidt, M.; Ratcliff, W.; Radaelli, P. G.; Refson, K.; Harrison, N. M.; Cheong, S.W. Spin singlet formation in MgTi2O4: evidence of a helical dimerization pattern. Phys. Rev. Lett. 2004, 92, 056402. (20) Talanov, V.M.; Shirokov, V.B.; Ivanov, V.V.; Talanov, M.V. Theory of structural phase transition in MgTi2O4. Crystallogr. Rep. 2015, 60, 101-110. (21) Okamoto, Y.; Nohara, M.; Aruga-Katori, H.; Takagi, H. Spin-liquid state in the hyperkagome antiferromagnet Na4Ir3O8. Phys. Rev. Lett. 2007, 99, 137207. (22) Talanov, V. M.; Shirokov, V. B.; Talanov, M. V. Unique atom hyper-kagome order in Na4Ir3O8 and in low-symmetry spinel modifications. Acta Crystallogr. A. 2015, 71, 301-318. (23) Takayama, T.; Yaresko, A.; Matsumoto, A.; Nuss, J.; Ishii, K.; Yoshida, M.; Mizuki, J.; Takagi, H. Spin-orbit coupling induced semi-metallic state in the 1/3 hole-doped hyper-kagome Na3Ir3O8. Sci. Rep. 2014, 4, 6818. (24) Talanov, M. V.; Talanov, V. M.; Shirokov, V. B. Unique hyper-kagome atomic order in geometrically frustrated iridium spinel-like structures. Russ. Chem. Bull. 2017, 66 (10), 17191727. (25) Goodenough, J.B.; Dutta, G.; Manthiram, A. Lattice instabilities near the critical VV separation for localized versus itinerant electrons in LiV1-yMyO2 (M=Cr or Ti) Li1-xVO2. Phys. Rev. B 1991, 43, 10170-10178. (26) Senn, M. S.; Wright, J. P.; Attfield, J. P. Charge order and three-site distortions in the Verwey structure of magnetite. Nature. 2012, 481, 173-176. (27) Senn, M. S.; Loa, I.; Wright, J. P.; Attfield, J. P. Electronic orders in the Verwey structure of magnetite. Phys. Rev. B. 2012, 85, 125119.

ACS Paragon Plus Environment

22

Page 23 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

(28) Senn, M. S.; Wright, J. P.; Cumby J.; Attfield, J. P. Charge localization in the Verwey structure of magnetite. Phys. Rev. B. 2015, 92, 024104. (29) Horibe, Y.; Shingu, M.; Kurushima, K.; Ishibashi, H.; Ikeda, N.; Kato, K.; Motome, Y.; Furukawa, N.; Mori, S.; Katsufuji, T. Spontaneous formation of vanadium “Molecules” in a geometrically frustrated crystal: AlV2O4. Phys. Rev. Lett. 2006, 96, 086406. (30) Browne, A.J.; Kimber, S.A.J.; Attfield, J.P. Persistent three- and four-atom orbital molecules in the spinel AlV2O4. Phys. Rev. Mater. 2017, 1, 052003(R). (31) Radaelli, P. G.; Horibe, Y.; Gutmann, M. J.; Ishibashi, H.; Chen, C. H.; Ibberson, R. M.; Koyama, Y.; Hor, Y.-S.; Kiryukhin, V.; Cheong, S.-W. Formation of isomorphic Ir3+ and Ir4+ octamers and spin dimerization in the spinel CuIr2S4. Nature 2002, 416, 155-158. (32) Talanov, V.M. Structural mechanism of the ordering of ions in tetrahedral holes in spinels. J. Struct. Chem. 1986, 27(2), 327–330. (33) Plumier, R.; Sougi, M.; Lecomte, M. Observation of an unusual short range magnetic ordering in spinel Cu12In12Cr2S4. Phys. Lett. A 1977, 60, 341–344. (34) Plumier, R.; Lecomte, M.; Miedan Gros, A.; Sougi, M. Observation of a first order macro to microdomain transition in chalcogenide spinel Ag12In12Cr2S4. Physica B 1977, 86–88, 1360– 1362. (35) Nauciel-Bloch, M.; Plumier, R. Magnetic structures in spinel b-site lattices with ordered diamagnetic ions at the A-sites. Solid State Commun. 1971, 9, 223–226. (36) Okamoto, Y.; Nilsen, G.J.; Attfield, J.P.; Hiroi, Z. Breathing pyrochlore lattice realized in A-site ordered spinel oxides LiGaCr4O8 and LiInCr4O8. Phys. Rev. Lett. 2013, 110, 097203.

ACS Paragon Plus Environment

23

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 32

(37) Talanov, M.V.; Talanov, V.M. Universal mechanism of tetrahedral metal cluster formation in structures with breathing pyrochlore sublattices. Nanosystems: Phys. Chem. Math. 2017, 8(5), 677–687. (38) Talanov, M.V.; Shirokov, V.B.; Talanov, V.M. Phenomenological thermodynamics and the structure formation mechanism of the CuTi2S4 rhombohedral phase. Phys. Chem. Chem. Phys. 2016, 18, 10600-6. (39) Soheilnia, N; Kleinke, K.M.; Dashjav, E.; Cuthbert, H.L.; Greedan, J.E.; Kleinke, H. Crystal structure and physical properties of a new CuTi2S4 modification in comparison to the thiospinel. Inorg. Chem. 2004, 43, 6473-6478. (40) Dong, Y.; McGuir, M. A.; Yun, H.; DiSalvo, F.J. Synthesis, crystal structure, and properties of the rhombohedral modification of the thiospinel CuZr1.86(1)S4. J. Solid State Chem. 2010, 183, 606-612. (41) Howard, C.J.; Stokes, H.T. Structures and phase transitions in perovskites–a grouptheoretical approach. Acta Crystallogr. A. 2005, 61, 93-111. (42) Campbell, B.J.; Stokes, H.T.; Tanner, D.E.; Hatch, D.M. ISODISPLACE: a web-based tool for exploring structural distortions. J. Appl. Cryst. 2006, 39, 607-614. (43) Orobengoa, D.; Capillas, C.; Aroyo, M.I.; Perez-Mato, J.M. AMPLIMODES: symmetrymode analysis on the Bilbao Crystallographic Server. J. Appl. Cryst. 2009, 42, 820-833. (44) Perez-Mato, J.M.; Orobengoa, D.; Aroyo, M.I. Mode crystallography of distorted structures. Acta Crystallogr. A. 2010, 66 , 558-590. (45) Kalavathi, S.; Amirthapandian, S.; Chandra, S.; Sahu, P.Ch.; Sahu, H.K. Valence state, hybridization and electronic band structure in the charge ordered AlV2O4. J. Phys.: Condens. Matter, 2014, 26, 015601.

ACS Paragon Plus Environment

24

Page 25 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

(46) Kalavathi, S.; Raju, S.V.; Chandra, S.; Williams, Q.; Sahu, P.Ch. Hydrostatic pressure dependence on the collapsing of heptamer clusters in the charge ordered spinel AlV2O4. Solid State Commun. 2017, 250, 23-27. (47) Kalavathi, S.; Raju, S.V.; Williams, Q.; Sahu, P.Ch.; Sastry, V.S.; Sahu, H.K. Pressureinduced frustration in charge ordered spinel AlV2O4. J. Phys.: Condens. Matter, 2013, 25, 292201. (48) Matsuno, K.; Katsufuji, T.; Mori, S.; Nohara, M.; Machida, A.; Moritomo, Y.; Kato, K.; Nishibori, E.; Takata, M.; Sakata, M.; Kitazawa, K.; Takagi H. Charge ordering and spin frustration in AlV2-xCrxO4. Phys. Rev. Lett. 2003, 90, 096404. (49) Matsuno, K.; Katsufuji, T.; Mori S.; Moritomo, Y.; Machida, A.; Nishibori, E.; Takata, M.; Sakata, M.; Yamamoto, N.; Takagi, H. Charge ordering in the geometrically frustrated spinel AlV2O4. J. Phys. Soc. Jpn. 2001, 70, 1456-1459. (50) Landau, L.D.; Lifshitz, E.M. Statistical Physics, Part 1. Pergamon Press: Oxford, 1980. (51) Kovalev, O. V. Representations of the Crystallographic Space Groups. Irreducible Representations, Induced Representations and Co-representations, Gordon and Breach: Amsterdam, 1993. (52) Miller, S. C.; Love, W. F. Tables of Irreducible Representations of Space Groups and Corepresentations of Magnetic Space Groups. Pruett: Boulder, 1967. (53) For k9(L): τ1 – L1+, τ2 – L1–, τ3 – L2+, τ4 – L2–, τ5 – L3+, τ6 – L3–; for k10(X): τ1 – X3, τ2 – X4, τ3 – X1, τ4 – X2; for k11(Γ): k11(τ1) (A1g) – GM1+, k11(τ2) (A1u) – GM1–, k11(τ3) (A2g) – GM2+, k11(τ4) (A2u) – GM2–, k11(τ5) (Eg) – GM3+, k11(τ6) (Eu) – GM3–, k11(τ7) (F2g) – GM5+, k11(τ8) (F2u) – GM5–, k11(τ9) (F1g) – GM4+, k11(τ10) (F1u) – GM4–.

ACS Paragon Plus Environment

25

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 32

(54) Sakhnenko, V.P.; Talanov, V.M.; Chechin, G.M. Group-theoretic analysis of the complete order parameters condensate that occurs at structural phase transitions . Phys. Met. Metallogr. 1986, 62(5), 10–18. (55) Streltsov, S.V.; Khomskii, D.I. Orbital physics in transition metal compounds: new trends. Physics–Uspekhi 2017,60(11), 1121-1146. (56) Rogers, D. B.; Arnott, R. J.; Wold, A.; Goodenough, J. B. J. The preparation and properties of some vanadium spinels. Phys. Chem. Solids. 1963, 24, 347-360. (57) Goodenough, J. B. Magnetism and the Chemical Bond. John Willey & Sons: New York, 1963. (58) Goodenough, J. B. Anomalous properties of the vanadium oxides. Annu. Rev. Mater. Sci. 1971, 1, 101-138. (59) Browne, A.J.; Lithgow, C.; Kimber, S.A.J.; Attfield, J.P. Orbital molecules in the new spinel GaV2O4. Inorg. Chem. 2018, 57, 2815-2822. (60) Talanov, M.V.;

Shirokov, V.B.; Avakyan, L.A.; Talanov, V.M.; Borlakov, Kh.Sh.

Vanadium clusters formation in geometrically frustrated spinel oxide AlV2O4. Acta Crystallogr. B. (under consideration). (61) Yamasaki, Y.; Miyasaka, S.; Kaneko, Y.; He, J.-P.; Arima, T.; Tokura, Y. Magnetic reversal of the ferroelectric polarization in a multiferroic spinel oxide. Phys. Rev. Lett. 2006, 96, 207204. (62) Ivanov, V.V.; Talanov, V.M.; Shirokov V.B.; Talanov, M.V. Crystal chemistry and formation mechanism of tetragonal MgTi2O4. Inorg. Mater. 2011. 47(9), 990-998. (63) Hanke, K.; Kupcik, V.; Lindqvist, O. The crystal structure of CuTe2O5. Acta Crystallogr.B 1973, 29, 963-970.

ACS Paragon Plus Environment

26

Page 27 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

(64) Ushakov, A.V.; Streltsov, S.V. Electronic and magnetic structure for the spin-gapped system CuTe2O5. J. Phys.: Condens. Matter. 2009, 21, 305501. (65) Imada, M.; Fujimori, A.; Tokura, Y. Metal-insulator transitions. Rev. Mod. Phys. 1998, 70, 1039-1263. (66) Hase, M.; Terasaki, I.; Uchinokura, K. Observation of the spin-Peierls transition in linear Cu2+ (spin-1/2) chains in an inorganic compound CuGeO3. Phys. Rev. Lett. 1993, 70, 3651. (67) Katayama, N.; Uchida, M.; Hashizume, D.; Niitaka, S.; Matsuno, J.; Matsumura, D.; Nishihata, Y.; Mizuki, J.; Takeshita, N.; Gauzzi, A.; Nohara, M.; Takagi, H. Anomalous metallic state in the vicinity of metal to valence-bond solid insulator transition in LiVS2. Phys. Rev. Lett. 2009, 103, 146405. (68) Kato, H.; Kato, M.; Yoshimura, K.; Kosuge, K.

23

Na-NMR study in NaV6O11. J. Phys.

Soc. Jpn. 2001, 70, 1404-1410. (69) Liu, G.; Greedan, J.E. Structure and magnetism in AB10O15+x (A= Sr, Ba; B= V, Cr; x= 0,∼ 1): evidence for geometric frustration. J. Solid State Chem. 1996, 122, 416-427. (70) Bridges, C.A.; Greedan, J.E.; Kleinke, H. Insights on the origin of the structural phase transition in BaV10O15 from electronic structure calculations and the effect of Ti-doping on its structure and electrical transport properties. J. Solid State Chem. 2004, 177, 4516-4527. (71) Bridges, A.; Greedan, J.E. Phase transitions and electrical transport in the mixed-valence V2+/V3+ oxide BaV10O15. J. Solid State Chem. 2004, 177, 1098-1110. (72) Kajita, T.; Kanzaki, T.; Suzuki, T.; Kim, J.E.; Kato, K.; Takata, M.; Katsufuji, T. Opening of a charge gap with V trimerization in BaV10O15. Phys. Rev. B., Condens. Mater. Phys. 2010, 81, 060405(R).

ACS Paragon Plus Environment

27

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 28 of 32

(73) Miyazaki, J.; Sonehara, T.; Akahoshi, D.; Kuwahara, H.; Kim, J.E.; Kato, K.; Takata, M.; Katsufuji, T. Impact of orbital degrees of freedom on geometrical frustration in the kagome-like magnet SrVxGa12-xO19. Phys. Rev. B., Condens. Mater. Phys. 2009, 79, 180410(R). (74) Ikeda, M.; Okuda, T.; Kato, K.; Takata, M.; Katsufuji, T. Competition between vanadium tetramerization and trimerization in Ba1-xSrxV13O18. Phys. Rev. B., Condens. Mater. Phys. 2011, 83, 134417. (75) Miyazaki, J.; Matsudaira, K.; Shimizu, Y.; Itoh, M.; Nagamine, Y.; Mori, S.; Kim, J.E.; Kato, K.; Takata, M.; Katsufuji, T. Formation of a three-dimensional network of V trimers in A2V13O22 (A = Ba, Sr). Phys. Rev. Lett. 2010, 104, 207201. (76) Streltsov, S.V.; Khomskii, D.I. Unconventional magnetism as a consequence of the charge disproportionation and the molecular orbital formation in Ba4Ru3O10. Phys. Rev. B., Condens. Mater. Phys. 2012, 86, 064429.

ACS Paragon Plus Environment

28

Page 29 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

For Table of Contents Use Only

Two different mechanisms of metal cluster formation in the rhombohedral spinel structures: AlV2O4 and CuZr1.86(1)S4

Mikhail V. Talanov

SYNOPSIS: The difference in structures of the AlV2O4 and CuZr1.86(1)S4 isosymmetric rhombohedral spinels is due to a critical influence of the order parameter value that describes the real (in the case of AlV2O4) and virtual (in the case of CuZr1.86(1)S4) structural phase transitions which are associated with the formation of various metal clusters: heptamer structural blocks and ‘‘bunch’’ of dimers, respectively.

ACS Paragon Plus Environment

29

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 1. The cubic spinel structure fragments. A- and B- cation tetrahedra in the spinel structure (a). A atoms and B atoms are shown by green and deep blue color respectively; X-anions are not shown. In the spinel structure, B-sublattices are formed by a network of corner-sharing B-tetrahedra (deep blue color). The B-sublattice can be visualized as alternating kagome (b) and triangular (b, c) layers of structure along a [111] direction as designated in (b). Kagome layers are shown by deep blue, triangular layers are shown by dark yellow color. 170x46mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Figure 2. The transformations of pyrochlore (a, b, c) and diamond (d, e) sublattices of the Fd-3m cubic spinel structure into sublattices of R-3m -rhombohedral spinel modification. As a result of the phase transition B-sublattice of spinels (pyrochlore sublattice) splits into layers that are mutually interconnected by B(1) and B(2)-atoms (only in AlV2O4) (a). In this case [B(1)B(3)6]-and [B(2)B(3)6]-heptamer structural blocks are formed too (only in AlV2O4). Kagome sublattices in low-symmetry R-3m-phase have two types of [B(3)2]-dimers (dark blue and red) and two types of [B(3)3]-trimers (dark blue and red) with different B(3) – B(3) interatomic distances (red-long, blue - short) (b). Triangular sublattices, formed the apical В-atoms in cubic spinel structure, are transformed into two types of triangular sublattices in rhombohedral structure (c). One type of sublattice is formed by B(2)-atoms (blue), and another type is formed by B(1)-atoms (deep blue). Triangular diamond sublattice, formed by tetrahedral А-atoms in cubic spinel structure, splits in the R3m-rhombohedral structure into two types of triangular sublattices too, namely one type of triangular sublattice is formed by A(1)-atoms (dark green) and another type is formed by A(2)-atoms (yellow) (d, e). The disposition of the layers formed by A(1)-and A(2)-atoms is shown in (e). 176x129mm (300 x 300 DPI)

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 3. Cations A(1) (green) and B(2) (blue) located between the kagome layers of cations B(3) (grey) (a). Anions are not shown. Kagome layers and the B(2) coordination polyhedra [B(2)X(3)6] along the [001] direction (b). The coordination polyhedra around B(2) and B(3) octahedra: O(1)/S(1) (red), O(2)/S(2) (purple), O(3)/S(3) (orange), O(4)/S(4) (dark purple) (c). 147x170mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 32 of 32