Logarithmic Relaxation of the Nonequilibrium State of the Charge Density Wave in TbTe3 and HoTe3 Compounds

The measurements of electronic transport including the dynamic properties of the charge density wave (CDW) in the quasi-two-dimensional compound HoTe3 have been performed. The effects of the slow relaxation of the nonequilibrium state of the CDW during isothermal exposure in the zero current mode, previously observed in TbTe3, have been discovered and studied. A significant increase in the exposure time made it possible to clearly demonstrate that the relaxation is logarithmic. Relaxation features were studied in different temperature and time ranges. The data obtained indicate the glassy behavior of the CDW pinning centers in rare-earth tritellurides.


INTRODUCTION
The charge density wave (CDW), which was studied for many years mainly in quasi-one-dimensional compounds [1], has recently found more and more new manifestations in quasi-two-dimensional systems [2]. The properties of the "two-dimensional" CDW, despite the similarities with one-dimensional systems, have quite a few differences that require deep studying and theoretical explanation [3]. Rare-earth tritellurides, RTe 3 (R: La, Ce, Pr, Nd, Gd, Tb, Dy, Er, Tm), are one of the few quasi-two-dimensional systems in which almost all the main CDW effects are observed. Compounds of this group attract considerable attention due to the existence of various collective quantum states in them: superconductivity [4,5], several CDW types [6,7], magnetic ordering [8,9], as well as their interaction/competition with each other [10,11].
The compounds of the RTe 3 family listed above have the same crystal structure and a close band filling. The lattice parameters increase monotonically with the atomic number of a rare-earth element, which leads to the monotonicity of a number of properties: CDW transition temperature and wave vectors, anisotropy in the layer plane and even environmental stability. Thus, replacing a rare-earth atom with a lighter/heavier one is equivalent to creating a so-called chemical pressure in the system [6,7] and allows to fine-tune its properties, including the parameters of the charge density wave.
In these compounds, the most interesting dynamic CDW effects are also clearly observed: collective charge transfer (CDW sliding) and rf-synchronization of the CDW (Shapiro steps) [12][13][14]. As in quasi-onedimensional systems, the CDW sliding effect manifests itself in a sharp increase in conductivity when an electric field is applied above a certain threshold value , which is determined by pinning, i.e., "coupling" of the charge density wave with defects.
Recently, a number of unusual effects of extremely slow relaxation (the relaxation time is tens of hours) associated with the CDW collective motion were discovered in the TbTe 3 compound: during isothermal exposure of the sample, an increase in the threshold field with saturation was observed [15,16]. The change in the value of the threshold field could be controlled using external factors: current and temperature [17,18]. The effects, apparently associated with the dynamics of pinning centers, have much in common with relaxation processes in glassy systems, but their nature is still unclear, especially at the microscopic level.
The next logical step in studying these relaxation processes is to search them in another compound of RTe 3 class, which was implemented in this work. HoTe 3 was chosen as such a compound. The choice is due to a number of factors:

CONDENSED MATTER
1. Holmium is one of the elements closest in properties to terbium: their atomic numbers are 67 and 65, respectively, their molar masses differ by several percent, atomic radii-by fractions of a percent. Accordingly, similar properties of the compounds are expected, which makes it possible to use the sample fabrication and measurement techniques that have been developed for TbTe 3 .
2. Also favorable is the fact that the CDW transition in the HoTe 3 compound occurs at K [7]-just below room temperature. This makes it possible to study samples that had never undergone a CDW transition prior to measurements. In the TbTe 3 ( K [7]), storage of crystals at room temperature resulted in the formation of a state with increased pinning, which was observed in the temperature dependence of the threshold field as a maximum at 295-300 K, and traces of which remained even after a long sample annealing.
In both TbTe 3 and HoTe 3 compounds, both highand low-temperature CDWs in orthogonal directions in the layer plane are observed: , . Transition temperatures of lowtemperature CDW are as follows: K [19] in TbTe 3 and 110 K in HoTe 3 [7].
EXPERIMENT Transport measurements (resistance, current-voltage characteristics) of HoTe 3 microbridges were carried out in the temperature range of 77-350 K. The samples were prepared as follows [17,20]: HoTe 3 single crystals were grown using the self-flux method [12] and attached to a sapphire substrate with a thin layer of Stycast epoxy, then submicron films were fabricated by micromechanical splitting. The quality of the crystals was controlled using X-ray diffraction, the smoothness and uniformity of the detached films were controlled using a Carl Zeiss Axio Imager A2 optical microscope in the differential interference contrast mode. The use of double-sided polished sapphire as a substrate turned out to be extremely successful for samples of this type, firstly, due to its high thermal conductivity, and secondly, because of the possibility to control optically the quality of the surface on both sides of the film. We also managed to reduce the thickness of the adhesive layer to fractions of a micron, as evidenced by the presence of a Newton's ring type structure when observing the sample surface from the side of the substrate.
Using focused ion beam (FIB) etching, bridge-like structures 20-250 μm long and 8-100 μm wide, oriented along the direction of the wave vector were milled. After this operation the bridge area was encapsulated with a thin layer of polymer. The electrical contacts were fabricated by laser ablation and indium cold soldering. HoTe 3 is less stable at room conditions than TbTe 3 : after a few hours, signs of surface degradation can be evidenced by optical imaging. The method of sample preparation described above using vacuum systems (FIB, laser deposition, desiccators) and encapsulation made it possible to reduce the time of exposure to air on the working area of the sample to tens of minutes.
Conductivity and current-voltage characteristics (CVCs) of the structures were measured by the fourterminal method using a Keithley 2400 precision current source and a Keysight 34420A nanovoltmeter. All measurements were carried out in an inert helium atmosphere.
Within the framework of this experiment, 11 bridge structures were fabricated and measured. The effects described below were observed on all samples. This paper presents the results obtained on a sample with the lowest values of resistivity and threshold field-a bridge with a length μm, a width 23.5 μm, and a thickness nm. In HoTe 3 , as in TbTe 3 , when the sample is exposed at a certain temperature below the CDW transition, the threshold field changes with time. In this work, series of CVCs were successively measured during long-time exposure at temperatures of 200, 220, 240, 260, and 280 K, similarly to [16]. Within one series, the sample was cooled from a temperature to a given exposure temperature and then isothermally kept at for a long time, up to 210 h. During exposure, the threshold CVCs were measured at certain time intervals, from 10 to 120 min. Between measurements of CVCs, to control the exposure process, the time dependence of the sample resistance at a current of μA was recorded. Passing through the sample a current much lower than the CDW sliding threshold current ( μA) does not affect the time evolution of the I-V characteristics [17]; thus, it is equivalent to exposure in the zero current mode. The value of the threshold field = , corresponding to the beginning of the CDW sliding, was determined from a sharp break in the dependence obtained by numerical differentiation.

RESULTS AND DISCUSSION
A typical change in the threshold CVCs at different temperatures before and after a long isothermal exposure is shown in Fig. 1a. The threshold field after exposure increased from two to five times, depending on temperature (Fig. 1b). During exposure, the resistance value corresponding to the static CDW regime remained unchanged (Fig. 1a). After the completion of the series, the sample was briefly heated to followed by cooling to the next exposure temperature . It is noteworthy that after short-term heating and cooling of the sample to the = 100 same temperature , the threshold CVCs were completely reproduced, regardless of the temperature and duration of the exposure process preceding heating. Restoration of the original CVC shape by short-term heating was observed in both TbTe 3 and HoTe 3 compounds and indicates that the growth of the threshold field is not a diffusion process. An increase in pinning due to the diffusion of mobile defects has recently been observed in the quasi-one-dimensional compound TaS 3 [21,22].
The evolution of the threshold field at different temperatures in the range of 200-280 K (Fig. 2a) has a form similar to that observed earlier [16]-the growth of the threshold field slows down with time. It was assumed in [16] that the dependence of the threshold field has the form E t (t) = E 0 --the value of tends asymptotically to the saturation value (Fig. 2b). For a more accurate determination of , the exposure time was After rearranging the plots on a X-axis logarithmic scale, the dependences become linear (Fig. 3), that is, the threshold field changes with time according to the law ). It can be seen in the inset to Fig. 2 that for TbTe 3 at 280 K the relaxation was very likely a logarithmic one, however, due to the limited exposure time, it was concluded that the decay of to the saturation value is exponential for all T. New results for HoTe 3 made it possible to clearly demonstrate that the common case is more complex.
The logarithmic relaxation law is often found in complex highly disordered systems: ordinary [23] and spin glasses [24], Anderson dielectrics [25], as well as in biological objects, for example, in DNA [26] or in corn roots [27]. Elements in such systems, as a rule, interact intensively with each other at the micro level, while the interaction strongly depends on the history of the system. Neurons in the human brain, which are known to be closely interconnected with each other, also have a logarithmic distribution [28].
In both compounds, the presence of two time intervals on the relaxation dependence, Figs. 3a, 3b, with different slope is traced. In HoTe 3 the transition is clearly visible for most temperatures (Fig. 3b), and the break point for different temperatures is practically the same: h. A change in the slope of the logarithmic dependence is often observed in various glassy systems and is usually associated with the phenomenon of "aging" [29]-the dependence of the relaxation behavior of the system on the time of being in a state that precedes the relaxation process. Moreover, the time , at which . Within the framework of this experiment, there are no characteristic times on a scale of 10 h: the process of cooling to a given temperature takes no more than 30 min, and the temperature stabilization time before the start of a series of measurements does not exceed several minutes at all. Probably, the occurrence of a break in the relaxation dependence is not directly related to the aging process. At the same time, the question of the existence of aging effects in our system remains open. Aging is inextricably linked with glass systems, so its experimental detection could clarify the nature of the nonequilibrium state in RTe 3 compounds.
Another feature that can be observed in Fig. 3 is that the growth of the threshold field reaches its maximum at intermediate temperatures. For all relaxation curves in Figs. 3a, 3b, the slope angles were determined for both time ranges. As expected, the dependences of the obtained values on the exposure temperature (Fig. 4a)  The appearance of the maximum in the temperature dependence of the threshold field growth rate, despite the fact that other parameters of the RTe 3 compounds (resistance, threshold field) behave monotonically in this temperature range, may seem unexpected. However, there are several possible reasons for the slowdown in the growth rate of as at high as well as at low temperatures.
During cooling, a slowdown of thermal processes occurs, including one in the CDW system. If we assume that the system of pinning centers is in a highly disordered glassy state, which, as is known, is nonequilibrium, and the system slowly relaxes to the equilibrium state, then a decrease in temperature leads to an increase in the transition time between metastable pinning states and, thereby, a strong deceleration of the relaxation process. The slowdown in the growth of  at high temperatures is probably due to the proximity to the CDW transition-fluctuations destroy the nonequilibrium state (similar phenomena were also observed in spin glasses).
In RTe 3 compounds, in contrast to one-dimensional CDW materials, the temperature dependence of the threshold field, , has a strict linear form (Fig. 1b) [13], despite the fact that the value of is nonequilibrium. At the moment there is no theoretical model that explains this linear relationship. It is noteworthy that the temperature dependence of the threshold field after a long exposure in Fig. 1b corresponding a strong decay of the relaxation process has a completely different form, resembling a typical temperature dependence of the CDW gap. We also note that if both dependences are interpolated to low temperatures, they can intersect at a temperature of K, near K.
Considering that in HoTe 3 the values of at 200 and 280 K differ by almost an order of magnitude, it was decided to analyze how not only the absolute value of the threshold field changes with temperature, but also its ratio to the initial nonequilibrium value, . Figure 4b shows the dependence of at different temperatures for both compounds. In the initial time interval, the growth rate of increases monotonically with increasing exposure temperature, while at close to , after h the growth of practically stops. Perhaps, at temperatures close to the CDW transition, the formation of a completely ordered structure of pinning centers is impossible due to fluctuations, as a result of which, during relaxation, a quasi-equilibrium state with pinning close to the maximum value for a given exposure temperature is formed. The dependences in Fig. 4b are similar for both compounds.
One of the main questions related to the discovery of a nonequilibrium CDW state in RTe 3 compounds is whether the evolution of is a change in the pinning parameters that determine value or a new type of pinning is formed. The results obtained show that the rate of growth of the threshold field during relaxation correlates weakly with the initial value of the threshold field , which is an indication of the formation and evolution of a special type of pinning in RTe 3 compounds. Since the conductivity of the samples and the CDW transition temperature do not change during isothermal exposure, the observed changes in the threshold field can be due to changes only in the CDW subsystem. Most likely, these may be the rearrangement of the CDW macrostructure, for example, the formation and evolution of a specific domain structure as suggested in [16]. Elucidation of the nature of the observed effect requires further research.