Suppression of stacking order with doping in 1T-TaS 2 − x Se x

In 1T-TaS 2 − 2 x Se 2 x , the charge density wave (CDW) state features a star of David lattice that expands across layers as the system becomes commensurate on cooling. The layers can also order along the c-axis and different stacking orders have been proposed. Using neutron scattering on powder samples, we compared the stacking order previously observed in 1T-TaS 2 as the system is doped with Se. While at low temperature, a 13c layer sequence stacking was observed in TaS 2 , this type of ordering was not evident with doping. Doping with Se results in a nearly commensurate state with the Mott state suppressed which may be linked to the absence of the layer stacking.


I. INTRODUCTION
Quasi-two dimensional (2D) in nature, transition metal dichalcogenides (TMDs) 1T-MX 2 (M = Ti, Ta and X = S, Se, Te) are prone to electronic instabilities [1].1T-Ta(S/Se) 2 exhibits an incredibly rich phase diagram with multiple charge density wave (CDW) transitions emerging as a function of temperature and upon doping.In 1T-TaS 2−2x Se 2x , macroscopic behaviors such as CDW and superconductivity [2][3][4][5] have been observed, and more recently, a quantum spin liquid (QSL) has been proposed in 1T-TaS 2 as well [6].In the typical Peierls model for CDW order [7], the instability of the coupled electronlattice system brings a structural phase transition that is driven by strong electron-phonon coupling [8][9][10][11].The CDW formation can bring electron localization where displacements along phonon modes lower the total electronic energy by opening up a gap at the Fermi level, E F [12,13].This scenario, although applicable to simple one-dimensional systems, does not fully describe the case of 1T-TaS 2−x Se x where the CDW behavior is intertwined with the opening of a Mott gap [14].The origin of the CDW has been highly debated in TMDs [15][16][17].The Fermi surface nesting scenario most often does not apply.Existing models of the CDW order are broadly classified into three types: in one, it involves an excitonic condensation mechanism; in two, it involves a Jahn-Teller-like distortion mechanism; and in three, it involves a hybrid model, a combination of Jahn-Teller and exciton condensation [17].
1T-TaS 2 exhibits a strong CDW instability and electronic localization that leads to several interesting effects.Upon cooling from high temperatures, three main phases form: the high temperature incommensurate CDW (ICDW), the intermediate temperature nearly commensurate CDW (NCCDW) and the low temperature commensurate CDW (CCDW) [18].The ICDW appears below 540 K on cooling from the high temperature metallic state, with a transition from the P 3m1 crystal symmetry shown in Fig. 1(a) to the P 3 structure shown in Fig. 1(b).This transition leads to displacements of Ta ions that gives rise to the well-known star of David motifs.Upon cooling from the normal, high temperature metallic state, systematic displacements of the transition metal Ta leads to a star of David formation consisting of 13 Ta ions, in-plane.Domains of these formations expand to a commensurate CDW phase on cooling.Distinct from other CDW systems, in 1T-TaS 2 , the commensurate CDW state is accompanied by a metal-insulator (MI) transition that has been proposed to arise either from Mott localization or from disorder induced Anderson localization.Important to the MI behavior are the orbital ordering and out of plane correlations, as well as layer stacking order.
In the ICDW, the stars have limited ordering in-plane.Further cooling leads to the ICDW becoming NCCDW at T=350 K, where the √ 13 • √ 13 structural modulation first appears with a 12 o tilt relative to the original ab-plane.An expansion of the star of David motifs occurs in-plane [19].Below 180 K, the √ 13 • √ 13 structural modulation persists with a rotation of 13.9 o relative to the plane while the CDW becomes commensurate.The steps in the CDW transitions coincide with the kinks observed in the transport [20] as the system goes from the metallic to the insulating state.On the other end of the phase diagram, in 1T-TaSe 2 , the CCDW sets in at T ≈ 430 K, and the system shows no MI transition.It remains metallic down to the lowest temperature according to transport data [21].Between these two ends, superconductivity emerges upon doping that coexists with a broad NCCDW region in the phase diagram [18].The coexistence of superconductivity with CDW domains has been observed in other TMDs such as in the 2H polytype and in other systems such as the cuprates [22][23][24][25].
The electronic bands appear to undergo a continuous change with decreasing temperature in going though the many transition steps [26][27][28].In the absence of high temperature angle resolved photoemission spectroscopy (ARPES) due to resolution, there is no apparent nesting of the Fermi surface and a CDW gap is not necessarily located at the Γ point.Measurements suggested that the gap appears elsewhere in k-space [7,29].The domainlike CDW structures of the NCCDW and ICDW states in 1T-TaS 2 are discommensurate and semi-metallic, but when the CDW becomes commensurate, the Fermi surface disappears and either a Mott-Hubbard localization or a disorder induced Anderson localization sets in [29].Across the NCCDW-CCDW boundary, the Fermi surface is continuously reduced.This effect is convoluted by delectron localization that opens up an energy gap.In the CCDW phase, the gap is fully present, leading to a semiconducting state with about a 200 meV bandgap [7].
In the normal phase above 540 K, the Ta 5d band at the Γ point should be above E F [30].As the system goes through the NCCDW phase, this band becomes visible with ARPES.Further cooling to the CCDW state, band folding is observed because of the smaller Brillouin zone between 180 and 160 K, and an abrupt energy shift occurs due to opening of the energy gap [4].The loss of the Fermi surface continues with further cooling while the CDW gap continues to grow.The first order transition seen in the transport at 180 K on cooling is most likely due to a Mott-Hubbard localization [31,32].On warming, a different behavior is observed where the resistivity exhibits a hysteresis, with its value dropping 280 K, marking the CCDW-NCCDW transition.This has been attributed to be due to changes in the c-axis stacking order [4,33].
We report on the nature of the layer stacking order with temperature and doping.Earlier, we observed that the c-axis expands in the CCDW phase of 1T-TaS 2 on warming but drops at the crossover between the CCDW-NCCDW transition [34].It has been suggested that the localization of the d-electrons that brings the gap in the electronic structure depends on the expansion of the caxis [31,35].This in turn is related to the c-axis stacking order where changes in the interlayer coupling might drive the Mott transition.Neutron diffraction measurements confirmed the presence of 13c stacking order that disappears on warming across the CCDW-NCDW transition in 1T-TaS 2 .The appearance of the 13c layer sequence is expected to drive the Mott localization.The 13c stacking sequence was previously suggested in Ref. [36] from X-ray diffraction data down to 80 K.This study extends the data down to 2 K.Moreover, from single crystal measurements, we previously identified a 3c layer stacking as well, that commences in the ICDW state and continues to grow through the NCCDW to CCDW crossover [37] It is possible that both the 3c and 13c coexist at low temperatures in 1T-TaS 2 .Our single crystal data only reached 150 K.With doping, the neutron diffraction data clearly indicate that the 13c structure is suppressed.Its signature diffraction peak around 0.6 Å is not observed with doping.At the same time, it is not clear what happens to the 3c stacking with doping.Further experiments using single crystals are underway to elucidate the doping dependence of the 3c order.

II. RESULTS AND DISCUSSION
Shown in Fig. 1(a) is the hexagonal crystal structure of the high temperature undistorted lattice.Layers of the transition metal are separated by the chalcogen ion creating a quasi-2D lattice where weak interlayer interactions are expected due to the van der Waals nature of the forces holding the layers together.However, orbitals play an important role in this TMD and out of plane electron correlations lead to layer ordering and a gap in the density of states.The out-of-plane coupling is important to understand the electronic characteristics of these materials where band structure calculations suggested that opening a gap at the Γ point depends on the orbital order and out-of-plane stacking [31,35].Also shown in Fig. 1(a) is the low temperature crystal structure in the CCDW phase where the high temperature cell has undergone a √ 13 • √ 13 structural expansion and a rotation of 13.9 o relative to the primary axis.The star formation is the result of Ta displacements towards the middle Ta ion.The extent to which the star lattice spreads in the ab-plane depends on temperature.The star clusters expand on cooling giving rise to large domains in the CCDW state that are highly ordered, but become disordered on warming, breaking up into domains with star formations separated by regions of undistorted lattice.Three samples were measured using neutron scattering and the diffraction data are shown in Figs.1(c) and 1(d).Powder samples of TaS 2 , TaSSe and TaSe 2 were measured as a function of temperature.At 300 K, several superlattice reflections are indicated that might be due to 3c stacking structure.Similarly, at 2 K, the same superlattice reflections are observed as indicated.However, the intensity of the peaks is too small to be discerned from the powder diffraction data and single crystal experiments will help elucidate their presence.
The reciprocal lattice vector Q was calculated using Q = ha * 0 + kb * 0 + lc * + m 1 q 1 + m 2 q 2 , following the formalism introduced in [19].The modulation wave vectors q 1 = σ 1 a * 0 + σ 2 b * 0 and q 2 = −σ 2 a * 0 + (σ 1 + σ 2 )b * 0 for the CCDW phase were obtained based on the commensurate wave vector parameters σ 1 and σ 2 .The reciprocal lattice commensurate wave vector can be described as q ccdw = (σ 1 a * 0 + σ 2 b * 0 ) and for √ 13 • √ 13 in-plane translation, σ 1 and σ 2 values are 0.2308 and 0.0769 respectively.In the NCCDW phase, σ 1 and σ 2 values are 0.2448 and 0.0681 respectively.Shown in Fig. 1(e) is a plot of the diffraction pattern at very small momentum transfers, Q.At 5 K, a superlattice reflection belonging to 13c ordering is observed in 1T-TaS 2 at low temperatures.This reflection only appears from the 13c stacking order and not from the 3c order, even though most of the higher order reflections overlap between the two stacking models as shown in Fig. 1(f).The calculated positions of the satellite peaks corresponding to 3 c stacking order (c * = c * 0 /3) and 13 c stacking order (c * = c * 0 /13) are shown.Also shown in Fig. 1(e) are data for TaSSe in the same region of momentum transfer.In the TaSSe data, the reflection ∼ 0.6 Å−1 is notably absent which indicates that there is no 13c ordering in the superconducting state.A similar measurement was carried out for 1T-TaSe 2 and even though the sample was not a single phase of 1T, no evidence for the 0.6 Å−1 was observed.This indicates that stacking order might only be present in 1T-TaS 2 .The temperature and composition dependence of the Ta and S/Se thermal factors, ⟨U ⟩ 2 , lattice constants and unit cell volume are plotted in Fig. 2. As a function of composition, superconducting TaSSe has the largest thermal factor for the Ta ion that continues to increase on warming.Shown in Fig. 2(c) are the thermal factors for S and Se.Fig. 2(b) is a plot of the c/a ratio.In the case of TaSSe and TaSe 2 , the c/a ratio is almost constant as a function of temperature which indicates that the unit cell expands uniformly in the a-and c-direction.However, in TaS 2 , the ratio drops between 200 and 300 K because of the contraction of the c-lattice constant as previously observed in our earlier study [34] and by others [38].The contraction of the c-axis corresponds to the transition from the commensurate to the nearly commensurate state.We observed that this transition is coupled to the disappearance of the 13c ordering.Shown in Fig. 2(d) is the unit cell volume for the three compositions as a function of temperature.
Fig. 3 shows the results from the pair density function (PDF) analysis of the Ta displacements from the local structure at 2 K.The local structure is obtained by

III. MATERIALS AND METHODS
Powders were prepared using solid-state reaction.The neutron powder diffraction measurements were performed to investigate the structure through the multiple CDW steps.The time-of-flight (TOF) neutron measurements were carried out at the Nanoscale Ordered Materials Diffractometer (NOMAD/BL-1B) and at SEQUOIA (BL-17), a direct geometry spectrometer, at the Spallation Neutron Source (SNS) of Oak Ridge National Laboratory (ORNL) at temperatures ranging from 1.8 to 500 K.The aluminium can was used for SEQUOIA measurements and the empty can data were subtracted from the data.The reason SEQUOIA was used is that it reaches very small momentum transfers, not accessible to NO-MAD.The diffraction data from NOMAD were analyzed using the Rietveld refinement to obtain the unit cell parameters characterizing the crystal structure [39], resulting in what is referred to as the average model.The pair density function (PDF) analysis [40,41] provides information on the local arrangement of atoms in real space without the assumption of periodicity.It was performed on the same neutron diffraction data as the ones used for the Rietveld refinement.NOMAD is a diffractometer with a large bandwidth of momentum transfer Q, and it provides the total structure function S(Q).The S(Q) was Fourier transformed into real-space to obtain the G(r) [42,43].The instrument background and empty sample container were subtracted from the S(Q) and the data were normalized by a vanadium rod.A maximum Q of 40 Å-1 was used.

IV. CONCLUSIONS
The formation of hetero Layer stacking can be engineered to enable new behaviors and new properties.For instance it has been theoretically proposed that stacking of the honeycomb ferromagnet CrI 3 has the potential to give rise to ferroelectricity [44].Moreover, stacking in moire superlattices can create polar domains because of local spontaneous polarization [45].Hexagonal boron nitride was shown to exhibit ferroelectric switching in bilayers, leading to new concepts for functional heterostructures [46].Similarly, ferromagnetic heterostructures were demonstrated by stacking non-magnetic WS 2 with antiferromagnetic FePS 3 .At the interface, the FePS 3 shows ferromagnetism [47].
Layer stacking in homostructure TMDs maybe similalry linked to the transport behavior.1T-TaS 2 has an insulating CDW in contrast to other CDW dichalcogenides that have metallic CDW's.The reason for this is linked to Mott-Hubbard electron-electron correlations.Every star of David contributes one 5d electron to a half filled narrow conduction band.In the 13c layer stacking, there is an odd number of electrons and in the presence of large Coulomb repulsion acting on the layers, the Mott-Hubbard transition occurs [29].Density functional theory (DFT) calculations [33] have shown that the insulating phase and the MI transition originate not from the 2D order of the stars of David but by the vertical order.Our results confirm the significance of interlayer coupling and the insulating property.Interlayer stacking order in the CCDW phase has been verified to be a 13c repeat unit cell.This result contradicts the notion that the stacking is partially disordered in the CCDW state.Hence it is less likely that Anderson localization drives the MI transition.This result also contradicts the bilayer stacking model.

V. ACKNOWLEDGEMENTS
A portion of this research used resources at the Spallation Neutron Source, a DOE Office of Science User Facility operated by Oak Ridge National Laboratory.We thank Dr. John Schneeloch (University of Virginia) for valuable inputs on the sample growth and Dr. Utpal Chatterjee for valuable discussions about the ARPES data.
FIG. 1.(a) The high-temperature and low temperature crystal structure of 1T-TaX2.The lattice symmetry is the trigonal P 3m1 at high temperature which becomes P 3 at low temperatures.The star is the result of the Ta displacements.(b) A plot of the diffraction pattern at low temperatures showing the presence of the 13c superlattice in 1T-TaS2 although absent in 1T-TaSSe.The calculated peak positions of 1T-TaSSe corresponding to 3c and 13c stacking order is shown at the bottom.(c) The neutron powder diffraction data collected at 300 K compared among 1T-TaS2, 1T-TaSSe and 1T-TaSe2.All data are fit well using the P 3 symmetry.The diffraction peaks shift to the left with doping because Se is nominally a larger ion than S. (d) The neutron powder diffraction data collected at 2 K are shown.The arrows mark the positions of the CDW superlattice reflections.(e)The diffraction data plotted at very low Q indicate a superlattice peak corresponding to the 13c stacking order present in 1T-TaS2.(f ) A plot of the expected positions of 13c and 3c stacking order.

,FIG. 2 .
FIG. 2. (a) The Ta atomic displacement ¡U¿ 2 is shown as a function of temperature for the three compositions.Of the three, superconducting TaSSe shows the largest thermal factors.Shown in (c) are the S and Se thermal factors for the three compositions.In (b) is a plot of the c/a ratio for the three compositions and in (d) is a plot of the unit cell volume.The c/a ratio in 1T-TaS2 shows a decline on warming past 200K.