Abstract

A novel PCF-based sensor has been presented in this paper to sense different chemicals and biocomponents. The proposed biochemical sensor (BCS) comprises a simple arrangement of rectangular holes. The competence of this BCS in detecting analytes is evaluated employing the full-vector finite element method (FEM). Performance metrics prove the competence of this BCS in sensing various analytes that have a refractive index in the range of 1.33–1.48. The proposed BCS shows ultralower values for both the bulk absorption and confinement loss. This BCS attains a maximum of about 95.82% relative sensitivity at 2.5 THz. Also, the dispersion for this sensor is only 0.12 ± 0.011 ps/THz/cm at the same point of interest. These results prove that the proposed BCS will play a key role in sensing toxic chemicals, illegal drugs, biocomponents, etc. Besides, the simple rectangle-based PCF structure ensures the feasibility of fabrication by practicing the existing fabrication strategies.

1. Introduction

The terahertz (THz) band is gaining high popularity as an apposite operating region for diverse application fields like spectroscopy, sensing, astrophysics, hereditary, biomedical, and telecommunication [18]. This band offers a broad frequency band extending from 0.1 to 10 THz. An entire THz system requires 3 basic elements that include the waveguide, source, and detector. The development of THz source and detector elements in recent years has accelerated THz-based sensing. Besides, numerous THz waveguides, for instance, the Bragg waveguide [9], metal wires [10], metallic tubes [11], etc., have been proposed recently. However, most of these waveguides fail to achieve higher relative sensitivity. To avoid this problem, several photonic crystal fibre (PCF)-based THz waveguides have been proposed in recent years.

Researchers have successfully modelled PCF-based chemical and biosensors functioning in the THz band. As a sensor, PCF is applied to sense multiple chemical and biological elements, for example, the detection of fuel adulteration [12], salinity [13], bane chemicals [14], blood components [5], and cholesterol [15]. This sensor resolves the snags of conventional fibres by maintaining higher diversity and standard values for optical parameters. Besides, this type of sensor is much faster than the existing detection systems. One of the splendid properties of this sensor is the tunability of the optical parameters through careful adjustment of the size, position, and structure of the core and cladding hole of the PCF. Also, this type of sensor offers a lot more advantageous features that includes higher relative sensitivity, lower material and confinement loss (CL), light weight, robustness, cheaper, etc. All these features have widened the potential application of PCF-based sensors in chemical and bio-sensing applications, which includes the study of RNA and DNA [1618], detection of cancerous cells [1921], classification and detection of alcohol [2226], detection of formalin [27, 28], and surface plasmon resonance-based biosensors [2931].

Since the last decades, several refractive index (RI)-based chemical and biosensors have been proposed. In [5], a PCF-based biosensor is proposed. The model is particularly aimed to sense blood components based on RI. The model maintains low CL, but it shows a high effective area (EA) of about 166000 μm2 at 1.5 THz. Besides, the relative sensitivity attained for the sensor is only 79.9% and the dispersion is higher. In [32], a PCF-based chemical sensor is proposed where the cladding holes are heptagonally positioned. The sensing performance of the sensor has been evaluated while detecting benzene, ethanol, and water. The sensor exhibits lower CL, but the relative sensitivity attained for the sensor is only 63.24% at 1.0 THz. Several key optical parameters including numerical aperture and effective material loss (EML) have not been determined for the sensor. Another PCF-based chemical sensor is presented in [14]. The sensor is modelled to detect bane chemicals, namely sarin, soman, and tabun. The sensor exposes 0.009 cm−1 EML and 1.71 × 10−14 cm−1 CL. The relative sensitivity attained for the sensor is 94.4% for tabun at 1.8 THz. However, the true sensing region and numerical aperture are ignored for the sensor. A PCF-based NaCl sensor is presented in [33], where the core of the sensor comprises tetra-rectangles. The sensor shows 91.5% relative sensitivity and 0.004 cm−1 EML, but it has a higher EA of 397340 μm2 at 1.8 THz. Recently, a PCF-based chemical sensor is proposed in [34] to detect ethanol, benzene, and water. The sensor shows 3.02 × 10−8 cm−1 CL, but the EML has not been estimated for this sensor. However, the sensor has achieved almost 78.1% relative sensitivity for benzene. In [35], a Ag-graphene-coated surface plasmon resonance (SPR) sensor is proposed to sense analytes with RI ranging from 1.33 to 1.41. This model shows 12.6 μm/RIU wavelength sensitivity and 79.4 × 10−5 RIU resolution, but this model suffers from high CL of about 1.75 × 102 dB/cm at 625 nm for water. Besides, this sensor comprises complex H-shaped PCF. In [36], another H-shaped PCF-based SPR sensor is presented to enhance the sensing range. This sensor can detect analytes with 1.33 to 1.49 RI. This model exhibits only 1.2 × 103 nm/RIU sensitivity and a high peak loss of 14.34 dB/cm for water.

Most of the existing hollow and porous core PCF sensors are particularly modelled to detect any specific analyte or a group of correlated analytes. Also, most of the available sensor models fail to attain higher relative sensitivity and lower losses at the same time, whereas SPR sensors have a RI limit as most of the SPR sensors fail to effectively sense the analytes with RI greater than 1.41 [35, 3740]. However, most of the existing sensors are modelled with a complex PCF structure that includes a compound arrangement of air holes. This complex structure will certainly increase fabrication constraints. Hence, there is still scope to model a biochemical sensor (BCS) based on fabrication-friendly PCF structure that will be capable to sense a wide variety of chemicals and biocomponents.

In this paper, a simple rectangle-based PCF (R-PCF) structure has been modelled as a BCS. The proposed R-PCF structure has been numerically explored in the THz band in the range of 1.5–4.0 THz. Optical parameters have been studied for the proposed BCS to examine the sensing performance. This study has been carried out for the analytes with the RI ranging from 1.33 to 1.48. In addition to getting the typical values for all the performance metrics, one of the key objectives of this study is to present a highly fabrication-friendly BCS, which is ensured by the presence of only rectangular holes in the PCF structure. R-PCF ensures heterodyne detection, and several sensors based on R-PCF have been proposed [4145]. Besides, asymmetrical PCF structures comprising rectangles have already been implemented [46, 47].

2. Sensor Design

This paper focuses on modelling a novel, uniform, and simple BCS based on PCF. The PCF-based sensor was modelled using FEM-based software. The x-y plane representation of the BCS is depicted in Figure 1. The core region of the sensor contains a single rectangle, which is denoted by Cr in the figure. The height of Cr is 400 μm, whereas the width of it is 410 μm. The cladding area is modelled through the simple arrangement of eight rectangular air holes. These eight rectangles can be grouped into two different classes based on their similarities in height and width. These two types of cladding rectangles are denoted by Cl1 and Cl2. The width and height of Cl1 are 200 μm and 1000 μm individually. On the other hand, the width and height of Cl2 are 200 μm and 400 μm exclusively. The radius of the sensor is 750 μm, including a 75 μm (10% of the radius) perfectly matched layer (PML) at the outer layer. The strut, the gap between two consecutive rectangles, is 10 μm for this PCF-based sensor. Multiple fibre materials are available that include silica, zeonex, topas, etc. [4850]. Various splendid properties of zeonex and topas have influenced the researchers to select them as the fibre material. Both of them have a fixed RI of 1.53, lower dispersion, and lower material loss in the THz band. However, zeonex has some additional features; for instance, higher chemical resistance and glass transition temperature [49, 50]. These two features are highly desired in the successful fabrication of any PCF. Hence, zeonex is chosen to be the fibre substance for the proposed BCS. The light propagation profile in the x and y polarization mode (PoM) are pictured in Figures 2(a) and 2(b), respectively, while sensing the analyte with 1.44 RI. This figure reveals the fact that propagating wave is highly confined within the core region.

The experimental setup of PCF-based models in sensing applications has been described in several papers [5154]. The model represented in Figure 1 has resulted from the trial and error method. Through modelling different structures and comparing their optical parameters, this model has been chosen due to its dominant results. We have denoted this model as the optimum (OPT) model. Keeping the strut value constant, we have tuned the OPT model by ±10% to represent a comparison among their optical parameters and justify the selection of the OPT model as the proposed BCS in this paper. The design parameters for the two models resulted from the tuning along with the OPT model are represented in Table 1 that includes the respective width (W) and height (H) of each rectangle.

3. Fabrication Feasibilities

In earlier stages, PCF fabrication was considered to be implemented using the stack and draw method [55], but this method fails to satisfactorily fabricate the PCF models with a complex and asymmetrical arrangement of holes [56]. To resolve these issues, several more methods have been advanced, which include drilling [57], sol-gel [58], extrusion [59], 3D-printing [46], etc. These methods are capable to fabricate complex PCF structures. The sol-gel method can potentially perform the fabrication of the PCF model with circular holes. On the other hand, both the 3D-printing and extrusion methods fabricate the asymmetrical PCF models, which include rectangular, square, and elliptical holes. The proposed BCS comprises a mono-rectangular core and a cladding with eight rectangular air holes. Hence, both 3D-printing and extrusion will be a suitable choice to fabricate the proposed R-PCF-based BCS. Besides, the fibre with rectangular holes guarantees the heterodyne detection [41, 42]. A few asymmetrical PCF models have been fabricated recently by the Max Planck Institute [46, 47]. The institute has successfully fabricated different PCF models comprising the elliptic, rectangle, etc. This clearly ensures the feasible fabrication of the proposed PCF structure comprising a simple arrangement of rectangles using the existing fabrication strategies.

The core of the BCS needs to be infiltrated with the target sample. There are several available techniques for the selective filling of PCF; for instance, selective filling by collapsing air holes [60], by splicing single-mode fibre [61], by femtosecond laser micromachining [62], etc. Among these strategies, we prefer the selective filling by the femtosecond laser micromachining method. This method comprises three steps. First, the PCF is spliced with a normal single-mode fibre. Then, the single-mode fibre is cut near to spliced junction using FL. Finally, selective laser drilling is performed to make a path to the core or target holes using FL. Through these three steps, the fibre will be prepared to be infiltrated with the target sample. The method has proved its higher efficiency, preciseness, and flexibility, and thereby, it has become a dependable technique of SF [63].

4. Mathematical Methods for Sensing Performance Evaluation

To assess the potentiality of the presented BCS in this paper, several optical parameters need to be calculated. These parameters are estimated separately for each analyte. Filling the core with the target analyte, light is injected into the core. Though the light should be bound into the core, a portion of light spreads outside the core region. The true sensing area, also known as the EA, is the measure of light spreading. The EA is evaluated using the following equation, where E defines the transverse electric field [5, 33].

The numerical aperture (NA) is a measure of the range of angles of the incident light that can enter and exit from a fibre. It is a unitless parameter, which is a function of the EA, the speed of light (c), and the operating frequency (f). Higher NA is desirable for sensing applications, and this desired value results from lower EA. The value for this parameter is calculated using the folllowing equation [64]:

The EML is the loss of optical power due to the absorption of the fibre materials. Zeonex, the fibre material of the proposed BCS, is responsible for this loss. Appropriate designing of the PCF model and lessening the background region will keep the EML at a negligible level. Following the perturbation theory [33, 65], equation (3) has been used to compute this key parameter for the proposed BCS.

The denominator of equation (3) represents the integral of the Poynting vector, where H stands for the magnetic field strength. ε0 and μ0 stand for the free space permittivity and permeability, respectively. η and α represent the RI and absorption coefficient, respectively.

The use of zeonex as a fibre material may introduce a negligible amount of material dispersion [66]. Hence, this dispersion can be avoided, and the only dispersion that needs to be considered is the waveguide dispersion. A PCF with flat dispersion facilitates the examination of the capability of the fibre to inherit multichannel signals. Therefore, the quality of signal propagation from the source to the sink can be evaluated through this parameter. The dispersion (D) for the proposed BCS is estimated using the following equation [67, 68]:where ω and ηeff represent the angular frequency and the effective RI, respectively.

Another key parameter that has an impact on the sensing performance is the CL. This loss results from improper PCF structure and the leaky nature of the mode. Hence, a proper PCF model can result in a reduced CL value. The amount of CL is quantified using the following equation:

The most important optical parameter that directly computes the sensing performance of any PCF-based sensor is the relative sensitivity (RS). This parameter is directly proportional to the power fraction (Pf), which quantifies the amount of light power present in the core compared with the entire injected power. The relative sensitivity for the proposed BCS is calculated usingwhere ηa denotes the RI of the analyte.

5. Simulation Results

Different types of analytes or target samples are injected into the core region of the sensor separately and once at a time. Then, the light propagation through the core region of the proposed BCS provides us with different values for the optical parameters for different types of analytes.

5.1. Optimum Model Selection

To validate the selection of the OPT model as the BCS in this paper, we have first carried out the values of the optical parameters for the OPT − 10%, OPT, and OPT + 10% models. We have performed this simulation for the analyte with 1.44 RI.

The values of EffA in x- and y-PoM for the three models in the THz band are shown in Figures 3(a) and 3(b), respectively. Since EffA is a measure of how light is spread inside the PCF, the lower values for it is desirable. The lowest value of it is attained for the OPT − 10% model, whereas the highest value is found for the OPT + 10% model. The curve for the EffA for the OPT model lies between the other two models. The OPT model exhibits 98975 μm2 and 99121 μm2 in x- and y-PoM, respectively, at 2.5 THz.

The NA indicates the efficiency of a PCF in collecting light. A higher value of NA is expected for every PCF-based chemical sensor. Figure 4 depicts the values of NA in the THz band for the three models. The OPT − 10% model shows the highest NA, while the OPT attains a NA of 0.2104 at 2.5 THz.

The EML in x- and y-PoM are presented in Figures 5(a) and 5(b) separately. In the case of x-polarization, the OPT − 10% model has the lowest EML of about 0.0058 cm−1 at 2.5 THz. Again, the OPT − 10% model shows the lowest EML of about 0.0057 cm−1 at 2.5 THz in y-PoM. However, the OPT model has an EML of only 0.0060 cm−1 and 0.0058 cm−1 in the x- and y-polarization mode, respectively, at the same point which is very close to the lowest value.

Figure 6 represents the dispersion for the three models in the THz band. All three models show low dispersion, but the OPT and OPT + 10% model show lower dispersion compared with the OPT − 10% model. However, all three models exhibit nearly flat dispersion in the frequency band ranging from 2.4 to 3.0 THz. The values of dispersions are 0.1467 ± 0.0129, 0.1206 ± 0.0116, and 0.1206 ± 0.0116 ps/THz/cm in these frequency range for OPT − 10%, OPT, and OPT + 10% model, respectively.

The CL in x- and y-PoM for the three models are presented in Figures 7(a) and 7(b), respectively. All the models show ultra-low CL in the THz band. The lowest values in both x- and y-PoM are attained for the OPT model at 2.5 THz. The numerical values of CL for the OPT model in x- and y-PoM are 3.20 × 10−12 and 9.00 × 10−12 cm−1 at 2.5 THz.

Figures 8(a) and 8(b) are the representation of the relative sensitivity for the three models in x- and y-PoM. The OPT + 10% model exhibits dominant relative sensitivity in both x- and y-PoM compared with the other two models since relative sensitivity is a function of the core power fraction, which is higher for this model. The values of relative sensitivity attained at 2.5 THz for the OPT + 10% model are 95.99 and 95.94%, respectively. On the other hand, the OPT model attains 95.20 and 95.13% relative sensitivity in x- and y-PoM, respectively, at the same frequency.

Table 2 represents the numerical values of the optical parameters for the three models. The optimum values are marked in italics. This italic marking shows that none of the three models is dominant for all the parameters. Considering a trade-off among the optical parameters, the OPT model has been chosen as the BCS in this paper. The OPT model maintains a balance between the losses and relative sensitivities. To attain higher power fraction, relative sensitivity, and NA as well as lower EML, EA, dispersion, and CL are the key objectives for any type of PCF-based sensors. As no particular model among the OPT − 10%, OPT, and OPT + 10% model has attained standard values for all the parameters at the same time, we have to select a model that has nearly attained typical values for the parameters. Hence, we have selected the OPT model as the BCS in this study.

5.2. Optimum Strut Selection

To evaluate the optimum value of strut for our BCS, we have varied the values of it and evaluated the optical parameters for each time. We have simulated the OPT model for three values of strut, i.e., 7.5 μm, 10 μm, and 12.5 μm. We have performed this simulation for the analyte with 1.44 RI.

Figures 9(a) and 9(b) show the value of EffA in the THz band in x- and y-PoM, respectively, for three dissimilar values of the strut. The 7.5 μm strut model has the lowest value for EffA at 2.5 THz in both x- and y-PoM. The values of EffA for the 10 μm strut model in x- and y-PoM are between the values of EffA for the other two models.

The values of the NA in the THz band for the 7.5 μm, 10 μm, and 12.5 μm strut models are shown in Figure 10. All three models show nearly equal values for the NA. However, the 7.5 μm strut model shows the highest value of 0.2136 for the NA at 2.5 THz. The value of NA for the 10 μm strut model is close to the maximum NA in the THz regime.

Figures 11(a) and 11(b) show the values of EML in x- and y-PoM, respectively, for the three models. The 12.5 μm strut model shows the lowest values for EML in both the x- and y-PoM. These values are 0.0059 and 0.00576 cm−1 at 2.5 THz in x- and y-PoM, respectively, whereas the 10 μm strut model has 0.0060 and 0.0058 cm−1 at the same frequency.

The amount of dispersion in the THz band for the three models with different strut values is shown in Figure 12. All these models show nearly flat dispersion from 2.4 to 3.0 THz. The 10 μm strut model exhibits the lowest dispersion at 2.5 THz. The values of dispersion are 0.1223 ± 0.0121, 0.1206 ± 0.0116, and 0.1206 ± 0.0116 ps/THz/cm for 7.5 μm, 10 μm, and 12.5 μm strut models, respectively.

The fluctuation of CL in the x- and y-PoM due to the variation in strut values is revealed in Figures 13(a) and 13(b), respectively. In both PoMs, the CL decreases sharply with the rise in frequency. Though the CL for the 12.5 μm strut model is the lowest initially, the 7.5 μm strut model has the lowest CL at 2.5 THz. The amount of the CL approach nearly zero immediately after 2.0 THz for the three models. However, the 10 μm strut model confirms a lower value of CL in both PoMs at 2.5 THz. The values of CL for the 7.5 μm strut model are 2.66 × 10−13 and 7.16 × 10−12 cm−1 in the x- and y-PoM, respectively, at 2.5 THz. Both these values are very close to the respective values of CL for the 10 μm strut model.

Figures 14(a) and 14(b) depict the values of relative sensitivity as a function of the THz band for the 7.5 μm, 10 μm, and 12.5 μm strut models in x- and y-PoM, respectively. The 7.5 μm strut model shows the highest relative sensitivity for both PoMs throughout the THz band. This model has attained 96.00% and 95.96% relative sensitivity at 2.5 THz for the x- and y-PoM, respectively. On the other hand, the 10 μm strut model shows higher relative sensitivity in both x- and y-PoM compared with the 12.5 μm strut model.

A numerical comparison among the values of optical parameters for the 7.5 μm, 10 μm, and 12.5 μm strut models is shown in Table 3. The standard values attained for each optical parameter are marked in italics. The table reflects that the 7.5 μm strut model has typical values for the optical parameters except for the EML and dispersion. On the other hand, the 10 μm strut model has attained nearly typical values for all the optical parameters maintaining lower losses and higher relative sensitivity. Accounting for this scenario, we have selected the 10 μm strut model in this paper as it balances among the optical parameters through maintaining comparatively lower loss along with higher relative sensitivity.

5.3. RI-Based Analyte Sensing

The optimum PCF structure and the optimum value for the strut have together modelled the proposed BCS of this paper. This sensor is then used to sense different chemicals with different RI. Injecting any particular chemical into the core region of the sensor, the optical parameters have been evaluated for that chemical using the finer mesh analysis. This strategy is repeated several times to evaluate the optical parameters for the chemicals with different RI values ranging from 1.33 to 1.48.

Figures 15(a) and 15(b) describe the effective sensing region in the THz band for the OPT model in x- and y-PoM, respectively. For both PoM, the value of EffA decreases with an increase in the light frequency. This property of this sensor indicated that the light intensity inside the core is higher for higher frequencies. The values of EffA are 104980, 104310, 103050, 101910, 100850, 99880, 98975, 98130, and 97341 μm2 for the chemicals with RI 1.33, 1.34, 1.36, 1.38, 1.40, 1.42, 1.44, 1.46, and 1.48, respectively, at 2.5 THz.

Figure 16 compares the NA of the proposed BCS for chemicals with different RI values. The proposed BCS maintains acceptable values for the NA throughout the THz band. The values of NA are 0.2045, 0.2052, 0.2064, 0.2075, 0.2085, 0.2095, 0.2104, 0.2112, and 0.2121 for the chemicals with RI 1.33, 1.34, 1.36, 1.38, 1.40, 1.42, 1.44, 1.46, and 1.48, respectively, at 2.5 THz.

The values of EML in the x and y-PoM for analytes with dissimilar RI are revealed in Figures 17(a) and 17(b) individually. The proposed model maintains low EML for both PoM. The model shows 0.0051, 0.0052, 0.0054, 0.0056, 0.0057, 0.0059, 0.0060, 0.0061, and 0.0062 cm−1 EML for the analytes with 1.33, 1.34, 1.36, 1.38, 1.40, 1.42, 1.44, 1.46, and 1.48 RI respectively at 2.5 THz in the x-PoM.

Figure 18 compares the dispersion of the proposed BCS for different analytes in the THz band. The model maintains ultra-low dispersion for all the analytes. Besides, the model exhibits nearly flat dispersion for the frequency band ranging from 2.4 to 3.0 THz. The values of dispersion found at 2.5 THz are 0.1185 ± 0.01225, 0.1185 ± 0.0123, 0.1201 ± 0.0127, 0.1212 ± 0.0129, 0.1215 ± 0.01325, 0.1184 ± 0.01215, 0.1206 ± 0.0116, 0.1176 ± 0.0112, 0.122 ± 0.0131 for the analytes with 1.33, 1.34, 1.36, 1.38, 1.40, 1.42, 1.44, 1.46, and 1.48 RI respectively.

The comparison of CL of the proposed BCS in the x- and y-PoM for different analytes are pictured in Figures 19(a) and 19(b), respectively. In both the propagation modes, the CL approaches absolute zero immediately after 2.0 THz. This indicates higher light confinement inside the core of the sensor model. The values of CL are 6.92 × 10−10, 4.28 × 10−10, 1.64 × 10−10, 6.30 × 10−11, 2.40 × 10−11, 8.90 × 10−12, 3.20 × 10−12, 1.10 × 10−12, and 3.00 × 10−13 cm−1 for the analytes with RI 1.33, 1.34, 1.36, 1.38, 1.40, 1.42, 1.44, 1.46, and 1.48, respectively, in the x-PoM at 2.5 THz.

The relative sensitivity profile in the x- and y-PoM are revealed in Figures 20(a) and 20(b), respectively. The proposed model shows higher relative sensitivity for the chemicals with higher RI. The relative sensitivities are 92.45, 92.79, 93.41, 93.95, 94.42, 94.83, 95.20, 95.53, and 95.82% for the analytes with RI 1.33, 1.34, 1.36, 1.38, 1.40, 1.42, 1.44, 1.46, and 1.48, respectively, in the x-PoM at 2.5 THz, whereas in the y-PoM, the proposed BCS has attained slightly lower relative sensitivities compared with the x-PoM.

The numerical comparison of the optical parameters for different analytes is presented in Table 4.

The performance indices of the proposed BCS are compared with the other sensors in Table 5. Compared with the other sensors, the proposed model has typical values for all the optical parameters.

6. Conclusions

This paper presents a simple rectangular-based BCS operating in the THz regime. This BCS can effectively be applied for both biocomponents and chemicals detection purposes. The sensing performance of this sensor is evaluated for the frequency band that ranges from 1.5 to 3.0 THz. Besides, the proposed BCS comprises a simple arrangement of rectangles, which will ensure the feasible fabrication of this sensor. The sensor attains a maximum of 95.82% relative sensitivity at 2.5 THz. Also, this sensor suffers from a negligible amount of EML and CL. Due to its ultra-lower losses and ultra-higher relative sensitivities, this BCS model can effectively and efficiently be applied in biochemical sensing applications for different samples with dissimilar RI ranging from 1.33 to 1.48.

Data Availability

This is a simulation-based study, and no external data have been used.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Abdullah Al-Mamun Bulbul, Abbas Z. Kouzani, M. A. Parvez Mahmud, and Abdullah-Al Nahid conceptualized the study. Formal analysis was performed by Abdullah Al-Mamun Bulbul. Funding acquisition was done by Abbas Z. Kouzani and M. A. Parvez Mahmud. Abdullah Al-Mamun Bulbul and Abdullah-Al Nahid investigated the study. Methodology was developed by Abdullah Al-Mamun Bulbul, M. A. Parvez Mahmud, and Abdullah-Al Nahid. Project administration was carried out by Abbas Z. Kouzani and M. A. Parvez Mahmud. Resources were obtained by Abbas Z. Kouzani and M. A. Parvez Mahmud. Supervision was done by Abbas Z. Kouzani, M. A. Parvez Mahmud, and Abdullah-Al Nahid. Abbas Z. Kouzani and Abdullah-Al Nahid validated the study results. Abdullah Al-Mamun Bulbul, M. A. Parvez Mahmud, and Abdullah-Al Nahid contributed to study visualization. Writing of original draft was performed by Abdullah Al-Mamun Bulbul and Abdullah-Al Nahid. Review and editing were done by Abbas Z. Kouzani and M. A. Parvez Mahmud.

Acknowledgments

The authors are grateful to Deakin University, Geelong, VIC 3216, for providing resources for this research.