{'fraction_non_alphanumeric': 0.033995323025374335, 'fraction_numerical': 0.021535024955498937, 'mean_word_length': 4.908847184986596, '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}}
Reversible switching between pressure-induced amorphization and thermal-driven recrystallization in VO2(B) nanosheets

Pressure-induced amorphization (PIA) and thermal-driven recrystallization have been observed in many crystalline materials. However, controllable switching between PIA and a metastable phase has not been described yet, due to the challenge to establish feasible switching methods to control the pressure and temperature precisely. Here, we demonstrate a reversible switching between PIA and thermally-driven recrystallization of VO2(B) nanosheets. Comprehensive in situ experiments are performed to establish the precise conditions of the reversible phase transformations, which are normally hindered but occur with stimuli beyond the energy barrier. Spectral evidence and theoretical calculations reveal the pressure–structure relationship and the role of flexible VOx polyhedra in the structural switching process. Anomalous resistivity evolution and the participation of spin in the reversible phase transition are observed for the first time. Our findings have significant implications for the design of phase switching devices and the exploration of hidden amorphous materials.

M aterials that contain no long-range structural order (for example, glass or amorphous phase) are fundamentally interesting to both the basic sciences and for practical industrial applications 1,2 . Glass, which has been considered a metastable and kinetically frozen disordered state, can be synthesized by rapidly quenching a melt from high temperature to prevent crystallization. Another frequently reported but distinct approach to amorphous phases or glasses is via pressure-induced amorphization (PIA) 3,4 . The disordering process in PIA differs significantly from the thermally-induced disorder in substance melting, where a much higher density is expected. PIA has been observed in many materials such as ice [5][6][7] , a-quartz SiO 2 (refs 8-10), AlPO 4 (ref. 11), R-Al 5 Li 3 Cu (ref. 12), and so on, and is widely accepted as an important condensed matter phenomenon. In some cases, mechanical instability caused by 'thermodynamic melting' and increased atomic coordination are considered to contribute most to PIA 8,13 . While in other cases, such as ZrW 2 O 8 , negative thermal expansion (NTE) may have some possible connection with the PIA [14][15][16][17] . However, determining the underlying mechanism still remains one of the most fascinating challenges and some recent evidence even shows that previously reported PIAs were due to the fragmentation of bulky particles into nanocrystals and strongly correlated pressurization environments 3 . Furthermore, new high-density PIA generated materials are expected to drive new theoretical approaches for modelling this phenomenon, which is essential for discovering their practical applications. Recently, pressureinduced crystallization revealed the topological ordering in metallic glasses at room temperature 18 . But upon heating, the recovered amorphous phase tends to return to its stable crystalline phase to minimize the system's energy, depending on the kinetic barrier it needs to overcome. Such examples include the main group compounds SnI 4 , LiKSO 4 , Ca(OH) 2 , clathrasils and berlinite AlPO 4 (refs 19-23). Molecular dynamics calculations on PIA-AlPO 4 (ref. 23) showed that the presence of non-deformable, fourfold coordinated PO 4 interlinking units in the crystal structure have a crucial role in the reversible amorphization and crystallization phase transition, where they acted as templates around which the original structure and even the original orientation were restored 23,24 . So far, reversible PIA has only been observed in the above-mentioned main group compounds, where only the lattice participates in the intriguing structure transformation. The study of more variables (that is, charge, spin and orbital) involved in PIA and the subsequent recrystallization processes could potentially allow a more comprehensive understanding of the order-disorder transformation mechanism and the electronic behaviour in highly disordered materials.
Discovering new external-stimuli responsive compounds with switchable ground states is a major objective in material science because these materials often lead to unusual phenomena or useful functionalities 25,26 . When considering PIA materials as metastable and dense phases distinct from their crystalline forms, it is important to achieve a controllable phase switching between either the PIA glass and their crystalline polymorphisms or even between two PIA phases.
In this work, we report a reversible structural switch between the metastable crystalline phase VO 2 (B) and its PIA glass, in the form of nanosheets. The phase switch between these two metastable phases is realized with low-pressure compression (B20 GPa) and relative low-temperature annealing (B200°C). The phase transformation and underlying mechanism is thoroughly studied using in situ synchrotron X-ray diffraction (XRD), infrared, Raman spectroscopy and theoretical calculations. High-pressure electronic transport properties and magnetic properties are also studied to provide direct evidence of the participation of charge and spin during the phase transformation, for the first time. Our findings may provide a new research platform for the exploration of novel amorphous materials and phenomena under controllable external stimuli.

Results
Material and crystal structure. There are several crystalline phases of VO 2 at ambient conditions: VO 2 (A), VO 2 (B), VO 2 (M), VO 2 (R) and other metastable phases [27][28][29][30][31] . Among them, VO 2 (M) is the thermodynamically stable phase but it undergoes an insulator-to-metal transition (ITM) under high pressure to a metallic phase VO 2 (R) around room temperature. The structures and physical properties of the phases differ significantly due to the distinct coordination environments of V atoms in the VO x polyhedra, the V 4 þ -V 4 þ interactions and the various cross linking manners ( Supplementary Fig. 1). VO 2 (B) is a thermodynamically metastable phase adopting an anisotropic layered structure, and is frequently used as a battery material 32 . Its unique structural features, such as well-embedded layers and hierarchical V-O bonding, also make VO 2 (B) an interesting candidate for structural stability studies under high pressure. In this work, VO 2 (B) nanosheets were synthesized from high pure V 2 O 5 raw material via a hydrothermal route with citric acid as the reducing agent 33 . The product consisted of well-shaped nanosheets that were several micrometres in length/width and tens of nanometres in thickness, and its single crystal nature is proven by the well-arranged spots in the selected-area electron diffraction pattern (Fig. 1a,b). The structure of VO 2 (B) is characterized by the layered feature in the bc plane with a twofold connected O1 atom between the layers (Fig. 1c). Less electron charge distribution between the layers makes the [100] direction more compressible, and there is also a distinct gap along the [001] direction. Overall, the VO 2 (B) structure has a hierarchical structure consisting of condensed face-sharing VO x polyhedra. More structural features of VO 2 (B) along different axes are shown in Supplementary Fig. 2.
PIA at room temperature. The in situ synchrotron XRD patterns of VO 2 (B) nanosheets, as a function of pressure without pressure transition medium, are shown in Fig. 2a. There was no structural transition until the onset of amorphization around 20 GPa. The XRD pattern at ambient conditions is wellindexed with the monoclinic space group C2/m and lattice constants a ¼ 12.054(3) Å, b ¼ 3.693(1) Å, c ¼ 6.424(2) Å and b ¼ 106.96(1)°, in good agreement with other reported values 12,13 . As discussed above, from a structural chemistry viewpoint, the VO 2 (B) nanosheet was expected to show an anisotropic compressibility under pressure. The lattice parameters at various pressures before PIA were obtained by Reitveld refinements of these XRD patterns, as shown in Supplementary  Fig. 3. During compression, the diffraction peaks became broader and weaker from about 17 GPa, and completely vanished around 20 GPa, indicating the loss of long-range ordering in the pressure-amorphized state. To check the shear and nonhydrostaticity effect on the PIA, we repeated the experiments using neon and helium as pressure-transmitting medium (PTM) for comparison. In all cases, the PIA process occurred and the amorphous PIA-VO 2 (B) was preserved to ambient conditions after releasing pressure. In contrast, no PIA was observed in the thermal-stable VO 2 (M) phase up to 55 GPa 34 .
Spectroscopic techniques can probe the short-range structural features of local atomic coordinates. We employed infrared and Raman spectroscopy to examine the local structural evolution of VO 2 (B) nanosheets during compression and decompression. Figure 2b displays the infrared spectrum of VO 2 (B) as a function of pressure up to 23.6 GPa. Characteristic peaks located around 1,000 cm À 1 (at ambient conditions) are assigned to the shortest V 4 þ ¼ O1 bonds pointing perpendicularly into the bc interlayer. The mode frequency barely shifted (o50 cm À 1 ) towards high wave numbers upon compression. However, the broadening of the vibrational mode increases with pressure, and finally merging with other broad bands down to 600 cm À 1 is a more profound change. These were associated with angle deformations of the VO x polyhedral upon compression and eventually a disordered state resulted in amorphization. This evidence clearly indicates the nature of the degenerated chemical environments of oxygen around the V 4 þ centres. As pressure increases, the enhanced oxygen coordination around vanadium atoms in VO 2 (B) is expected to lead to the dynamical lattice instability, which triggers the PIA. It is interesting that the preserved PIA-VO 2 (B) has similar local vibrational modes to the pristine VO 2 (B) nanosheets, indicating that the short-range structure features of VO 2 (B) are preserved, as indicated in Fig. 2d. Moreover, Raman spectroscopy was employed to evaluate the contribution of the electron-phonon interactions to the PIA phase transition of VO 2 (B) (Fig. 2c). At low pressure, only five bands near 198, 260, 570, 780 and 1,020 cm À 1 were observed with moderate intensity in the VO 2 (B) nanosheets. Typically, the peak located at 1,020 cm À 1 can be assigned to the V 4 þ ¼ O1 stretching modes of terminal oxygen atoms. Upon compression above 10 GPa, all the bands weakened and finally vanished. Meanwhile, broad bands between 600 and 1,000 cm À 1 were observed corresponding with the PIA process. As anticipated, similarly to the infrared results, all of the atomic location information was preserved with a more disordered state as indicated by the broadening of the vibrational bands in Fig. 2e, which closely resembles the Raman changes during the ITM process from VO 2 (M) to metallic VO 2 (R).
Thermal-driven recrystallization and the phase diagram. The PIA-VO 2 (B) phase returns to the pristine VO 2 (B) structure upon annealing at relatively low temperature (B200°C by in situ annealing experiment) for a short time (B5 min). Figure 3 shows the XRD pattern of the starting VO 2 (B) nanosheets and PIA-VO 2 (B) before and after annealing at 250°C (50°C above the critical temperature to guarantee a complete recrystallization). To ensure that the PIA process was achieved, a higher pressure of 31.5 GPa (far beyond the PIA starting point) was applied. As discussed, the loss of the X-ray diffraction peaks in the recovered PIA-VO 2 (B) sample suggests the absence of long-range ordering. After annealing at 250°C in a vacuum environment for 5 min, surprisingly, we noticed that the diffraction pattern from this heat-treated sample showed the same powder characteristic peaks as the VO 2 (B) phase with a monoclinic space group C2/m, lattice Supplementary Fig. 4) and a little expanded cell volume (273.5 Å 3 versus 272.2 Å 3 of the pristine sample). This structure switching phenomenon is distinct from the previously reported reversible PIA phenomena with the following characteristics: 23,24 (1) In the cases of some so-called 'memory glasses', such as AlPO 4 , the PIA phases can be restored to their initial crystalline structures spontaneously at room temperature once pressure is released. In contrast, PIA-VO 2 (B) can exist as a metastable, high density, intermediate phase, and the small active energy barrier between PIA-VO 2 (B) and crystalline VO 2 (B) enables feasible control of the phase switch; (2) The reversible PIA and recrystallization processes are guaranteed to occur between VO 2 (B) and PIA-VO 2 (B) (except some thermodynamically stable phases such as VO 2 (M)), due to the relatively low operating temperature (compression at room temperature and annealing at B200°C). The illustration that dynamically low temperature can hinder both the pressure-and temperature-induced traditional phase transitions indicates that The latter will be discussed in detail in the electron transport property section. Recent investigations by more powerful and precise techniques show that PIA is highly related to the pressurization environments and crystallinity of the starting materials. Some of the phenomena previously reported as PIA were likely due to either the formation of multiple polymorphic phases, or even that the PIA process did not occur if single crystal samples were adopted or more hydrostatic compression was applied [8][9][10]35 .
To check this issue in the VO 2 (B) system and obtain the exact conditions where the phase transformations occurred, in situ powder XRD experiments were conducted. Firstly, the PIA process of VO 2 (B) with different PTMs were evaluated. Figure 4a displays the integrated (110) peak intensity of VO 2 (B) as a function of pressure without a PTM or with Ne, He as the PTMs. VO 2 (B) nanosheets under all of these three conditions showed the onset of PIA around 10 GPa. Without PTM, the PIA process accomplished concluded around 20 GPa or a little higher when a better hydrostatic pressure condition was given (30 and 35 GPa for neon and helium as PTMs, respectively). This indicates that the PIA of VO 2 (B) is somewhat related to the deviatoric stress 36 , which is reasonable considering the hierarchical structural features of VO 2 (B). Fortunately, the value of 20 GPa was low enough to obtain bulk samples using a large volume press apparatus for routine magnetism measurements. Figure 4b shows the (110) peak intensity of PIA-VO 2 (B) as a function of the annealing temperature. The recrystallization starts from as low as 100°C and completes around 200°C. The relatively low annealing temperature makes the switch between PIA-VO 2 (B) and crystalline VO 2 (B) feasible. The pressure-and    temperature-induced phase transformations of VO 2 (B) nanosheets observed in our work are summarized schematically in Fig. 4c with other known VO 2 crystalline phases. Under compression at room temperature, the expected phase transformation from metastable VO 2 (B) to thermodynamically stable VO 2 (M) does not occur. This indicates a higher kinetic barrier (E barrier 2) of atomic diffusion, and the large surface energy from the nanosheets at room temperature. The formed high-density PIA-VO 2 (B) returned to the pristine metastable VO 2 (B) structure after short-time annealing at relatively low temperature (250-300°C, 10 min), instead of transforming to the thermodynamically stable VO 2 (M) phase. This indicates a relatively low kinetic barrier transforming to metastable VO 2 (B) compared with the high crystallization energy associated with the VO 2 (M) phase and demonstrate that hidden structural switching can be realized using a compression and heating route, as shown in Fig. 4c. Such an intriguing phase transition is the first ever example observed in a 3d metal involved material. Similar structure switching behaviour was observed in PbTe nanoparticles 37 , but in our case the switching between two thermally metastable phases was highly dependent on the temperature parameter. We believe that the key factor enabling this reversible phase transformation is the temperature, which is kinetically low enough for sufficient atomic mobility, which thus causes the elastic deformation of the structure. However, when the annealing temperature was high enough to surpass the VO 2 (B) to VO 2 (M) phase transition temperature (B550°C) 33 , the PIA-VO 2 (B) can also transfer to VO 2 (M) via the intermediate VO 2 (B) phase.
First-principles calculation of the pressure-phase relationship.
The experimental results evidently show the pressure-and thermal-controlled structure switching of VO 2 (B) nanosheets. To provide further thermodynamic and kinetic insight into the pressure-phase relationship during the phase transformation, we performed density functional theory (DFT) simulations using CASTEP code 38  under ambient pressure at room temperature, and that VO 2 (M) and VO 2 (R) have a higher density than VO 2 (A) and VO 2 (B) (Supplementary Fig. 5). Upon compression up to 100 GPa, VO 2 (B) was thermodynamically metastable compared with VO 2 (M) and VO 2 (R). However, no crystal-to-crystal phase transition from VO 2 (B) to VO 2 (M) or VO 2 (R) was observed, due to the high-kinetic energy barrier between them at room temperature. Figure 5 shows both the calculated pressurevolume equations of states of the four crystalline VO 2 phases and the experimental data of VO 2 (B) nanosheets upon compression. The experimental data of the VO 2 (B) nanosheets exhibit a larger unit cell volume than those theoretically predicted, which is reasonable considering the nanosheet surface defects. This large cell volume falls into the traditional VO 2 glass region. Upon compressing to around 20 GPa, the VO 2 (B) nanosheets passed through the ideal pressure-volume boundary of VO 2 (B) and became amorphous with a higher density. The stable region of PIA-VO 2 (B) is indicated in Fig. 5, which passed between VO 2 (B) and VO 2 (M) under high pressure. Several possibilities have been proposed for the underlying mechanism of a PIA process, such as the mechanical instability beyond the Born stability conditions and 'thermodynamic melting' 1,9 . In the case of VO 2 (B), we propose that the combination of kinetic hindrance of the phase transformation to thermodynamically stable VO 2 (M) and the increased atomic coordination are the driving forces of PIA. Upon moderate temperature annealing, the preference of VO x polyhedra in thermodynamically metastable PIA-VO 2 (B) drives the structure to return to the pristine VO 2 (B).
In the P-V phase diagram, one can clearly see the difference between PIA-VO 2 (B) and VO 2 glass. It is important to note that other hidden phases of VO 2 may exist, either in the low-density glass range or the high-density pressure-induced amorphous range. In this study, we observed profound structure switching behaviour of PIA-VO 2 (B) at kinetically moderate temperatures, which demonstrates an interesting phenomenon of how temperature has a critical role in governing phase transformations under high pressure. The exploration of yet more interesting physical properties related to high density VO 2 glass with a short-range order may greatly boost the development of condensed matter science.
Discussion VO 2 phases as the famous strongly correlated family always shows interesting magnetic and electric properties. To gain more insight into the participation of charge and spin interactions during the PIA, in situ electrical resistance within a diamond anvil cell (DAC) and ex situ magnetic susceptibility measurements (using an amorphized sample synthesized by a large volume press apparatus) were performed. Figure 6 shows the resistance change during compression and decompression, and those at 7.5 and 26 GPa (before and after amorphization) as a function of temperature. At the first stage (Po10 GPa), the resistance drop can be associated with the broadening and partial overlapping of the valence and conduction bands, caused by the normal pressure-induced shortening of V-V distances and bending of the V-O bonds. A steep increase in the resistance, which is caused by the onset and completion of the PIA, dominates the profile between 10 and 18 GPa. During the PIA process, point defects accompany the breaking of the long-range order, and thus the electron transport is suppressed by these increased defects, which serve as scattering centres. After the PIA process, band overlapping proceeds, as revealed by the continuous decrease of the electrical resistance, and finally drops to 30 GPa; nearly the same level as the crystalline VO 2 (B) before PIA. However, after decompression, the resistance increases by two orders of magnitude than the pristine VO 2 (B) sample (Fig. 6b), which indicates a reinforced electronic localization in PIA-VO 2 (B). The temperature dependence of the electric resistance before and after the PIA process (7.5 and 26.0 GPa) proves that PIA-VO 2 (B) remains a semiconductor but has poorer electrical conductivity (Fig. 6c). We also fabricated PIA-VO 2 (B) samples using a large volume press apparatus at room temperature and 20 GPa. The decrease of the effective Bohr magnetic moment in PIA-VO 2 (B) from 1.35m B to 1.15m B per V 4 þ , as derived from the ex situ magnetic susceptibility measurement, suggests that the electron localization may be due to the formation of more V-V pairs/dimers ( Supplementary  Fig. 6). More investigations are required to further probe the role of charge and spin in the PIA phenomenon and structure memory of VO 2 (B) nanosheets.
In conclusion, we report a precise control of the structure switching between crystalline and pressure-induced amorphous phases in VO 2 (B) nanosheets for the first time. The PIA-VO 2 (B) has no long-range periodicity, but preserves the inherited VO x polyhedral from VO 2 (B) with highly distorted local ordering. Upon moderate temperature annealing, these VO x polyhedra release excess energy and restore the pristine long-range periodicity. Spectral evidence and DFT calculations provide additional insight into the phase relationship between crystalline VO 2 phases and PIA-VO 2 (B) glass. The dynamically moderate annealing temperature provides a key mechanism of energetic control on the hidden reversible amorphization of metastable VO 2 (B). Preliminary evidence for the participation of charge and spin interactions during PIA were also observed for the first time.
The robust control of the phase transition between the metastable crystalline VO 2 (B) phase and the PIA-VO 2 (B) amorphous phase, via compression and low-temperature annealing, reveals the structure switching behaviour within a tetravalent vanadium-based material for the first time. This highlights exploration of thermodynamically hindered phase transformations with pressure and temperature tuning. Further investigations on the physical mechanism behind the phenomenon are expected, including orbital, electron and magnetic interactions.

Methods
Sample preparation. VO 2 (B) nanosheets were synthesized from high pure V 2 O 5 raw material via a hydrothermal route, with citric acid as the reducing agent. Briefly, 0.182 g V 2 O 5 powder (1 mmol) and 0.288 g citric acid (1.5 mmol) were added into 20 ml distilled water and continuously stirred to obtain the precursor. The precursor was then transferred into a Teflon-lined autoclave (25 ml capacity, 80% filling). The autoclave was heated to 200°C at a rate of 3°C min À 1 and maintained at 200°C for 5 h, followed by air cooling to room temperature by switching the power off. The resulting precipitates were washed with deionized water and dried at 80°C overnight. Transmission electron microscopy (TEM) techniques were applied to check the morphology of the as-synthesized nanosheets at ambient pressure.
In situ high-pressure characterizations. A symmetric DAC was employed to generate high pressure. A stainless steel gasket was pre-indented to a 40 mm thickness, followed by laser-drilling the central part to form a 120 mm diameter hole to serve as the sample chamber. Pre-compressed VO 2 (B) powder pellets and a small ruby ball were loaded in the chamber. Helium or neon was used as the pressure-transmitting medium and the pressures were determined by the ruby fluorescence method 39 . The in situ high-pressure angular-dispersive XRD experiments were carried out at the 16BM-D station of the High-Pressure Collaborative Access Team (HPCAT) at the Advanced Phonon Source (APS), Argonne National Laboratory (ANL). A focused monochromatic X-ray beam about 5 mm in diameter (FWHM) and with wavelengths of 0.4246 Å was used for the diffraction experiments. The diffraction data were recorded by a MAR345 image plate and processed with the Fit2D programme. High-pressure infrared spectroscopic experiments were performed at the U2A beamline of National Synchrotron Light Source (NSLS), Brookhaven National Laboratory and the mid-IR beamline of the Canadian Light Source. The infrared spectra were collected in transmission mode by a Bruker FTIR spectrometer using a nitrogen-cooled mid-band MCT (MCT-A) detector. The recorded frequencies were in the range of 600-8,000 cm À 1 with a resolution of 4 cm À 1 . High-pressure Raman spectra were measured by a Raman spectrometer with a 532.1 nm excitation laser at HPCAT. For the ex situ annealing experiments, the PIA-VO 2 (B) samples were removed with the gaskets after depressurization. They were then annealed at different temperatures (200, 250, 300, 350, 400, 450, 500, 550°C) each for 5 min within a vacuum furnace, and the resulting phases were checked with synchrotron XRD. While in the in situ annealing experiments, the bulky PIA-VO 2 (B) sample was made using a large volume press apparatus and sealed in a glass tube. Then XRD patterns were taken at 30, 60, 90, 120, 160, 190 and 220°C with electric heating, respectively. The holding time at each temperature was 5 min. Structure refinements were performed with the FULLPROF programme 40 .
Resistance and magnetic susceptibility measurements. In situ electrical resistance was measured by a four-probe resistance measurement system consisting of a Keithley 6,221 current source, a 2182A nanovoltmeter and a 7,001 voltage/current switch system. A DAC device was used to generate pressures up to 30 GPa, and a cubic boron nitride layer was inserted between the steel gasket and diamond anvil to provide electrical insulation between the electrical leads and gasket. Four gold wires were arranged to contact the sample in the chamber for resistance measurements. For the magnetism measurement, 0.1 g sample was made using a large volume press apparatus at 20 GPa for 1 h at room temperature. The DC magnetic susceptibility was measured using a SQUID magnetometer (Quantum Design) with an applied magnetic field of 500 Oe.
DFT calculations. Total energy and pressure-volume equations of state for the four VO 2 crystalline phases were calculated at the atomic level using the first-principles plane-wave pseudopotential method 41 based on the density functional theory, with the CASTEP package 38 . The exchange-correlation function is described by the local density approximation 42 . The ion-electron interactions were modelled with ultrasoft pseudopotentials 43 for all constituent elements, where the V 4s 2 3d 3 and O 2s 2 2p 4 electrons were treated as the valence electrons, respectively. The kinetic energy cut-off of 380 eV and Monkhorst-Pack k-point meshes 44 spanning o0.04 Å À 3 in the Brillouin zone were chosen. The starting structure models were obtained from the Inorganic Crystal Structure Database (ICSD) for VO 2 (A), VO 2 (B), VO 2 (M) and VO 2 (R). The cell parameters and atomic positions within the unit cell of these four phases, under hydrostatic pressure ranging between 0 and 100 GPa (with the interval of 5 GPa), were fully optimized using the quasi-Newton method 45 . The convergence thresholds between the optimization cycles for energy change, maximum force, maximum stress, and maximum displacement were set as 5.0 Â 10 -6 eV per atom, 0.01 eV Å À 1 , 0.02 GPa and 5.0 Â 10 -4 Å, respectively. The optimization terminates when all of these criteria are satisfied. All these computational parameters were tested to ensure the sufficient accuracy for the present purposes.
Data availability. The data that support the findings of this study are available from the corresponding author upon request.