{'fraction_non_alphanumeric': 0.02804084276768362, 'fraction_numerical': 0.01887898324952957, 'mean_word_length': 4.959191176470588, '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}}
Characterization of thermal and mechanical properties of stanene nanoribbons: a molecular dynamics study

Stanene, a buckled honeycomb structure of monolayer tin, has several intriguing electrical and thermoelectrical applications that closely depend on its thermal, mechanical, and electrical properties. However, thermal and mechanical characterizations of stanene nanoribbons (STNRs) have not yet been comprehensively investigated. In this study, we have performed an equilibrium molecular dynamics simulation to characterize the thermal and mechanical properties of STNRs using the modified embedded-atom method potential. The room temperature thermal conductivities of pristine 10 nm 3 nm zigzag and armchair stanene nanoribbon were estimated to be 0.95 0.024 W m 1 K 1 and 0.89 0.026 W m 1 K , respectively. We also studied the thermal conductivity as a function of temperature and width of the ribbon. The thermal conductivity was found to decrease with increasing temperature, whereas it tends to increase with increasing width for both configurations. In all cases, the zigzag STNR exhibited a higher thermal conductivity than its armchair counterpart did. Furthermore, our study includes an investigation of the thermal transport in defected STNRs. For a defect concentration of 1.5%, the thermal conductivity of defected stanene nanoribbon experiences a reduction of approximately 30–50%, whereas a 70–90% reduction was observed at a vacancy concentration of 5% for various types of defects. Finally, the stress–strain behavior of STNRs with varying width was analyzed using uniaxial loading. Zigzag STNRs were found to have higher fracture strength than their armchair counterparts. Moreover, with increasing width, both fracture strain and fracture stress of armchair STNRs were found to show small variations compared with their zigzag counterparts. This study provides insights for tuning the thermo-mechanical characteristics of stanene-based nanostructures for thermal management and possible applications as thermoelectrics.


Introduction
Following the recent technological advancements 1 triggered by graphene 1-3 because of its combination of electrical properties, 4,5 high thermal conductivity 6,7 and mechanical strength, 8 there has been ongoing increasing interest in the research of two dimensional (2D) nanomaterials. Apart from graphene, researchers have performed experiments with other 2D materials, such as silicene, 9 germanene, 10 hexagonal boron-nitride (h-BN) 11 and molybdenum disulde (MoS 2 ), 12 due to their novel properties [13][14][15] and broad application prospects that range from next generation energy efficient nanoelectronic devices 16 to aerospace structures.
Recently, a new 2D group IV material, stanene, which is a buckled honeycomb structure of a 2D hexagonal tin lm, has incited noteworthy interest due to its intriguing prospect as a quantum Hall insulator, 17,18 topological insulator 19 and topological superconductor. 20 Additionally, the spin orbiting coupling (SOC) 19 -induced bandgap of $0.1 eV for free-standing stanene makes it a promising candidate for nanoelectronic applications. 18,21 Moreover, heterobilayers such as stanene/graphene 22 and stanene/MoS 2 (ref. 23) are reported to show nite band gaps with high carrier mobility which would further promote the development of Sn-based nanoelectronics. In addition, the stable low buckled form of stanene, 21,24,25 which was synthesized on a Bi 2 Te 3 (111) substrate by Zhu et al., 26 was theoretically investigated and found to support a large-gap 2D quantum spin Hall (QSH) state at room temperature, thereby enabling superior electric conduction with zero dissipation. 18,27 Such 2D QSH systems, such as stanene and stanene with dumbbell units, 19,28 can lead to a signicant improvement in the thermoelectric gure of merit (zT). 29,30 Further improvement in the zT of stanene can be achieved by proper optimization and tuning of the sample size, 31 which could enhance the contribution of the gapless edge states, thereby eventually optimizing the electron transport in stanene.
The high thermoelectric performance of a system depends simultaneously on the good electron conductivity and low transport of phonons. 32,33 Hence, apart from the electrical properties, the prospect of stanene in thermoelectric applications urges the optimization of the thermal transport characteristics of stanene and its nanostructures. However, there have been markedly fewer investigations of the thermal transport characteristics of stanene and its nanostructures compared with the signicant research concentration on its electrical properties. By combining rst-principle calculations and the phonon Boltzmann transport equation, Peng et al. 34 calculated the lattice thermal conductivity of stanene and obtained a much lower thermal conductivity in stanene lattice in comparison with other 2D materials. The results of Peng et al. 34 predict that stanene is a suitable candidate for next-generation thermoelectric devices with its high thermoelectric efficiency. Nissimagoudar et al. 35 studied the diffusive nature of the thermal conductivity of a stanene sheet using a phonon Boltzmann transport formalism and density functional theory calculations. For a stanene sheet with a sample size of 1 mm Â 1 mm, the calculated lattice thermal conductivities along the zigzag and armchair directions are 10.83 W m À1 K À1 and 9.2 W m À1 K À1 respectively, at room temperature. 35 Recently, Cherukara et al. 36 proposed a bond order potential for stanene sheet and found the thermal conductivity to be substantially lower than graphene suggesting that stanene can be a promising thermoelectric material. They further expected that, their results would encourage further characterization of stanene, such as the mechanical and thermal properties of stanene ribbons and nanotubes. However, the thermal transport characterization of either pristine or defected stanene nanoribbons (STNRs) has yet to be explored.
On the other hand, the electronic properties of a QSH insulator have been reported to change under deformation and mechanical loading in gallium arsenide (GaAs) lms. 37 In this context, Mojumder et al. 38 investigated the effect of the strain rate and variation of temperature under uniaxial and biaxial loading in both the armchair and zigzag direction of a 20.5 nm Â 20.5 nm stanene sheet using molecular dynamics (MD) simulations. However, further investigation is required to study the effect of size on the fracture strain and strength under uniaxial tensile loading of STNRs for proper characterization of the mechanical properties of stanene nanoribbons.
Hence, in this study, we performed molecular dynamics simulation using modied embedded-atom method (MEAM) potential 39 parameters to characterize the thermal and mechanical properties of stanene nanoribbons (STNRs). During recent decades, the MEAM potential has been well established for studying material properties, 40,41 such as thermal conductivity, thermal resistance 40,42 and mechanical stress-strain. 38,43 In fact, using MEAM potential parameters, our obtained equilibrium buckling height (0.88Å) and Sn-Sn bond length (2.83Å) of stanene match reasonably well with the values obtained experimentally as well as with prior density functional theory (DFT) simulations. 44 Using an equilibrium molecular dynamics (EMD) simulation, we computed the thermal conductivity of both armchair and zigzag-congured $10 nm Â 3 nm STNRs with a variation of temperature in the range of 100 K to 600 K. To study the size dependence of thermal conductivity in STNRs, we also computed the room temperature thermal conductivity of STNRs with varying widths in the range of $2 nm to 6 nm while the length was kept xed at 10 nm. Next, considering the inevitable presence of some form of structural defects during the synthesis, integration and processing of nanoribbons, we studied the impact of three types of vacancies, point vacancy, edge vacancy and bi-vacancy, on the thermal conductivity of STNRs. Finally, we report the stress-strain behavior of both armchair-and zigzag-congured STNRs with increasing widths under uniaxial tensile loading.

Methods
In this study, we conducted molecular dynamics simulation using Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) 45 to investigate the thermal and mechanical properties of stanene nanoribbons. The Modied Embedded-Atom Method (MEAM) potential parameters for stanene 38,39 were employed to perform the MD simulations.
To compute the thermal conductivity of both armchair and zigzag STNRs, an EMD simulation based on the linear response theorem was used in this study. In this method, the thermal conductivity computation is related to the ensemble average of the heat current auto correlation function (HCACF) by following the Green-Kubo formulation: here, K x is the thermal conductivity in the x direction, K B is the Boltzmann constant, T is the system temperature, V is the system volume dened as the area of STNR multiplied by the van der Waals thickness including the buckling height (4.5Å), 36 and s is the correlation time required for the reasonable decay of the HCACF. hJ x (t)$J x (0)i is the ensemble averaging term, which can be calculated from the heat current equation: where 3 i is the total site energy, v i is the time-dependent velocity of atom i, r ij ¼ r i À r j (r i is the time-dependent position of atom i), and F ij is the force exerted by atom j on atom i. A periodic boundary condition along the zigzag and armchair direction was applied. A xed time step of 0.5 fs was used throughout all of the simulations. The velocity-Verlet integrator 46 was used to integrate Newton's classical equations of motion. The structures were equilibrated using a Nose-Hoover thermostat for 3 Â 10 5 time steps followed by constant volume, constant energy ensemble (NVE) for 10 5 time steps, whereas the energy minimization was performed using the steepest descent algorithm. In-plane heat current data were recorded every 5 steps by performing constant energy simulations, and heat ux autocorrelation values were calculated by averaging 10 obtained HCACFs. Then, a 10 ps correlation time was used for all simulations to allow the HCACFs to decay to zero. Finally, using eqn (1), a converged value of the average thermal conductivity was calculated using the computed thermal conductivity values from 5 independent microcanonical ensembles (constant number of particles, volume and energy), i.e., NVE ensembles, each starting with a different initial velocity of the molecules.
For the calculation of the phonon density of states (PDOS), the FixPhonon command 47 of LAMMPS was employed to directly obtain the dynamical matrices from the MD simulation based on uctuation-dissipation theory. Using dynamical matrices, an auxiliary post-processing code, 'phana', was then used for the evaluation of the PDOS of the STNRs. In this study, alongside the uniform generation of q (wave vector) points, we used a tricubic 48 interpolation method to compute the phonon density of states for both pristine and defected STNRs.
For the specic characterization of the mechanical properties to observe the strain-stress behavior of the stanene nanoribbon, we performed MD simulations in LAMMPS using MEAM potential parameters. 38,39 Equations of motions were integrated with a time step of 1 fs. The steepest descent algorithm was used to relax the initial geometry. The initial structure was equilibrated using a constant pressure and constant temperature (NPT) ensemble for 30 ps followed by an NVE ensemble for 10 ps to stabilize the system temperature and pressure. Then, a uniaxial tensile loading was applied along the zigzag or armchair direction at a constant strain rate of 10 8 s À1 and at a temperature of 100 K. Stresses in LAMMPS were calculated using the following formula: 38,49 where s is the virial stress, V is the volume of the stanene nanoribbon, and the summation is over all the atoms occupying the total volume. m i is the mass of atom i, ui is the velocity of the atoms, r ij is the position vector of atom, 5 is the cross product, and f ij is the interatomic force applied on atom i by atom j. The thickness of the stanene nanoribbon was taken to be 4.5Å. To minimize the stress uctuations caused by thermal uctuations, the atomic stresses were averaged over a time interval of every 100 time steps. The electronic properties of the proposed heterobilayer structure are investigated using the density-functional theory (DFT) with a plane-wave basis set using the Ab initio code PWSCF package of Quantum Espresso. 50 The electron-ion interactions are accounted using Troullier Martin pseudopotentials. To describe the electron exchange correlation energy, the generalized gradient approximation (GGA) with Perdew-Burke-Ernzerhof (PBE) exchange correlation functional is implemented. For the structural optimization, the plane wave basis cut off is set at 550 eV (40 Ry) with a convergence threshold on force of 10 À4 Ry a.u À1 . The rst Brillouine zone of the unit cell is sampled with a 12 Â 12 Â 1 Monkhorst-Pack grid for the geometry optimization and a 15 Â 15 Â 1 grid for the subsequent calculations. Fig. 1 is a schematic representation of the atomic structure of a low buckled phase stanene nanoribbon with a Sn-Sn bond length of 2.83Å. Structural optimization i.e., energy minimization, was performed using the steepest descent algorithm. The buckling height of the relaxed structure is 0.88Å. This hexagonal structure with low buckling is geometrically similar to silicene, germanene and blue phosphorene 10,51-53 but different from the planar geometric conguration of graphene nanoribbons (GNRs) and hexagonal boron nitride nanoribbons (h-BNNRs). 54, 55 We considered 10 nm Â 3 nm armchair stanene nanoribbon structures with 3 types of vacancies as represented in Fig. 2. The structure with point vacancy results from the random removal of a single atom from the lattice is depicted in Fig. 2(a). As shown in Fig. 2(b), edge vacancy can be considered as a special form of point vacancy and involves removing atoms solely from  the lattice boundary. Fig. 2(c) presents the stanene nanoribbon structure with bi-vacancy, which originates from either the lumping of two consecutive point vacancies or from the random elimination of a pair of bonding atoms.

Results and discussions
Our computed average thermal conductivity for a 10 nm Â 3 nm zigzag stanene nanoribbon using MEAM potential parameters at room temperature (300 K) is 0.95 AE 0.024 W m À1 K À1 , whereas the estimated value of the average thermal conductivity is 0.89 AE 0.026 W m À1 K À1 for armchair stanene nanoribbon as shown in Fig. 3. Cherukura et al. 36 reported a room temperature thermal conductivity value of $3 W m À1 K À1 for both zigzag and armchair congured 80 nm Â 80 nm stanene sheets, and they expected these values to decrease even further with nanostructuring, which is in accordance with our reported result for a stanene nanoribbon in both congurations.
Additionally, using the phonon Boltzmann transport formalism and density functional theory calculations, Nissimagoudar et al. 35 reported a thermal conductivity value of 1.6 W m À1 K À1 for a 100 nm Â 100 nm stanene sheet. The authors 35 also mentioned the possibility of a further reduction in the thermal conductivity of stanene with a reduction in the dimensionality of the sample, i.e., nanoribbons, which is in line with the computed value of the thermal conductivity in our study. The calculated values of thermal conductivity for the stanene nanoribbon in this study are signicantly smaller compared with other 2D hexagonal materials and nanostructures such as GNRs ($3000 W m À1 K À1 ) 56,57 and h-BNNRs ($400-600 W m À1 K À1 ). 55,58 This can rst be attributed to the lower Debye temperature of $72 K for stanene 34   higher phonon scattering rates, 34 lower the thermal conductivity of STNRs. Fig. 3 shows that with an increase of temperature from 100 K to 600 K, the thermal conductivities of the STNRs decays for both armchair and zigzag congurations at slightly different rates. This decaying trend above the Debye temperature ($72 K for stanene) is similar to the study of Peng et al. 34 for a stanene lattice structure, which used the single mode relaxation approximation method and to a molecular dynamics study by Cherukara et al., 36 which used an optimized Tersoff potential for 80 nm Â 80 nm stanene sheets. This is also in agreement with the reported trend for other 2D hexagonal materials, such as graphene nanoribbon, 54,56 BNNRs 55,58 and monolayer MoS 2 . 63 Moreover, as shown in Fig. 3, at lower temperatures, zigzag nanoribbons exhibit a higher thermal conductivity than the armchair congurations, which is in agreement with the study of stanene nanotubes 36 as well as other 2D nanostructures, such as GNRs, 64 carbon nanotubes, 65 and h-BNNRs. 58 However, at an elevated temperature, the armchair and zigzag congurations exhibit a very small anharmonicity, i.e., a small difference in thermal conductivity values. This is due to the comparatively higher group velocities of the armchair-congured stanene nanoribbon compared with the zigzag conguration at higher temperature. 35    . 4(a) shows the total energies during the simulation for our considered temperature range. As shown in the gure, the energy uctuations were negligible, thereby reecting the stability of the stanene nanoribbon structure. The dotted lines for each temperature further depict the convergence of the total energies for different temperatures. Fig. 4(b) depicts the reasonable decay of the heat current autocorrelation function required for the computation of the thermal conductivity using the Green-Kubo method. Furthermore, Fig. 4(b) can be taken into consideration to explain the decrease in thermal conductivity with the increasing temperature for the STNRs. With increasing temperature, phonon-phonon scattering, specically Umklapp scattering, increases due to the increasing number of phonons. [34][35][36] As a result, with the increase in temperature, the HCACF prole decays to zero within a quicker time, 66 which results in the computation of a decreasing thermal conductivity for STNRs.
The inuence of width of the sample on the thermal conductivity of zigzag and armchair STNRs is shown in Fig. 5. A reasonable decay for HCACF was maintained, which ensured the good convergence of the Green-Kubo integral as shown in the inset of Fig. 5(a). As depicted in Fig. 5(a), for both STNR congurations, the room temperature thermal conductivity increased with an increase in width ranging from $2 nm to $6 nm. This increasing trend for the thermal conductivity with width is in agreement with the study by Peng et al. 34 for stanene nanowires. Furthermore, our study shows that thermal conductivity increases rapidly initially with increasing width, and further increases in width results in smaller increments of thermal conductivity. Two-dimensional hexagonal GNRs and h-BNNRs have been reported to show a similar variation in the literature. 54,57,58 The width dependence of thermal conductivity is inuenced by both the edge localized phonon boundary scattering effect and anharmonic phonon-phonon scattering effects. Each of these has negative effects on the thermal conductivity value of stanene. For the smaller sample sizes considered in this study, phonon boundary scattering mostly dominates the thermal conductivity of STNRs. 35 In this case, with an increasing width, the edge-localized phonon scattering effect becomes suppressed, resulting in an increased thermal conductivity. 35,57,67 In fact, for a 2D stanene sheet, using phonon Boltzmann transport  formalism, Nissimagoudar et al. 35 have reported that, the inuence of phonon boundary scattering increases with the decrease in width, which reduces the thermal conductivity. This is in good agreement with our study on the width dependence of thermal conductivity in stanene nanoribbons.
To provide further insight into the width dependence of the thermal conductivity of STNRs, Fig. 5(b) can be taken into consideration. The phonon density of states for varying widths of STNRs is shown in Fig. 5(b) with a xed length of 10 nm. As shown in the gure, the magnitude of the dominant peaks is higher in case of 5 nm width. In fact, more low frequency i.e. long wavelength acoustic phonons are available with the increasing width for STNRs. Recently, Peng et al. 34 showed that because of higher group velocities and a lower scattering rate, low frequency phonons provide the major contribution to the thermal conductivity of stanene. As a result, with the increase in width, thermal conductivity increases.
Furthermore, the phonon density of states (PDOS) of stanene is limited in the low frequency region of the phonon spectrum, whereas the range of PDOS for graphene is up to 50 to 60 THz. 66 In fact, the presence of high-frequency peaks ($50 THz) in the phonon power spectrum of graphene as opposed to the lowfrequency peaks in stanene provides further insight into the signicantly smaller thermal conductivity of STNRs compared with to GNRs. 68 The experimental synthesis, integration and fabrication process of stanene is inevitably associated with some forms of structural defects and vacancies. Such defects have a signicant impact on the chemical, physical and thermal properties of materials. 69 Thus, from a practical perspective, a study on how these defects inuence the thermal conductivity of STNRs is essential for its proper thermal transport characterization. Moreover, it would enhance the possibility of further tuning the thermal conductivity of STNRs for a better thermoelectric performance. Hence, from this point of view, by varying the defect concentrations, we next performed EMD simulations using the MEAM potential to investigate the thermal conductivity of armchair STNRs with different types of vacancies, namely, point vacancy, edge vacancy and bi-vacancy. In this study, defect or vacancy concentrations are dened as the percentage ratio of the number of random vacancies to the total number of tin atoms considered. Fig. 6(a) suggests that, for all types of vacancies, the thermal conductivity of stanene nanoribbon decreases monotonically with an increase in the defect percentage. As observed from the gure, edge vacancy and point vacancies impose a signicant reduction in thermal conductivity, whereas bi-vacancy was found to be less severe in terms of reduction in thermal conductivity. For GNRs and h-BNNRs, a similar exponential decay of the thermal conductivity in the presence of various types of defects has been reported in the literature. 66,[70][71][72] Again, for all three types of vacancies, the thermal conductivity prole shows a fast decrease at low defect concentrations (0.5% to 3%); then, the decay rate ceases for higher defect concentrations up to 5%.
The percentage reduction in thermal conductivity of STNRs with different types of vacancies compared with the same-sized pristine stanene nanoribbon is presented in Fig. 6(b). For a defect concentration of $1.5%, the thermal conductivity of defected stanene nanoribbon experiences a reduction of approximately 30-50%, whereas a $70-90% reduction was observed at a vacancy concentration of $5% for various types of defects.
The PDOS shown in the inset of Fig. 6(b) can be taken into consideration for the reduction of the thermal conductivity in the presence of a vacancy in STNRs. The peaks in the phonon spectrum in the low frequency region are damped out in the case of STNRs with vacancy, which is in accordance with the reported literature for defected graphene and GNRs with vacancies. 66,[70][71][72] These low frequency phonons are found to be the main contributor to the thermal conductivity of stanene. 34 Hence, the soening of phonon peaks in the low frequency region for the defected structure as opposed to the pristine stanene nanoribbon can be attributed to the observed percentage decrease in thermal conductivity in Fig. 6(b). Furthermore, signicant phonon scattering at and around vacancy centers causes a reduction in the phonon mean free path, 66,73 thereby reducing the thermal conductivity of STNRs with several types of vacancies. A similar phenomenon is expected to be observed in the zigzag conguration with vacancies.
The low thermal conductivity of pristine STNRs and further lower thermal conductivity of the defected STNRs with several types of vacancies enhance the applicability of STNRs as thermoelectrics. However, the high thermoelectric performance of a system depends simultaneously on the low transport of Fig. 9 Fracture pattern of a $10 nm Â 6 nm zigzag stanene nanoribbon under uniaxial tension at a strain rate of 10 8 s À1 and at 100 K. The atomic stress in GPa is shown in the color bar.
phonons and the good electron conductivity. Hence, using DFT, we have computed the electronic band structures of pristine stanene and stanene with a point vacancy as shown in Fig. 7(a) and (b), respectively. As can be seen from Fig. 7, in the presence of spin orbital coupling, a nite band gap is opened in the band structures for both the pristine and defected stanene structures while the linear Dirac dispersion relation is nearly preserved at the K point. Using the conduction band in the obtained band structure, the electron effective mass can be calculated from the following formula: where, m* is the particle effective mass, E(k) is the dispersion relation, k is the wave vector and ħ is the reduced Planck constant.
The calculated electron effective mass for pristine and defected stanene are $0.49 Â 10 6 ms À1 and $0.2 Â 10 6 ms À1 respectively. Here, m o is the mass of an electron. The small electron effective mass and the high Fermi velocity indicate the high carrier mobility in our studied structure, thereby ensuring its applicability as thermoelectrics. Our computed electron transport characteristics of pristine as well as defected stanene are in good agreement with the reported literature. 21,74 In this section, using a molecular dynamics simulation with MEAM potential parameters, we present a comprehensive study of the size and chirality effects on the mechanical properties of STNRs under uniaxial tensile strain. Nominal strain and virial stress were used in this study to analyze the strain stress behavior. Fig. 8(a) represents the stress-strain curve under uniaxial tensile loading in the zigzag direction for a 20 nm Â 20 nm stanene sheet (shown at inset) at 100 K with a strain rate of 10 8 s À1 . As shown in the gure, a maximum stress of $3 GPa was found at a fracture strain of $24%, which is in good agreement with the study by Mojumder et al. 38 Stress-strain curves of $10 nm Â 6 nm STNRs under zigzag and armchair uniaxial tensile loading are shown in Fig. 8(b). The zigzag stanene nanoribbon has a higher stress as well as a higher fracture strain compared with the armchair stanene nanoribbon with a sudden drop in the stress prole aer a certain stress for both congurations. This trend is in agreement with the study of a stanene sheet by Mojumder et al. 38 as well as studies on GNRs. 75,76 For $10 nm Â 6 nm STNRs, our computed value for the fracture strain along the armchair direction is $16%, whereas along the zigzag direction, it is $20%. The breaking stresses for armchair and zigzag stanene nanoribbons are $0.2 GPa and $0.4 GPa, respectively. The higher strength in the zigzag direction compared with the armchair direction can be attributed to the difference in bond orientation along the direction of loading. 38 The fracture pattern for uniaxial zigzag loading for stanene nanoribbon at a strain rate of 10 8 s À1 and at 100 K is shown in Fig. 9. As the gure shows, with an increase in strain, a crack is initiated and propagated with the breaking of Sn-Sn atomic bonds, and there is a subsequent decrease in the cross-sectional area, leading to the ultimate fracture of the stanene nanoribbon. Similar fracture behavior has been observed in other 2D materials, such as silicene and GNRs. 75 Fig. 10(a) and (b) interprets the strain-stress dependence of STNRs for various widths with zigzag and armchair edges, respectively. Fig. 10(c) shows the variation of fracture strain of both armchair and zigzag STNRs with increasing width, whereas Fig. 10(d) represents the failure stress for stanene nanoribbons of varying width in both armchair and zigzag directions. As the gures show, both the fracture strain and fracture stress of zigzag stanene nanoribbon are always higher compared with those of armchair stanene nanoribbons. Moreover, the fracture strain and maximum stress or fracture stress of armchair stanene nanoribbons have small variations with an increasing width compared with their zigzag counterparts. For similarly sized GNRs, this trend is in agreement with the study by Chu et al. 75 This phenomenon can be attributed to the difference in failure mechanism for the armchair and zigzag directions. Nucleation of fracture starts from the edge for the armchair direction, whereas for the zigzag direction, it begins from the interior region. 75 Therefore, the initial bond breaking of Sn-Sn in the zigzag direction depends on the width. As a result, the variation in the fracture stress and the fracture strain with increasing width is higher for zigzag STNRs compared with their armchair counterparts.

Conclusion
In summary, we studied the thermo-mechanical properties of stanene nanoribbons using MD simulations with MEAM potential parameters. In this study, we investigated the thermal conductivity of both zigzag and armchair STNRs using EMD simulations based on the Greek-Kubo method. At room temperature, our computed thermal conductivity for 10 nm Â 3 nm STNRs in the zigzag and armchair directions are 0.95 W m À1 K À1 and 0.89 W m À1 K À1 , respectively. The computed values of thermal conductivity are signicantly smaller in STNRs due to lower Debye temperature of stanene accompanied by lower phonon group velocities and a higher scattering rate compared with other similarly sized 2D hexagonal nanoribbons, such as GNRs and h-BNNRs. Our study shows that the thermal conductivity of STNRs follows an inverse relation with temperature for both armchair and zigzag congurations. We also studied the width dependence of the thermal conductivity, and the thermal conductivity was found to increase with increasing width for both congurations, which is in agreement with our computed PDOS for STNRs. In all cases, the zigzag STNR exhibited a higher thermal conductivity than its armchair counterpart did. Signicant phonon scattering at and around vacancy centers causes a reduction of the thermal conductivity with vacancies, thereby enhancing the applicability of stanene nanoribbons in thermoelectrics. Edge vacancy and point vacancies were found to be more destructive than bi-vacancy in terms of reduction in thermal conductivity. Finally, the effects of size and chirality on the mechanical properties of STNRs under a uniaxial tensile strain at 100 K and at a strain rate of 10 8 s À1 were studied using MD simulations. The stress-strain behavior of zigzag-and armchair-congured STNRs showed that zigzag STNRs have a comparatively higher strength and higher fracture strain compared with that of armchair STNRs. Furthermore, with increasing width, both fracture strain and fracture stress of armchair STNRs were found to show small variations compared with their zigzag counterparts. Overall, apart from providing an insight into the thermo-mechanical characteristics of stanene nanoribbons, this study investigated how to tune these characteristics for efficient practical design of thermoelectrics. Thus, our study will aid in the design of thermoelectrics as well as nanoelectronic devices based on stanene nanoribbons and encourages further investigation on the proper optimization of thermo-mechanical characteristics of stanene nanostructures.

Author contributions
A. I. K. conceived the work; A. I. K. and R. P. performed the simulations under the supervision of S. S. All authors discussed and analyzed the data and interpretation and contributed during the writing of the manuscript.

Conflicts of interest
There are no conicts of interest to declare.