{'fraction_non_alphanumeric': 0.0303881368308712, 'fraction_numerical': 0.021653653878559793, 'mean_word_length': 4.722757955641273, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Antiferromagnetism in the van der Waals layered spin-lozenge semiconductor CrTe3

The crystallographic, magnetic, and transport properties of the van der Waals bonded, layered compound CrTe3 have been investigated on single crystal and polycrystalline materials. The crystal structure contains layers made up of lozenge shaped Cr_4 tetramers. Electrical resistivity measurements show the crystals to be semiconducting, with a temperature dependence consistent with a band gap of 0.3eV. The magnetic susceptibility exhibits a broad maximum near 300K characteristic of low dimensional magnetic systems. Weak anomalies are observed in the susceptibility and heat capacity near 55K, and single crystal neutron diffraction reveals the onset of long range antiferromagnetic order at this temperature. Strongly dispersive spin-waves are observed in the ordered state. Significant magneto-elastic coupling is indicated by the anomalous temperature dependence of the lattice parameters and is apparent in structural optimization in van der Waals density functional theory calculations for different magnetic configurations. The cleavability of the compound is apparent from its handling and is confirmed by first principles calculations, which predict a cleavage energy 0.5J/m^2, similar to graphite. Based on these results, CrTe3 is identified as a promising compound for studies of low dimensional magnetism in bulk crystals as well as magnetic order in monolayer materials and van der Waals heterostructures.


I. INTRODUCTION
Electronic peculiarities and proposed applications of graphene [1,2] ignited intense research into two-dimensional materials composed of layers held together by van der Waals bonding. Efforts to extend this research to other materials led, for example, to the discovery of interesting and potentially useful optoelectronic behaviors in ultra-thin transition metal dichalcogenides [3][4][5][6], providing further demonstrations of new physics emerging as bulk materials are reduced to monolayers. Interest has also risen in exploring van der Waals heterostructures [7], which requires a library of layered materials that can be mechanically exfoliated and recombined in customized stacks. The incorporation of materials with complementary functionalities, including magnetism, holds promise for the development of new physics and devices that can be controlled by external fields. While the set of twodimensional materials beyond graphene and transition metal dichalcogenides continues to expand [8][9][10], there are still relatively few studies of easily cleavable magnetic compounds and interest in such materials is growing.
The compound CrTe 3 was first reported by Klepp and Ipser [30,31], and received little subsequent attention. It adopts the unique layered crystal structure shown in Fig. 1a that contains tetramers of Cr atoms arranged in a rhombus or lozenge shape. This structure is described in detail in the literature [31][32][33] and in the discussion below. Canadell et al. [33] identified CrTe 3 as a Mott-Hubbard insulator based on electronic structure calculations and the observation of semiconducting behavior in electrical resistance measurements, but were unable to characterize the magnetic properties due to the presence of strongly magnetic impurities in powder samples.
We have grown single crystals of CrTe 3 and have studied the structural, thermal, magnetic, and transport properties of the material using x-ray and neutron diffraction, physical property measurements, inelastic neutron scattering, and first principles calculations. We find very anisotropic and unusual thermal expansion below room temperature, with evidence of strong magneto-elastic coupling as in-plane magnetic correlations develop. The compound is easily cleavable, which is clear from handling the crystals and from density functional theory calculations incorporating magnetism and the van der Waals interactions. These calculations predict a band gap of 0.26 eV for a magnetically ordered state, close to the values determined from electrical resistivity measurements near room temperature. The magnetic susceptibility follows qualitatively the behavior expected for a two-dimensional antiferromagnet. Alternatively, the data can be described by adapting a spin-lozenge model previously developed for other materials containing spin-3/2 tetramers like those found in CrTe 3 . Low dimensional magnetic correlations develop above room temperature, and magnetic order is observed in neutron diffrac-tion measurements below T N = 55 K. Moments are arranged in an up-down-up-down pattern around the perimeters of each lozenge-shaped tetramers, and the long-range magnetic order consists of stripe-like antiferromagnetic arrangement of the tetramers within each layer, with layers stacked antiferromagnetically. Weak anomalies are observed in the heat capacity and magnetic susceptibility at T N . Inelastic neutron scattering reveals dispersive magnetic excitations at low temperature, supporting the presence of non-negligible inter-tetramer interactions. Together the experimental and theoretical results presented here suggest CrTe 3 to be (1) an interesting material in which to study low dimensional magnetism in bulk crystals, (2) an easily cleavable material for investigating magnetic monolayers or few layer systems, and (3) a candidate semiconducting material for incorporating magnetism into van der Waals heterostructures.

II. PROCEDURES
CrTe 3 powder was produced by reacting Cr powder and Te shot in evacuated and sealed silica ampoules. The elements were mixed in a 1:3 molar ratio and heated to 700 • C and held for several hours, then cooled quickly to 425 • C and annealed for about 24 hours. The resulting sample was single phase CrTe 3 by powder x-ray diffraction, but some ferromagnetic Cr 5 Te 8 was detected by magnetization measurements. This impurity could be reduced in concentration by annealing the powder at similar temperatures in the presence of Te vapor produced by excess Te shot in the same ampoule but physically separated from the sample, or by starting with a Te-rich composition like CrTe 3.1 .
Single crystals were grown by two methods. The first employed a halide flux consisting of a mixture of KCl and AlCl 3 in a 45:55 molar ratio (m.p. 220 • C). Working in a helium filled glovebox, about 100 mg of CrTe 3 powder produced as described above was added to the mixture of dry salts in a silica tube (9 mm I.D., 1.5 mm wall thickness) that was sealed under vacuum without exposure to air. The tube was heated in a horizontal furnace so that one end was held at 450 • C and the other at 425 • C. After several days to one week the furnace was turned off and the tube cooled to room temperature. Platelike crystals were recovered by washing the flux away with water. This method produced mm-sized thin crystals like those show in Fig. 1b. The second crystal growth method used excess Te to grow CrTe 3 directly from the melt. Mixtures of Cr and Te containing 1-3 atomic % Cr were placed in an alumina crucible covered with an alumina frit [34], and heated to 1050 • C and held for 16 h, cooled to 600 • C and held for 24 h, then cooled at a rate of 1.5 • C per hour to 455 • C. At this temperature the excess Te was removed by centrifugation. Crystal produced by this method are shown in Fig. 1c and were typically several mm on a side and often grew in thick, blocklike forms with thicknesses up to several mm, but are easily delaminated into thin sheets. Increasing the Cr concentration in the melt to 4% produced single crystals of Cr 5 Te 8 instead of CrTe 3 , consistent with the published phase diagram [35].
Powder and single crystal x-ray diffraction were performed using a PANalytical X'Pert Pro MPD powder diffractometer (Cu-K α1 radiation) with an Oxford PheniX cryostat and a Bruker SMART APEX CCD single crystal diffractometer (Mo K α radiation). Powder data was analyzed using FullProf [36]. Single crystal data were analyzed using SADABS for absorption corrections, XPREP for symmetry analysis, and the ShelX suite for structural solution and refinement. Mag-netization measurements were conducted using a Magnetic Property Measurement System (Quantum Design). Heat capacity and electrical resistivity measurements were performed with a Physical Property Measurement System (Quantum Design). Electrical contacts were made using silver paste and platinum wires. Heat capacity was measured on single crystal samples, with measurements down to 380 mK carried out using the Helium-3 Option. Neutron powder diffraction measurements were performed using the POWGEN instrument at the Spallation Neutron Source (SNS) at Oak Ridge National Laboratory. For these measurements, the powder sample was loaded into 8 mm diameter vanadium can, and cooled down to 10 K using a closed-cycle-refrigerator. Two frames with center neutronwavelengths of 1.333Å and 2.665Å were used for collecting data over a sufficiently large d-spacing range (0.5Å to 10Å) at two different temperatures: 10 K and 300 K. A temperature dependence study was carried out by collecting data continuously using the 2.665Å frame while cooling the sample between 300 K and 10 K at a nominal rate of 0.5 K/min. The data were rebinned into datasets every 10 K with a nominal counting time of 20 minutes. The structural data was refined using FullProf [36]. Neutron diffraction study of a single crystal specimen (identical to that used for the magnetization study) was performed using the HB3A 4-circle diffractometer at the High Flux Isotope Reactor. Data were collected using the neutron wavelength of 1.546Å. Inelastic neutron scattering (INS) measurements were performed using the hybrid spectrometer HYSPEC at the SNS. HYSPEC is a highly versatile direct geometry spectrometer that combines time-of-flight spectroscopy with the focusing Bragg optics [37]. The incident neutron beam is monochromated using a Fermi chopper and is then vertically focused by Bragg scattering onto the sample position by highly oriented pyrolytic graphite. For the INS measurements, approximately 5 g of CrTe 3 powder was loaded in an aluminum can and placed in an Orange cryostat capable of reaching 1.5 K. Data was collected with incident energies E i = 50 meV and 35 meV, and Fermi chopper frequency of 300 Hz.
Density functional theory calculations were performed in order to understand the atomic, magnetic, and electronic properties of CrTe 3 . The cleavage energy is calculated using the spin polarized vdW-DF-C09 functional [38][39][40][41] with norm conserving psedopotential using the Quantum Espresso package v5.4.0 [42], in which the spin is included in the non-local part of the exchange-correlation functional. For these calculations a model containing six layers per unit cell was constructed that included 192 atoms. The plane wave cut off was set to 80 Ry and the Brillouin zone sampling used for the supercell was 1 × 3 × 3. This approach does not, however, allow automated relaxation and energy minimization, so for the structural optimizations the Vienna Ab-initio Simulation Package (VASP version 5.3.5) [43] with the projectoraugmented wave (PAW) [44] potentials was used. In both codes, 3p 6 4s 1 3d 5 for Cr and 5s 2 5p 4 for Te were explicitly included as valence electrons in the pseudopotentials. Bulk CrTe 3 with various magnetic configurations were relaxed to get the optimized geometry. The exchange-correlation was approximated with the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) functionals [45] as well as the vdW-DF-optB86b functional [38][39][40]46] in order to incorporate van der Waals interaction during structural optimization. The Monkhorst-Pack [47] k-point sampling method in the Brillouin zone with a 4x4x4 and a 8x8x8 meshes for ionic and electronic optimization, respectively. The energy cutoff was 500 eV and the criteria for energy and force convergence are set to be 1×10 −4 eV and 0.01 eV/Å, respectively.

A. Crystal structure and thermal expansion
Diffraction data collected in this study are consistent with the published crystal structure of CrTe 3 [31], and crystallographic parameters determined at several temperatures can be found in the Supplemental Material [48]. The layered nature of the compound is apparent in Fig. 1a, where the van der Waals gap separating layers of composition CrTe 3 is visible. A plan view of a layer is also shown. The layers are made up of Cr 4 Te 16 units (Fig. 1a) comprising four edge-sharing CrTe 8 octahedra. The octahedra share corners with those in neighboring Cr 4 Te 16 units to form the two dimensional layers. The four Cr atoms in each unit are arranged in a rhombus or lozenge shape. The intra-lozenge Cr−Cr distances labeled by J 1 and J 2 in Fig. 1a are 3.62−3.64Å and 3.89Å, respectively, so no Cr-Cr chemical bonding is expected (the Cr−Cr dis-tance in Cr metal is 2.5Å). There are, however, short distances between Te atoms on neighboring Cr 4 Te 16 units, shown by the dark lines in the plan view in Fig. 1a, forming Te 2 dimers and Te 3 trimers [31]. Such bonding between anions is relatively common in telluride compounds [32]. Simple electron counting gives oxidation states of zero for the Te atoms with two bonds (centers of the trimers), 1-for the Te atoms with one bond (dimers and ends of trimers), and 2-for the Te atoms bonded only to Cr. This results in the expected oxidation state of 3+ for Cr. First principles calculations indicate strong covalency in this material, with nearly equal contribution from Cr and Te in states near the Fermi level, as shown in the Supplemental Material [48].
As noted above, crystals of CrTe 3 are easily separated into thin sheets due to the van der Waals bonding between the layers, suggesting these materials could be of some interest in forming van der Waals heterostructures. To further support and quantify this observation, first principles calculations were performed to estimate the cleavage energy, defined as the energy per unit area required to separate a crystal at the van der Waals gap into two. For the purpose of these calculations, a structural model was built that contained six CrTe 3 layers per unit cell and the distance (d) between two of the layers was gradually increased from the equilibrium value (d 0 = 3.5Å), as shown in Fig. 1d. The figure shows the calculated energy relative to the equilibrium energy as a function of (d − d 0 )/d 0 . The curve approaches a saturation value, the cleavage energy, of about 0.5 J/m 2 . This is comparable to well known cleavable materials used for studies of monolayer materials and construction of van der Waals heterostructures, including graphite 0. 43 [19,51].
The temperature dependence of the monoclinic lattice parameters, unit cell volume, and layer spacing determined from neutron powder diffraction measurements are shown in Fig.  2 with more details in the Supplemental Material [48]. The layer spacing is defined as the distance between neighboring layers measured perpendicular to the planes of the layers (along the a* reciprocal lattice direction) and is equal to a·sin(β). Very anisotropic thermal expansion is observed, as is often the case for layered, van der Waals bonded materials, in which the bonding between the layers is weak compared to the intralayer bonding. However, in CrTe 3 , the thermal expansion is largest along the b-axis, which lies within the layers, and the a-axis and layer separation have only a weak temperature dependence. This is seen more clearly in the lower panel of Fig. 2.
To investigate the origin of this unusual thermal expansion behavior, full single crystal diffraction datasets were collected at 250, 173, and 100 K, and refined to determine accurately the atomic positions and interatomic distances. Results of the refinements are included in the Supplemental Material [48]. Knowledge of the coordinates of each Te atom allows measurement of the van der Waals gap, defined as the perpendicular distance between Te planes on each side of the gap, and the layer thickness, defined as the perpendicular distance between Te planes on each side of a Cr plane (Fig. 1a). The layer spacing as defined above is the sum of these two distances. This analysis shows that the van der Waals gap distance does indeed decrease upon cooling, by 0.2% between 250 and 100 K, but this is compensated by an increase of nearly the same amount in the layer thickness, resulting in little change in the layer spacing. It is expected that the increase in thickness of the layers upon cooling is driven by the very strong contraction of the in-plane b lattice parameter seen in Fig. 2a. Within the Cr sublattice, the in-plane contraction arises mostly from the contraction of the tetramers themselves, which actually become further separated from neighboring tetramers upon cooling. The Te-Te dimers and Te-Te-Te trimers are quite rigid, with intra-dimer and intra-trimer distances changing by about 0.1% between 250 and 100 K, while distances between Te atoms not bound into these polymeric units decrease by about 1%.
In addition to the unusual relative thermal expansion among the different crystallographic directions, anomalous behavior in the temperature dependence of individual lattice parameters is also seen (Fig. 2). This is most apparent in the layer spacing, but the in-plane lattice parameter b and unit cell volume V also behave in an unusual way. In particular, negative curvature is observed in these data from room temperature down to about 150 K. These behaviors are expected to be related to the onset of magnetic correlations, since first principles calculations discussed below show notable variation in the optimized lattice parameters when different magnetic configurations are used. Thus the present results indicate strong magneto-elastic coupling is present in CrTe 3 .

B. Electrical resistivity
The semiconducting nature of CrTe 3 is illustrated by the electrical resistivity (ρ), which was measured with the current in the bc-plane of two crystals grown using a tellurium flux (Fig. 3). The same qualitative behavior is seen in both samples, but the absolute values differ by about a factor of 2 over the entire temperature range. This is likely attributable to uncertainty in the actual thickness through which the current flows in this easily delaminated material. Comparison of data collected in zero field and in a 50 kOe applied magnetic field indicates negligible magnetoresistance. The transport in CrTe 3 is clearly activated, with a change in activation energy occurring near room temperature. Linear fits to log(ρ) vs 1/T at high and low temperatures are shown for both samples in Fig. 3. From these fits, activation energies of 0.07−0.08 eV and 0.14−0.15 eV are determined at low and high temperatures, respectively. The low temperature behavior likely arises from activation of defects or impurities acting as donors or acceptors. Intrinsic semiconducting behavior may be responsible for the behavior above room temperature. If that is the case, a band gap of 0.28−0.30 eV can be estimated. An activation energy of 0.35 eV has been reported based on resistivity measurements on a cold-pressed polycrystalline sample of CrTe 3 , but that determination was made using data only below 300 K [33]. Interestingly, first principles calculations (see [48]) give a similar gap of 0.26 eV in the magnetically ordered state.

C. Heat capacity
The specific heat capacity (c P ) of a CrTe 3 single crystal is shown in Fig. 4. The data reach 24.4 J/mol-at./K at 200 K (Fig. 4a), close to the Dulong-Petit limit of 24.9 J/molat./K. No large anomalies typically associated with long range magnetic order are seen; however, close inspection reveals a weak anomaly that likely indicates a broad lambda-like peak centered at about 53 K (Fig. 4a, inset). Data were collected through this transition in zero field and in an applied field of 50 kOe, and no field effect is observed.
Lower temperature data on a single crystal was also collected to look for additional anomalies. The results are shown in Fig. 4b. Debye-like behavior is seen to persist down to 380 mK. A linear fit to c P /T vs T 2 for the data collected in zero field is shown. The y-axis intercept of the fit is nearly zero, as expected for a semiconductor with no T-linear electronic contribution to the heat capacity. Assuming the heat capacity in this temperature range is entirely due to phonons, a Debye temperature of θ D = 170 K can be determined (although antiferromagnetic magnons would also contribute to the T 3 term). This is consistent with the high temperature behavior of c P , which reaches 95% of the 3R at this temperature. Data collected with a 50 kOe magnetic field applied perpendicular to the bc-plane is shown in Fig. 4 as well. The magnetic field is seen to have a negligible effect.  Figure 5a shows the magnetic susceptibility (χ) measured on a single crystal between 2 and 380 K with the field along the layer stacking direction (out of plane) and with the field in the plane of the CrTe 3 layers (in plane). Higher temperature measurements were performed on a sample comprising four single crystals. Figure 5b is a composite of that high temperature data (300−720 K) and the powder average from the single crystal measurements shown in Fig. 5a. (2−380 K). At higher temperature χ shows an unusual temperature dependence, with a broad hump at room temperature. A Curie-Weiss model can only describe the data at the highest temper- Results of fitting the high temperature magnetic susceptibility data using models described in the text. atures (Fig. 5b). Fitting of data above 475 K gives an effective moment of 4.0 µ B per Cr (expect 3.87 µ B for Cr 3+ ) and a Weiss temperature of -285 K, indicating relatively strong antiferromagnetic interactions are present.

D. Magnetic Properties
Inspection of χ near the temperature at which the weak heat capacity anomaly was noted above (see Fig. 4a) reveals an anomaly in the magnetic susceptibility as well (Fig. 5a). Upon cooling below about 55 K anisotropy develops in the magnetic susceptibility due to a change in slope of the inplane data. This is highlighted in the inset of Fig. 5a. This behavior is interpreted as the onset of magnetic order, as supported by the neutron diffraction results below. A similarly weak anomaly is observed in the magnetic susceptibility of the quasi-1D magnetic material CrSb 2 near its Neel temperature [52,53]. In analogy with CrSb 2 , which was shown to have strongly one dimensional magnetic interactions [54], it is expected that magnetic correlations in CrTe 3 developing well above the long range ordering temperature are responsible for the weakness of the heat capacity and magnetic susceptibility anomalies near 55 K.
The broad hump in χ is suggestive of low dimensional magnetic correlations. For example the susceptibility data shown in Fig. 5a are reminiscent of the quasi-2D antiferromagnet K 2 V 3 O 8 , which shows a broad maximum near 10 K and long range magnetic order below 4 K [55,56], and layered thiophosphates like MnPS 3 , which shows a broad maximum near 120 K and long range order below 78 K [57,58]. The interpretation of the magnetic data as arising from 2D correlations is consistent with the layered nature of the crystal structure of CrTe 3 ; however, since the CrTe 3 layers themselves are made up of Cr 4 Te 16 units containing tetramers of Cr, a zero-dimensional or molecular magnetism viewpoint can also be considered. Magnetism associated with lozenge shaped tetramers formed by spin- 3 2 magnetic ions like those found in CrTe 3 has been studied in several systems [59][60][61][62], by assuming magnetically isolated tetramers with intra-tetramer interactions described by J 1 and J 2 defined in Fig. 1a. Results of fitting the the magnetic susceptibility of CrTe 3 using the model of Drillon et al. [61], assuming g = 2 and a spin of 3/2 per Cr, are shown in Fig. 5b. The model works reasonably well down to 350 K, and gives J 1 /k B = -39.9(1) K and J 2 /k B = 35(2) K. The model does not work well at lower temperatures, likely due to the other magnetic interactions that are not taken into account. Since long range magnetic order is observed, it is assumed that non-negligible inter-tetramer interactions are present in CrTe 3 . The same was noted for Na 3 RuO 4 by the authors of Ref. 63, and it was concluded that the noninteracting dimer or tetramer approximation is expected to be violated as a result. However, it is interesting to note that the negative J 1 and positive J 2 determined by the fit would favor an "up-down-up-down" spin configuration around the perimeter of the lozenges.
Neutron powder diffraction was used to look for further evidence of magnetic ordering in CrTe 3 . Comparison of data collected at 300 and 10 K reveals the presence of additional scattering at 10 K that is absent at higher temperatures (see Supplemental Material [48]). No additional scattering was observed in x-ray diffraction between 300 and 20 K, indicating the feature noted in Fig. 6a is magnetic in nature. The observation of this magnetic Bragg peak, which is centered near Q = 0.945Å −1 (d = 6.65Å), indicates the onset of magnetic order.
Further investigation using a single crystal sample at the HB3A diffractometer allowed identification of the magnetic ordering vector (wavevector) as k=(0.5, 0.5, 0.5) and the ordering temperature T N = 56(1) K, which is consistent with magnetization and heat capacity results. Figure 6a shows the order parameter measurement of the (0.5, 1.5, -0.5) magnetic Bragg peak as a function of temperature. The temperature dependence of the peak intensity below T N was fit to a power law I(T) ∝ (1-T/T N ) 2β , where β is the critical exponent. This is shown on the lower left inset of Fig. 6a. The data obeys this relationship reasonably well down to about 16 K below T N . A linear fit on the log-log plot in this temperature range gives a value of β = 0.37 (2), which is similar to the theoretical prediction for the 3D Heisenberg model. To evaluate the possible spin configurations, a magnetic symmetry analysis has been carried out using the tools available at the Bilbao Crystallographic Server [64]. The best fit of the 43 independent magnetic reflections has been obtained using a magnetic structure model defined by the magnetic space group P S 1 , shown in Fig 6b. In this model, each lozengeshaped tetramer has an up-down-up-down spin arrangement, with the alternating nonequivalent Cr sites carrying antiparallel spins. The successive tetramers alternate their spins' directions along all three crystallographic axes, leading to a doubling of the magnetic unit cell with respect to the nuclear cell in all three directions. The nearest neighbor tetramers can be arranged inside the bc plane in two equivalent configurations (magnetic domains) that produce stripes aligned along either the [0,1,1] or [0,-1,1] directions. The two magnetic domains are displayed in Figure 6c and d. Our neutron data indicates an almost equal domain population (0.502(7)% and 0.498(7)%). The refined magnetic moment is 2.0(2) µ B , with the following components along the crystallographic directions: m a = 2.1(1), m b = 0.7 (1), and m c = 1.1(1)µ B . The amplitude of the ordered moment is significantly smaller than the 3 µ B moment expected for trivalent Cr (S = 3/2). Strong covalency between Cr and Te or direct Cr-Cr chemical bonding could reduce the Cr moment; however, electronic structure calculations described below give the full 3.0 µ B moment for Cr. Thus the reduced moment observed by neutron diffraction may be ascribed to incomplete magnetic ordering between the spin-lozenges, with residual disorder most likely along the layer stacking direction. This is supported by the neutron powder diffraction data, in which the magnetic (0.5, 1.5, -0.5) is broader and more asymmetric than the nearby nuclear 100 reflection [48]. The intra-lozenge order is expected to be complete, as evidenced by the small value of χ and lack of a Curie tail at low temperatures (Fig. 5).
Stability of possible magnetic structures were compared using DFT calculations. Non-magnetic (NM), ferromagnetic (FM), and several antiferromagnetic (AFM) structures were considered. Table I shows the optimized lattice parameters from DFT calculations performed using different magnetic configurations and calculational methods. All calculations used a 2x1x1 supercell model. As expected, neglecting the long range van der Waals interaction (PBE) strongly overestimates the a-axis length. Of the magnetic models used, the AFM2 structure is the closest approximation to the experimental structure, which would require a 2x2x2 supercell to represent exactly. In AFM2 (and AFM3) every lozenge within one layer has the same up-down-up-down configuration, while in the experimental structure each layer contains stripes of lozenges with up-down-up-down and down-up-down-up configurations. All of the magnetic structures considered are at least 330 meV more stable than the NM state. Among the different magnetic configurations, AFM2 and AFM3 are more stable than FM and AFM1 by about 30 meV. Although this is near the limit of precision for the calculations, it does support the experimentally determined magnetic structure. The calculations give a magnetic moment of 3.0 µ B per Cr. The associated density of states in AFM2 state is included in the Supplemental Material [48].
Some variation in the relaxed lattice constants is observed when different magnetic configuration are used in the calculations (Table I). In particular, the in-plane lattice parameters (b and c) are are significantly affected by the inclusion of magnetism. This is identified as an indication of magneto-elastic coupling in CrTe 3 , which was also noted above in the temperature dependence of the lattice parameters (Fig. 2).
Results of inelastic neutron scattering experiments are shown in Fig. 7. At 1.6 K clear evidence of dispersive mag-  netic excitations are seen that could be associated with the magnetic ordering, which, according to the diffraction data, arises below about 56 K. The spin-wave-like excitations propagate out of the ordering wave-vector, at Q≈0.945Å −1 and extend up to approximatively 25 meV. The inelastic scattering intensity decreases with increasing wave-vector, following the decay of the magnetic form factor, which confirms the magnetic nature of these excitations. At 200 K , the intensity of inelastic scattering becomes much weaker and broadens in both energy and momentum space, consistent with an evolution to a paramagnetic state. Phonon scattering becomes more prominent at larger wavevectors above 4Å −1 . For an isolated tetramer scenario, dispersionless excitations would be expected with a Q dependence arising from the dynamic structure factor dependent on intra-tetramer atomic separations. Using the exchange values determined from the magnetic susceptibility, and the calculations for the allowed excitation en-ergies for a spin-3/2 lozenge model given in Eq. 3 of Ref. 63, one would expect to observe two localized magnetic modes at 3.4 meV and 6.0 meV. Instead, dispersive magnetic excitations are observed indicating the presence of significant intertetramer couplings and consistent with the long-range ordered ground state. Preliminary analysis indicates that the excitation spectrum can be reproduced using a spin wave model based on three exchange couplings, which includes the two intra-tetramer couplings (J 1 and J 2 ) and one antiferromagnetic inter-tetramer interaction (J 3 ). Similar dispersive magnetic excitations have been observed in Na 3 RuO 4 [63], and copper tellurate compounds Cu 2 Te 2 O 5 Cl 2 and Cu 2 Te 2 O 5 Br 2 [65], that were previously proposed to consist of coupled spin tetrahedra and exhibiting magnetic long-range order at low temperatures.

IV. SUMMARY AND CONCLUSIONS
CrTe 3 is an easily cleavable magnetic material with a unique layered crystal structure. The results presented here suggest significant magneto-elastic coupling in CrTe 3 and reveal low dimensional magnetic correlations extending to above room temperature and long range magnetic order below T N = 55 K. The structure contains lozenge shaped tetramers of Cr that adopt an up-down-up-down spin arrangement around their perimeters in the magnetically ordered state. There is stripe-like ordering among tetramers within planes and antiferromagnetic stacking between planes. Transport properties show semiconducting behavior consistent with a small band gap well above T N . This may arise from a Mott-Hubbard mechanism since Cr 3+ has half-filled t 2g orbitals, although strong hybridization between Cr and Te is indicated by electronic structure calculations. Based upon the observed electronic and magnetic properties, bulk CrTe 3 crystals should be of interest for the study of low dimensional magnetism, and the cleavability of the material makes it an exciting candidate for exploring magnetism in monolayer materials and for incorporation into van der Waals heterostructures. and Engineering Division, and Scientific User Facilities Division (neutron scattering). The authors thanks Ben S. Conner and Andrew F. May for helpful discussions through the course of the work, and Radu Custelcean for use of and assistance with the single crystal diffractometer.

SUPPLEMENTARY MATERIAL
Magnetic susceptibility vs temperature for two polycrystalline samples of CrTe 3 is shown in Fig. S1. Sample 1 is relatively high purity and the data shown are from measurements carried out in a field of 50 kOe. Sample 2 is less pure, with Cr 5 Te 8 providing a ferromagnetic component below about 280 K. For this sample the data shown were determined from linear fits to M vs H data between 30 and 50 kOe.
Rietveld fits to neutron powder diffraction data collected at 10 and 300 K are shown in Figure S2. The fitted crystallographic parameters are collected in Table S1. Low Q portions of the data are shown in Fig. S3, in which the magnetic Bragg peak (0.5, 1.5, -0.5) can be seen. The limited statistics of this data make determining the ordering temperature difficult, but the contour plot of intensity shown in the inset appears to be consistent with the single crystal data, which showed onset of magnetic scattering near 56 K. The relative breadth and apparent asymmetry of the magnetic peak suggests that there may be some residual magnetic disorder along the stacking direc-tion at 10 K.
The temperature dependence of the lattice parameters determined from neutron and x-ray powder diffraction measurements are shown in Figure S4. The largest discrepancy is seen in a, which varies by about 0.1% between the two measurements. This is expected to be due primarily to small sample to sample variations, and it may indicate some degree of stoichiometric variation in the CrTe 3 phase. Values for a ranging from 7.884 to 7.899Å were observed at room temperature based on powder x-ray diffraction measurements.
Results of single crystal x-ray diffraction refinements are collected in Table S2, which lists lattice parameters and agreement indices, as well as atomic coordinates and equivalent isotropic displacement parameters.
The calculated total and projected densities of states for CrTe 3 are shown in Fig. S5. This calculation was performed using the AFM2 configuration.
The positions and moments on the Cr atoms determined by single crystal neutron diffraction refinement are listed in Table  S3.   x-ray neutron a (Å) x-ray neutron b (Å) x-ray