Abstract

This paper presents a systematic study in reviewing the application of finite element method for the analysis of correction mechanism of spine deformity due to scoliosis. The study is aimed at systematically (1) reviewing the use of finite element analysis in spine deformity case, (2) reviewing the modelling of pedicle screw and rod system in scoliosis surgery, and (3) analysing and discussing gap between the studies. Using the restricted key phrases, the review gathered studies from 2001 to 2021 from various electronic databases (Scopus, ScienceDirect, PubMed, Medline, and WorldCAT). Studies were included if they reported a finite element study on spine deformity. Studies that did not fully disclose their methodology and results had significant discrepancies, not published in English or not yet published were all disqualified. Regardless of inconsistencies in the methodological design of the studies, the quality of all papers was above the acceptable level. A total of fifteen manuscripts were considered for inclusion and were given a comprehensive review. This study indicates that analysing the forces acting on the spine, as well as the interrelationship between the force, stress, and degree of correction (which measured as the Cobb angle), could help to improve the corrective mechanism procedure of spine deformity. Pedicle screws and its placement strategies are also important as it influence the corrective forces for scoliosis treatment. Hence, the findings of this study could potentially be used as a guidance to develop a reliable finite element analysis that can predict the biomechanics responses during the corrective spine deformity treatment.

1. Introduction

Scoliosis is a three-dimensional (3D) spinal deformity characterised by axial vertebral rotation. The Cobb angle value is used to determine the severity of a scoliotic deformity. The Cobb angle is the maximum angle made in the frontal plane by two lines drawn parallel to the endplates of scoliotic vertebrae. Surgical treatment with implant fixation is used when the Cobb angle is more than 50° [1]. There are several surgical procedures that have evolved to be more advanced in applying the three-dimensional corrective forces for the correction of severe scoliotic deformity. Cotrel and Dubousset’s (CD) rod derotation technique, ventral derotation spondylodesis (VDS), Halm–Zielke instrumentation (HZI), simultaneous dual rod rotation method (SDRRT), direct incremental segmental translation (DIST), and others are examples of these procedures [1, 2]. These surgical therapies for severe scoliosis need the use of surgical methods to secure implanted devices such as rods, screws, hooks, and wires. For example, the anterior single rod correction procedure in Figure 1 involves removing the malformed intervertebral discs, implanting material to stimulate intervertebral joint fusion, and attaching metal rods to the spinal vertebrae [3] as spine deformity correction procedure.

The usage of computer simulations of the spine has skyrocketed in the last two decades. Computational techniques, notably finite element method (FEM), have previously been demonstrated to be effective in analysing the mechanics of the scoliotic spine during surgery [3]. For instance, Aubin et al. [4] created a model of the spine that had rigid bodies that represented the thoracic and lumbar vertebrae seen on intraoperative radiographs, as well as flexible parts that represented the intervertebral structures. Developing an optimal system of corrective force for an individual scoliosis patient via trial and error during surgery is actually unrealistic. Hence, FEM has been used to simulate surgical process changes and predicts the three-dimensional outcome in terms of deformity treatment and build flexibility. The FEM used in scoliosis research was divided into four groups based on model complexity: reflective simple variants of beam element-based models, representative complicated versions of beam element-based models, representative segmental volumetric models, and representative extensive volumetric models, according to the previous study [5].

By developing FEM, the relationships between the magnitude of corrective forces, number of screws, screw placement configuration, and degree of correction can be further elucidated [2]. Actually, the well-developed FEM of the spine allows for more complete assessments of internal stress distribution [6]. The studies of internal stress distributions and hypothetical scenarios cannot be properly examined without the help of real patient records, which is a possible drawback of rigid body modelling and patient-based FEM design.

Scoliosis correction aims to distort and restore the scoliotic spine to its original shape without inflicting injury or neurological complications. This could be accomplished by using implant rods and screws to impart appropriate correction forces to the spine. To obtain the optimum correction, the corrective forces required to rectify the deformity must be adequate [6]. Despite the fact that the number of patients was small, there was a growing correlation between the applied force and the degree of correction. Nevertheless, because of the rotating device is only linked to the implant rod, the correction forces acting at every screw were quite hard to be measured.

2. Methods

2.1. Search Strategy

Database search through the internet performed in November 2021 was restricted to the last twenty years of publication in Scopus (2001-2021), ScienceDirect (2001-2021), WorldCAT (2001-2021), PubMed (2001-2021), and Medline (2001-2021). The key MESH terms included “correction,” “spine,” and “deformity.” A thorough search was also conducted using the additional keyword search query: “finite element.” This extra manual search was implemented by manual screening conducted for relevant articles based on the reference lists of the retrieved articles. To rule out the likelihood of those items being overlooked, an additional search was conducted. All of the articles were retrieved and checked to ensure that the database search results were accurate and related to other articles. The ultimate decision was drawn in order to keep our findings within the scope of the papers.

2.2. Inclusion/Exclusion Criteria

From the electronic databases, only full-text articles in English were chosen. If there was a disagreement about an article during the screening process, it was debated to reach a consensus. Articles on studies focusing on corrective deformed scoliosis were included in the titles and abstracts screening procedure. The following criteria were used to evaluate articles: (1) rod and pedicle screw system, (2) Cobb angle, and (3) finite element analysis, full nonclinical articles. There were no restrictions on the participants’ age, gender, BMI, or medical history. Articles by the same author were double-checked to ensure there were no duplicates.

2.3. Review Process

Two reviewers (K.G. and K.S.B.) screened the search results according to the inclusion criteria. After the screening procedure, the final articles were extracted and separated from duplicated articles in several databases. Duplicate articles were eliminated from various databases. The eligibility criteria were initially applied to the titles and abstracts of the papers that were chosen. A full-text review was undertaken if the title and abstract did not give enough information in the article screening procedure. To avoid misinterpretation, rejected items were rescreened.

2.4. Assessment of Methodological Quality

The papers were reviewed and analysed using a systematic quality approach that helped to evaluate the quality of the articles identified, as well as extracting the most relevant information from those publications. There are no standardised methods for evaluating the reliability and credibility of each of the articles examined aside than the processes provided in this study. To evaluate credibility, 13 questions were adapted from Azizan et al. [7] and Ku et al. [8]. Few questions were had to drop as they are not justifying to spine deformity and its correction mechanism. At the same time, some of the questions were later altered according to FEA. Each question was given a score of “2” if the answer satisfied the standard questions, and a score of “1” if the information was limited. If no information was provided, the questions were marked with “0” or “no,” while questions that were not applicable were marked with “NA.” The 13 questions are as follows: (1)Is the study’s objective stated in a clear and concise manner?(2)Is the study design outlined in detail?(3)Are the patients’/models’ characteristics and details clearly provided?(4)Is the process of geometrical model development clearly explained?(5)Is the convergence test in the study clearly stated?(6)Is the boundary condition clearly described?(7)Are the appropriate mathematical models used to calculate parameters?(8)Are the mechanical properties of the model distinctly defined?(9)Are appropriate numerical methods used in data analysis clearly defined?(10)Is the predicted numerical value appropriately verified or validated?(11)Is the key outcome measure mentioned clearly?(12)Are the study’s limitations disclosed clearly?(13)Is the study’s conclusion conveyed in a clear and concise manner?

3. Results

3.1. Primary Search Results

Since the number of findings was limited in terms of quantity of information and the amount of materials available, the authors then conducted full-text reviews of articles. After a thorough screening process, twenty-two retrieved articles were finalised. Figure 2 shows the systematic review process of the present study. A total of 967 items were found after the database screening procedure. 81 of these items, however, were found duplicates and were removed. The relevancy of the studies undertaken was determined by looking at the titles and abstracts, and then, 602 publications were eliminated. Additional screening was carried out by reading the rest of the article’s content in order to ascertain the study’s goal based on the standard parameters that was assessed. After removing another 42 articles, there were 15 articles that meet all the criteria and were further reviewed.

3.2. Analysed Data Quality

Table 1 shows the quality ratings of the 15 articles that were assessed. The reviewed papers have a quality score ranging from 80% to 96%. Articles with a score of more than 85% are regarded good because they provide answer for all thirteen questions. These publications offered detailed information about their objectives, study design, key findings, and conclusion [1, 3, 9]. The remaining eleven studies fulfilled at minimum 70% to 90% of the questions [6, 10]–[11]. Other elements that could have aided in the comprehension of the questions were not considered in this review.

3.3. Participant Characteristics

Table 2 shows a list of physical characteristics and anthropometric factors from 15 articles. The majority of the participants were young people, with only a few adults of average age. Nine studies involved adolescences (aged between 10 and 19) [1, 3, 9, 11, 12, 14, 15, 17], two articles assessed middle-aged individuals (age range from 30–59) [10, 19], and none of the articles involved elderly persons (aged 60 and above). The number of participants in the evaluated publications varied, with the highest number being 20 peoples and seven articles keeping one as patient data for their investigation [6, 9, 11, 12, 16, 18, 19]. Four articles provide no anthropometric data of tested participants as those studies used spine finite element models [6, 16, 18, 20]. The participants were also categorised into presurgical [1, 2] and postsurgical patients [3, 12].

3.4. Finite Element Modelling

By conducting an exhaustive search for all published papers, a systematic review aims to minimise the incidence of bias. The modelling parameters in FEA, such as loads and boundary conditions, element types and sizes, geometrical model, type of material, and mechanical properties, have a significant influence on prognostic accuracy. These variables have an impact on the simulation’s predictive accuracy and should be considered when evaluating the results and conclusions. Table 3 shows finite element modelling variables that were used by the reviewed articles on the corrective mechanism of spine deformity due to scoliosis. These data can provide additional information on the simulations, allowing for replication and comparison.

Ten out of fifteen studies used ANSYS software for FEA studies [1, 2, 6, 9, 15, 18, 19]. Whereas Abolaeha et al., Little et al., and Guan et al. used ABAQUS software [3, 12, 18]. Wang et al. mentioned only radiographic software used for his FEA studies [20]. Special mention to Chen et al. [19] where Solidworks was the primary software for FEA studies where the pedicle screws created based on the imported blueprint then into Hypermesh and assembled into scoliotic spine FE model. Elements are the basic building block of FEA. Several authors used 10 node tetrahedral solid elements [1, 2, 9, 11, 14]. However, Wang et al. [15] and Guan et al. [18] used hexahedron elements. Moreover, beam elements for vertebral and pelvic sections, tension-only cable elements for ligaments, surface contact elements for articular facets, and modified beam elements are all included in FE model by Dumas et al. [10].

In order to produce an accurate reconstruction, the image quality is critical. To develop a rebuilt model as accurate as possible, it can be seen that the models in most of the studies were created from computed tomography (CT) scan images that was widely used by the researchers [1, 3, 6, 912, 15, 16, 18], either directly or via a previous model’s scaling. Meanwhile, only Wang et al. [13] used ADAMS 2005 software (Mechanical Dynamics Inc., Ann Arbor, MI) and the ADAMS Software Development Kit (SDK). Finite analysis has been shown to be useful in understanding the aetiology of scoliosis in these models. To ensure that the model performs realistically, the produced FEM must be tested against existing scientific results on motion and material properties.

As for the material, titanium alloy was widely used for both rod and screws compared to cobalt chromium which also has elastoplastic properties [1, 3, 12, 15, 18, 19]. Abolaeha et al. [12] used stainless steel that also has elastoplastic material behaviour as aluminium alloy. Some authors used both aluminium alloy and cobalt for screws and rod, respectively [3, 17]. As for the mechanical properties, the range of elastic modulus for rod is between 105 GPa to 213 GPa [1, 3, 6, 15, 17, 20].

Modern segmental spinal instrumentation systems are used to execute a variety of deformity correction techniques, including vertebral translation, rod derotation, direct vertebral derotation, and in situ rod contouring, all of which have been thoroughly studied. Dumas et al. [10] introduced situ contouring technique in his studies. Whereas Abolaeha et al. [12] used spinal growth rod instrumentation, an early onset scoliosis management method, and emerging technology that treating scoliosis without fusion hold the exciting prospect of a new paradigm in spinal deformity care. Salmingo et al. in their both studies and Abe et al. [1, 2, 14] used simultaneous double rod rotation technique (SDRRT) surgical technique with rods and polyaxial pedicle screws. Wang et al. [13] studied a posterior instrumentation using monoaxial pedicle screws, whereas Chen et al. [19] presented a procedure that includes rod derotation on the concave side and rod implantation on the convex side for strengthening.

Setting boundary condition is an important procedure that has to be done while set up simulation process of the finite element analysis. Dumas et al. [10] set displacement to the model as one of the boundary conditions. Displacement by means, the 3D motion of T1 with regards to the pelvis was computed between the bending test measurements and the standing position measurements. Force that could be directly applied was also set as the boundary conditions of the surgical manoeuvre.

There are several studies that have set the boundary conditions as constraining motions of particular vertebra. Balamurugan et al. [16] set the L5 vertebra as constrained from all degrees of freedom, and a compressive force of 50 N along the -axis was applied on T1 vertebra for the analysis. On the other hand, in the Salmingo et al. works [1, 2], the boundary condition was set considering the manner of rod fixation during the surgical treatment. The screws’ coordinates were reoriented such that the most superior screw coincides with the -axis (located on top of the most inferior screw) because each patient has different implant rod orientation and fixation levels. The most inferior screw at the end of the rod was fixed in all translations but free to rotate. The most superior screw was also fixed except that it was free to move along the superior direction only. The same practice was applied by Guan et al. [17] in their research where the boundary conditions were included as a fixed pelvis in rotation translation, and was limited to transverse plan movements. Zhang et al. [9] also set constraint on the displacement and rotation of all nodes on the base of the L5 vertebral body in all directions. However, there is no thoracic regions vertebra was included to set as boundary conditions.

Apart from abovementioned boundary conditions, the contact between rod and screw surface was also introduced by Little et al. [3] in their research. A “no separation” normal contact and frictionless tangential contact definition were defined between the screw head and the surface along the rod during surgery.

Researchers, on the other hand, have developed a variety of approaches to model the loading circumstances and limitations that are relevant to the corrective spine deformity process. Commonly, to achieve the desired correction, the force required to rectify the deformity must be significant. After a set of iterations with the force optimization method, the corrective forces acting on the implant rod were obtained [1]. Salmingo et al. [1, 2] analysed forces of screws set with initial values of zero before surgery on each screw’s matching point on the rod geometry. According to Wang et al. [15], Balamurugan et al. [17], Abolaeha et al. [12], and He et al. [11] studies, load was applied to the upper region of the vertebrae for observation of stress distribution. This was applied as boundary condition.

3.5. Variability in Measured Parameters

The output and findings of the parameters are summarised in Table 4. The focus of this data analysis is to look at the relevant biomechanical criteria that are often used to identify FEA and/or have clinical value.

Most of the articles focused on the influence of Cobb angles which is to indicate magnitude of spine deformity except for two articles which are Balamurugan et al. [17] and He et al. [11] that concentrated on the effect of surgery on deformity treatment in a scoliotic spine. Preoperative and postoperative main curves were described in four different investigations after follow-up period, and the degree of correction varies from 14° to 70° [1, 3, 12].

In this review, scoliosis in several planes’ views such as lateral, sagittal, axial, frontal, transverse, and coronal can be observed. Most of the authors studied scoliosis deformity in both sagittal [2, 9, 12, 13, 18] and coronal plane [3, 9, 11, 15, 16, 18]. However, Salmingo et al. and Zhang et al. focused on frontal plane [1, 9]. Figure 3 shows an example of a lateral displacement of the spine from the midline in the coronal (frontal) plane, decreased curvature in the thoracic region in the sagittal (side) plane and rotation in the axial plane.

Most of the authors carried out force analysis during the treatment of the spinal deformity and growth periods, on the rods and the spine for their studies [1, 3, 6, 12, 14, 16, 19]. Few of them created stress profile to understand the stress concentration profile on the vertebra and disc under different loads [1, 9, 11, 17]. Meanwhile, Dumas et al. [10] focused on rod rotation analysis on lateral, sagittal, and axial plane, and Wang et al. [15] demonstrated ranges of motion for L2 to L5 under various loading scenarios.

4. Discussion

4.1. Quality of Search

The aim of this systematic review was to analyse the biomechanical characteristics and parameters that are typically used in finite element analysis to investigate the corrective mechanism of scoliosis-related spine deformity. Understanding the corrective mechanism requires a comprehensive analysis on the parameters used in each investigation. In the present study, fifteen articles were included for the extensive review. In the reviewed studies, participants’ characteristics, Cobb angle, pedicle screw systems, and biomechanical responses can be further discussed.

4.2. Effect of Deformity Angles on Spine Corrective Forces and Stresses

The simulated corrected Cobb angle is usually attributed to the clinically established postoperative Cobb angle in the immediate postoperative period. This could provide the models’ accuracy in forecasting the change in coronal deformity following surgery. The Cobb angle which is used for comparison of deformity level is the maximum angle that can be projected between the upper and lower endplates of the scoliotic curve. To treat and prevent worsening deformity in severe cases of scoliosis (Cobb angle more than 45°-50°), surgical instrumentation or, in some cases, spinal fusion is sometimes utilised. Developing an optimum method of corrective force and predicting surgically imposed contact stresses between adjacent vertebral endplates for scoliosis patient during surgery through practical experiment is quite difficult. Hence, FEA can be used to model different surgical procedures and anticipate the three-dimensional results in the form of deformity correction and construct flexibility. Table 3 covers the variables involved in FEA studies.

Understanding and analysing the forces acting on the spine, as well as the interrelationship between both the force and the Cobb angle, will enable us to advance with improved systems [12]. According to Abolaeha et al., the resultant Cobb angles are inversely proportional to the progression of growth, rod lengths, and distraction force during a two-year period. To assess the forces required to treat scoliosis, Salmingo et al. [1] created an elastoplastic FEM. Based on differences in implanted geometry before and after surgery, the three-dimensional forces required to deform a rod were calculated. The instrumented spine’s at the lowest level experienced the highest forces.

Using the same FE model, Salmingo et al. continued to study the relationship between the magnitude of corrective forces and the degree of correction, which they measured changes of Cobb angle [2]. Actually, these values can be obtained by calculating the difference between preoperative and postoperative Cobb angle. They claimed that the degree of correction and the corrective forces operating on the rod were unrelated too. They also suggested that other factors, such as screw implantation arrangement and spine rigidity, may be linked to scoliosis repair.

However, from the study carried out by Little et al. [3] revealed that increasing the simulated intraoperative forces caused the anticipated corrected Cobb angle to decrease. Force, geometry (human anatomical), and tissue stress are involved in coronal plane deformity treatment. These are the most significant considerations in getting the best possible correction for a patient with the least amount of risk of high stresses on the spinal tissues, which could lead to implant-related problems.

After adding a growing rod, Abolaeha et al. [12] created a scoliotic spine FEM to model the spine growth over a two-year period. Based on the analysis, the pressures required to induce the correct Cobb angle changes are identical to those seen in patients. To distribute the load, the rod was linked to both vertebrae in the pair, which is identical to the present surgical insertion process. At the time of the original operation, it was expected that the deformity angle would be rectified by 50%. Following the initial operation, invasive lengthening treatments (similarly referred as distractions) were conducted every six months over a year to keep up with the growth of the spine.

Meanwhile, Guan et al. [18] concluded that whenever 3D correction forces rose, the thoracolumbar segment’s Cobb angle steadily reduced, and the vertebral body’s rotation angle lowered as well. The combined force correction effects were even higher. When correction forces were applied, stress at intervertebral discs in the distorted region changed drastically. Essentially, during scoliosis surgery, corrective force cannot be applied to the spinal implant beyond the anchor holding strength limit. If the corrective force exceeds the anchor’s strength, the implant may break or the bone may fracture, resulting in “screw ploughing.” Destabilization of the spinal segment by releasing soft tissue or the facet joint could be more critical than using an excessive correction manoeuvre with rigid implant to avoid implant fracture or pedicle rupture during a more severe curve correction procedure [14].

4.3. Analysis of Pedicle Screws and Implant Rod System due to Spine Deformity

For the treatment of spinal malformations, pedicle screw fixation has become a common surgical instrumentation approach. The better bone-implant connection allows surgeons to diagnose more corrective movements and employ larger correction pressures when translating and derotating the deformed spine. Hence, pedicle screws and its placement strategies are important as it has a minor influence on the curve correction scoliosis treatment [14]. In recent years, fewer screws have been used in scoliosis surgery for cost considerations, and correction rates have been reported as being similar between the less density group and the high-density group. For instance, Salmingo et al. [2] carried out a study by increase in absolute number of implant screws which resulted in reduction of the magnitude of corrective forces and did not give a greater degree of correction, and it was hypothesised that additional screws might prolong the surgery and result in more blood loss for the patients.

On the contrary, Clin et al. [16] discovered that reducing the number of screws raised the postoperative stresses that each screw could withstand, but that the influence on potential problems has to be investigated further. In their investigation, they found that independent of screw type, both high-density and low-density implant designs achieved comparable coronal correction and shared corrective forces equally well. Increased degrees of freedom of the screw head were also discovered to decrease the potential to cure coronal deformity while generating reduced bone-screw forces.

Theoretically, a greater number of implantations might give higher correction forces, resulting in better coronal and sagittal plane correction rates. A screw–rod connection which provides degrees of freedom, on the other hand, may make it difficult to perform the desired manoeuvres [21]. In addition, other variables such as curve flexibility, surgeon-specific objectives, and procedures may also play a role in the contradictory findings [19]. Wang et al. [20] used three forms of screws namely monoaxial, polyaxial, and dorsoaxial pedicle screws for their study. At each step, external forces must be raised until the rods can easily lock into the screw head saddles, which is linked to minimum “true corrective forces” (TCF) and little to no “Extra Forces” (EF) available to deliver the desired correction. The results showed that the dorsoaxial screws allowed for the least amount of EF to be created while forcing in order to make it certain appropriate rod seating and locking at all pedicle screws for corrective deformity correction. Clin et al. [16] claimed that lowering implant density by 30% permitted almost same degree of coronal correction as a fully instrumented construct irrespective of pedicle screw type, but that the influence on potential complications has to be investigated further.

He et al. [11] also claimed as pedicle-screw-rod system (PSRS) has always been regarded as the gold standard for the scoliosis treatment even though it has its own limitations. PSRS has several advantages, including rigorous fixation of deformities, increased osseous fusion, and a lower rate of pseudarthrosis. Front and mid columns are protected by rigid fixing, which counteracts eccentric stress. In addition, the fusion rate of stiff fixing is higher than that of semirigid fixing or no fixing. High stiffness, on the other hand, promotes fast scoliosis stability and minimises the physiological stress on the deformed vertebra.

Despite the fact that in recent years, the quality of studies in this field has improved, this review underlines the present literatures lack of regular use of standardised measures of end results and methodologies for preoperative and postoperative assessment. This standard should be broadened to include procedures for classifying and reporting complications. For example, past research has shown that excessive correction forces might could result in implant or bone fractures, which could result in screw extraction from the vertebra. As a result, gaining a better knowledge of scoliosis correction biomechanics necessitates an examination of corrective forces acting on the deformed rod [1].

Once the implant has been installed, the stress is centred on the two ends of the vertebral body, the rod, and the pedicle screw, resulting in a stress shielding effect. The stress shielding reduces the pressure on the intermediate vertebral body, and it just may result in bone loss and osteoporosis. Another limitation is that the cephalic and posterior sides of fixed segments have their rotation centres shifted. Because of these disadvantages, some orthopaedics professionals and researchers recommend using biodegradable or internal fixing materials with a low Young’s modulus [11].

Another notable highlight is lack of studies by previous researchers which focus only one specific region either it is lumbar or thoracolumbar region. Most of the researchers either involves whole spine region or mostly thoracolumbar region for their studies ([1, 2, 9, 11, 12, 15, 17, 20]. Actually, between thoracic and lumbar vertebrae, there are changes in transverse process bone mass and anatomical structure [11]. As a result, more research involving patients’ thoracic regions is required. The search method was confined to English-language articles, which is a limitation of the study. To discover articles, only five databases were used, and it is possible that some articles were overlooked. For the missing relevant articles, a manual search was conducted.

5. Conclusions

The present review highlights fifteen articles related to corrective mechanism of spine deformity that is published from 2001 until 2021. The correction mechanism, pedicle screw, rod system, Cobb angle, and other variable characteristics related to scoliosis surgery on patients’ bodies were the subject of this review. The collected data were able to furnish basic details about the simulations as well as some variables that may affect the predictive accuracy of the simulation. However, insufficient information in certain aspects prevents the analysis of related measured variables. There are various aspects that associated to scoliotic patients such as muscular activation, spine rigidity, deformity severity, the amount of stress that an internal fixation could withstand, and inter-individual differences have yet to be investigated. Insufficient information prevents the analysis of related measured variables. Hence, to improve and provide a better knowledge of the finite element approach for the analysis of correction mechanisms of spine deformity due to scoliosis, further research is needed in the areas stated above.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they do have no conflicting interests which might affect the results of the study.

Acknowledgments

The authors would like to acknowledge the support from the Fundamental Research Grant Scheme (FRGS) under a grant number of FRGS/1/2020/TK0/UNIMAP/02/20.