ORIGINAL RESEARCH article

Front. Chem., 17 June 2026

Sec. Theoretical and Computational Chemistry

Volume 14 - 2026 | https://doi.org/10.3389/fchem.2026.1823507

Face reverse degree topological analysis of TP-COFs, existence of isentropic COFs and spectral characteristics

  • 1. Department of Mathematics, Loyola College, University of Madras, Chennai, India

  • 2. Department of Mathematics, Loyola College, Chennai, India

  • 3. Department of Mathematics, Women’s Christian College, Chennai, India

  • 4. Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia

  • 5. Department of Biology, College of Science, Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia

Abstract

The recent synthesis of triple-pore covalent organic frameworks (TP-COFs), specifically TP-COF-DAB and TP-COF-BZ, represents a significant milestone in framework chemistry, as replicating such multi-pore architectures remains highly challenging. Despite their structural novelty, the mathematical perspective of their topology and associated properties is still limited. In this study, closed-form expressions are derived for reverse degree topological indices of these TP-COFs based on a bitrapezium topological arrangements. Scaled face reverse degree indices are employed to explicitly characterize pore geometry in triple-pore architectures. Furthermore, reverse degree based entropy measures are computed to quantify structural orderness and to explore isentropic configurations within the TP-COF family. Using these structural parameters, quantitative structure-property models are constructed for spectral energy, demonstrating strong predictive capability. In addition, the spectral diameter and HOMO–LUMO gap are analyzed using graph spectral techniques, providing insights into the spectral characteristics of TP-COFs.

1 Introduction

Covalent organic frameworks have emerged as highly promising materials over the past decades due to their remarkable crystallinity, extensive surface areas, and wide-ranging applications, including catalysis, sensing, energy storage, drug delivery, optoelectronic devices, separation processes, and gas adsorption and storage (Bai et al., 2016; Ding et al., 2016; Kandambeth et al., 2017; Liao et al., 2016; Liu et al., 2026). As an emerging class of crystalline porous organic materials, the properties and performance of COFs strongly depend on the characteristics of their pores, which are determined by the topological structures of the networks and their sizes (Carrington et al., 2022; Liang et al., 2020; Yusran et al., 2024). In particular, appropriate pore architectures and functionalization enable COFs to act as efficient adsorbents for the capture of hazardous metal ions, organic and biological pollutants, as well as greenhouse gases (Ghazi et al., 2018; Li et al., 2021; Xiao et al., 2021; Xin et al., 2021). Traditionally, COFs have been constructed through the principles of reticular chemistry, which enable the precise prediction of their structures based on the symmetry and geometry of the building blocks used for condensation reactions (Abuzeid et al., 2021; Yaghi, 2016; Yu et al., 2022). The design strategies for COFs have therefore focused mainly on combining building units with compatible symmetries. Using this approach, various COFs with tetragonal, hexagonal, or triangular pores have been designed and synthesized since the first two COFs were reported in 2005 (Chen et al., 2025; Côté et al., 2005; Gao et al., 2018; Geng et al., 2020; Jing et al., 2023). In these traditional COFs, there is usually only one kind of pore in a given framework, which limits their structural complexity and functional diversity.

Recent advances have introduced a paradigm shift with the development of heteropore COFs that incorporate two or three distinct types of pores within a single continuous network (Dalapati et al., 2016; Du et al., 2016; Han et al., 2020). These structures enrich the structural diversity of the COF family and open new avenues for multifunctionality. A pioneering design strategy combining vertex truncation and multiple linking site approaches has enabled the fabrication of triple-pore COFs (Qian et al., 2017). In this design, the building block (TPTCA) is shown in Figure 1a, while the linear linkers 1,4-diaminobenzene (DAB) and benzidine (BZ) are shown in Figures 1b,c. Condensation of TPTCA with these linkers yields the structural units for TP-COF-DAB and TP-COF-BZ, presented in Figure 2. These COFs exhibit an fxt topology, where the V-shaped geometry of the TPTCA unit is crucial for generating the connectivity that gives rise to three distinct pore types (Bhambri et al., 2022). This advancement reflects a broader challenge in polymer and materials science, moving beyond primary connectivity toward precisely controlled higher-order structures. While biological polymers achieve remarkable hierarchical organization, replicating such structural precision in synthetic systems remains a significant challenge (Kricheldorf, 2006). The design of multi-pore COFs such as TP-COF-DAB and TP-COF-BZ represents a meaningful step toward this goal, embodying the idea of tailor-made frameworks in materials design. Analyzing their intricate architectures is therefore not only of theoretical interest but also essential for understanding and exploiting their full functional potential.

FIGURE 1

FIGURE 2

Since the pore architecture is the primary factor governing the physical and chemical properties of COFs, a topological analysis that quantitatively characterizes the pore environment is essential. Topological indices, which are numerical descriptors derived from graph theory, have been widely employed to characterize the structure of molecular and extended frameworks (Arockiaraj et al., 2025a; Junias and Clement, 2024; Kalaam et al., 2024; Mondal et al., 2026; Tang et al., 2025). In the context of COFs, several studies have applied degree-based and related topological indices to correlate structural features with properties (Arockiaraj et al., 2025b; Augustine and Roy, 2022; Kurian et al., 2025; Tu et al., 2025; Zhang X. et al., 2025). These indices provide a powerful way to understand how the structure of the framework influences its physicochemical properties (Arockiaraj et al., 2025c; Manuel and Angamuthu, 2025; Paul et al., 2025; Raza et al., 2024). However, conventional degree-based topological indices mainly describe atom connectivity, which may not fully capture the complexity of the pores themselves. To address this limitation, a face-degree-based topological framework is employed, in which indices are defined in terms of face degrees, thereby providing a more direct and accurate representation of the pore structure. Notably, scaled face reverse-degree based indices have been shown to be effective for benzenoid hydrocarbons compared with conventional degree-based indices (Arockiaraj et al., 2026a). In parallel, reverse-degree topological indices offer an alternative perspective for capturing structural information by emphasizing complementary connectivity patterns within the network (Ahmad et al., 2023; Kalaam and Greeni, 2024; Kalaam and Greeni, 2025; Rao et al., 2024; Youssef et al., 2024). Consequently, the reverse-degree framework provides additional structural information that can be used for subsequent QSPR modeling (Aq et al., 2026; Ravi, 2024). Thus, in this study, we analyze TP-COFs with reverse-degree indices, and scaled face reverse-degree indices. We further investigate the existence of isentropic structures by means of reverse-degree based graph entropies, which have attracted considerable attention in the literature as measures of structural complexity and information content (Kavitha et al., 2021; Naeem et al., 2024; Zhang G. et al., 2025). In addition, entropy measures are relevant to both thermodynamics and information theory, which provide valuable insights into the behaviour of complex systems across scientific and engineering domains (Abraham et al., 2022; Hussain et al., 2025; Jawahar and Clement, 2026; Peter et al., 2025; Yu et al., 2024). Building upon these analyses, we further investigate spectral descriptors and develop predictive models for spectral graph energy using the computed topological indices.

2 Computational techniques

The triple-pore covalent organic frameworks investigated in this study are modeled as hydrogen-suppressed graph structures , where represents the set of vertices corresponding to atoms, and represents the set of edges corresponding to bonds. In this work, we aim to explore the degree parameter of each vertex , defined as the number of edges incident on the vertex, to its full potential. Accordingly, a modified reverse-degree parameter is employed to reveal different degree combinations and to represent vertices with lower connectivity. The reverse degree is defined as

Here, denotes the maximum degree in , and the reversing parameter takes integer values from 1 to . For , the degree values repeat cyclically. The reverse-degree parameter defined in Equation 1 follows existing formulations in the literature, where reverse-degree concepts are employed to capture complementary connectivity patterns in complex networks. Such a definition has been effectively used in recent studies to improve degree-based indices (Ahmad et al., 2023; Kalaam and Greeni, 2024; Kalaam and Greeni, 2025; Rao et al., 2024; Youssef et al., 2024). In this work, the same formulation is adopted in the context of TP-COFs to appropriately reflect their connectivity variations. The degree-based topological indices are defined as . In order to incorporate the reverse-degree parameters, we define the associated topological indices as

Furthermore, to integrate the effect of reverse-degree modification with the structural influence of pores, we employ the face-reverse-degree indices for each face , defined aswhere denotes the set of boundary edges of the face , and represents the set of all faces of . The set includes both internal faces , which are completely bounded cycles, and the external face , which outlines the outer boundary of the structure. During the boundary traversal of a face, pendant edges are counted twice (Arockiaraj et al., 2026b). The total number of faces is given by the cardinality . Aggregating the face contributions derived from Equation 3, the scaled face-reverse-degree index of is defined as

To emphasize the structural influence of these reverse-degree modifications and their pore-level extension, we consider the first Zagreb , second Zagreb , hyper-Zagreb , forgotten , arithmetic , bi-Zagreb , tri-Zagreb , bi-Zagreb-arithmetic , tri-Zagreb-harmonic , and tri-Zagreb-arithmetic , which are defined respectively as

Among these topological indices, the first Zagreb index was originally introduced in relation to the total -electron energy of molecular graphs, while the forgotten index was later shown to provide improved structural sensitivity in subsequent investigations (Furtula and Gutman, 2015; Gutman and Trinajstić, 1972). In recent studies, these indices have been extended through reverse-degree and scaled face-based formulations to better represent structural variations (Arockiaraj et al., 2026a; Kalaam and Greeni, 2025). These extensions have been shown to provide additional structural information compared to existing indices and support their applicability in structural analysis and potential QSPR studies (Ahmad et al., 2023; Aq et al., 2026; Ravi, 2024).

3 Results and discussion

In TP-COFs, the structural units can be arranged in various configurations, such as linear, rectangular, hexagonal, parallelogram, and bi-trapezium, depending on the topology and connectivity of the building units. These configurations influence the overall geometry, porosity, and functional characteristics of the materials. In the present study, the bi-trapezium configuration is used to model the arrangement of structural units as graph structures. The representative bi-trapezium configurations considered in this work are illustrated in Figures 3, 4. These figures depict the structural arrangement of the TP-COFs, where vertices represent atoms and edges represent bonds. In particular, the figures highlight variations in connectivity and the arrangement of faces, which form the basis for computing the reverse-degree and scaled face reverse-degree topological indices. This configuration includes two important structural forms, namely linear and hexagonal, as subcases, and therefore provides a general framework suitable for representing a wide range of geometric variations observed in TP-COFs. Let and represent the bi-trapezium graph structures corresponding to the two TP-COFs, namely TP-COF-DAB and TP-COF-BZ, respectively. Here, and denote the structural parameters that determine the extent of the bi-trapezium arrangement. The vertex and edge set counts of these COFs are given, respectively, by and .

FIGURE 3

FIGURE 4

3.1 Reverse degree topological indices

To streamline the computation of the reverse-degree topological indices, we first partition the bonds into bond classes using the bond-partition technique, which groups edges according to the degrees of their end vertices. Mathematically, this is defined asLet denote the set of bond classes of the TP-COFs. Then, the frameworks and share the same bond classes, . However, their cardinalities differ, and the corresponding values are presented in Table 1.

TABLE 1

Bond classes

Bond partitions of and .

From Table 1, the degree set of is , with maximum degree 3. Accordingly, the reversing parameter ranges from 1 to 3, and the corresponding modifications are

Thus, the modified bond partitions for are given by . For , the partitions are and for , . The corresponding cardinalities remain unchanged as per order given in Table 1. Then, the general formula for reverse-degree topological indices, obtained by substituting into Equation 2, is given by

By substituting the topological indices into Equations 5 and 6, the corresponding closed-form expressions for the reverse-degree topological indices are obtained, and the results are presented in the form in Results 1 and 2.

Result 1

For the bi-trapezium configurationof dimension, the reverse-degree indices are given by

Result 2

For the bi-trapezium configurationof dimension, the reverse-degree indices are given by

Based on Results 1–2, the reverse-degree topological indices of and for are computed, and the values are presented in Tables 24.

TABLE 2

24001848792042241200424860722600102243688
4040310013296709620207140101964382.6667171246209.3333
652049882140811432326011508164207081.33332751610022.6667
56804352186729968284010032143206165.3333240248730.6667
98407512322561723249201735224744106964140015128
7320560424048128403660129241844479483092411252
131601003643104230326580231963306814310.66675528420233.3333
173201319656688302968660305164349218841.33337266026630.6667
896068562942415712448015816225689730.66673782413773.3333
164801256053952288328240290404139217925.33336916825338.6667
22320169927300839024111603931256016242889352834320
2648020152865924628813240466326644028818.666711090440717.3333
319224601051255921596565280523476135724900
5360412017616937626809480134965842.6667227048229.3333
863266202832015080431615252217009417.33333644413254.6667
752857802472013160376413308189408209.33333183611558.6667
130089960426242270465042296832664142005479219976
96967440318241694448481713624384105764096814888
173841330056928303288692306844362818982.66677314026697.3333
2286417480748323987211432403445735224973.33339609635114.6667
11864910038928207285932209642982812942.66675010018217.3333
2176016640712323795210880384005459223765.33339148833418.6667
294482250096336513361472451948738363217212366045228
34928266801142406088017464616088756038162.666714661653645.3333

Reverse-degree topological indices of and for varying .

TABLE 3

38884992204961051219448880155043379.2412085884.8
655283963448017688327614948260845698.4692129917.6
10584135405561628536529224124420769210.411148416021.6
9216118004846424864460821016366648017.69721613950.4
159842042483904430567992364086348013915.216802424196.8
118801520462448320405940270844724410336.812522017983.2
2138427308112192575761069248692848841862022456432372
281523593214763275768140766408411170024517.629537242618.4
14544186087643239216727233152578241265615322422016
267843419214048072096133926097610628823324.828110440547.2
36288462961902249763218144825841439283160838044854936
4305654920225664115824215289797617074437505.645125665182.4
518466482726413968259211832206164516.8547447843.2
8712111564576023448435619868346047594.49177213181.6
14040179567366437752702031996557081224414758021244
12240156646425632928612027904485921067212880018520
21168270481109765688010584482168392818465.622216832030.4
157682017282752424087884359406258013749.616582823858.4
282963614014828876008141486443611214824687.229675642816.8
372244752419500899960186128474814748432480.839012456327.2
1929624680101248518889648439767656816827.220285629196.8
354244523218560095136177128065614036830908.837134453603.2
47952612002511361287362397610915218993641846.450227272561.6
5688072584297856152688284401294642252724964059564086072

 Reverse-degree topological indices of and for varying .

TABLE 4

264019921008060961320463280882544162484224
44803400171841038422407880137844304278087168
7280555228032169283640128322248069764551211648
6320480824288146723160111281948060643936810112
110408448426242572855201948834176105606936017664
8160621631392189604080143762517678245092813056
1480011344572163452874002614445872141449320823680
19520149847555245584976034504605681864012320031232
10000762438496232485000176243087295846248816000
18560142407180843328928032800575681772811705629696
252001936897632588961260044568782642404815933640320
2992023008115968699521496052928929602854418932847872
3576282013968832817886396111483444231245700
604047802366414104302010820188845804392689628
9776776038400228804888175363064093766384815584
8504674033360198804252152442662081645541213556
1478411760581763465673922654446416141609686423568
109688700430562565654841966834356105247155617484
19792157607795246432989635552621921894412988031552
2607220780102768612081303646852819882494017133241564
1343210660527523143267162409242092128848770021412
248001976097728582081240044560779682372816289639536
33624268201326247898416812604441058043214822122053604
39904318401574409376019952717441256003814426267263616

 Reverse-degree topological indices of and for varying .

3.2 Scaled face reverse degree topological indices

In the case of scaled face reverse degree topological indices, the TP-COFs are partitioned into faces. For both and , the face partitions consist of five internal faces and one external face , where the internal face partitions are illustrated in Figures 5, 6. In particular, these figures illustrate the classification of faces using different colours, where each colour corresponds to a distinct face class. Higher-dimensional structures follow a similar pattern. The cardinalities of the faces , and are the same for both TP-COFs. However, the edge composition of each face, represented by , varies. An exception occurs for the face , whose cardinality varies while its edge composition remains the same. The corresponding cardinalities and edge compositions are presented in Table 5. Moreover, the total number of faces in and are given by and , respectively. Let . Then and . On substituting and the modified reverse-degrees into Equation 4, the general formulas for scaled face reverse-degree topological indices are presented below.

FIGURE 5

FIGURE 6

TABLE 5

1

Cardinalities and edge compositions of the face partitions of and .

By substituting the topological indices into Equations 7 and 8, the closed-form expressions for the scaled face reverse-degree-based topological indices are presented as in Results 3 and 4.

Result 3

For the bi-trapezium configurationof dimension, the scaled face reverse-degree indices are given by

Result 4

For the bi-trapezium configurationof dimension, the scaled face reverse-degree indices are given by

The scaled face reverse-degree topological indices of and are computed using Results 3 and 4, and the corresponding values are summarized in Tables 68.

TABLE 6

2402.33111835.15777901.30884230.99341201.16564237.48896066.15112595.719310150.44313696.9430
4045.17633078.117313267.76347111.52882022.58817123.293610189.64614375.372417000.46826226.9801
6530.36754953.127121369.356911463.10263265.183811483.494716416.22977069.709127322.314210055.0260
5688.09684321.203918634.69829992.29032844.048410009.300714313.49436155.068023851.36888757.1256
9857.98427460.346432206.645317285.95254928.992117318.330624746.298910678.782141117.054815181.1862
7331.04575564.337924001.813012873.13713665.522912895.383618437.47507934.779330702.596111287.3121
13185.67449967.701743044.426523109.02316592.837223153.376133076.724814287.900654912.721820307.4482
17355.782713107.256356617.860630403.34808677.891330463.039043510.604318811.752172181.211826731.8132
8974.00856807.494929369.015315754.02554487.004315781.503422561.52049714.498337553.982013817.5187
16513.395112475.112953882.411128932.18548256.697628988.508041407.298217897.037668708.770025433.7526
22368.327516879.045572927.073639168.982611184.163839247.373156048.028124250.348592922.764534454.3065
26538.493120018.706286500.891646463.479313269.246546557.199366482.185428774.2354110191.976740878.7508
3239.55852506.082710746.84105734.67561619.77935745.64128240.75833491.360913919.83394979.7562
5440.84834196.908218012.36129618.54482720.42419637.756513815.45305868.302723286.72178365.3938
8763.43016743.340328960.970715474.29024381.715015506.770322217.63059458.446837381.186613476.4134
7642.19015887.822525278.215413502.57043821.095113530.012619390.39298245.274232654.220311751.1061
13207.367210145.503543593.108723302.10176603.683623352.870833447.605214261.834556204.072120312.8999
9843.55197578.770332544.196017386.65554921.775917422.322124965.425810622.256842021.949315136.8469
17651.359813547.769358225.618031130.07958825.679931199.129144677.848719065.256475027.655327149.4632
23216.705317805.763576541.642840930.115811608.352641022.468858735.879325081.921298579.638935711.4894
12044.92339269.734339810.238521270.76996022.461721314.657630540.504212999.244851389.790618522.6019
22095.375516950.077272858.280938958.126411047.687739045.452755908.203723868.692493851.526633986.0586
29903.418622919.514598541.284752702.255714951.709352822.933175621.770232311.8382126860.214945998.9990
35468.808227177.5904116857.605662502.424717734.404162646.398789680.015238328.5302150412.755054561.0863

Scaled face reverse-degree topological indices of and for varying .

TABLE 7

3890.33114981.488920486.633310523.65561945.16568871.820015505.14453377.097841111.04755891.5645
6557.17638379.293634472.468517713.88143278.588114936.469826093.17495695.058769063.05269931.2938
10594.367513515.494755618.827028587.83765297.183824109.862242103.33239205.4989111275.233116047.2361
9224.096811777.300748459.085324904.48394612.048421001.397536681.78468013.060197017.345813971.1334
16001.984220390.330683926.581943145.92088000.992136392.314763536.251413908.4661167750.250424239.5022
11891.045715175.383662445.995832095.22855945.522927066.429347270.612210331.0767124972.496318011.0147
21409.674427265.3761112235.123957704.371810704.837248675.050484969.747918611.4752224227.620932431.8735
28187.782735879.0390147704.991275946.913314093.891364066.8216111825.952224507.0078294967.628942700.5575
14558.008518573.503476433.049339286.04267279.004333131.511957859.546012649.1007152928.063722050.9164
26817.395134140.5080140543.991672262.975613408.697660957.9031106403.483623314.5016280705.962940624.2886
36336.327546231.3731190336.383797873.637618168.163882567.7006144105.010731595.0746379970.876955045.5804
43114.493154845.1993225806.8640116116.465421557.246597959.6924170961.664737490.6399450712.718465314.3463
5231.55856741.641227689.075114205.79262615.779311973.199720947.43384535.528255846.86707919.5888
8792.848311313.756546479.754223852.24134396.424120106.604735165.99787625.951393632.533613311.7452
14171.430118210.770374830.691038409.15037085.715032382.200456619.920712294.9541150587.708421455.9061
12354.190115886.012665270.975933498.95066177.095128240.202849384.963310716.4025131419.791918703.9777
21367.367227432.8708112742.577657876.836110683.683648800.238085309.706918542.5742226714.492932352.1603
15915.551920458.322184062.403443145.75927957.775936373.873963604.081313806.8644169207.652424096.2393
28563.359836655.1291150655.057277344.799114281.679965218.4889113999.928224790.2258302843.050143248.4938
37576.705348202.4688198128.4639101723.526418788.352685779.1740149925.995132616.4938398143.156156896.9168
19476.923325030.6576102853.931752792.61659738.461744507.580977823.274116897.3315206995.807429488.5151
35759.375545877.4527188567.782896812.877417879.687781636.8282142690.330131037.8905378972.341354144.8604
48407.418662074.9331255162.9590131013.092824203.7093110482.3516193088.025942021.3866512615.295873297.4505
57420.808273622.3987302636.8385155392.041228710.4041131043.2069229014.439949847.6797607916.815686945.9367

 Scaled face reverse-degree topological indices of and for varying .

TABLE 8

2635.33781917.31989911.99066077.35101317.66894552.65757994.67082525.329915634.24024233.3456
4469.64753260.706516864.002710342.58982234.82377730.353913603.29634269.176626657.41777182.1183
7259.26505308.833327462.726416845.05983629.632512568.098222153.89316915.208343494.093711667.3216
6303.80654603.848023814.921814607.22583151.903210907.654419211.07386012.962037678.106510130.6509
11004.03178061.544241707.215225584.12685502.015919065.575933645.671010463.386166142.464617688.6773
8137.90865946.894630765.423618871.63444068.954314084.803224818.52907756.723748697.837013079.0935
14748.651310814.063255950.731534322.60517374.325625562.714545136.668314011.515888788.785123709.7868
19448.434714271.198773840.136145297.73889724.217433719.633459568.937518461.7997117247.110931267.0697
9971.98307289.894437715.720623135.93194985.991517261.877330425.82639500.473659717.095316027.4924
18493.209813566.500370193.839743060.83909246.604932059.710156627.339417559.6251111434.236829730.7944
25103.345018432.907095375.193758509.379812551.672543536.251976942.286823814.2267151516.963040360.4632
29803.013821889.8939113263.843169484.055314901.506951692.907791373.949228264.4735179973.699447917.5542
3480.88302528.745513005.02287947.53181740.44156009.628510476.27733371.186420461.05445582.5795
5878.30354279.150322015.514413457.21392939.151710157.453717736.36425678.787634695.26429429.8194
9513.13996939.060335706.680121828.55954756.569916452.200228767.61989170.765156360.472215263.5147
8275.61986029.388831025.256718966.47904137.809914305.008624995.86787986.347248927.779813276.8923
14385.265610508.342054077.746233061.06237192.632824893.607643569.404313847.085585455.165423083.4456
10672.89637779.562540034.710324475.58525336.448218452.458832255.147810293.890663159.637317123.9020
19257.280414077.478972448.079244293.12159628.640233334.759358370.600418523.3697114548.312430903.1910
25366.589518556.001395498.360358386.357812683.294743922.590776942.359124384.0003151071.073140709.1786
13070.15349529.704049044.021629984.61356535.076722599.857439514.317512601.426077391.168920970.8808
24129.249017646.554290818.104555524.996212064.624541775.803273171.550323199.6385143640.804838722.8596
32713.162923943.8724123228.396375340.651516356.581456657.035399284.523931428.9681195023.049152501.3576
38822.383528422.2794146278.092989433.534119411.191867244.6629117855.813537289.5699231544.577562307.1972

 Scaled face reverse-degree topological indices of and for varying .

3.3 Existence of isentropic structures

The term isentropic, meaning constant entropy, has been of long-standing interest in thermodynamics and has led to diverse concepts such as isentropic analysis, isentropic efficiency, and isentropic flow. In the context of topological analysis, isentropic structures are defined as systems that attain equal entropy values even in the presence of variations in topology, geometry, or composition. The recognition of such structures is significant, as it reveals underlying structural similarity, indicates equivalent information content, and enables the classification of networks into entropy-equivalent families. These measures provide an effective framework for quantifying structural complexity and information distribution in molecular and networked systems (Anand and Bianconi, 2009; Bonchev, 1983; Mowshowitz, 1968). They have been successfully applied to both organic and inorganic materials, where entropy-based approaches contribute to the analysis of framework organization and support rational material design (Jacob et al., 2023; Rahul et al., 2022). These entropy formulations are constructed through additive and multiplicative self-powered transformations of the underlying indices, thereby incorporating both frequency and magnitude information associated with bond partitions. Recent advancements have incorporated reverse-degree modifications into entropy measures. As a result reverse-degree-based entropy of takes the formwhere, is the multiplicative self-powered transformation defined as . Explicit mathematical expressions for the entropy measures can be derived from the corresponding topological and self-powered indices. However, the resulting formulations are complex and difficult to present concisely. Therefore, we focus on the numerical evaluation of the entropy values for and . The computed values are presented in Tables 911.

TABLE 9

9.19128.920810.37879.75198.49689.758810.11379.274610.62819.6216
9.75789.485610.943710.31729.063810.324410.67889.842011.191810.1882
9.89199.619511.077710.45129.197910.458410.81279.976111.325710.3223
9.89199.619511.077710.45129.197910.458410.81279.976111.325710.3223
9.09738.829110.28669.65938.40279.665810.02149.179810.53759.5276
10.18299.909611.367910.74169.489210.749011.103010.267611.615310.6133
10.864510.590012.048511.422410.171011.429911.783610.949612.295011.2949
10.864510.590012.048511.422410.171011.429911.783610.949612.295011.2949
11.512211.236712.695312.069410.818812.077112.430511.597712.941211.9426
11.512211.236712.695312.069410.818812.077112.430511.597712.941211.9426
11.999611.723613.182212.556411.306312.564312.917512.085313.427712.4300
11.999611.723613.182212.556411.306312.564312.917512.085313.427712.4300
9.47109.20368.776710.658110.028310.039810.39209.558710.90909.9002
10.03599.767211.221810.59229.342110.604010.955810.124211.471810.4651
10.16989.90109.476011.355610.726110.737811.089610.258111.605610.5990
10.16989.90109.476011.355610.726110.737811.089610.258111.605610.5990
9.37879.11318.684410.56729.93709.948310.30109.465710.81939.8079
10.460110.190511.645211.01599.766411.027711.379310.548711.894710.8893
11.140710.870210.447212.325011.695811.707812.059111.229612.573911.5699
11.140710.870210.447212.325011.695811.707812.059111.229612.573911.5699
11.787511.516411.094212.971312.342212.354312.705411.876813.219612.2168
11.787511.516411.094212.971312.342212.354312.705411.876813.219612.2168
12.274512.002911.581213.457912.828912.841113.192012.363913.705912.7038
12.274512.002911.581213.457912.828912.841113.192012.363913.705912.7038

Reverse-degree based entropy measures of and for varying .

TABLE 10

9.67669.921811.334210.66748.982910.499711.05549.537912.028110.0913
10.243810.487811.900611.23389.550311.066311.621810.105512.593810.6585
10.377810.621812.034711.36789.684311.200311.755810.239512.727910.7925
10.377810.621812.034711.36789.684311.200311.755810.239512.727910.7925
9.58209.828411.240610.57388.888310.405810.96189.442811.93559.9966
10.669210.912612.325611.65889.975811.491412.046810.531113.018411.0839
11.351111.593813.007012.340110.657812.172912.728111.213313.699311.7658
11.351111.593813.007012.340110.657812.172912.728111.213313.699311.7658
11.999012.241213.654512.987711.305812.820613.375611.861414.346312.4138
11.999012.241213.654512.987711.305812.820613.375611.861414.346312.4138
12.486612.728414.141713.475011.793413.307913.862912.349114.833412.9014
12.486612.728414.141713.475011.793413.307913.862912.349114.833412.9014
9.958110.20339.264511.614610.946610.781311.33549.821512.308310.3724
10.523510.767812.179411.51139.830011.346211.900110.387112.872610.9377
10.657410.90169.964012.313311.645211.480112.034010.521013.006511.0716
10.657410.90169.964012.313311.645211.480112.034010.521013.006511.0716
9.865410.11159.171711.522610.854510.689011.24349.728412.217110.2796
10.947911.191712.603411.935310.254511.770312.324210.811713.296311.3621
11.628711.872010.935413.283812.615712.450913.004611.492713.976412.0430
11.628711.872010.935413.283812.615712.450913.004611.492713.976412.0430
12.275812.518611.582513.930613.262513.097713.651312.139914.622812.6901
12.275812.518611.582513.930613.262513.097713.651312.139914.622812.6901
12.762813.005412.069614.417413.749413.584614.138212.627115.109413.1771
12.762813.005412.069614.417413.749413.584614.138212.627115.109413.1771

 Reverse-degree based entropy measures of and for varying .

TABLE 11

9.30609.036710.655510.15168.61189.873610.43509.261511.14049.7762
9.87729.611811.229910.72529.183410.446611.00919.831111.717810.3474
10.01159.746411.364410.85969.317810.581011.14369.965311.852510.4817
10.01159.746411.364410.85969.317810.581011.14369.965311.852510.4817
9.20678.933610.553210.04998.51259.772710.33299.164311.03539.6770
10.305010.041611.659311.15429.611310.875211.438310.257912.148810.7751
10.989310.728212.345411.839810.295911.560512.124310.941212.836511.4594
10.989310.728212.345411.839810.295911.560512.124310.941212.836511.4594
11.639411.380012.996812.490910.946012.211212.775511.590413.489112.1094
11.639411.380012.996812.490910.946012.211212.775511.590413.489112.1094
12.128111.869713.486312.980311.434812.700413.265012.078613.979312.5981
12.128111.869713.486312.980311.434812.700413.265012.078613.979312.5981
9.59879.36888.904810.967810.450410.183610.74239.555611.476310.0653
10.16719.939611.538311.02049.473410.753111.312510.122712.048610.6336
10.301210.07389.607611.672511.154610.887211.446710.256712.183010.7677
10.301210.07389.607611.672511.154610.887211.446710.256712.183010.7677
9.50269.27028.808610.869610.352510.086410.64429.461011.37649.9691
10.593210.367011.965511.44749.899711.179711.739710.548112.476911.0597
11.275911.051110.582512.649312.131011.863012.423411.230013.161711.7423
11.275911.051110.582512.649312.131011.863012.423411.230013.161711.7423
11.924511.700811.231213.298912.780412.512113.072911.877913.812112.3909
11.924511.700811.231213.298912.780412.512113.072911.877913.812112.3909
12.412412.189411.719213.787313.268813.000313.561312.365514.301012.8789
12.412412.189411.719213.787313.268813.000313.561312.365514.301012.8789

 Reverse-degree based entropy measures of and for varying .

From Tables 911, we clearly observe that the configurations pair , , , , and so on, yield similar entropy values, leading to the existence of isentropic structures. Furthermore, since the isentropic dimensional patterns are consistent for both and frameworks, this behavior can be expressed in the general forms . Selected examples illustrating the isentropic correspondence between and are listed in Table 12. It should be noted that the term “isentropic” used in this study is based on graph-theoretic entropy and does not refer to thermodynamic entropy in the classical materials science sense. The entropy considered here is a structural descriptor derived from topological indices, quantifying the distribution of connectivity within the network. Accordingly, two TP-COF structures are termed isentropic if they exhibit equal values of this graph-based entropy, even if their topology differs. Thus, the term indicates structural equivalence in terms of information content rather than thermodynamic behavior.

TABLE 12

Isentropic structures
Bond classes
1476205238405352
446456161185614880
2525406961452

 Bond partitions of isentropic structures.

4 Spectral descriptors and QSPR analysis

In this section, the spectral characteristics of the considered structures are investigated to gain insight into their structural properties and connectivity patterns. However, direct computation through density functional theory (DFT) is computationally demanding, time-consuming, and often challenging for large-scale frameworks. Therefore, we adopt a graph-theoretical approach based on the eigenvalues of the adjacency matrix. In particular, graph energy, HOMO-LUMO gap, and spectral diameter are computed using Python programming through adjacency matrices. The graph energy of a comprising vertices is determined from its spectrum obtained via the adjacency matrix representation. Let denote the eigenvalues of this matrix. The corresponding graph energy, is defined as

Furthermore, the HOMO-LUMO gap and spectral diameter are computed from the ordered eigenvalue set. Specifically, the HOMO-LUMO gap is evaluated as the difference between the smallest positive eigenvalue and the largest negative eigenvalue, whereas the spectral diameter is defined as the difference between the maximum and minimum eigenvalues of the spectrum. The computed spectral descriptors presented in Table 13 demonstrate distinct behavioural patterns across the investigated structures. In particular, the HOMO–LUMO gap remains constant for each framework, attaining fixed values of for and for irrespective of the dimensional parameters . This invariance indicates structural control of the spectral gap within each class. A similar near-constant trend is observed for the spectral diameter, which exhibits only marginal variation with increasing structural size. In contrast, graph energy displays a clear monotonic growth as the parameters increase, reflecting its sensitivity to structural expansion and complexity. Owing to this discriminative behaviour, graph energy is selected as the principal descriptor for subsequent QSPR modelling.

TABLE 13

(1,1)888.81254.84240.72951190.78334.84230.7492
(2,1)1498.87524.84290.72952002.15994.84260.7492
(2,2)2422.72974.84330.72953227.98524.84280.7492
(3,1)2108.93794.84310.72952813.53644.84270.7492
(3,2)3660.37594.84350.72954868.25934.84290.7492
(4,1)2719.00054.84310.72953624.91304.84270.7492
(4,2)4898.02224.84350.72956508.53354.84300.7492
(4,3)6449.46034.84360.72958563.25644.84300.7492
(5,1)3329.06324.84310.72954436.28954.84270.7492
(5,2)6135.66854.84360.72958148.80764.84300.7492
(5,3)8314.69024.84360.729511032.42814.84300.7492
(5,4)9866.12834.84370.729513087.15104.84300.7492

Computed spectral descriptors of for varying .

Both univariate and multivariate regression analyses are employed for effective model construction. The univariate regression model is expressed as where denotes the slope and represents the intercept. The multivariate regression model with two independent variables is given by where and are the regression coefficients associated with the indices and , respectively, and is the intercept. The selection of indices for both univariate and multivariate modelling is initially guided by correlation analysis. For the univariate case, the correlation results are illustrated through heatmaps shown in Figures 7, 8. Based on the best-performing index, univariate regression models are constructed. Similarly, for the multivariate case, all possible pairs of indices are correlated with graph energy using Algorithm 1 for multivariate regression analysis. The most significant index pairs are subsequently employed to develop multivariate regression models. The results of both univariate and multivariate regression analyses are presented in Table 14, where the coefficient of determination () and root mean square error (RMSE) are reported to nine and six decimal places, respectively.

FIGURE 7

FIGURE 8

TABLE 14

EquationR2RMSE
Univariate10.029072
10.00822
0.9999999990.062665
10.010652
Multivariate10.000006
10.000006
10.013664
10.010339

 Results of univariate and multivariate regression modelling for and .

Algorithm 1

  • Input: Dataset containing reverse degree indices , face reverse degree indices , and

  •    graph energy ;

  • Output: Best reverse index pair, regression model, and error measures for ;

  • Preprocessing:

  • Load dataset and remove non-numerical identifiers;

  • Separate reverse degree indices ;

  • Separate scaled face-reverse degree indices ;

  • Model construction (reverse-degree indices):

  • Initialize ;

  • foreachpairdo

  •   Fit regression model ;

  •   Compute ;

  • ifthen

  •     Update reverse-degree based pair and model

  • Model construction (face-reverse degree indices)

  • Initialize ;

  • foreachpairdo

  •   Fit regression model ;

  •   Compute ;

  • ifthen

  •     Update best scaled face-reverse degree based pair and model;

  • Error evaluation:

  • Compute ;

  • Output:

  • Report best reverse degree index pair;

  • Report best face reverse degree index pair;

  • Report regression equations;

  • Report and RMSE;

Table 14 clearly shows that for both and , the association between graph energy and the topological indices follows an almost identical pattern in its correlation behaviour with the indices. In the univariate analysis, the same index consistently exhibits the strongest correlation with energy for both structures. A slight variation is observed only in the case of face reverse-degree-based indices when moving to the multivariate setting. This overall consistency highlights the strong structural similarity between and with respect to their spectral characteristics. In the univariate models, the index dominates the predictive performance. However, in the multivariate framework, the models based on and indices emerge as the most effective. The multivariate regression models exhibit comparatively low RMSE values, indicating improved predictive accuracy. Among these, the reverse degree-based indices perform particularly well. Specifically, for , the reverse degree descriptor with modification parameter provides the best performance, whereas for , the descriptor with performs better. The resulting regression equations corresponding to the best-performing models, along with the standard error (.) and -value, are given below.

Substituting the index values into these reverse-degree regression equations yields

The low and high -value indicate the good fit of the model. We examine the two best-performing multivariate models, statistical significance and stability using the p-values, leave-one-out cross-validation coefficient , and the variance inflation factor (VIF). Both models yielded , confirming statistical significance at the 1 level. The cross-validation results produced , indicating perfect internal predictive agreement within the dataset. However, the VIF values are very high , suggesting the presence of multicollinearity between the and predictors. The elevated VIF values indicate a strong linear association between the and descriptors, suggesting potential redundancy. This behaviour is not unexpected, as both descriptors are degree-based and derived from closely related degree distributions. Consequently, partial linear dependence may arise, leading to possible coefficient instability in the multivariate regression model. Therefore, while the combined model provides superior fitting performance, caution is required when interpreting the individual regression coefficients, as multicollinearity may inflate statistical measures.

To confirm that the observed perfect predictive ability is not merely a consequence of linear dependence between and , external validation was performed. The regression model was applied to higher-dimensional structures that were not involved in the model fitting process, and the results are presented in Table 15. The sustained predictive accuracy observed in this independent dataset confirms that the excellent performance is not an artefact of multicollinearity but instead reflects a stable and structurally consistent relationship between the selected indices and the energy. To demonstrate the performance of the multivariate models, first-degree polynomial surface fitting were carried out, and the resulting visualizations are depicted in Figures 9, 10.

TABLE 15

ActualPredictedActualPredicted
(5,5)10789.9827710789.9827714312.9763414312.97634
(6,3)10179.9201110179.9201113501.5997813501.59978
(7,4)14851.7553514851.7553419683.2895419683.28954
(8,2)9848.607349848.6073313069.6299813069.62998
(9,5)23271.5713423271.5713430820.8437230820.8437

Actual vs. predicted energy values of and .

FIGURE 9

FIGURE 10

Figures 9, 10 clearly illustrate the effectiveness of the polynomial surface fitting for the considered dataset. The fitted surface closely aligns with the observed data points, indicating a strong agreement between the predicted and computed values. The smooth planar nature of the surface suggests a stable and linear relationship between the spectral descriptor and the selected topological descriptors. This confirms the suitability of the proposed multivariate model for predicting graph energy.

4.1 Y-scrambling

To further validate the robustness and predictive reliability of the proposed QSPR models, a Y-scrambling test is performed. In this procedure, the dependent variable values are randomly shuffled while keeping the independent variables unchanged, and new regression models are constructed for each permutation. The RMSE obtained from each scrambled model is then compared with that of the original model. In this study, the dependent variable is randomly permuted 500 times. For each permutation, the model is recalibrated and the corresponding RMSE is computed. The distribution of RMSE values obtained from the scrambled models is subsequently compared with the RMSE of the original model to assess the presence of chance correlations. The comparative results of the calibrated RMSE () and the recalibrated RMSE () are graphically illustrated in Figure 11.

FIGURE 11

From Figure 11, a substantial difference between and is observed. The significantly lower compared to the values confirms that the developed QSPR models are not influenced by random correlation and demonstrates their statistical robustness.

5 Conclusion

In this work, we carried out a comprehensive structural and spectral investigation of triple-pore covalent organic frameworks using reverse-degree, and scaled face reverse-degree indices. The study provides a systematic mathematical characterization of these frameworks and establishes clear relationships between their structure and spectral behavior. The entropy analysis revealed that the considered TP-COF structures exhibit isentropic behavior, indicating structural regularity and uniform information distribution. From the spectral perspective, graph energy increases monotonically with structural growth, while the HOMO–LUMO gap remains invariant and the spectral diameter shows only marginal variation. The invariance of HOMO–LUMO indicates structural control of the gap within the series in the graph-theoretical model, whereas the monotonic growth of graph energy highlights its suitability as a descriptor for QSPR modelling. Regression analysis further strengthened these observations. Both univariate and multivariate models were constructed to examine the predictive capability of the proposed indices. The multivariate models, particularly those involving the first and second Zagreb indices, produced extremely low RMSE (0.000006) with a leave-one-out cross-validation coefficient close to one. The statistical significance was confirmed through very low p-values, and external validation demonstrated high predictive accuracy. Furthermore, the regression models were validated using a Y-scrambling test. These results provide a mathematical characterization of TP-COFs and establish their structural, spectral, and predictive properties within a unified framework.

Statements

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

TR: Visualization, Validation, Conceptualization, Methodology, Writing – original draft, Software. MA: Writing – review and editing, Methodology, Formal Analysis, Supervision, Conceptualization. AS: Investigation, Writing – review and editing, Data curation, Conceptualization, Methodology. HA: Funding acquisition, Writing – review and editing, Visualization, Validation, Methodology. NA-H: Formal Analysis, Writing – review and editing, Data curation, Investigation, Funding acquisition.

Funding

The author(s) declared that financial support was received for this work and/or its publication. Authors are thankful to Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2026R299), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

The author MA declared that they were an editorial board member of Frontiers at the time of submission. This had no impact on the peer review process and the final decision.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AbrahamJ.ArockiarajM.JencyJ.KavithaS. R. J.BalasubramanianK. (2022). Graph entropies, enumeration of circuits, walks and topological properties of three classes of isoreticular metal organic frameworks. J. Math. Chem.60, 695732. 10.1007/s10910-021-01321-8

  • 2

    AbuzeidH. R.EL-MahdyA. F. M.KuoS.-W. (2021). Covalent organic frameworks: design principles, synthetic strategies, and diverse applications. Giant6, 100054. 10.1016/j.giant.2021.100054

  • 3

    AhmadA.KoamA. N. A.AzeemM. (2023). Reverse-degree-based topological indices of fullerene cage networks. Mol. Phys.121 (14), e2212533. 10.1080/00268976.2023.2212533

  • 4

    AnandK.BianconiG. (2009). Entropy measures for networks: toward an information theory of complex topologies. Phys. Rev. E Stat. Nonlin. Soft Matter Phys.80, 045102. 10.1103/PhysRevE.80.045102

  • 5

    AqibM.MalikM. A.AliY.JuttS.ImranM. (2026). Predicting optimal drugs for fin rot by reverse topological QSPR analysis. Comput. Biol. Chem.121, 108847. 10.1016/j.compbiolchem.2025.108847

  • 6

    ArockiarajM.Jeni GodlinJ. J.RadhaS.AzizT.Al-HarbiM. (2025a). Comparative study of degree, neighborhood and reverse degree based indices for drugs used in lung cancer treatment through QSPR analysis. Sci. Rep.15, 3639. 10.1038/s41598-025-88044-x

  • 7

    ArockiarajM.MaaranA.DossC. I. A. (2025b). Topological characterization and predictive modeling of graph energy in ionic covalent organic frameworks. Comput. Mater. Contin.85 (1), 637655. 10.32604/cmc.2025.065674

  • 8

    ArockiarajM.JencyJ.ShaliniA. J.BalasubramanianK. (2025c). Degree-based topological insights and graph entropies of kagome lattice covalent organic frameworks. Chem. Pap.79, 135157. 10.1007/s11696-024-03743-5

  • 9

    ArockiarajM.RajendranT.BalasubramanianK. (2026a). Face-degree and reverse-degree topological indices for quantitative structure property predictions of benzenoid hydrocarbons. Chem. Pap.80, 377402. 10.1007/s11696-025-04399-5

  • 10

    ArockiarajM.RajendranT.ShaliniA. J.BalasubramanianK. (2026b). Graph theoretical characterization and QSPR modeling of aluminum and magnesium hydroxide networks through integrated face degree topological indices. Eur. Phys. J. E49, 14. 10.1140/epje/s10189-025-00554-8

  • 11

    AugustineT.RoyS. (2022). Topological study on triazine-based covalent-organic frameworks. Symmetry14 (8), 1590. 10.3390/sym14081590

  • 12

    BaiL.PhuaS. Z. F.LimW. Q.JanaA.LuoZ.ThamH. P.et al (2016). Nanoscale covalent organic frameworks as smart carriers for drug delivery. Chem. Commun.52 (22), 41284131. 10.1039/c6cc00853d

  • 13

    BhambriH.KhullarS.MandalS. K. (2022). Nitrogen-rich covalent organic frameworks: a promising class of sensory materials. Mater. Adv.3 (1), 19124. 10.1039/d1ma00506e

  • 14

    BonchevD. G. (1983). Information theoretic indices for characterization of chemical structures. Chichester, England: John Wiley Sons.

  • 15

    CarringtonM. E.RampalN.MaddenD. G.O’NolanD.CasatiN. P. M.DivitiniG.et al (2022). Sol-gel processing of a covalent organic framework for the generation of hierarchically porous monolithic adsorbents. Chem8 (11), 29612977. 10.1016/j.chempr.2022.07.013

  • 16

    ChenF.ZhengH.YusranY.LiH.QiuS.FangQ. (2025). Exploring high-connectivity three-dimensional covalent organic frameworks: topologies, structures, and emerging applications. Chem. Soc. Rev.54 (1), 484514. 10.1039/d4cs00703d

  • 17

    CôtéA. P.BeninA. I.OckwigN. W.O’KeeffeM.MatzgerA. J.YaghiO. M. (2005). Porous, crystalline, covalent organic frameworks. Science310 (5751), 11661170. 10.1126/science.1120411

  • 18

    DalapatiS.JinE.AddicoatM.HeineT.JiangD. (2016). Highly emissive covalent organic frameworks. J. Am. Chem. Soc.138 (18), 57975800. 10.1021/jacs.6b02700

  • 19

    DingS.-Y.DongM.WangY.-W.ChenY.-T.WangH.-Z.SuC.-Y.et al (2016). Thioether-based fluorescent covalent organic framework for selective detection and facile removal of mercury(II). J. Am. Chem. Soc.138 (9), 30313037. 10.1021/jacs.5b10754

  • 20

    DuY.YangH.WhiteleyJ. M.WanS.JinY.LeeS.-H.et al (2016). Ionic covalent organic frameworks with spiroborate linkage. Angew. Chem. Int. Ed. Engl.55 (5), 17371741. 10.1002/anie.201509014

  • 21

    FurtulaB.GutmanI. (2015). A forgotten topological index. J. Math. Chem.53, 11841190. 10.1007/s10910-015-0480-z

  • 22

    GaoQ.LiX.NingG.-H.XuH.-S.LiuC.TianB.et al (2018). Covalent organic framework with frustrated bonding network for enhanced carbon dioxide storage. Chem. Mater.30 (5), 17621768. 10.1021/acs.chemmater.8b00117

  • 23

    GengK.HeT.LiuR.DalapatiS.TanK. T.LiZ.et al (2020). Covalent organic frameworks: design, synthesis, and functions. Chem. Rev.120 (16), 88148933. 10.1021/acs.chemrev.9b00550

  • 24

    GhaziZ. A.KhattakA. M.IqbalR.AhmadR.KhanA. A.UsmanM.et al (2018). Adsorptive removal of Cd2+ from aqueous solutions by a highly stable covalent triazine-based framework. New J. Chem.42 (12), 1023410242. 10.1039/c8nj01778f

  • 25

    GutmanI.TrinajstićN. (1972). Graph theory and molecular orbitals. Total pi-electron energy of alternant hydrocarbons. Chem. Phys. Lett.17, 535538. 10.1016/0009-2614(72)85099-1

  • 26

    HanX.-H.QiQ.-Y.ZhouZ.-B.ZhaoX. (2020). Designed synthesis of a two-dimensional covalent organic framework with three-level hierarchical porosity. Chin. J. Chem.38 (12), 16761680. 10.1002/cjoc.202000317

  • 27

    HussainN.HussainM.RiazA.HaidarA. (2025). Statistical modeling between entropy- and degree-based topological indices of two dimension carbon nitride monolayers. Chem. Pap.79, 85618573. 10.1007/s11696-025-04335-7

  • 28

    JacobK.ClementJ.ArockiarajM.PaulD.BalasubramanianK. (2023). Topological characterization and entropy measures of tetragonal zeolite merlinoites. J. Mol. Struct.1277, 134786. 10.1016/j.molstruc.2022.134786

  • 29

    JawaharK.ClementJ. (2026). Exploring entropy-energy relationships in PHI zeolite through topological indices. J. Chem. Inf. Model.66, 28202829. 10.1021/acs.jcim.5c03095

  • 30

    JingX.ZhangM.MuZ.ShaoP.ZhuY.LiJ.et al (2023). Gradient channel segmentation in covalent organic framework membranes with highly oriented nanochannels. J. Am. Chem. Soc.145 (38), 2107721085. 10.1021/jacs.3c07393

  • 31

    JuniasJ. S.ClementJ. (2024). Predictive analytics of conductance and HOMO-LUMO gaps with topological descriptors of porphyrin nanosheets. Phys. Scr.99, 015208. 10.1088/1402-4896/ad0c1b

  • 32

    KalaamA. R. A.GreeniA. B. (2024). Comparative analysis of modified reverse degree topological indices for certain carbon nanosheets using entropy measures and multi criteria decision-making analysis. Int. J. Quantum Chem.124 (1), e27326. 10.1002/qua.27326

  • 33

    KalaamA. R. A.GreeniA. B. (2025). Comparative analysis of kekulene tessellation patterns using generalized reverse degree-sum descriptors combined with graph entropy and energy properties. Int. J. Quantum Chem.125 (8), e70043. 10.1002/qua.70043

  • 34

    KalaamA. R. A.GreeniA. B.ArockiarajM. (2024). Modified reverse degree descriptors for combined topological and entropy characterizations of 2D metal organic frameworks: applications in graph energy prediction. Front. Chem.12, 1470231. 10.3389/fchem.2024.1470231

  • 35

    KandambethS.BiswalB. P.ChaudhariH. D.RoutK. C.Kunjattu HS.MitraS.et al (2017). Selective molecular sieving in self-standing porous covalent-organic-framework membranes. Adv. Mater.29 (2), 1603945. 10.1002/adma.201603945

  • 36

    KavithaS. R. J.AbrahamJ.ArockiarajM.JencyJ.BalasubramanianK. (2021). Topological characterization and graph entropies of tessellations of kekulene structures: existence of isentropic structures and applications to thermochemistry, nuclear magnetic resonance, and electron spin resonance. J. Phys. Chem. A125 (36), 81408158. 10.1021/acs.jpca.1c06264

  • 37

    KricheldorfH. R. (2006). Polypeptides and 100 years of chemistry of alpha-amino acid N-carboxyanhydrides. Angew. Chem. Int. Ed. Engl.45 (35), 57525784. 10.1002/anie.200600693

  • 38

    KurianS.RoyS.GayathriK. B.JyothishK. (2025). Entropy comparison of benzothiadiazole-based covalent-organic frameworks. Front. Chem.13, 1704165. 10.3389/fchem.2025.1704165

  • 39

    LiX.CuiY.-Y.YangC.-X. (2021). Covalent coupling fabrication of microporous organic network bonded capillary columns for gas chromatographic separation. Talanta224, 121914. 10.1016/j.talanta.2020.121914

  • 40

    LiangR.-R.JiangS.-Y.AR.-H.ZhaoX. (2020). Two-dimensional covalent organic frameworks with hierarchical porosity. Chem. Soc. Rev.49 (12), 39203951. 10.1039/d0cs00049c

  • 41

    LiaoH.WangH.DingH.MengX.XuH.WangB.et al (2016). A 2D porous porphyrin-based covalent organic framework for sulfur storage in lithium–sulfur batteries. J. Mater. Chem. A4 (19), 74167421. 10.1039/c6ta00483k

  • 42

    LiuR.JiaY.XiaY.WangS. (2026). Metalloporphyrin-based covalent organic frameworks: design, construction, and photocatalytic applications. Catalysts16 (1), 76. 10.3390/catal16010076

  • 43

    ManuelM.AngamuthuP. (2025). Mathematical modeling of H1-antihistamines: a QSPR approach using topological indices. ACS Omega10 (41), 4901949034. 10.1021/acsomega.5c07577

  • 44

    MondalS.DasP.RazaZ.PalA.GhorbaniM. (2026). Graph spectrum of neighbourhood sombor matrix and structure-property modelling. Theor. Comput. Sci.1068, 115758. 10.1016/j.tcs.2026.115758

  • 45

    MowshowitzA. (1968). Entropy and the complexity of graphs. I. An index of the relative complexity of a graph. Bull. Math. Biophys.30, 175204. 10.1007/BF02476948

  • 46

    NaeemM.RaufA.AkhtarM. S.IqbalZ. (2024). QSPR modeling with curvilinear regression on the reverse entropy indices for the prediction of physicochemical properties of benzene derivatives. Polycycl. Aromat. Compd.44 (3), 14521469. 10.1080/10406638.2023.2196429

  • 47

    PaulD.ArockiarajM.EmiletD. A.GreeniA. B.KalaamA. R. A. (2025). Molecular descriptor-based QSPR analysis of physicochemical properties in neuromuscular drugs. Mod. Phys. Lett. B39 (28), 2550155. 10.1142/s0217984925501556

  • 48

    PeterP.ClementJ.ArockiarajM.JacobK. (2025). Predictive modeling of molecular interaction energies using topological and spectral entropies of zeolite AWW. Front. Chem.13, 1543588. 10.3389/fchem.2025.1543588

  • 49

    QianC.QiQ.-Y.JiangG.-F.CuiF.-Z.TianY.ZhaoX. (2017). Toward covalent organic frameworks bearing three different kinds of pores: the strategy for construction and COF-to-COF transformation via heterogeneous linker exchange. J. Am. Chem. Soc.139 (19), 67366743. 10.1021/jacs.7b02303

  • 50

    RahulM. P.ClementJ.JuniasJ. S.ArockiarajM.BalasubramanianK. (2022). Degree-based entropies of graphene, graphyne and graphdiyne using Shannon’s approach. J. Mol. Struct.1260, 132797. 10.1016/j.molstruc.2022.132797

  • 51

    RaoY.ChenR.AhmadH.AhmadU. (2024). Reverse Zagreb indices and their application in the evaluation of physiochemical properties of anticancer/antibacterial drugs. ACS Omega9 (28), 3105631080. 10.1021/acsomega.4c04409

  • 52

    RaviV. (2024). QSPR analysis of drugs used for treatment of hepatitis via reduced reverse degree-based topological descriptors. Phys. Scr.99, 105236. 10.1088/1402-4896/ad729d

  • 53

    RazaZ.ArockiarajM.MaaranA.ShaliniA. J. (2024). A comparative study of topological entropy characterization and graph energy prediction for marta variants of covalent organic frameworks. Front. Chem.12, 1511678. 10.3389/fchem.2024.1511678

  • 54

    TangJ.-H.HamrazA.UllahA.HamedY. S.ZamanS.BelayM. B. (2025). Chemical applicability and predictive potential of certain graphical indices for determining structure-property relationships in polycrystalline acid magenta (C20H17N3Na2O9S3). Sci. Rep.15, 13886. 10.1038/s41598-025-98024-w

  • 55

    TuY.YousafS.TariqS.TawfiqF. M.AslamA.TolaK. A. (2025). Computational analysis of metal organic framework and covalent organic framework using degree based topological indices with QSPR validation. Sci. Rep.15, 37311. 10.1038/s41598-025-21240-x

  • 56

    XiaoY.MaC.JinZ.WangJ.HeL.MuX.et al (2021). Functional covalent organic framework for exceptional Fe2+, Co2+ and Ni2+ removal: an upcycling strategy to achieve water decontamination and reutilization as smoke suppressant and flame retardant simultaneously. Chem. Eng. J.421, 127837. 10.1016/j.cej.2020.127837

  • 57

    XinJ.ZhouY.WangX.XuG.XieM.LiuL.et al (2021). Room-temperature synthesis of magnetic covalent organic frameworks for analyzing trace benzoylurea insecticide residue in tea beverages. Food Chem.347, 129075. 10.1016/j.foodchem.2021.129075

  • 58

    YaghiO. M. (2016). Reticular chemistry–construction, properties, and precision reactions of frameworks. J. Am. Chem. Soc.138 (48), 1550715509. 10.1021/jacs.6b11821

  • 59

    YoussefM. Z.Al-DayelI.HanifM. F.SiddiquiM. K.AhmedH.TolasaF. T. (2024). On statistical evaluation of reverse degree based topological indices for iron telluride networks. Sci. Rep.14, 20533. 10.1038/s41598-024-71645-3

  • 60

    YuX.ChenZ.LiX. (2022). Reticular chemistry in action: 3D porphyrinic COFs with scu topology. Chem8 (11), 28992901. 10.1016/j.chempr.2022.10.013

  • 61

    YuG.SiddiquiM. K.HussainM.HussainN.SaddiqueZ.PetrosF. B. (2024). On topological indices and entropy measures of beryllonitrene network via logarithmic regression model. Sci. Rep.14, 7187. 10.1038/s41598-024-57601-1

  • 62

    YusranY.MiaoB.QiuS.FangQ. (2024). Functional covalent organic frameworks: design principles to potential applications. Acc. Mater. Res.5 (10), 12631278. 10.1021/accountsmr.4c00195

  • 63

    ZhangX.ArockiarajM.MaaranA.DossC. I. A. (2025). Topological information entropy and spectral energy analysis in aminal-linked covalent organic frameworks. Sci. Rep.15, 19861. 10.1038/s41598-025-01210-z

  • 64

    ZhangG.LiY.RaufA.AfzalM. A.AliP.AslamA. (2025). Estimating the prediction ability of reverse degree-based entropy indices for the physicochemical properties of lymes disease drugs. Front. Phys.13, 1536995. 10.3389/fphy.2025.1536995

Summary

Keywords

face reverse degree indices, HOMO‐LUMO gap, isentropic structures, QSPR modeling, spectral energy, triple-pore covalent organic frameworks

Citation

Rajendran T, Arockiaraj M, Shalini AJ, Alshanbari HM and Al-Hoshani N (2026) Face reverse degree topological analysis of TP-COFs, existence of isentropic COFs and spectral characteristics. Front. Chem. 14:1823507. doi: 10.3389/fchem.2026.1823507

Received

05 March 2026

Revised

31 March 2026

Accepted

20 April 2026

Published

17 June 2026

Volume

14 - 2026

Edited by

Xinguo Ren, Chinese Academy of Sciences (CAS), China

Reviewed by

Sreejith Shankar, Council of Scientific and Industrial Research (CSIR), India

Janani Ezhilan, Joy University, Vadakkankulam, India

Updates

Copyright

*Correspondence: Micheal Arockiaraj,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics