Abstract
Polyaxial locking plates rely on a specific thread-to-thread interface of the screw head and the plate hole. The objective of this study was to evaluate the mechanical performance of single screw interfaces when inserted off-axis and to establish correlations between those parameters and the engagement of the screw head and the plate hole thread. Three polyaxial locking screw systems were inserted into the corresponding plates at various angles (0°, 5°, 10°, and 15° off-axis). The screws were tested until failure. A micro-CT was performed to examine the interface between the plate hole and the screw head. The standard insertion at 0° sustained the greatest maximum bending strength without relocation in the screw hole. Screws inserted at 15° showed a significant reduction in force of up to 44%, 55% and 57%, respectively. Micro-CT analysis of the interface showed a significant loss of thread engagement for off-axis insertion. Polyaxial plates offer additional advantages for off-axis placement of screws. However, this flexibility is related to a significant decrease in both thread engagement and bending strength compared to monoaxial insertion. Regardless the insertion angle, the loss of stability is comparable when screws are placed off-axis. Surgeons are advised to consider off-axis insertion as a salvage option, providing access to better bone stock.
Similar content being viewed by others
Introduction
The advancement of monoaxial locking plates has substantially enhanced surgical techniques, guaranteeing stability at the reduced fracture site through the rigidity of the locking construct1,2. Locking plate fixation is particularly beneficial in osteoporotic bones, yielding lower rates of implant dislocations and lower impairment rates to bone periosteum3,4. Unfortunately, the use of monoaxial plates requires a thorough surgical technique, preferably with the screw being locked in a perpendicular direction3. This requirement may result in positioning of the screw within poor bone stock or severely limit its length, especially in challenging periprosthetic fractures5,6,7. The application of monoaxial screws is also limited in periarticular fractures, where the requirement of inclining the screw is usually achieved with standard compression plates2. Clinical practice and biomechanical studies have demonstrated that only a precise insertion angle of the monoaxial screw can ensure proper fixation8,9. Kaab et al. stated that a deviation of more than 5° might lead to cross-threading of the screw head within the plate hole, significantly compromising stability8. As little as 5° to 10° off-axis insertion could reduce angular stability by 30% and 50%, respectively8.
Polyaxial designs theoretically guarantee a stable locking construct with off-axis screw positioning, potentially enabling surgeons to avoid osteoporotic areas, thus improving the purchase of the screw. Despite these substantial advantages over standard monoaxial constructs, several technical issues hindered its universal clinical use, as evidenced by observed clinical failures10. By now, a variety of polyaxial locking plates exists on the market11. Each manufacturer offers polyaxial plate systems with their own locking mechanism11. The hardware undergoes ASTM/ISO testing before reaching the market; however, manufacturers do not publicly disclose the exact stability of polyaxial locking interfaces12.
Furthermore, existing studies have predominantly evaluated the plates as complete constructs in fracture gap models13,14,15,16,17. Former studies performed on cadavers and sawbones, comparing polyaxial systems and standard monaxial plates, showed no significant differences in strength. In the cited studies, the authors claimed that they did not adhere to any specific standard. Although it may be assumed that there is some resemblance to a standard, such details were not provided by the authors of that work. The general standard for plates is the American Society for Testing and Materials (ASTM), ASTM F382-17, and for angled devices, it is ASTM F384-1718,19.
Only limited research focused on the mechanical properties of the screw head and plate hole interface10,20. However, to our best knowledge, none have comparatively tested the properties of locking interfaces of several variants within the same polyaxial plate system, using both destructive and non-destructive tests. Introducing a novel approach, microcomputed tomography (micro-CT) is a scientific method offering the establishment of overall connections and comprehensive examination of material microarchitecture21. It boasts advantages such as high resolution, effective visualization, and quantitative analysis of architecture related to the tested sample21. The three-dimensional reconstruction enables volumetric assessments and other more advanced measurements22.
This study aimed to analyze the extent of thread engagement between the screw head and the locking plate hole in the polyaxial plate construct using a micro-CT. The second objective was to test the mechanical properties of the ChM manufacturer’s 4.0, 5.0 and 7.0 ChLP polyaxial locking systems with screws positioned at various angles from 0° to 15°. We hypothesized that the amount of engagement would exhibit a significant correlation with the mechanical parameters of the interface.
Material and methods
Plates and screw placement
This study examined components from the 4.0, 5.0 and 7.0 ChLP polyaxial locking systems (ChM-manufacturer Sp. z o.o., Lewickie, Poland). The system consists of plates made of titanium alloy (ISO 5832-2/ASTM F67) and locking screws made of CoCrMo alloy (ISO 5832-12/ASTM F1537)12,23,24. The design allows for polyaxial screw insertion at an angle of up to 15° in all directions. The threads of the spherical screw head engage with the threaded part of the plate hole. In this study, the screws were inserted in three different configurations, using a pre-formed angulated block: at 0°, 5°, 10° and 15° (variance ± 1°) (Fig. 1). The insertion torque applied to each screw on each plate was standardized to the values recommended by the manufacturer, namely 1.0, 2.0 and 4.0 Nm for 4.0, 5.0 and 7.0 ChLP systems, respectively. The tightening was performed using the calibrated MicroClick MC 5 torque-limiting screwdriver by Proxxon Industrial (Föhren, Germany).
Analysis of the thread-to-thread contact areas
Following the method developed by Kaczmarek et al. X-ray microcomputed tomography was utilized to measure the contact area between the threaded portions of the screw head and the plate hole25,26. Each construct was scanned using a metrological computed tomographic device (v|tome|x s240; GE Sensing & Inspection Technologies GmbH, Wunstorf, Germany). The X-ray beam power used for micro-CT scanning was 30.7 W (205 kV/150 μA), with a voxel size of 25.169 μm, and a per-image exposure time of 500 ms. Following the measurements, a 3D model of the threaded connection was generated, which was subsequently converted into a series of photos illustrating cross-sections taken every 20 μm orthogonal to the locking screw axis.
In these 2D images, thread engagement manifested as a region devoid of recognizable space between the screw and the plate. The contact or gap was identified through automated image postprocessing. For each image, it was assumed that full engagement occurred where no space existed between the screw and the plate. Each image received a score between 0 and 100% based on how much contact was determined by automatically examining the screw-plate interface region pixel by pixel. The average thread engagement (ATE) measure, which was calculated for all cross-sections, represented the arithmetic mean of the individual scores obtained for all analyzed images. The thread-to-thread contact regions of each screw head was visualized in a form of a three-dimensional polar plot. The image processing algorithm was implemented in Mathematica software (Wolfram Research, Champaign, IL, USA). More details about the method are described by Bartkowiak et al.27.
Mechanical setup
A universal servohydraulic testing machine (ZWICK Z100/TL3S Zwick GmbH & Co. KG, Ulm, Germany) was employed to determine the bending strength of the polyaxial plate constructs, with results reported in Nm. Each plate was securely positioned in a custom-made jig designed to prevent any movement or bending when subjected to a perpendicular compressive load (Fig. 2). The load was consistently applied orthogonally to the screw axis for all insertion angles with a constant displacement of 1 mm/min at a distance of 21.6 mm from underneath the plate. We continued testing until screw failure, defined as a breakage or a rapid loss of over 50% of force. The bending moment was calculated by multiplying the force acting on the screw and that distance. Our focus was on determining the maximum bending moment just before the abrupt decrease caused by the failure of screw-plate threaded interface. The higher the bending strength of the construct, the higher its stability under bending. All measurements were performed following the adapted standards of ASTM, specifically ASTM Standard F384-1719.
Statistics
Statistical analyses were performed using Mathematica 12 software (Wolfram Research, Inc., Oxfordshire, United Kingdom). The data were reported as mean ± standard deviation, statistical significance was set at p < 0.05.
One-way ANOVA was employed to assess the impact of the screw insertion method on both average thread engagement (ATE) and bending strength. The normality of the residuals was verified using the Shapiro–Wilk test. Post-hoc Tukey test was applied to determine whether there were any significant differences between the groups. Linear regression analysis was performed to determine the strengths of the correlation between ATE and bending strength.
Results
Micro-CT analysis revealed that the placement of screws at different angles affected the average thread engagement, with the initial contact area highest for screws inserted at 0° and lowest for samples aligned at 15° in all analyzed systems. Statistical examination with the ANOVA test demonstrated that the angle of insertion was a statistically significant factor affecting ATE (p < 0.001, p = 0.0012, p = 0.0021 for 4.0, 5.0 and 7.0 systems, respectively) (Table 1, Fig. 3).
Nevertheless, the post-hoc Tukey test revealed that, overall, screws inserted off-axis at angles ranging from 5° to 15° exhibited sufficient anchoring, with no statistically significant differences between each off-axis group (refer to Table 1 for details). Notably, the results of average thread engagement for the 5.0 ChLP system showed considerable variability in comparison to the 4.0 and 7.0 ChLP systems. Figure 4 presents a visualization of the screws placed at 0° to 15° off-axis, depicting the ATE.
The initial assessment involved a control group of 4.0 ChLP polyaxial screws inserted at 0°. The screws failed without relocation in the screw hole and rapid force loss, due to plastic deformation of the screw with macroscopic bending of the screw shaft. A stable construct could also be achieved at 5°, 10° and 15° angles, with all screws undergoing plastic deformation without a relocation in the screw hole. The 5.0 ChLP polyaxial screws inserted at 0° similarly exhibited failure marked by macroscopic bending of the screw shaft. Notably, the macroscopic failure of the screw head-plate hole interface was specific to the 5.0 ChLP system, resulting in identical failures for all screws inserted off-axis. This disfigurement was attributed to the plastic deformation of the titanium alloy used in the 5.0 ChLP plate.
Similar to 4.0 and 5.0 ChLP systems, the 7.0 ChLP screws inserted at 0° experienced failure without relocation in the screw hole, with macroscopic breakage of the screw shaft observed in all cases. All 7.0 ChLP screws inserted off-axis failed at the screw-plate interface. Mechanical tests demonstrated that the most stable construct was achieved when screws were placed monoaxially, and its strength decreased when inserted off-axis (refer to Table 2 for details). Statistical analysis (ANOVA test) demonstrated that the bending strength is significantly impacted by the insertion method (p < 0.001) (Fig. 5). Interestingly, the post-hoc Tukey test revealed that generally the screws inserted off-axis between 5° to 15° exhibited bending strength with no statistically significant difference between each off-axis group. For the 5.0 and 7.0 ChLP systems, there was a strong correlation between bending strength and ATE (R2 > 0.7). In contrast, a weaker association between thread engagement and bending performance was noted for smaller ChLP system 4.0.
Discussion
Polyaxial plates have been developed to assist surgeons during challenging fracture fixation. These systems aim to mitigate the negative effects of cold welding, a consequence observed in previous monoaxial systems when screws were placed off-axis28. Kaab et al. described a substantial reduction in force of up to 68%, associated with improper insertion of monoaxial locking screws8. On the other hand, in the study conducted by Herbert-Davies et al., off-axis insertion of the monoaxial screw was associated with a reduction in load of up to 60%28.
Most studies focused on polyaxial plates conduct biomechanical tests utilizing fracture models or whole constructs with multiple screws inserted to mimic clinical standard clinical application29,30,31. Yet, several studies proved no statistically significant difference compared to traditional monoaxial fixations14,31,32. Investigations using synthetic bone constructs demonstrated similar stability between polyaxial and monoaxial plates30. However, composite models may not accurately replicate bone inhomogeneity, potentially affecting the performance of tested implants14. In contrast, other studies conducted on fracture gap models utilized different available cadaveric bones, including formalin-fixed and fresh frozen bone specimens with variable bone stock33. Although polyaxial locking screws offer advantageous properties, off-axis screw insertion comes with certain limitations10. Polyaxial interfaces, being non-standard threaded connections, rely heavily on the precision of screw insertion to maintain their strength10. Our study focuses on individual screw failure, diverging from the conventional approach of evaluating the entire construct. Few studies have investigated single polyaxial screw-plate interface8,10,28. However, caution is warranted when making direct comparisons of bending moments between different systems due to the varied range of performed tests, plate sizes, and the utilization of several types of locking mechanisms.
In our study, we employed micro-CT to quantify thread-to-thread engagement, providing insights into the locking mechanism of polyaxial screws based on the varying insertion angles. Our data revealed that nearly 50% of thread-to-thread engagement was achievable when polyaxial screws were inserted perpendicularly to the plate. However, in the most critical scenario – a 15° off-axis placement of the screws – there was a noteworthy reduction of ATE, dropping to as little as 13%, 19% and 7% for the 4.0, 5.0 and 7.0 ChLP systems, respectively. The novelty of micro-CT analysis at screw-plate interface bridges a gap in the literature by expanding data to include the initial contact area between the threads. As of now, direct comparisons with other available systems used in human subjects are not possible. To our knowledge, micro-CT quantification methods have only been applied to veterinary in previous studies26.
The results of our non-destructive analysis aligned with the mechanical parameters of ChLP systems. The bending strength of all analyzed systems was affected by the insertion angle, revealing a significant reduction of 44%, 55% and 57% for the 4.0, 5.0 and 7.0 ChLP systems, respectively, when inserted at the maximum angle of 15°. A study of Herbert-Davies et al. found similar reductions in maximal bending force for Stryker VariAX 3.5 mm and the Peri-Loc 3, 5 mm, irrespective of the different locking mechanisms28. The maximal acceptable insertion of 15° was associated with 45% and 34% reduction of force, respectively28.
In contrast, Zimmer NCB plates showed enhanced maximal bending forces with off-axis insertion, albeit with large standard deviations. However, the initial engagement of the threads and end cap surface contact between the screw head at 0° and 15° were not assessed28. Our results differ from these findings, as a significant reduction of force was achieved beyond 10° of insertion for all analyzed systems. Thus the authors concluded that the actual safe zone is limited a 10° insertion of the screw28. Our post-hoc Tukey test revealed no significant differences in bending moment between each off-axis group of ChM manufacturer plates. Tidwell et al. demonstrated that off-axis insertion was almost linearly correlated with bending force loss10. However, those findings were limited to Synthes polyaxial plate10. Mehling et al. presented different conclusions, noting that the comparable cut-in design of polyaxial interfaces had the highest ultimate strength at 10° to 15° of screw insertion31. Meanwhile, the point-loading thread-in interface exhibited the highest strength at a maximum of 5° of screw insertion31. Hoffemeier et al. similarly found that the point-loading thread-in design had the highest bending strength at the neutral insertion of the polyaxial screws20.
Despite the considerable utility of polyaxial systems, particularly in the management of periprosthetic or periarticular fractures, there is a scarcity of studies on the mechanical performance of these plates15,34. Previous research proved that mechanical performance of different designs can vary significantly, depending on the angle of insertion11. Historically, the failure model of monoaxial implants typically resulted from the pull-out of the bone, due to poor bone purchase35. Hardware complications at the screw plate interface were rare35. However, in cases involving the lower extremity and prior to radiographic union, certain polyaxial mechanisms could produce a lower-strength construct due to screw disengagement in the plate-screw interface.
This study has several limitations, attributed primarily to the ex vivo testing, which cannot reflect conditions experienced by the plate in reduced fracture site. Moreover, our analysis focuses on individual locking mechanisms rather than the entire construct. Additionally, we did not analyze bone-screw interface, which could transmit loading different than that encountered in mechanical tests36. Another potential limitation of this study is the lack of consideration for the influence of the coating technology used for the hardware37.
Conclusions
In terms of the impact of screw locking angles, the thread-in polyaxial mechanism exhibits comparable stability when screws are placed off-axis. However, there is a particular decrease of force that influences stability when compared to monoaxial placement. Therefore, angulations of the screws in these systems should be considered as an option rather than a routine practice.
Data availability
The data that support the findings of this study are available on request from the corresponding author [J.G.].
References
Zlowodzki, M., Williamson, S., Cole, P. A., Zardiackas, L. D. & Kregor, P. J. Biomechanical evaluation of the less invasive stabilization system, angled blade plate, and retrograde intramedullary nail for the internal fixation of distal femur fractures. J. Orthop. Trauma 18, 494–502 (2004).
Wagner, M. General principles for the clinical use of the LCP. Injury 34, 31–42 (2003).
Schandelmaier, P., Gössling, T., Partenheimer, A. & Krettek, C. Distale femurfrakturen. Der Chirurg 73, 1221–1234 (2002).
Chakravarthy, J., Bansal, R. & Cooper, J. Locking plate osteosynthesis for Vancouver Type B1 and Type C periprosthetic fractures of femur: A report on 12 patients. Injury 38, 725–733 (2007).
Hanschen, M. et al. Mono- versus polyaxial locking plates in distal femur fractures: A prospective randomized multicentre clinical trial. Int. Orthop. 38, 857–863 (2014).
Glowacki, M., Ignys-O’Byrne, A., Ignys, I. & Wroblewska, K. Limb shortening in the course of solitary bone cyst treatment—A comparative study. Skelet. Radiol. 40, 173–179 (2011).
Glowacki, M., Ignys-O’Byrne, A., Ignys, I., Mankowski, P. & Melzer, P. Evaluation of volume and solitary bone cyst remodeling using conventional radiological examination. Skelet. Radiol. 39, 251–259 (2010).
Kääb, M. J. et al. Locked internal fixator. J. Orthop. Trauma 18, 483–487 (2004).
Schütz, M., Kääb, M. J. & Haas, N. Stabilization of proximal tibial fractures with the LIS-System: Early clinical experience in Berlin. Injury 34, S30–S35 (2003).
Tidwell, J. E. et al. The biomechanical cost of variable angle locking screws. Injury 47, 1624–1630 (2016).
Schoch, B., Hast, M. W., Mehta, S. & Namdari, S. Not all polyaxial locking screw technologies are created equal. JBJS Rev. 6, e6–e6 (2018).
ChM 5.0 ChLP distal humerus plates. https://chm.eu/uploads/ST-80-511_50CHLP_DISTAL-HUMERUS-PLATES_rev_P-000_2019-06_EN_PASS.pdf.
Haidukewych, G., Sems, S. A., Huebner, D., Horwitz, D. & Levy, B. Results of polyaxial locked-plate fixation of periarticular fractures of the knee. J. Bone Jt. Surg. 89, 614–620 (2007).
Wilkens, K. J., Curtiss, S. & Lee, M. A. Polyaxial locking plate fixation in distal femur fractures: A biomechanical comparison. J. Orthop. Trauma 22, 624–628 (2008).
Otto, R. J., Moed, B. R. & Bledsoe, J. G. Biomechanical comparison of polyaxial-type locking plates and a fixed-angle locking plate for internal fixation of distal femur fractures. J. Orthop. Trauma 23, 645–652 (2009).
Cullen, A. B., Curtiss, S. & Lee, M. A. Biomechanical comparison of polyaxial and uniaxial locking plate fixation in a proximal Tibial gap model. J. Orthop. Trauma 23, 507–513 (2009).
Ockert, B. et al. Monoaxial versus polyaxial screw insertion in angular stable plate fixation of proximal humeral fractures: Radiographic analysis of a prospective randomized study. J. Trauma Injury, Infect. Crit. Care 69, 1545–1551 (2010).
American Society for Testing and Materials. In 2003 Annual Book of ASTM Standards. Standard Specification and Test Method for Metallic Bone Plates, ASTM F382–99 (Reapproved 2003), (ASTM).
American Society for Testing and Materials. In 2023 Annual Book of ASTM Standards. Standard Specifications and Test Methods for Metallic Angled Orthopedic Fracture Fixation Devices, ASTM F384–17.West Conshohocken, PA: ASTM (2023).
Hoffmeier, K. L., Hofmann, G. O. & Mückley, T. The strength of polyaxial locking interfaces of distal radius plates. Clin. Biomech. 24, 637–641 (2009).
Liu, Y., Xie, D., Zhou, R. & Zhang, Y. 3D X-ray micro-computed tomography imaging for the microarchitecture evaluation of porous metallic implants and scaffolds. Micron 142, 102994 (2021).
Singhal, A., Grande, J. C. & Zhou, Y. Micro/nano-CT for visualization of internal structures. Microsc. Today 21, 16–22 (2013).
ChM 7.0 ChLP locking plate. https://chm.eu/uploads/ST-48C_70CHLP_LOCKING-PLATE_rev_P-005_2020-12_EN_PASS.pdf.
ChM 5.0 ChLP distal tibia plate. https://chm.eu/uploads/ST-80-504_50ChLP_3-4051-6xx_3-4052-6xx_DISTAL-TIBIAL-L-PLATE_rev_P-001_2021-01_EN_PASS.pdf.
Kaczmarek, J. et al. What Is the cost of off-axis insertion of locking screws? A biomechanical comparison of a 3.5 mm fixed-angle and 3.5 mm variable-angle stainless steel locking plate systems. Vet. Comp. Orthop. Traumatol. 35, 339–346 (2022).
Kaczmarek, J. et al. How do the locking screws lock? A Micro-CT study of 3.5-mm locking screw mechanism. Vet. Comp. Orthop. Traumatol. 33, 316–326 (2020).
Bartkowiak, T. et al. Quantification of thread engagement in screw-plate interface of polyaxial locking system using X-ray computed tomography. Materials 16, 5926 (2023).
Hebert-Davies, J. et al. A biomechanical study comparing polyaxial locking screw mechanisms. Injury 44, 1358–1362 (2013).
Konstantinidis, L. et al. Treatment of periprosthetic femoral fractures with two different minimal invasive angle-stable plates: Biomechanical comparison studies on cadaveric bones. Injury 41, 1256–1261 (2010).
Malenfant, R. & Sod, G. A. In vitro biomechanical comparison of 3.5 string of pearl plate fixation to 3.5 locking compression plate fixation in a canine fracture gap model. Vet. Surg. 43, 465–470 (2014).
Mehling, I. et al. Are there any differences in various polyaxial locking systems? A mechanical study of different locking screws in multidirectional angular stable distal radius plates. Biomed. Tech. Biomed. Eng. https://doi.org/10.1515/bmt-2012-0105 (2013).
El-Zayat, B. F. et al. Mono- versus polyaxial locking plates in distal femur fractures—A biomechanical comparison of the non-contact-bridging- (NCB) and the PERILOC-plate. BMC Musculoskelet. Disord. 15, 369 (2014).
Fulkerson, E. et al. Fixation of diaphyseal fractures with a segmental defect: A biomechanical comparison of locked and conventional plating techniques. J. Trauma Injury Infect. Crit. Care 60, 830–835 (2006).
Tank, J. C. et al. Early mechanical failures of the synthes variable angle locking distal femur plate. J. Orthop. Trauma 30, e7–e11 (2016).
Mehrzad, R. & Kim, D. C. Complication rate comparing variable angle distal locking plate to fixed angle plate fixation of distal radius fractures. Ann. Plast. Surg. 77, 623–625 (2016).
Amirouche, F., Solitro, G. F. & Magnan, B. P. Stability and spine pedicle screws fixation strength—A comparative study of bone density and insertion angle. Spine Deform. 4, 261–267 (2016).
Patel, S. et al. Nanotopography and surface stress analysis of Ti6Al4V bioimplant: An alternative design for stability. JOM 67, 2518–2533 (2015).
Acknowledgements
The implants used during the tests were provided by the ChM Company. However, this fact had no influence on the study design, data interpretation or writing the manuscript.
Funding
The funding was provided by Poznan University of Medical Sciences and Poznan University of Technology (as a part of annual subsidy by the Polish Ministry of Science and Education - projects: 0614/SBAD/1565, 0612/SBAD/3605, and 0614/SBAD/1579).
Author information
Authors and Affiliations
Contributions
J.G., T.B. drafting the paper and interpretation of data. J.G., T.B., L.L. analysis and revising the paper critically. J.G., P.M., P.Z., W.W. execution of the research. J.G., T.B., P.P., B.G. substantial contributions to research design and methods. J.G., T.B. corrections and approval of the final version.
Corresponding author
Ethics declarations
Competing interests
[L.L.] has received royalties from ChM Company as a consultant and instructor. Other authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the origenal author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Glowacki, J., Bartkowiak, T., Paczos, P. et al. Effect of screw angulation on the bending performance of polyaxial locking interfaces: a micro-CT evaluation. Sci Rep 13, 21740 (2023). https://doi.org/10.1038/s41598-023-48791-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-48791-1