, where E is the total energy of the deformed system, E 0 is the equilibrium energy, and S 0 is equilibrium area. Strain is applied in both directions (xy plain), ranging from -2% to +2%. The computed elastic constants C 2D in units of Jm -2 are 438.01 (459.26), 430.03 (452.18), and 423.60 (452.56) along the x (y) direction for qHP's of C 60 , C 58 B 1 N 1 , and C 54 B 3 N 3 , respectively. For comparison, the experimental and theoretical elastic constants (C β ) are 260 J/m² for BN, 350 Jm -2 for graphene, and 385 J/m² for qTP C 60 , indicating their strain-induced stability . The conductivity of the given systems was determined using semiclassical Boltzmann theory within the constant relaxation time approximation, employing the BoltzTrap code . Using the Boltz-Trap code, conductivity and mobility can be obtained only under the relaxation time approximation. In this approach, the conductivity is calculated using the equation