DFT Study of WS2-Based Nanotubes Electronic Properties under Torsion Deformations

In this study, the influence of torsional deformations on the properties of chiral WS2-based nanotubes was investigated. All calculations presented in this study were performed using the density functional theory (DFT) and atomic gaussian type orbitals basis set. Nanotubes with chirality indices (8, 2), (12, 3), (24, 6) and (36, 9) corresponding to diameters of 10.68 Å, 14.90 Å, 28.26 Å and 41.90 Å, respectively, are examined. Our results reveal that for nanotubes with smaller diameters, the structure obtained through rolling from a slab is not optimal and undergoes spontaneous deformation. Furthermore, this study demonstrates that the nanotube torsion deformation leads to a reduction in the band gap. This observation suggests the potential for utilizing such torsional deformations to enhance the photocatalytic activity of the nanotubes.


Introduction
Tungsten disulfide (WS 2 ) is a promising material for nanotechnology due to its layer structure and hexagonal lattice, which allow the formation of nanosheets and nanotubes.The first WS 2 nanotubes were synthesized in 1992 by Tenne [1], and since then, extensive research has been conducted to explore their unique properties and applications.Numerous scientific articles are devoted to experimental studies on WS 2 nanotubes.Modern synthesis methods involve a heterogeneous reaction between powdered WO 3−x and a gaseous mixture of H 2 S and H 2 /N 2 at temperatures ranging from 800-950 • C. Using this method, nanotubes with lengths of 1-10 µm and diameters of 20-150 nm can be obtained [2].In 2017, Chen et al. [3] successfully measured the precise diameters and chirality indices by analyzing electron diffraction patterns of the corresponding nanotubes.
Moreover, WS 2 -based nanotubes have been comprehensively investigated theoretically using quantum mechanics, molecular dynamics, and molecular mechanics methods.In [4,5], the electronic structure of WS 2 single-walled (SWNT) and multi-walled nanotubes with various chirality indices was calculated using DFT-LCAO and DFT-LACW methods.This allows obtaining the desired value of the band gap by varying the diameter of NTs [4].It was also shown that dopants and structural defects can also influence the electronic properties of NT [5].Studies [6,7] indicated that one-dimensional transition metal structures, like MoSSe-WSSe double-walled Janus nanotubes and MoS 2 NT/PbS x Se 1−x quantum dots heterostructures, demonstrated the possibility of ultrafast charge transfer, which is remarkable for their application as solar cells.Alternatively, one can explore the mechanical deformation of nanotubes with fixed chirality, rather than studying nanotubes with different chirality indices.
The geometry of nanotubes can be completely characterized by chirality indices (n 1 , n 2 ).These indices provide information about the nanotubes' symmetry, length of the translation vector, number of atoms in the unit cell, and diameter [8][9][10].Nanotubes with (n, n) and (n, 0) indices are achiral, possessing inversion symmetry, while other nanotubes are chiral, exhibiting left and right screw axes [10].Different chirality indices correspond to nanotubes with distinct electronic properties, including the electronic band gap energy.
It is well-known that the electronic properties of nanotubes, such as the band gap, the positions of the valence band edge (E VB ) and conduction band edge (E CB ), depend on the nanotube's deformation.Several studies have been dedicated to the investigation of the axial deformation of MoS 2 and WS 2 nanotubes [11][12][13].These studies have demonstrated that MoS 2 nanotubes can be stretched up to 16%, and WS 2 nanotubes up to 12% before rupture.The electronic band structure of these nanotubes can be tuned by applying tensile stress, making them potential candidates for applications in nanoelectronics as switches and piezoresistive strain sensors.
Another type of deformation is torsional deformation, which has been studied theoretically in carbon and silicon nanotubes [14][15][16].In particular, D'yachkov investigated [13] carbon and silicon nanotubes using the linearized augmented cylindrical waves (LACW) method.He demonstrated that adjusting the torsion angle ω affects the band gap as well as the positions of the valence band top (E VB ) and conduction band bottom (E CB ).The dependence of the band structure on ω differs between chiral and achiral nanotubes.However, D'yachkov's study limited the range of torsion angle ω to [−2  ] degrees to prevent potential rupture of the nanotubes.It is important to note that within this approach, quantum chemical calculations are performed without geometry optimization, which can potentially impact the accuracy of the obtained results.
In [17], torsional deformations of MoS 2 nanotubes were investigated using molecular dynamics, and focusing solely on mechanical properties.In a subsequent work [18], ab initio DFT calculations were performed on TMD nanotubes with different metals, but only achiral nanotubes were studied, and only axial deformations were considered.Therefore, there is a need for detailed first-principle calculations and the characterization of chiral WS 2 -based NTs, taking into account the dependence of the electronic properties on the torsion angle.
Several studies have been conducted to experimentally investigate torsion deformations of WS 2 -based nanotubes.In [19,20], a method is described where a WS 2 nanotube is fixed at its edges, and an AFM (Atomic Force Microscopy) needle pushes a lever fixed at the middle of the NT, causing torsional deformation of the structure.This technique ensures accurate measurement of mechanical and electrical properties, such as the shear modulus and conductivity of the material.In addition, theoretical studies have been conducted using semi-empirical DFTB (Density Functional based Tight Binding) method to examine the band structure changes of achiral WS 2 -based nanotubes under torsional deformations [21].The results have demonstrated a correlation between experimental and theoretical data, indicating that conductivity increases as the torsion angle increases.However, since most synthesized nanotubes are chiral, it would be more accurate to study this type of nanotubes.
This study focused on the effect of torsional deformations on the properties of WS 2 nanotubes with chirality indices n(4, 1) where n = 2, 3, 6, 9 (see Figure 1).Torsion deformations were modeled using the line groups theory approach.By taking symmetry into account at each stage of the calculations, the complexity of the quantum chemical calculations was significantly reduced.We obtained dependencies of the relative energies, the electronic band gap energies, and the energies of the top of the valence band and the bottom of the conduction band on the torsion angle.This technique has already been used to study a variety of carbon nanostructures, including polytwistane [22], nanohelicenes [23], and polyacetylene [24].However, distinct from [17][18][19], ab initio methods are explored for the first time in this study to investigate the properties of nanotubes as a function of torsion angle.Currently, methods for determining the chirality indices of experimentally synthesized nanotubes are known [3,25].Furthermore, the recent work by An et al. [26] proposed a method for synthesizing nanotubes with a fixed chiral angle.However, multi-walled nanotubes with large diameters (approximately 100-200 Å) are typically synthesized.One approach to modeling such nanotubes is using molecular mechanics [27]; however, even this cannot provide an accurate representation of experimentally obtained nanotubes.Further research is required to develop ab initio simulation techniques capable of accurate reproducing nanotubes with diameters closely matching those observed in experiments.

Computational Details
Ab initio calculations in this work were performed with CRYSTAL17 [28,29] program using a hybrid version of DFT exchange-correlation functional HSE06.The relativistic effective core pseudopotentials and the corresponding basis sets were taken the same as in [4].To provide a balanced summation in both direct and reciprocal lattices, the reciprocal space integration was performed by sampling the Brillouin zone (BZ) with the 18 × 1 × 1 Monkhorst-Pack mesh [30], which results in a total of 10 k-points evenly distributed over BZ.Calculations were considered as converged only when the total energy differs by less than 10 −9 a.u. in two successive cycles of the self-consistent field procedure.Strict accuracy criteria (10 −8 10 −8 10 −8 10 −8 10 −18 ) were used in the calculations of the Coulomb and exchange series.These tolerances mean that during the direct lattice summations, the one-electron integrals and two-electron Coulomb integrals less than 10 −8 are estimated by the multipolar expansion, and two-electron exchange integrals less than 10 −16 are ignored.For more comprehensive information, please refer to the CRYSTAL17 User's Manual [31].Both the translational period and atomic positions have been optimized.The gradient threshold was set equal to 0.00005 Ha•Bohr −1 .Convergence criterion on the atomic displacement was set equal to 0.0001 Bohr.
By incorporating roto-helical symmetry throughout the calculation, both computational and memory costs were significantly reduced.To validate the computational scheme, atomic and electronic structure calculations were performed for 2H-WS2 and a single 0001 monolayer composed of three atomic planes (S-W-S).Experimental measurements revealed a fundamental band gap of 2.4 eV for the WS2(0001) nanosheet with 1 Currently, methods for determining the chirality indices of experimentally synthesized nanotubes are known [3,25].Furthermore, the recent work by An et al. [26] proposed a method for synthesizing nanotubes with a fixed chiral angle.However, multi-walled nanotubes with large diameters (approximately 100-200 Å) are typically synthesized.One approach to modeling such nanotubes is using molecular mechanics [27]; however, even this cannot provide an accurate representation of experimentally obtained nanotubes.Further research is required to develop ab initio simulation techniques capable of accurate reproducing nanotubes with diameters closely matching those observed in experiments.

Computational Details
Ab initio calculations in this work were performed with CRYSTAL17 [28,29] program using a hybrid version of DFT exchange-correlation functional HSE06.The relativistic effective core pseudopotentials and the corresponding basis sets were taken the same as in [4].To provide a balanced summation in both direct and reciprocal lattices, the reciprocal space integration was performed by sampling the Brillouin zone (BZ) with the 18 × 1 × 1 Monkhorst-Pack mesh [30], which results in a total of 10 k-points evenly distributed over BZ.Calculations were considered as converged only when the total energy differs by less than 10 −9 a.u. in two successive cycles of the self-consistent field procedure.Strict accuracy criteria (10 −8 10 −8 10 −8 10 −8 10 −18 ) were used in the calculations of the Coulomb and exchange series.These tolerances mean that during the direct lattice summations, the oneelectron integrals and two-electron Coulomb integrals less than 10 −8 are estimated by the multipolar expansion, and two-electron exchange integrals less than 10 −16 are ignored.For more comprehensive information, please refer to the CRYSTAL17 User's Manual [31].Both the translational period and atomic positions have been optimized.The gradient threshold was set equal to 0.00005 Ha•Bohr −1 .Convergence criterion on the atomic displacement was set equal to 0.0001 Bohr.
By incorporating roto-helical symmetry throughout the calculation, both computational and memory costs were significantly reduced.To validate the computational scheme, atomic and electronic structure calculations were performed for 2H-WS 2 and a single 0001 monolayer composed of three atomic planes (S-W-S).Experimental measurements revealed a fundamental band gap of 2.4 eV for the WS 2 (0001) nanosheet with 1 monolayer thickness when supported by Si [32] or α-Al 2 O 3 [33] substrates, which is considerably larger than the experimental band gap of bulk WS 2 (1.4

eV).
The results obtained by us from the quantum chemical simulations demonstrate good agreement with the experimental crystal lattice parameters and reasonable agreement with the electronic band gap, see Table 1.

Application of Line Group Theory to Nanotubes with Hexagonal Morphology
Let us consider a brief overview of the theoretical modeling of nanotubes with hexagonal morphology using the theory of line groups.For more comprehensive information, please refer to the monographs [10,36] and articles [8,9,37].Article [38] is devoted to the study of the symmetry of nanotubes based on transition metal dichalcogenides.
Structures of tubes with a hexagonal morphology, obtained by rolling a hexagonal Here, → T is a translation vector and → R is called a chiral vector, which becomes a circumference of nanotube.Furthermore, where n 1 and n 2 -are the chirality indices of the nanotube.This can also be rewritten as n( ñ1 , ñ2 ), where n is the greatest common divisor for n 1 and n 2 .The chirality indices also determine the components of the translation vector a 2 by the formula: Another important characteristic of nanotube structure is the spiral vector The indices h 1 and h 2 can be determined from the formula: Figure 2 shows the chirality vectors → R for nanotubes with the different chirality indices and hexagonal morphology.It is seen that the chirality vectors belong to one ray if ( ñ1 , ñ2 ) is fixed and n = 1, 2, 3, . . .monolayer thickness when supported by Si [32] or α-Al2O3 [33] substrates, which is considerably larger than the experimental band gap of bulk WS2 (1.4 eV).
The results obtained by us from the quantum chemical simulations demonstrate good agreement with the experimental crystal lattice parameters and reasonable agreement with the electronic band gap, see Table 1.

Application of Line Group Theory to Nanotubes with Hexagonal Morphology
Let us consider a brief overview of the theoretical modeling of nanotubes with hexagonal morphology using the theory of line groups.For more comprehensive information, please refer to the monographs [10,36] and articles [8,9,37].Article [38] is devoted to the study of the symmetry of nanotubes based on transition metal dichalcogenides.
Structures of tubes with a hexagonal morphology, obtained by rolling a hexagonal layer (a Here,  ⃗ is a translation vector and  ⃗ is called a chiral vector, which becomes a circumference of nanotube.Furthermore,  ⃗ =   ⃗ +   ⃗ where  and  -are the chirality indices of the nanotube.This can also be rewritten as ( ,  ), where  is the greatest common divisor for  and  .The chirality indices also determine the components of the translation vector  ⃗ =   ⃗ + t  ⃗ by the formula: Another important characteristic of nanotube structure is the spiral vector  ⃗ = ℎ  ⃗ + ℎ  ⃗ .The indices ℎ and ℎ can be determined from the formula: Figure 2 shows the chirality vectors  ⃗ for nanotubes with the different chirality indices and hexagonal morphology.It is seen that the chirality vectors belong to one ray if ( ,  ) is fixed and n = 1, 2, 3, … When folding a WS2 layer into a chiral nanotube, the layer's symmetry group is replaced by the line group  =   , which belongs to the first family of line groups. is When folding a WS 2 layer into a chiral nanotube, the layer's symmetry group is replaced by the line group L = T r q C n , which belongs to the first family of line groups.T r q is referred to as the generalized translation group and is represented by the operation Z = T C Q f , where T is nanotube translational group if Q is rational, i.e., Q = q r = n q r ≥ 1 represents the order of the helical axis (rotation angle of the helical axis ϕ = 360 • Q ), and f = t q is a partial translation (t = | → T |).C n represents the symmetry group of the monomer, consisting of 3•n atoms.This form of representing line group is known as the "polymer" or "helical" factorization.Conversion to the crystallographic notation is based on the calculation of the p index, required for this type of notation L = Lq p [36].All parameters of the symmetry group can be easily determined by knowing the vectors → R, → T , and Here, l is the positive integer or zero.
For instance, let us consider a nanotube with chirality indices (10,5).By finding the greatest common divisor (GCD) between n 1 and n 2 , one can simplify the notation to 5(2, 1), where n = 5.Using Formulas (1) and ( 2), one can determine the components of the vectors → T and With this information, one can calculate the parameters for the symmetry group of the selected nanotube using Formulas (3)- (6).It can be found that q = n ∼ q = 70, r = 9, Q = 70 9 , and p = 55.Therefore, the symmetry group of the (10,5) nanotube in the helical factorization is written as T 9  70 C 5 = ZP, where Z ≡ T 9 70 = T C 9 70 f , and P ≡ C 5 .

Torsion Deformation of NTs with Hexagonal Morphology
From a practical standpoint, it is convenient to introduce the variable Q = q r , which determines the parameters of the helical axis, and the angle ϕ = 360 • Q .The use of Q is advantageous because, for any set of tubes specified by the formula n( ñ1 , ñ2 ), the value of Q remains unchanged even if Q undergoes multiple changes.
In order to identify groups that describe deformed nanotube (NT) structures, it is necessary to find for each of them Qx = qx r x such that the torsion angle ω = ϕ x − ϕ falls within a certain range.This range is determined by the specific structure and composition of the nanotube.Table 2 presents the torsion angles and their corresponding symmetry groups, indicated in both crystallographic (Ln• qn• ∼ p ) and helical notations (T r q C n ).By utilizing the Formulas (3) and ( 4), the value of Q for the initial nanotube can be easily determined.With the aid of an enumeration method, it is possible to identify a set of q and r values for which the torsion angle ω falls within the specified range.It is important to impose a maximum limit on the q value to avoid excessive computational complexity.Increasing q < qmax leads to a corresponding increase in the number of formula units in the elementary cell according to the formula q = n q, thereby complicating the quantum chemical calculations.

Results of Calculations for Non-Deformed Nanotubes
In the present paper, we consider torsion deformation of WS 2 -based chiral nanotubes n(4, 1) for n = 2, 3, 6, 9.The results of our calculations are provided in Table 3.The crystallographic factorization Lq p of the corresponding line groups is provided in the second line of Table 3.

Bulk
Mono ( The input values of NT diameters D = a π n ñ2 1 + ñ1 ñ2 + ñ2 2 are increasing from 9.22 Å to 41.49 Å.It is seen that due to the geometry optimization the change of the diameter ∆D decreases with the input diameter increasing as for the monolayer, this change becomes formally zero.For all NTs in consideration, the translation vector → T length t is the same as it is defined by the Z subgroup of line group L in the polymer factorization L = ZP (P = C 2 , C 3 , C 6 , C 9 ; It is seen from Table 3 that the translation vector length t changes due to the geometry optimization decreasing with NT diameter increasing, and moves to value 0.01 Å-the difference of the bulk and monolayer hexagonal lattice vector length a.
Table 3 shows that the primitive unit cell energy E value per formula unit decreases with NT diameter increasing, thus moving to that for monolayer.
The formation E f orm and strain E str energies of NTs are defined by relations: where N, N bulk and N mono are the numbers of formula units in the primitive cell for nanotube, bulk crystal, and monolayer.The corresponding energies per primitive unit cell are E, E bulk and E mono .It is evident that the difference (E f orm − E str ) = 22.56 kJ/mol is the same for all the nanotubes as it is the monolayer formation energy from bulk crystal.Table 3 shows the forbidden energy gap (E gap ) dependence on the nanotube diameter.As the energy gap for bulk is smaller than that for the slab (1.68 eV instead 2.53 eV), the band gap enlarges as the diameter of the nanotube (NT) increases.Moreover, it shifts toward the value of the band gap observed in a monolayer configuration.In the latter two lines of Table 3, we provide the one-electron energies E v and E c for the top of the valence band and the bottom of the conduction band, respectively.It is known that visible-light-driven water-splitting catalysts should meet following energy requirements: E v < −5.67 eV and E c > −4.44 eV (values of oxidation and reduction potentials of water).Our results show that this condition is fulfilled only for NTs (24,6) and (36,9).This conclusion is in agreement with that made in [4]: NTs with a diameter larger than 19 Å show the potential ability to serve as photocatalysts for water-splitting reactions.

Energy Minima and Atomic Structures of Torsional Deformed Nanotubes
For each nanotube analyzed, the impact of torsional deformation was studied within the range of ω ∈ [−2.857 • , 3.025 • ].This specific torsion angle range was selected because it allows obtaining a single energy minimum for each nanotube in consideration.Further increasing the absolute value of torsion angle beyond this range would lead to an increase in the relative energy, indicating structural damage of the nanotube.The atomic structure of deformed NTs was optimized for each chosen torsion angle (see Table 2), and subsequent calculations were performed to determine NTs properties.
Figure 3 shows the dependence of the relative torsion energy per formula unit, which was calculated using the following formula: where relative energy difference, denoted as ∆E rel (ω), is calculated as the difference be- tween the total electronic energies of the original nanotube structure, E(0), and the deformed nanotube structure with a torsion angle ω, E(ω).The number of formula units in the original nanotube is represented in (9) by q 0 , while q ω corresponds to the number of formula units in the nanotube structure with the torsion angle ω.
Figure 3 shows the dependence ∆E rel on torsion angle.It is seen that when the nanotube diameter decreases, the position of the energy minimum deviates more significantly from ω = 0. Furthermore, the preferred torsion direction varies based on the nanotube's diameter.In the case of the smallest nanotube (8, 2), torsion towards positive ω is less favorable compared to larger diameter NTs.Conversely, within the range of negative ω, torsion of the largest nanotube is less preferable.
In Table 4, the characteristics of the energy minima obtained through ∆E rel (ω) inter- polation are presented.An observed trend indicates that with an increase in nanotube diameter, the band gap increases while the shear modulus decreases.Additionally, the energy minimum shifts towards that of the undeformed nanotube structure.

Figure 3.
The dependence of the Δ on the torsion angle.Δ of structures with irrational Q are obtained here by interpolation [23] using SciPy [39] implementation of Akima interpolation [40].Matplotlib [41] and Seaborn [42] were used for data visualizing.
In Table 4, the characteristics of the energy minima obtained through Δ (ω) interpolation are presented.An observed trend indicates that with an increase in nanotube diameter, the band gap increases while the shear modulus decreases.Additionally, the energy minimum shifts towards that of the undeformed nanotube structure.(36,9) 41.90 −0.053 −0.029 2.38 0.17 a -The diameter values are measured relate to the positions of the W atom centers.b -Electronic band gap; c -Shear module, calculation technique is described in works [22,23]; *-The values correspond to the results obtained from the quantum chemical calculation for the torsion angle that is closest to the energy minimum.
Figure 4 illustrates the dependence on the torsion angle of structural parameters of the investigated nanotubes.In this context, an increase in the parameter partial (Figure 4a) translation  = corresponds to the stretching of the nanotube.As the nanotube diameter increases, the partial translation demonstrates an upward trend and moves to a constant value.Overall, the dependence of the partial translation on torsion angle remains, approximately, the same across all nanotubes; however, for the smallest nanotube (8,2), the increase becomes slightly stronger when undergoing torsion within the range of positive torsion angles.[23] using SciPy [39] implementation of Akima interpolation [40].Matplotlib [41] and Seaborn [42] were used for data visualizing.[22,23]; *-The values correspond to the results obtained from the quantum chemical calculation for the torsion angle that is closest to the energy minimum.
Figure 4 illustrates the dependence on the torsion angle of structural parameters of the investigated nanotubes.In this context, an increase in the parameter partial (Figure 4a) translation f = t q corresponds to the stretching of the nanotube.As the nanotube diameter increases, the partial translation demonstrates an upward trend and moves to a constant value.Overall, the dependence of the partial translation on torsion angle remains, approximately, the same across all nanotubes; however, for the smallest nanotube (8, 2), the increase becomes slightly stronger when undergoing torsion within the range of positive torsion angles.
It is evident from Figure 4b that the diameter of the examined nanotubes displays slight fluctuations within the torsion angle range of [−2 • , 2 • ].Across all of the nanotubes assessed, the most significant deformation occurs at negative torsion angles, particularly with the nanotube (36,9) displaying the most significant diameter changes.A comparable pattern is noticeable when deforming moves into the positive torsion angle region, however, the change in diameter is notably smaller.It is evident from Figure 4b that the diameter of the examined nanotubes displays slight fluctuations within the torsion angle range of −2 ∘ , 2 ∘ .Across all of the nanotubes assessed, the most significant deformation occurs at negative torsion angles, particularly with the nanotube (36,9) displaying the most significant diameter changes.A comparable pattern is noticeable when deforming moves into the positive torsion angle region, however, the change in diameter is notably smaller.
Notably, it can be observed that the physical properties of nanotubes undergo more essential changes upon torsion as the initial nanotube diameter increases.The asymmetric band gap values relative to the torsion angle ω = 0 are a consequence of the chiral nature of the nanotubes under investigation.
Figure 5 shows that the band gap exhibits slight variations for nanotubes with the smallest diameter (8, 2) and (12, 3) within the range of torsional deformation, with changes of no more than 0.15 eV and 0.25 eV, respectively.
For the nanotubes (24,6) and (36,9), the band gap decreases with torsional deformation, a trend that aligns with the findings of previous studies [19,21], demonstrating an increase in the electrical conductivity of WS2-based nanotubes due to torsion-induced deformation.While the observed trends in band gap for (24,6) and (36,9) are comparable, it is reasonable to hypothesize that nanotubes (NTs) with larger diameters exhibit similar tendencies in electronic properties.
Figure 6 illustrates the dependencies of the energies at the top of the valence band and the bottom of the conduction band for the nanotubes under investigation.It is evident that nanotubes (24,6) and (36,9) show potential applicability for water splitting, aligning with the established notion that nanotubes with smaller diameters are less suitable for this purpose.Notably, it can be observed that the physical properties of nanotubes undergo more essential changes upon torsion as the initial nanotube diameter increases.The asymmetric band gap values relative to the torsion angle ω = 0 are a consequence of the chiral nature of the nanotubes under investigation.
Figure 5 shows that the band gap exhibits slight variations for nanotubes with the smallest diameter (8, 2) and (12, 3) within the range of torsional deformation, with changes of no more than 0.15 eV and 0.25 eV, respectively.For the nanotubes (24,6) and (36,9), the band gap decreases with torsional deformation, a trend that aligns with the findings of previous studies [19,21], demonstrating an increase in the electrical conductivity of WS 2 -based nanotubes due to torsion-induced deformation.While the observed trends in band gap for (24,6) and (36,9) are comparable, it is reasonable to hypothesize that nanotubes (NTs) with larger diameters exhibit similar tendencies in electronic properties.
Figure 6 illustrates the dependencies of the energies at the top of the valence band and the bottom of the conduction band for the nanotubes under investigation.It is evident that nanotubes (24,6) and (36,9) show potential applicability for water splitting, aligning with the established notion that nanotubes with smaller diameters are less suitable for this purpose.Furthermore, it is apparent that as the diameter increases, the contribution of the top of the valence band to the change in E gap diminishes.Additionally, significant torsional deformation results in E c becoming lower than the reduction potentials of water.
Table 5 provides data on the behavior of the band gap under torsional deformations.Notably, for nanotube (8,2), there is a transition from an indirect to a direct band gap in one case, whereas no similar transition occurs for the other nanotubes examined.Currently, we are unable to determine if this change is specific to nanotubes with small diameters or if it can be observed with stronger torsion.Nevertheless, this phenomenon warrants further investigation.
Based on the obtained data, it is evident that the nanotubes under investigation are not suitable for photocatalytic water-splitting, as they exhibit indirect bandgap semiconductor behavior.

Conclusions
In this study, we investigated the impact of torsional deformation on the properties of WS 2 -based chiral nanotubes using first principles analysis.We utilized line groups theory to reduce computational costs and applied density functional theory (DFT) for quantum chemical calculations.
Our findings reveal that nanotubes with smaller diameters exhibit a more significant deviation in the position of the energy minimum on the torsion curve.In this context, the smaller the diameter, the more substantial the deviation.A decrease in diameter also causes diminishing of relative diameter changes with partial translation increases, aligning with the nanotube's stretching.
Additionally, our study shows that torsion has a minor impact on nanotube diameter within the range of −2 • to 2 • , but leads to more significant changes under stronger torsion, especially in the case of the nanotube (36,9).
We observed that torsional deformation has a limited effect on the band gap in nanotubes with small diameters but becomes more noticeable as the nanotube diameter increases, aligning with previous research [21].This property holds promise for applications in nanoelectronics based on WS 2 nanotubes.
The study showed that the studied nanotubes cannot be used for photocatalytic watersplitting due to the indirect band gap.However, in the case of tube (8,2), strong torsional deformation changed the band gap from indirect to direct, a phenomenon that deserves further study.Plots of band zones can be found in the supporting information, specifically in Figures S1-S4.

Figure 2 .
Figure 2. Chirality vectors of the nanotubes rolled up from a layer with hexagonal morphology.

Figure 2 .
Figure 2. Chirality vectors of the nanotubes rolled up from a layer with hexagonal morphology.

Figure 3 .
Figure 3.The dependence of the ∆E rel on the torsion angle.∆E rel of structures with irrational Q are obtained here by interpolation[23] using SciPy[39] implementation of Akima interpolation[40].Matplotlib [41]  and Seaborn[42] were used for data visualizing.

Figure 4 .
Figure 4.The torsion angle dependence on the values of (a) partial translation, (b) diameter changes.

Figure 4 .
Figure 4.The torsion angle dependence on the values of (a) partial translation, (b) diameter changes.

Nanomaterials 2023 , 13 Figure 5 .
Figure 5.The dependencies of the values of band gap as the functions of torsion angle.Figure 5.The dependencies of the values of band gap as the functions of torsion angle.

Figure 5 .
Figure 5.The dependencies of the values of band gap as the functions of torsion angle.Figure 5.The dependencies of the values of band gap as the functions of torsion angle.

Figure 5 .
The dependencies of the values of band gap as the functions of torsion angle.

Figure 6 .
Figure 6.Energies of the top of the valence band and the bottom of the conduction band as the functions of torsion angle.

Figure 6 .
Figure 6.Energies of the top of the valence band and the bottom of the conduction band as the functions of torsion angle.

Table 1 .
Comparison of calculated properties of bulk 2H-WS 2 and monolayer with experimental data.

Table 1 .
Comparison of calculated properties of bulk 2H-WS2 and monolayer with experimental data.

Table 2 .
The torsion angles for the nanotubes under investigation.The corresponding symmetry groups are provided for each torsion angle.

Table 4 .
Properties of NTs structures corresponding to obtained energy minima.

Table 4 .
Properties of NTs structures corresponding to obtained energy minima.The diameter values are measured relate to the positions of the W atom centers.
ab -Electronic band gap; c -Shear module, calculation technique is described in works

Table 5 .
Band gap behavior of investigated nanotubes under torsional deformation.
*-Cannot be obtained due to technical reasons.