A novel perspective for M-polynomials to compute molecular descriptors of borophene nanosheet

Nanomaterials feature exceptional, one-of-a-kind qualities that might be used in electronics, medicine, and other industries. Two-dimensional nanomaterials called borophene have a variety of intriguing characteristics, which helped them to leave an indelible impression in the fields of chemistry, material science, nanotechnology, and condensed matter physics. The concept of modelling the structure of a molecule or chemical network to a chemical graph and then quantitatively analysing them with the aid of topological descriptors was a major advance in the fields of mathematics and chemistry, with a wide range of applications. M-polynomial approach is a very versatile and quick method for computing the degree-based descriptors of chemical graphs or networks. The degree-based descriptors of the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta _{12}$$\end{document}β12-Borophene nanosheet are established in this study utilising the M-polynomial technique. A program code that enables to generate the M-polynomial of any chemical structure was developed in Java platform and the same is displayed. At the conclusion, the numerical and graphical comparison based on the identified analytic expressions is also provided. Additionally, the QSPR analysis was also carried out and the outcoms are presented therein.

The discovery of the two dimensional nanomaterial, graphene bought a stupendous research zest in the field of two dimensional materials.It has led to a remarkable response in the material field and appreciably improved the applications of two dimensional materials in various fields.The unique and superlative properties, the outstanding performances and potential applications of these materials have grabbed major attention from the scientific experts.Two dimensional nanosheets are ultrathin compared to bulk materials and have majority of atoms exposed to the surface, so that it can have greater surface areas, higher chemical and physical activity and quantum confinement effects 1 .This nature enrich them with special photonic, electronic, catalytic, and magnetic properties.Also these nanomaterials have found great application in bio-like materials, drug carriers, biosensors, electronic devices, etc.
Among the two-dimensional allotropes of boron, borophene offers a variety of intriguing characteristics.Due to its minimal weight, borophene has exceptional mechanical stiffness and superconducting characteristics.It also possesses a number of exceptional therapeutic features, which expands the scope of its application in the medical and technological sectors.Owing to its remarkable structural endurance and high reliability, borophene has the potential to be used as a channel for drug distribution in immunotherapy and is predicted to have strong immunological function and catalytic activity 1 .These nanoparticles, with their exceptional physical and chemical attributes are ideal for a wide range of scientific and technological purposes, including higher energy fuel, elevated temperature devices, coatings, atomic engineering, and atmospheric emissions.In the areas of chemistry, bioengineering, nanoelectronics, and quantum physics, borophene has steadily gained significance.In contrast to graphene, borophene exhibits a triangular framework with hexagonal cavities (HH), and the distribution of these cavities changes over the entire structure.β 12 , χ 3 , 2 − Pmmn and honeycomb are the most researched borophene nanostructures.While the β 12 and χ 3 variants of borophene have planar configurations devoid of vertical oscillations and various patterns of systematic boron voids, the 2 − Pmmn phase of borophene has a corrugated structure, while the honeycomb borophene seems to have a structural resemblance to graphene 2 .
Graph theory, a sub-branch of discrete mathematics deals with the study of connection between things.This study was extended to chemical structures and networks paving way to the concept of chemical graph theory.Thus combining graph theory and chemistry, another sub-branch, chemical graph theory emerged, which utilised the concept of graph theory to characterize molecular structures 3 .Here the idea is to model the structure of a chemical compound or networks to chemical graph and is further analyzed quantitatively with the help of topological descriptors.Topological descriptors have found application in the field of chemical graph theory and are crucial in study of structural properties with many applications in structural chemistry.The IUPAC 4 defines the topological index as "A topological index is a numerical value which is associated with the chemical constitution and is used for correlation of chemical structure with various physical properties, chemical reactivity, or biological activity".Topological descriptors help in a very appealing way to determe the physical, chemical, biological or pharmacological properties of a chemical structure.It has application in vast areas of chemistry, informatics, biology, Quantitative Structure-Property Relationships (QSPR) and Quantitative Structure-Activity Relationships (QSAR), etc. 3,[5][6][7][8][9] .Wiener index is the first known structural descriptor, introduced by H. Wiener in 1947 5 .Using the concept, Wiener defined the boiling points for alkanes.Subsequently, in 1972, Gutman and Trinajstić 10 introduced the first and second Zagreb indices.Further, many modified forms and variants of Zagreb indices were introduced, such as modified Zagreb index 11 , Augmented Zagreb index 12 etc.In the study of heat formation of heptanes and octanes, the augmented Zagreb index is useful.In 1975, Randić proposed the Randić index 8 , which has application in the field of drug design 13 .A generalized version of Randić index 14 was also defined.Another variation of the Randić index, which is known as Harmonic index was first introduced and defined by Siemion Fajtlowicz 15 .In chemical graph theory, there are many more topological descriptors, especially many degree based descriptors.The application of these structural descriptors in various fields are remarkable 16 and researchers working on this concept are coming up with new descriptors that correlates with the structural properties of a chemical structure more precisely 17,18 .
Traditionally, the structural descriptors were computed directly using the definitions, which are time consuming.Renowned researchers established several techniques to compute the structural descriptors [19][20][21] , inorder to save computational time.One among them was the polynomial representation of structural descriptors, which acquired broad attention and acceptance in the literature 22 .To compute wiener index, the Hosoya polynomial, which is also known as Wiener polynomial 23,24 was used.It was also useful in the computation of Hosoya index, Hyper Wiener index etc.The wiener index is obtained by evaluating the derivative of Hosoya polynomial by equating the variable to 1. Later, in 2015, the concept of M-polynomial was proposed by Deutsch and Klavžar 25 .The concept of M-polynomial is to provide a general polynomial with the help of which, one can determine various degree based structural descriptors.
In this work, various degree based descriptors of β 12 -Borophene nanosheet 2 are computed using M-polyno- mial approach.The structural descriptors of borophene nanosheets and borophene nanotubes of other variations, namely, 2 − Pmmn , χ 3 and honeycomb were studied by researchers [26][27][28] .In this paper, the β 12 -Borophene nanosheet is denoted as β 12 − BN(r, c) , where r is the count of rows and c is the count of columns that hexagons with centre vertex are repeated.The molecular graph of β 12 − BN(r, c) with r taking value 3 and c taking value 4 is can be denoted as β 12 − BN(3, 4) , and the same is presented in Fig. 1.

Mathematical concepts
This section comprises of the notations and the concepts that are employed for the study.Consider a simple graph , which is connected.The vertex set is denoted by V (�) and E(�) denotes the edge set of .The degree of a vertex µ ∈ V (�) is denoted as d µ , which is the number of edges incident to the vertex.The first Zagreb index and second Zagreb index 10 are defined as, These indices are recogonised to provide quantitative measures of molecular branching 10 and even used in the study of complexity of chemical compounds 16 .
Researchers have developed modified versions Zagreb indices and among them two are the second modified Zagreb index 11 , M m 2 (�) and augmented Zagreb index 12 , A(�).
In the line of Zagreb indices, the eminent scientist Randić proposed the Randić index 8 , Later a generalized version of the Randić index was introduced, replacing (− 1 2 ) by a fixed number κ ∈ R , known as the general Randić index 14 given by, R κ (�) The general reciprocal Randić index 29 is defined as Harmonic index 15 is an another version of Randić index given by, The hyper Zagreb index 30 and the forgotten index 31 (also known as F-index) are defined as, One of the major application of F-index is that it is found to be appropriate for testing the properties of drug related molecular structures 31 .
The sigma index 32 is given by, To improve the Quantitative Structure-Activity-Property Relationships studies, Vukičević and Gašperov 33 established a new class of structural descriptors, called as "discrete Adriatic indices" and it comprises of one hundred and forty-eight descriptors.The symmetric division index 34 is one among them and is a significant predictor of physicochemical properties of chemical compounds.It is also noted that, among all the existing structural descriptors, the symmetric division index has the strongest correlation for predicting the total surface area of polychlorobiphenyls 33 .It is given as, SDD(�) = µη∈E(�) The inverse sum index 35 , is defined as The atom bond connectivity index is utilized in theoretical chemistry for designing chemical compounds with specific physicochemical properties or given pharmacological and biological activities.The geometric arithmetic index and sum-connectivity index are also correlated with the physical and chemical properties of chemical compounds and used in QSPR and QSAR studies.The atom bond connectivity index 36 , ABC(�) , geometric arithmetic index 37 , GA(�) and sum-connectivity index 38 , SC(�) are defined as, ABC(�) For be a graph and m ij (�), i, j ≥ 1 be the number of edges µη ∈ E(�) such that (d µ , d η ) = (i, j) , the M-pol- ynomial is depicted as 25 , The following operators are required while computing the degree based descriptors using the M-polynomial concept 25 .
It can be also observed that By definition of M-polynomial, Therefore, Table 1.The formulas for deriving the degree based descriptors from M-polynomials.

Degree based descriptor t(g, h)
Derivation from M(�; g, h) Second modified Zagreb Symmetric division Inverse sum gh g+h Vol.:(0123456789) M-Polynomial of a chemical graph.In this section, a program code using Java programming language was developed, that enables us to generate the M-polynomial of any chemical graph.The program code is as follows.

Output Window:
Once the program is runned, the output window as given in Fig. 3 is obtained.In this window, the edge partition data of the chemical structure is needed to be inputed inorder to obtain the M-polynomial.Suppose there are x number of edges with end vertex degree (d µ , d η ) , then in 'Input X' column, we have to input the value of d µ , in 'Input Y' column, we have to input the value of d η and in 'Input Equation' column, we have to input the value of x.Next, the inputed data is added by clicking the ' Add' button.Likewise, each data is added and as a final outcome, the required M-polynomial equation is obtained.

Results based on edge partition technique.
This section comprises of theorems based on the computation of degree based descriptors of β 12 -borophene nanosheet, β 12 − BN(r, c); r, c ≥ 1 .In this section, the degree counting method and edge partition technique is applied and the expressions are derived from the definition of the corresponding degree based descriptors.
Proof Let � = β 12 − BN(r, c); r, c ≥ 1 .has 5rc + 2r vertices and (12r − 1)c + r edges.The edge set parti- tions based on the vertex degree are (3,3) , (3,4) , (3,5) , (3,6) , (4,4) , (4,5) , (4,6) , (5,5) and (5,6) for G, where, (i,j) = {µη ∈ E(G) :  www.nature.com/scientificreports/properties.The process begins with a suitable molecular topological descriptor and ends with some inference, hypothesis, and prediction about the molecule's behaviour, properties and characteristics.The availability of highquality experimental data, as well as its correctness, are vital in this investigation.The selection of the optimal topological descriptor for modelling in analysis is a crucial phase in the QSPR/QSAR process.Because there is no consensus on the best molecular description, all feasible descriptors are determined.A major application of QSPR/QSAR models is that the properties, activities, behavior, etc. of a newly designed or untested chemical compound can be inferred from the molecular structure of similar compounds whose properties, activities,  Linear regression model.As mentioned, 2 − Pmmn , β 12 , χ 3 are the most researched borophene nanostructures.In this section, the linear regression models for few properties of these borophene nanostructures with their degree based descriptors are obtained using least square procedure.
The following regression model is considered for the analysis.
where P and SD stands for property and structural descriptor respectively.For the study, the properties such as shear modulus (G), Young's modulus (Y) and Poisson's ratio (ν) are considered.The numerical values of various topological descriptors and the properties 39 for corresponding borophene nanostuctures are given in Tables 5 and 6 respectively.
Using the Microsoft Excel analytical tools, the regression analysis was carried out and the correlation coefficients ( R ) generated are given in Tables 7 and 8.
Based on the correlation coefficient ( R ), the regression model for Forgotten index, (F) is found to be best fitted for the shear modulus, and the regression equation can be given as For Young's modulus, the regression model for both second Zagreb, (M 2 ) and augmented Zagreb index, (A) is found to have same correlation coefficent.
The regression equation with respect to first Zagreb index can be given as P = x(SD) + y, P = 0.0154(F) − 8.0449.www.nature.com/scientificreports/For Poisson's ratio, the regression model for first Zagreb index, (M 1 ) is found to be best fitted and the regres- sion equation can be given as The scatter diagrams that show relation between properties and their best fitted descriptors are given in Fig. 8. Thus it can be concluded that, the above mentioned properties of different variations of borophene nanostructures can be predicted with the help of regression equations using the computed topological descriptors.

Conclusion
In this work, the degree based structural descriptors of β 12 -borophene nanosheet were computed using the M-polynomial approach.That degree based structural descriptors can be routinely computed from the M-polynomial.This helps in reducing the problem of determining the expression in each particular case to a single problem of determining the M-polynomial, thus making the process general, fast and efficient.Our work is reduced to determining the edge partitions due to the developed program code's relative ease of use.The descriptors is having application in disparate fields and thus the study is relevant in present scenario.The analytic expressions derived and the numerical and graphical comparison together with the Quantitative Structure-Activity-Property Relationship model may help researchers and scientists in future advancements in the field of nanomaterials and to study the nature and behaviour of the β 12 -borophene nanosheet and other borophene nanostructures without any laboratory experimentations.

Table 2 .
Computed numerical values of degree based descriptors for β 12 − BN(r, c).

Table 3 .
Computed numerical values of degree based descriptors for β 12 − BN(r, c).

Table 4 .
Computed numerical values of degree based descriptors for β 12 − BN(r, c).

Table 5 .
Numerical values of degree based descriptors for the borophene structures.

Table 6 .
Values of the properties for the borophene structures.

Table 7 .
Correlation coefficients ( R ) between the properties and degree based descriptors of the borophene structures.

Table 8 .
Correlation coefficients ( R ) between the properties and degree based descriptors of the borophene structures.