Switching VO2 Single Crystals and Related Phenomena: Sliding Domains and Crack Formation

VO2 is the prototype material for insulator–metal transition (IMT). Its transition at TIMT = 340 K is fast and consists of a large resistance jump (up to approximately five orders of magnitude), a large change in its optical properties in the visible range, and symmetry change from monoclinic to tetragonal (expansion by 1% along the tetragonal c-axis and 0.5% contraction in the perpendicular direction). It is a candidate for potential applications such as smart windows, fast optoelectronic switches, and field-effect transistors. The change in optical properties at the IMT allows distinguishing between the insulating and the metallic phases in the mixed state. Static or dynamic domain patterns in the mixed-state of self-heated single crystals during electric-field induced switching are in strong contrast with the percolative nature of the mixed state in switching VO2 films. The most impressive effect—so far unique to VO2—is the sliding of narrow semiconducting domains within a metallic background in the positive sense of the electric current. Here we show images from videos obtained using optical microscopy for sliding domains along VO2 needles and confirm a relation suggested in the past for their velocity. We also show images for the disturbing damage induced by the structural changes in switching VO2 crystals obtained for only a few current–voltage cycles.


Introduction
Many decades after its discovery [1][2][3], the steep and fast insulator-metal transition (IMT) of VO 2 at T IMT = 340 K remains one of the most attractive topics of research in condensed matter. To emphasize the importance of VO 2 and its potential applications, the reader is directed to a rather incomplete collection of recent articles that address this topic [4][5][6][7][8][9][10][11]. Here we focus on the I-M domain structure in the mixed state of VO 2 crystals during switching and to the damage caused by the structural changes [12]-both visible under the microscope due to the optical changes [13]. Twin boundaries of the insulating state are visible under the microscope due to polarization of the reflected light [14]. These twin boundaries disappear in the metallic state. Metal-insulator (semiconductor in the old notation) static or sliding domain patterns in the mixed state induced by heating [15] or self-heating [16] (upon passage of an electric current through the sample) have been recorded in movies. The most interesting dynamic pattern recorded in the seventies by one of us (Bertina Fisher) on a single crystal under constant current was of narrow semiconducting domains that nucleate at the positive edge of a static metallic domain and propagate towards the negative edge of the sample (in the sense of the applied current) [16]. The proposed interpretation-based on Reference [17]-suggested that sliding of the semiconducting (S) domains is driven by the Peltier effect at the two oppositely biased metal-semiconductor (MS) and semiconductor-metal (SM) interfaces [18]. Sliding in the conventional positive sense of the current is due to VO 2 being an n-type semiconductor (its Seebeck and thus its Peltier coefficient are negative). According to this model, Bertina Fisher suggested that the maximal velocity of the semiconducting domains (u max ) is obtained when the Peltier heat is fully exchanged with the latent heat; then: where Π SM is the Peltier coefficient of the VO 2 insulator (semiconductor)-metal couple, J is the current density, and L is the latent heat per unit volume for the IMT [18]. The measured velocity obtained there for one given current was smaller but close to the calculated u max . The sliding domains were accompanied by voltage oscillations [19]; each creation/annihilation of an insulating (semiconducting) domain caused an upwards/downwards jump of the voltage along the sample. Other interesting metal-insulator domain patterns consisted of chains of alternating metallic (dark) and insulating (bright) triangles [18]; these were often static or mobile with much lower velocity. A unique pattern found in a single sample was of stripes along the c-axis [20]. The topic of sliding domains remained dormant for several decades and was awoken in the nano-era [21][22][23]. In the respective reports, the drift of the M-I domain walls in pure or doped VO 2 current-carrying nano-wires were again attributed to the Peltier effect, on the basis of the sense of the drift which was consistent with the sign of Π SM . Relaxation-oscillations of frequency (f ) increasing linearly with current were obtained in a W-doped VO 2 nano-wire shunted by a capacitor; f as high as 5 MHz was obtained in a 1 µm nano-wire which sustained a single sliding metal-insulator boundary [21]. The frequency depends on the domains' velocities as well on the number of domains along the sample. Thus, the early single result of u(J) [18] remained solitary and called for additional attempts to confirm the validity of the old Equation (1). Growth runs that produce VO 2 single crystals-single domains in the insulating state-with shiny, perfect surfaces and dimensions suitable for studying sliding IM boundaries using optical microscopy are rare, but possible. Along the way, many needles were found in which Joule heating-induced switching caused cracks and breaks coinciding with rather erratic propagation of M-I walls and the formation of hot filaments. The damage caused by heating or self-heating in VO 2 due to the large symmetry and volume changes is well known. This acute problem sheds doubt even on the metallic resistivity of freshly-grown crystals of VO 2 because these crystals underwent at least one cooling-heating cycle through the transition [24]. We therefore decided to show some examples of the formation of cracks and breaks along with the impressive sliding domains. Three essentially identical growth runs (see Section 4) contained crystals with very different geometries and morphologies; D1 contained many thin, high-quality free-standing fibers; D2 contained platelets, many of them imperfect; V were of poor geometry, with multi-domain insulating patterns. Although sliding domains seem to be a common property of many of the switching VO 2 single crystals investigated, their velocity, u, as function of current-density, J, could be measured in only a few from batches D1 and D2. The results confirm Equation (1). Crack formation was recorded during switching of crystals from batch V. Figure 1a,b shows R(1/T) for two single crystals of VO 2 : P4 from batch V and P2 from batch D1 in the four-probe and two-probe configurations, respectively. The dimensions of P4 were: cross-section area A = 0.026 × 0.037 cm 2 and the distance between voltage probes d = 0.12 cm; the dimensions of P2 were A ≈ 0.005 × 0.005 cm 2 and Length l = 0.6 cm. The resistance jump at T IMT is close to five orders of magnitude for P4 and more than four orders of magnitude for P2. The calculated resistivities (ρ) of P4 and P2 at 300 K are 120 and 220 Ω cm, respectively. The resistivity of P4 just below the transition (320 µΩ cm) is comparable with ρ M (IMT) = 220 µΩ cm documented in Reference [24]. The difference between the values of ρ (300 K) of P4 and P2 may be due in part to the uncertainty in the geometry of P2 (its small A) and to some damage; while P4 was a freshly grown single crystal, P2 underwent at least five I-V loops with switching events before the R(T) run. The activation energies of conduction (∆)-which do not depend on geometry-are the same for the two samples, and they are in the range of the largest found in the literature, implying that conduction is intrinsic ( [25] and the references within).

R(T)
Materials 2017, 10, 554 3 of 11 in the range of the largest found in the literature, implying that conduction is intrinsic ( [25] and the references within).
(a) (b) The I-V characteristics were measured in the two-probe configuration using amalgam dots for contacts. With these contacts, the samples are free to move, being held only by surface tension. The large resistance jump in Figure 1b indicates that the resistance of the amalgam dots is low enough, even for metallic VO2, and certainly for its mixed state.

I-V Characteristics and Sliding Domains
I-V measurements were carried out at ambient temperature (RT) in the two-probe configuration. The measurements were carried out over many orders of magnitude of currents by varying the applied voltage and the resistance in series with the sample (the load resistance-RL). Upon increasing current, I(V) becomes non-linear due to self-heating, the voltage V reaches a maximum, and a current-controlled negative differential resistance (NDR) regime is obtained (see Appendix A). TIMT(VO2) is reached in a portion of the sample within the NDR regime, and the sample enters the mixed state.
In search of sliding domains, we started with crystals with close-to-perfect morphology. The best-looking sample from series D1 was labeled D1(1). The I-V loop shown at the top of Figure 2 for crystal D1(1) (dimensions 0.005 × 0.005 × 0.5 cm 3 ) is typical for switching of a high quality VO2 single crystal. The resistance in series with the sample (the load resistance-RL) for this loop is 100 K. The initial resistance of the sample at ambient temperature (20 °C) is 4.5 M. With increasing voltage, the I(V) trace becomes nonlinear due to self-heating; onset of the NDR regime occurs around 113 V. Upward switching between two steady states occurs along the upper dashed line, which represents the load line; the slope of the line is ≈ 2 × 10 −6 A/V (inverse of RL). Above the end of this short segment of the load line, I(V) enters a perfectly reproducible range. At the highest current, the sample's resistance is about two orders of magnitude lower than the R(RT) and less than one order of magnitude below RI(TIMT)-the sample is in a mixed state with a fraction of its volume in the insulating state. Backward switching occurs along the lower dashed line (parallel to the upper load line). Images of D1(1) were recorded on video during cycling through the I-V loops. Upward switching is marked by shuddering and twisting of the sample that often drive it out of focus, while downward switching is marked by a sudden freezing of motion. Sliding domains were observed for the whole range of the I-V loop between upward to downward switching. The sliding velocity increases with increasing current. The best microscope images were obtained for the slow domains (i.e., close to upward or downward switching); the contrast between the images of sliding domains and background becomes fuzzy when sliding is fast. Shuddering of the sample close to upward switching also disturbs the time-dependent The I-V characteristics were measured in the two-probe configuration using amalgam dots for contacts. With these contacts, the samples are free to move, being held only by surface tension. The large resistance jump in Figure 1b indicates that the resistance of the amalgam dots is low enough, even for metallic VO 2 , and certainly for its mixed state.

I-V Characteristics and Sliding Domains
I-V measurements were carried out at ambient temperature (RT) in the two-probe configuration. The measurements were carried out over many orders of magnitude of currents by varying the applied voltage and the resistance in series with the sample (the load resistance-R L ). Upon increasing current, I(V) becomes non-linear due to self-heating, the voltage V reaches a maximum, and a current-controlled negative differential resistance (NDR) regime is obtained (see Appendix A). T IMT (VO 2 ) is reached in a portion of the sample within the NDR regime, and the sample enters the mixed state.
In search of sliding domains, we started with crystals with close-to-perfect morphology. The best-looking sample from series D1 was labeled D1(1). The I-V loop shown at the top of Figure 2 for crystal D1(1) (dimensions 0.005 × 0.005 × 0.5 cm 3 ) is typical for switching of a high quality VO 2 single crystal. The resistance in series with the sample (the load resistance-R L ) for this loop is 100 KΩ. The initial resistance of the sample at ambient temperature (20 • C) is 4.5 MΩ. With increasing voltage, the I(V) trace becomes nonlinear due to self-heating; onset of the NDR regime occurs around 113 V. Upward switching between two steady states occurs along the upper dashed line, which represents the load line; the slope of the line is ≈2 × 10 −6 A/V (inverse of R L ). Above the end of this short segment of the load line, I(V) enters a perfectly reproducible range. At the highest current, the sample's resistance is about two orders of magnitude lower than the R(RT) and less than one order of magnitude below R I (T IMT )-the sample is in a mixed state with a fraction of its volume in the insulating state. Backward switching occurs along the lower dashed line (parallel to the upper load line). Images of D1(1) were recorded on video during cycling through the I-V loops. Upward switching is marked by shuddering and twisting of the sample that often drive it out of focus, while downward switching is marked by a sudden freezing of motion. Sliding domains were observed for the whole range of the I-V loop between upward to downward switching. The sliding velocity increases with increasing current. The best microscope images were obtained for the slow domains (i.e., close to upward or downward switching); the contrast between the images of sliding domains and background becomes fuzzy when sliding is fast. Shuddering of the sample close to upward switching also disturbs the time-dependent image. The four frames shown in Figure 2 were snipped from the video at the scene that corresponds to the point on the loop shown above the frames. The current at that point was 0.4 mA. The images correspond to four instants 0.3 s apart (see times at right of images). Bright (insulating) domains move towards the right (in the sense of the current) in a blue (metallic) background. The domain at the left of the upper image covers a distance of~200 µm (see right of the lowest image) over 0.9 s; the corresponding speed is u ≈ 0.022 cm/s. The current density J = 0.4/0.005 2 = 16 A/cm 2 . The latent heat of the IMT of 1020 cal/mole [26] corresponds to L = 240 J/cm 3 . If these data fit Equation (1), then uL/J = 0.33 V ≤ Π IM (T IMT ) = T IMT S IM (T IMT ), where S IM (T IMT ) is the relative thermopower of the metal-insulator couple (S I -S M ) at the transition. The above value of Π IM corresponds to a rather high value of S IM ; while S M << S I , it is safe to take S IM ≈ S I . S I (340 K) of −950 µV/K is rare in the literature, but comparable values were already measured in the past [18,26,27]. Π IM (T IMT ) was measured by two methods in Reference [17]; the value found there was 0.37 V. image. The four frames shown in Figure 2 were snipped from the video at the scene that corresponds to the point on the loop shown above the frames. The current at that point was 0.4 mA. The images correspond to four instants 0.3 s apart (see times at right of images). Bright (insulating) domains move towards the right (in the sense of the current) in a blue (metallic) background. The domain at the left of the upper image covers a distance of ~200 µm (see right of the lowest image) over 0.9 s; the corresponding speed is u ≈ 0.022 cm/s. The current density J = 0.4/0.005 2 = 16 A/cm 2 . The latent heat of the IMT of 1020 cal/mole [26] corresponds to L = 240 J/cm 3 . If these data fit Equation (1), then uL/J = 0.33 V ≤ IM(TIMT) = TIMT SIM(TIMT), where SIM(TIMT) is the relative thermopower of the metal-insulator couple (SI-SM) at the transition. The above value of IM corresponds to a rather high value of SIM; while SM << SI, it is safe to take SIM ≈ SI. SI(340 K) of −950 µV/K is rare in the literature, but comparable values were already measured in the past [18,26,27]. IM (TIMT) was measured by two methods in Reference [17]; the value found there was 0.37 V. Sliding domains were found in additional samples from the same series (D1), but no simple sequence could be obtained for measuring velocity. The samples were thinner than D1(1), and the pictures were blurred due to shuddering. In one case (sample D1(6) of dimensions 0.005 × 0.002 × 0.1 cm 3 ), sliding domains were observed in a wide range of currents and in a few snipped pictures were sharp enough to distinguish shapes. No narrow domains were observed in this sample, and the patterns consisting of wide often of triangular shapes changed frequently, divided, or coalesced as function of time. This scenario repeated itself in a few additional samples from this batch. Figure 6c in [23] illustrates the development of triangular semiconducting domains from the metallic state near the compressed edge of a nanoplatelet.
Sample D2 (1)   Sliding domains were found in additional samples from the same series (D1), but no simple sequence could be obtained for measuring velocity. The samples were thinner than D1(1), and the pictures were blurred due to shuddering. In one case (sample D1(6) of dimensions 0.005 × 0.002 × 0.1 cm 3 ), sliding domains were observed in a wide range of currents and in a few snipped pictures were sharp enough to distinguish shapes. No narrow domains were observed in this sample, and the patterns consisting of wide often of triangular shapes changed frequently, divided, or coalesced as function of time. This scenario repeated itself in a few additional samples from this batch. Figure 6c in [23] illustrates the development of triangular semiconducting domains from the metallic state near the compressed edge of a nanoplatelet.
Sample D2 (1) Figure 7. I(V) (y-axis at left) and the corresponding log R(V) (y-axis at right) obtained in a range of much higher currents is shown at the right of this Figure. Steady state and reproducibility are reached at very low voltages when most of the sample is metallic. At the highest currents, the slope of I(V) is slightly positive (i.e., metallic-like). Upon backward switching, the resistance reaches its initial value; this unusually robust sample underwent this loop undamaged. The maximal change in the R(V) was more than four orders of magnitude. No sliding domains were observed in the range of currents of the loop at right. In a following attempt to increase the range of currents (by decreasing R L ) the sample's resistance became infinite-the sample collapsed to thermal runaway.     Figure. This sample exhibits features similar to those of D1(6) (i.e., patterns that change with time). The width of a thin, horizontal, barely visible semiconducting domain in the first frame increases with time, the domain changes inclination, and after 2 s, becomes a thick domain at the expense of the neighboring domain. It is notable that the I-V characteristic shown at the right is very similar to that at the right of Figure 3.
Several I-V loops were recorded for sample D2(8) (dimensions 0.019 × 0.009 × 0.15 cm 3 ). Three loops (1, 2, and 3) for relatively high currents were practically identical with (V max (1) = 30.8 V, V max (2) = 30 V, and V max (3) = 29 V), respectively; for clarity, only the first loop is shown in Figure 5. Pairs of frames snipped from the videos at the scenes corresponding to the point marked on the first I-V characteristic, for the instants shown at the right of the frames are shown in the Figure. The M-I patterns of the three loops at the same current are very different. They have in common the appearance of a thin, slowly drifting semiconducting domain.   Sample D2(9) (dimensions 0.0057 × 0.0045 × 0.2 cm 3 ) with its rough matte surface was quite surprising. Sliding domains were observed over a wide range of currents above the switching threshold. On the loop shown in Figure 6, we stopped the current four times to take videos (see the I-V loop in Figure 6). Additional data were obtained for a second loop (not shown) which was similar but not identical. The domains were hardly visible within the rough portion of the surface; their rapid movement helped distinguish them from the noisy background. They became visible in a small portion of the sample that was smooth enough. A sequence of three snipped images, 0.1 s apart, are shown at the right of the I-V characteristic; these show a domain travelling in the smooth part of the sample. The domain shown here is perpendicular to the samples' edges-a rather rare case. The inclinations of other domains observed in this sample are of the more frequent ±50 • . Sample D2(9) (dimensions 0.0057 × 0.0045 × 0.2 cm 3 ) with its rough matte surface was quite surprising. Sliding domains were observed over a wide range of currents above the switching threshold. On the loop shown in Figure 6, we stopped the current four times to take videos (see the I-V loop in Figure 6). Additional data were obtained for a second loop (not shown) which was similar but not identical. The domains were hardly visible within the rough portion of the surface; their rapid movement helped distinguish them from the noisy background. They became visible in a small portion of the sample that was smooth enough. A sequence of three snipped images, 0.1 s apart, are shown at the right of the I-V characteristic; these show a domain travelling in the smooth part of the sample. The domain shown here is perpendicular to the samples' edges-a rather rare case. The inclinations of other domains observed in this sample are of the more frequent 50°.
The old and new data for the sliding velocity u(J) of narrow semiconducting domains are plotted in Figure 7. A linear trend-line was fitted to the data for the sample of 1974 [14], D1(1) (Figure 2), and D2(1) (Figure 3). The slope of the fitted line ( SM/L according to Equation (1)) is 0.0012 A −1 ·s −1 ·cm 3 . For L = 240 J/cm 3 , SM  0.29 V, in good agreement with the values of SM for pure VO2. Note that u(J) of D2(1) at low currents lie below the fitted line. The data of u for the narrow domains seen in sample D2(8) (Figure 5) lie far below the line (the lowest point is for Loop 1) of the series. The set of u(J) data for sample D2(9) lie in between. The slope of the fitted straight line to these data is 0.0005 A −1 ·s − 1·cm 3 . This is probably related to the poor morphology of the sample.
Supplementary material S1 contains video clips for samples "1974", D1(1), D2(1), D2(6), D2(8), and D2(9). The old and new data for the sliding velocity u(J) of narrow semiconducting domains are plotted in Figure 7. A linear trend-line was fitted to the data for the sample of 1974 [14], D1(1) (Figure 2), and D2(1) (Figure 3). The slope of the fitted line (≤ Π SM /L according to Equation (1)) is 0.0012 A −1 ·s −1 ·cm 3 . For L = 240 J/cm 3 , Π SM ≥ 0.29 V, in good agreement with the values of Π SM for pure VO 2 . Note that u(J) of D2(1) at low currents lie below the fitted line. The data of u for the narrow domains seen in sample D2(8) (Figure 5) lie far below the line (the lowest point is for Loop 1) of the series. The set of u(J) data for sample D2(9) lie in between. The slope of the fitted straight line to these data is 0.0005 A −1 ·s − 1·cm 3 . This is probably related to the poor morphology of the sample.
Supplementary material S1 contains video clips for samples "1974", D1(1), D2(1), D2(6), D2(8), and D2 (9). Sample D2(9) (dimensions 0.0057 × 0.0045 × 0.2 cm 3 ) with its rough matte surface was quite surprising. Sliding domains were observed over a wide range of currents above the switching threshold. On the loop shown in Figure 6, we stopped the current four times to take videos (see the I-V loop in Figure 6). Additional data were obtained for a second loop (not shown) which was similar but not identical. The domains were hardly visible within the rough portion of the surface; their rapid movement helped distinguish them from the noisy background. They became visible in a small portion of the sample that was smooth enough. A sequence of three snipped images, 0.1 s apart, are shown at the right of the I-V characteristic; these show a domain travelling in the smooth part of the sample. The domain shown here is perpendicular to the samples' edges-a rather rare case. The inclinations of other domains observed in this sample are of the more frequent 50°.
The old and new data for the sliding velocity u(J) of narrow semiconducting domains are plotted in Figure 7. A linear trend-line was fitted to the data for the sample of 1974 [14], D1(1) (Figure 2), and D2(1) (Figure 3). The slope of the fitted line ( SM/L according to Equation (1)

Crack Formation in Switching VO 2 Crystals
Crack formation during switching is vivid and colorful in videos. These become non-interesting and even boring when turned into sets of still pictures. Nevertheless, the problem of damage in switching VO 2 samples is too serious to be ignored. From a variety of videos, only a few instants from the troubled history of two samples were chosen. The two crystals are V1 (dimensions 0.010 × 0.0034 × 0.1 cm 3 ) and V2 (dimensions 0.013 × 0.0078 × 0.152 cm 3 ), and the chosen images are shown in Figure 8a,b. Figure 8a shows a portion of sample V1 before I-V cycling and during switching in the fifth I-V loop (left, when the cracks first became visible), and during the seventh loop. The dark portions of the sample are metallic. The maximal current density through the sample was 60 A/cm 2 , about one order of magnitude lower than the current densities in that left sample D2(1) intact. Figure 8b shows a portion of sample V2 at three instants during switching in its eighth I-V loop (left). The maximal current density in this sample was J max ≈ m10 A/cm 2 , but the crack appeared at lower J. A few seconds later, the crack expanded to the right.

Crack Formation in Switching VO2 Crystals
Crack formation during switching is vivid and colorful in videos. These become non-interesting and even boring when turned into sets of still pictures. Nevertheless, the problem of damage in switching VO2 samples is too serious to be ignored. From a variety of videos, only a few instants from the troubled history of two samples were chosen. The two crystals are V1 (dimensions 0.010 × 0.0034 × 0.1 cm 3 ) and V2 (dimensions 0.013 × 0.0078 × 0.152 cm 3 ), and the chosen images are shown in Figure 8a,b. Figure 8a shows a portion of sample V1 before I-V cycling and during switching in the fifth I-V loop (left, when the cracks first became visible), and during the seventh loop. The dark portions of the sample are metallic. The maximal current density through the sample was 60 A/cm 2 , about one order of magnitude lower than the current densities in that left sample D2(1) intact. Figure 8b shows a portion of sample V2 at three instants during switching in its eighth I-V loop (left). The maximal current density in this sample was Jmax  m10 A/cm 2 , but the crack appeared at lower J. A few seconds later, the crack expanded to the right.
Supplementary material S2 contains video clips for samples V1 and V2.

Discussion
The detailed properties of materials are best studied in single crystals. In this work, we investigated VO2 crystals prepared according to the same protocol. Samples from three different batches had the same TIMT-a property regarded as very sensitive to the purity of the material-but Supplementary material S2 contains video clips for samples V1 and V2.

Discussion
The detailed properties of materials are best studied in single crystals. In this work, we investigated VO 2 crystals prepared according to the same protocol. Samples from three different batches had the same T IMT -a property regarded as very sensitive to the purity of the material-but behaved very differently upon switching to the mixed M-I state. This shows that very subtle details of the crystal growth which were not yet under control govern the quality of the products.
The main property investigated in this work is so far unique to switching VO 2 single crystals: sliding semiconducting (SC) domains within a metallic background, in the conventional positive sense of the electric current. Moving boundaries were detected in all the samples of VO 2 at currents above switching to the mixed M-I state. Narrow SC domains of constant shape and inclination sliding with constant velocity were scarce. Nevertheless, we were able to add data that confirm Equation (1). We suggest that the material property that governs sliding domains in VO 2 samples is homogeneity at all scales (implying an absence of twin boundaries in the insulating state, or of any other defects). One cannot imagine a pair of parallel M-I and I-M barriers sliding as a solid entity if the current density is inhomogeneous.
Very large inhomogeneity of the current density is probably the cause of crack formation in crystals. This is also a warning signal for (polycrystalline) films of VO 2 used in applications; films are more robust than the single crystals, but they may not be immune to the damage accumulated after cycling countless times through the transition. A significant effort is presently invested in developing strategies to avoid the physical destruction of VO 2 . Highly repetitive (>10 4 ) cycling through the phase transition without structural damage was demonstrated in a nanocomposite of VO 2 [28].
Thermally-, Joule-heating-or photo-induced resistive switching in materials exhibiting IMT above RT have many potential applications [29], but such materials are scarce and the high temperature at which the devices function affect their stability. One such material is Ca 2 RuO 4 ; T IMT of pure single crystals is 357 K, the resistance jump is less than one order of magnitude, and the magnitude of the crystal axes change is comparable to that of VO 2 [30]. Single crystals shatter violently upon cycling through the transition, which makes the resistivity measurements a difficult task. From the large family of vanadium oxides exhibiting IMT, only VO 2 and V 3 O 5 have T IMT > RT. T IMT (V 3 O 5 ) = 428 K, the resistance jump is approximately two orders of magnitude, and the changes in symmetry and lattice parameters are small. Homogeneous heating of V 3 O 5 samples up to 450 K-even for long durations-do not affect R(T). Self-heating by high currents induced significant changes in repeated I-V characteristics due to changes in oxygen stoichiometry [31]; the most dramatic effect was the reversible transformation of a small portion of the sample into VO 2 and back to V 3 O 5 . These undesired changes were prevented by limiting the range of currents.
The static or dynamic natural M-I boundaries in the mixed state of VO 2 obtained upon applying a thermal gradient on the sample around T IMT or by passing an electric current through the sample provide an interesting and timely topic of investigation in the intense modern studies of surfaces and interfaces [32].
A final question is in order: are the sliding domains a property unique to VO 2 ? Can such domains be found in a system with IMT below RT, such as V 2 O 3 ?

Materials and Methods
Single crystals of VO 2 of various shapes, needles or platelets of various sizes (the longer, the thinner), have been grown by self-flux-evaporation [16] from Aldrich 99.99% V 2 O 5 powder with the addition of small quantities of VO 2 crystallites. Three essentially-identical growth runs (D1, D2, and V) contained crystals with very different microstructures; D1 contained many thin high-quality free-standing fibers; D2 contained platelets, many of them imperfect; V were of poor geometry, with multi-domain insulating patterns. Nevertheless, R(T) of all three types of crystals proved that they were pure VO 2 . R(T) of the VO 2 crystals was measured in the four-probe or two-probe configurations (with connected adjacent voltage and current probes). The four-probe configuration was employed for determination of the activation energy of conduction in the semiconductor close to IMT and the magnitude of the resistance jump. The contacts were indium-amalgam dots. These are ohmic contacts of low resistance for semiconducting VO 2 . The magnitude of the activation energies and the resistance jump at the IMT of two-probe and four-probe devices from these batches are usually comparable, the latter showing that the resistance of the amalgam dots is low enough even for metallic VO 2 samples and certainly for samples in the mixed state. With these contacts, the samples are free to move, being held only by surface tension. I-V measurements were carried out at ambient temperature in the two-probe configurations. These measurements were carried out over many orders of magnitude of currents by varying the applied voltage and the load resistance. The I-V loops were recorded on a YEW type 3036 X-Y recorder while the samples were viewed under the microscope. Photos and videos were taken with the camera of an iPhone 5 attached to a Zeiss stereo microscope. Dynamic phenomena were processed using the Movie Maker Editor.