Activated Amorphous Carbon With High-Porosity Derived From Camellia Pollen Grains as Anode Materials for Lithium/Sodium Ion Batteries.

Carbonaceous anode materials are commonly utilized in the energy storage systems, while their unsatisfied electrochemical performances hardly meet the increasing requirements for advanced anode materials. Here, activated amorphous carbon (AAC) is synthesized by carbonizing renewable camellia pollen grains with naturally hierarchical structure, which not only maintains abundant micro- and mesopores with surprising specific surface area (660 m2 g-1), but also enlarges the interlayer spacing from 0.352 to 0.4 nm, effectively facilitating ions transport, intercalation, and adsorption. Benefiting from such unique characteristic, AAC exhibits 691.7 mAh g-1 after 1200 cycles at 2 A g-1, and achieves 459.7, 335.4, 288.7, 251.7, and 213.5 mAh g-1 at 0.1, 0.5, 1, 2, 5 A g-1 in rate response for lithium-ion batteries (LIBs). Additionally, reversible capacities of 324.8, 321.6, 312.1, 298.9, 282.3, 272.4 mAh g-1 at various rates of 0.1, 0.2, 0.5, 1, 2, 5 A g-1 are preserved for sodium-ion batteries (SIBs). The results reveal that the AAC anode derived from camellia pollen grains can display excellent cyclic life and superior rate performances, endowing the infinite potential to extend its applications in LIBs and SIBs.

Introduction

The emergence and rapid development of portable electronic devices, electric vehicles and renewable energy industries urge preeminent large-scale energy storage and conversion systems, which exhibit superior energy density and specific capacity as well as phenomenal life-span with trustworthy safety (Hannan et al., 2017; Zhang et al., 2017). Among prevalent energy storage systems, lithium ion batteries (LIBs) and sodium ion batteries (SIBs) attract most attentions owing to their high working potential and extraordinary storage capacity (Shen et al., 2018). It is well-known that in order to promote practical applications of LIBs and SIBs, developing high-performing electrode materials are inevitable (Zheng et al., 2015; Luo et al., 2016). With increasing improvements in cathode materials, commercial graphite as the mainstream anode material gradually becomes dissatisfactory due to its limited theoretical capacity (372 mAh g−1 for LIBs) (Hou et al., 2015; Zhang et al., 2016b). Besides, the highly regular distance between graphene layer (0.335 nm) can hardly accommodate the size of Na+, prohibiting the possibility to form stable Na-intercalation compounds (Irisarri et al., 2015; Yang et al., 2017). Though graphite hinders the improvement of LIBs and SIBs, carbonaceous materials still consumingly appeal to researchers, thanks to their significant advantages of low lithium/sodium ion insertion/desertion potential, highly structural stability, cost-efficiency, and sufficient resources (Hou et al., 2017a,b).

Amorphous carbon, compared to graphite, can reach higher specific capacity with enlarged interlayer distance and shortened ion transportation paths. Owing to environment amiability, sustainability and cost efficiency, ta large amount of diverse selectable biomass materials, such as protein (Li et al., 2013), peanut shell (Lv et al., 2015), bee pollen grain (Tang et al., 2016), rice husk (Zhang et al., 2016a), cotton (Li et al., 2016; Xiong et al., 2018; Yang et al., 2018), garlic skin (Zhang et al., 2018c), have been investigated in depth of their electrochemical capabilities. Generally, biomass-derived amorphous carbon can be synthesized by directly pyrolysis (Li et al., 2016)−(Wang et al., 2015) or hydrothermal treatment (Li et al., 2015). It is noted that most of the biomass-derived carbon can successful maintain the uniform sphere shape, inherent textile and hollow tubes. Despite of advantages of biomass-derived carbon, it is convinced that its electrochemical potential has not been ultimately exploited. To address this challenge, researchers have exerted their efforts on constructing inner microstructures and modifying surface properties of biomass-derived carbon by utilizing the activation methods (Gu et al., 2014; Wu et al., 2016a). Noticeably, chemical activation by mixing carbon materials and chemical reagents, KOH for instance, is less energy-consuming than physical activation resulting from the relatively low anneal temperature (Gao et al., 2016).

Pollen grains, distinctively, presenting uniformly discrete particle distribution and divergent particle shapes from various kinds of flowers, are usually employed in synthesis of various materials as templates (Chen et al., 2016). When regarded as anode for LIBs, cattail pollen carbon after physically activating exhibits 382 mAh g−1 at 25°C at rate of 37.2 mA g−1, distinguishing itself as a competitive candidate in anode materials beyond graphite (Tang et al., 2016). Impractical test current density is reason of the slightly better capacity, which still need further improvements. Therefore, further exploration of beneficiated activations for fabricating pollen grain's microstructure is required.

In this work, taking the advantage of pollen grains, we put forward a homogeneous and tunable synthetic method for activated amorphous carbon (AAC) derived from camellia pollen grains. Pollen grains were activated through chemical activation using KOH as activating reagent (Su et al., 2017), followed by a two-stage carbonization process. The eventual product is identified as poly-porous amorphous carbon with enhanced specific surface area around 660 m2 g−1 and enlarged layer distance from 0.35 to 0.4 nm, which contributes to advanced specific capacities and cyclic performances with excellent rate when acting as anode materials for LIBs and SIBs. The outstanding characteristics provide camellia-pollen-grain-derived carbon with the opportunity to meet the increasing demand of advanced storage devices.

Experimental

Material preparation

The activated amorphous carbon was synthesized by annealing camellia pollen grain. After rinsing with acetone and deionized water, 3 g raw pollen grains were mixed with 80 mL KOH solution (0.025 g mL−1) and stirred for 6 h at 40°C, following by drying process at 80°C. The obtained pollen grains were then transferred into a tube furnace and underwent the two-stage heat treatment under N2 atmosphere, which suggested maintaining temperature at 400°C for 1 h at the first stage and 800°C for 3 h at the second stage with the fixed heating rate of 5°C/min for the entire calcination. Finally, the sample was etched by 3M HCL and thoroughly washed with deionized water and desiccated at 80°C for 24 h. The as-prepared activated amorphous carbon was denoted as AAC.

For comparison, another 3 g raw pollen grains were directly annealed through two-stage heat treatment aforementioned to acquire pristine amorphous carbon, denoting as AC.

Material characterizations

Scanning electron microscope (SEM, FEIQuanta 200 FEG) and transmission electron microscope (TEM, Tecnai G2 F20 S-TWIN, Japan) were performed to unveil the morphology and microscopic structure of both samples. To gain the intrinsic features of the samples, X-ray diffraction (XRD, Bruker D8 Advance, Germany) (Cu, Kα, λ = 1.5405 Å), Raman spectrometer (JOBIN-Yvon HR800) and X-ray photoelectron spectroscopy (XPS, Thermo/ESCALAB 250XI) were conducted. N2 adsorption/desorption isotherms for calculation of specific surface area (SSA) and accumulative pore volume based on Brunauer-Emmett-Teller and density functional theory methods was carried out on ASAP 2020 Micromeritics.

Electrochemical measurements

The working electrode was fabricated using copper foil covered with slurry comprised of 70 wt % of AAC or AC, 20 wt % of PVDF binder and 10 wt % carbon black. Subsequently the electrode is dried out with 0.5 mg cm−2 average active material loading to assemble CR2032 coin cells together with PE/PP films or glass-fiber papers as separators, metal lithium or sodium foils as counter electrodes and corresponding electrolyte (1 M LiPF6 or 1 M NaCF3SO3 in ethylene carbonate, dimethyl carbonate and ethyl methyl carbonate = 1:1:1(v: v: v) for lithium ion batteries and sodium ion batteries, respectfully). CV and EIS results were attributed to the employment of CHI660A electrochemical workstations while galvanostatic charge/discharge profiles were resulted from the measurements of LAND CT2001A battery-testing instruments.

Results and discussions

The overall synthesis of AAC is shown in Figure 1. The KOH activation followed by carbonization and acid etching process are involved in this strategy to introduce poly-pores into the interior structures of camellia pollen grains (details in experimentally section). The microscopic morphologies of the original pollen grains, AC and AAC are analyzed via SEM, TEM and HRTEM images (Figure 2 and Figure S1). Camellia pollen grains are consisted of irregular spheres with rough surfaces, and the uniform diameter is about 30 μm (Figures S1A,B). In comparison with AC (Figures S1C,D), Figures 2A,B visualizes that, in morphology, AAC possesses smooth surface, vast amount of mesopores and micropores, while primitive bulk shape is well-maintained, proving the effective function of KOH agent as impetus aiming at high porosity. Consistency can be found in corresponding TEM images (Figure 2C), in which explicit micro- and mesopores morphology. Moreover, it is noticeable in HRTEM (the inset figure of Figure 2D) that AAC is comprised of turbostratic nanodomains with ~0.4 nm interlayer distance in average graphene interlayers. The evidence of dispersed diffraction ring in selected area electron diffraction pattern (SAED) can be observed, consistent with the HRTEM results, confirming the amorphous characteristic of AAC.

Figure 1

Schematic illustration of the synthesis process for the activated amorphous carbon derived from camellia pollen grains.

Figure 2

(A, B) SEM images and (C, D) TEM images of AAC (the inner figures are SAED and HRTEM).

To gain a deeper insight of the composition and crystal structure of both samples, XRD patterns, and Raman spectra are investigated as displayed in Figures 3A–C. The broad peaks in XRD pattern referring to (002) facet are detected at 2θ of 21.8° and 25° for AAC and AC, respectively (Figure 2A), confirming the amorphous property of both samples. It is noted that the left shift of 2θ for AAC is ascribed to the enlarged interlayer. According to Bragg equation (Ou et al., 2016), interlayer distance between graphene layer can be calculated to be 0.409 and 0.352 nm for AAC and AC, respectively, both larger than 0.335 nm for graphite (Hou et al., 2017a). Such appealing expansion for AAC should be ascribed to KOH activation mechanism, which is explained by three steps as followed: (i) Firstly potassium-containing compounds commence etching carbon substrate to form the initial porous framework, which is called chemical activation. (ii) Then H2O and CO2 escape to impel the increment of porosity for physical activation. (iii) Metallic K permeates into carbon matrices, achieving the larger graphene interlayer distance, resulting in the lattices expansion (Romanos et al., 2012; Wang and Kaskel, 2012). The enlarged interlayer distance has favorable influences on the electrochemical performance of AAC, which enables reversible intercalation for Li+/Na+ ion and well maintain the structural integrity at the same time (Wen et al., 2014; Kim et al., 2017; Zou et al., 2017; Yang et al., 2018; Zhu et al., 2018). Subsequently, D-band representing defects induction at ~1334 cm−1 and G-band representing in-plane vibration at ~1577 cm−1 are located in Raman spectra (Figure 3B), whereby the ID/IG can be calculated as 0.84 for AAC and 1.11 for AC. Furthermore, based on the ID/IG ratio, the AAC is of higher graphitization degree, providing a crucial foundation for improving electrical conductivity (Lv et al., 2015; Zhang et al., 2018c).

Figure 3

(A) XRD patterns, (B) Raman spectra, and (C) N2 sorption isotherms (inner figure is pore size distributions) of AAC and AC; (D) XPS survey scan spectrum of AAC; (E,F) high resolution XPS spectra of (E) C 1s and (F) O 1s.

In addition, N2 adsorption/desorption isotherms measurement and differential pore volume distributions are tested for both sample in Figure 3C. The specific surface area (SSA) of AAC dramatically increases from 21.32 to 660.05 m2 g−1. While concurrently the accumulative pore volume of AAC scales up to 0.29 cm3 g−1, overwhelming larger than the AC of merely 0.01253 cm3 g−1. Unlike AC showing negligible few pores, AAC contains ample micropores and mesopores concentrating at ~6 and 11 nm. The improved SSA and porosity, resulted from synergic effect between the chemical conversion of KOH and formation of H2O and CO2, offers not only extra active sites for Li+/Na+ ion but also shortened ion transportation path with reduced inner resistance (Yu et al., 2016). To understand the existing chemical states of all elements in AAC, XPS measurement is characterized. Judging from XPS curve as displayed in Figure 3D, the peaks at 285.08 and 534.08 eV are attributed to C and O elements, yielding an occupancy of 94.56 and rest 5.44%, respectively. Individually, as displayed in Figures 3E,F, four predominant peaks are specified as C-C (284.78 eV), C-OH (285.48 eV), C=O (286.88 eV), and C-OOH (290.08 eV) in deconvoluted C 1s spectrum. Besides, three apparent peaks in deconvoluted O 1s spectrum are assigned to the representative C=O (531.88 eV), C-O-H (533.28 eV), and C-O-C (534.98 eV) (Lin et al., 2016; Li et al., 2017a).

Figure 4 illustrates the electrochemical performances of AAC and AC as anode materials for LIBs. Electrochemical behaviors of AAC and AC during the initial three cycles are examined through cyclic voltammetry (CV) test at a scan rate of 0.1 mV s−1 within voltage range between 0.01 and 3.0 V, as presented at Figure 4A and Figure S3A, respectively. Like AC sample, three identifiable cathodic peaks for AAC, which center at range of 1.8–0.4 V in the first CV curves and vanish in the following cycles are ascribed to the formation of irreversible solid electrolyte interface (SEI) and unanticipated side reactions due to the consumption of oxygenic functional groups (Lv et al., 2015; Li et al., 2017b; Huang et al., 2018) Meanwhile, a reversible cathodic peak appears at potential window of 0.5–0.01 V, in concert with the anodic peak at 0.2 V, is originated from the reversible lithiation/delithiation process. It is observed that the 2nd and 3rd cycle of CV curves are scarcely discriminated, suggesting that SEI films formed at the first cycles securely protect the structural integrity of AAC an AC, and enduring the rapid ion insert/extraction. Correspondently, the first galvanostatic charge/discharge portrait at 0.1 A g−1 for AAC (Figure S2A) shows charge and discharge plateaus at around 1.6 and 1.2 V, respectively, which are generated from unavoidable SEI fabrication and electrolyte decomposition, leading to an irreversible capacity loss. Therefore, the charge and discharge capacities of AAC for the first cycle are 777.8 mAh g−1 and 1803.8 mAh g−1, with initial coulombic efficiency (ICE) about 43.12%, which is exceeding the AC anode.

Figure 4

Electrochemical performances of AAC and AC as anode materials for LIBs: (A) cyclic voltammetry (CV) curves at 0.1 mV s−1 of AAC; (B) charge/discharge curves of AAC at different current densities; (C) rate performances; (D,E) long-term cyclic performances at 0.5 A g−1 (D) and 2 A g−1 (E).

Evaluation of rate performance for AAC anode is conducted by gradient current densities, as depicted in Figures 4B,C. AAC anode achieves the charge capacities of 459.7, 335.4, 288.7, 251.7, and 213.5 mAh g−1 at rate of 0.1, 0.5, 1, 2, and 5 A g−1, respectively, whereas AC delivers an overall inferior performance with unsatisfactory capacity decline at high current density of 5 A g−1 (Figure S3B). When the current density decreases from 5 to 0.1 A g−1 by stages, the charge capacity of AAC are restored to nearly 100% compared with the pristine state, reflecting the outstanding rate tolerance and stable structural maintenance of AAC. In galvanostatic charge/discharge test (Figures 4D,E), the charge capacity of AAC attractively reaches 479.2 mAh g−1 after 150 cycles at a current density of 0.5 A g−1 with 96% retention. Furthermore, it can achieve reversible capacity of 691.7 mAh g−1 after 1200 cycles at a high rate of 2 A g−1 with an ascending tendency in the first few hundred cycles, which is stemmed from gradual electrochemical activation of AAC and is similar to most carbonaceous anode materials (Qie et al., 2012; Lv et al., 2015; Xiong et al., 2018; Yang et al., 2018; Zhang et al., 2018b). On the contrary, ordinary AC merely delivers 261.6 mAh g−1 at 0.5 A g−1 and 229.1 mAh g−1 at 2 A g−1 under the same cycling condition. Conspicuous improvements in electrochemical performance for AAC anode can be explained by synergistic effort of highly porous structure and expanded interlayer of d002 plane due to KOH activation. Specifically, the interface between electrode and electrolyte is modified by multiple pores for wide wettability, facilitating a much more efficient accesses for ions transportation through minimizing inherent impedance and shorten the diffusion paths. Besides, the larger interlayer distance favors the reversibility of ion intercalation and extraction. Furthermore, taking the advantages of this unique structure, additional active sites for Li-ions adsorption are provided, and integral structure to mitigate volume variation are guaranteed, making AAC into a favorable anode material for LIBs (Qin et al., 2014; Xiong et al., 2016, 2017).

For further evidences, CV curves acquired at different scan rate varied from 0.1 to 10.0 mV s−1 with similar feature are shown in Figure 5A. Capacitive behavior of AAC is estimated via the following equation: or: where i is the responsive current at various scan rates (v) at randomly chosen potential (V), and k (both k1 and k2) are constants. k1 can be attained by calculating the slope of the liner relationship between i(V)/ν1/2 and ν1/2 at every fixed potential, evaluating the proportion of charge storage from capacitance and calculating capacity contributions from pseudocapacitance and diffusion (Chen et al., 2017; Wang et al., 2018; Zhang et al., 2018a). As scan rate increases, capacitive distribution increases in the same time and achieves 42.7% at a scan rate of 0.2 mV s−1 (Figures 5B,C). When tested at a scan rate of 10 mV s−1, the capacitive distribution can reach as high as 84.6%. The Li+ ions are easy to make absorption or reactions in macro-micropores, resulting in the pseudocapacitive behavior at edges or on the surface of amorphous carbon. The higher ratio of pseudocapacitance allows charges to transfer in an effective way. With rational design and modification of AAC anode, excellent rate capability can be accomplished by reallocating the proportion of capacitive contribution to a more dominant occupancy. In order to confirm the smaller resistance of AAC, the electrochemical impedance spectroscopies (EIS) for AAC and AC anodes are carried out in Figure 5D. Obviously, the charge-transfer resistance (Rct) of AAC located in high frequency is much smaller than AC, reducing from 319.5 to 89.07 Ω, which results in the considerably improved electronic conductivity and reaction kinetics (Wu et al., 2016b, 2017). In accordance with capacity distribution, the effective rate response of AAC can be ascribed to the excellent reaction dynamics (Wang et al., 2018).

Figure 5

(A) CV curves measured between 0.01 and 3.0 V at different sweep rate from 0.1 to 10.0 mV s−1; (B) CV curve with a blue closed area representing pseudocapacitance at a sweep rate of 0.2 mV s−1; (C) bar chart of pseudocapacitive contribution at corresponding sweep rate; (D) electrochemical impedance spectra in form of Nyquist plots of AAC and AC before cycling with equivalent circuit diagram as the inset.

The increasing price for Li-containing cathodes incents researchers to seek for alternative energy storage system, among which SIB is suitable and promising. Therefore, with the instruction of the employment in LIBs, the AAC can be also applied for SIBs (Figure 6). Figure 6A and Figure S3C present the CV curves of AAC and AC performed at a sweep rate of 0.1 mV s−1 within 0.01-3.0 V, respectively. Analyzing the CV curve of AAC anode (Figure 6A), it can observe a broad hump with two protuberant peaks at 1.0 and 0.4 V, which are disappeared in next two cycles. These irreversible peaks should be assigned to the formation of protective SEI film and undesirable reactions of oxygen-containing functional groups (Zhu et al., 2018), and this analogous phenomenon is also demonstrated in CV profiles of AC anode (Figure S3D) While a redox peak at 0.03 and 0.06 V is observed, referring to Na+-ion insertion/extraction. The overlapped CV profiles after the first cycle indicate a fully covered SEI film on AAC, suggesting the highly reversible redox reaction and stable structure. As shown in Figure S2B, AAC anode can exhibit an initial charge capacity of 328.7 mAh g−1 with the coulombic efficiency about 70%, while it maintains almost 100% columbic efficiency in the succeeding charge/discharge cycles. It is worth noticing that the initial coulombic efficiency of AAC is surprisingly higher than AC (Figure S2B), which can be reasonably elucidated to the activation of AAC. Based on adsorption-insertion sodium storage mechanism of amorphous carbon (Hou et al., 2017a; Li et al., 2017c), increasing amount of active sites (pores and defects) can attract more Na+ ion and the enlarged interlayer spacing is found suitable for Na+ intercalation/extraction between interlayers.

Figure 6

Electrochemical measurements of AAC and AC as anode materials for SIBs: (A) cyclic voltammetry (CV) curves at 0.1 mV s−1 of AAC; (B) electrochemical impedance spectra in form of Nyquist plots of AAC and AC before cycling with detailed view in high frequency area as the inset; (C) rate performances of AAC and AC and (D) charge/discharge curves of AAC at different current densities; (E) long-term cyclic performances at 0.5 A g−1; (F) bar chart of pseudocapacitive contribution at various sweep rate; (G) long-term cyclic performances at 2 A g−1.

The comparison of AAC and AC anodes is measured sequentially at different current rates from 0.1 to 5 A g−1 (Figures 6C,D). The reversible capacities of AAC anode at every rate are much higher than AC anode, reaching to 324.8, 321.6, 312.1, 298.9, 282.3, and 272.4 mAh g−1 at 0.1, 0.2, 0.5, 1.0, 2.0, and 5.0 A g−1, respectively. Surprisingly, it can recover to 343.5 mAh g−1 when current density rolls back to 0.5 A g−1, indicating the superior rate performance. Besides, the long-term galvanostatic charge/discharge property of AAC and AC anodes are presented in Figures 6E,G. AAC anode obtains a charge capacity of 361.7 mAh g−1 at 0.5 A g−1 after 200 cycles and maintains 189.1 mAh g−1 at high rate of 2 A g−1 beyond 1000 cycles, surpassing the AC anode, which possesses awfully low charge capacity of only 113.1 mAh g−1 and 41.1 mAh g−1 at 0.5 A g−1 and 2 A g−1, respectively, with continuous degradation. The above enhanced cycling performance is originated from the three-dimension hierarchical porous structure equipped with well-fabricated SEI film as well as intrinsic shortened transport routes and facile intercalation for Na+ (Ou et al., 2017).

EIS measurements (Figure 6B) validate the excellent cyclic and rate capability of AAC by simulating the equivalent circuits. Figure 6B displays that AAC exhibits satisfactory low impedance Rct of about 26.2 Ω, which is much smaller than that of AC (160.3 Ω), proving the low charge transfer resistance and fluent reaction process in AAC. Furthermore, the calculation of pseudocapacitive contributions at various sweep rates are also carried out in Figure 6F and Figure S4. Capacitive behavior of AAC is accountable for 78.9% of total capacity at sweep rate of 0.2 mV s−1, signifying that the critical capacitive behavior induced by larger specific surface area ensures the outstanding cyclic and rate performance of AAC. Compared with other published carbon anodes (Tables S1 and S2), it confirms that AAC exhibits attracting electrochemical performance for both lithium and sodium storage.

Conclusion

In general, the activated amorphous carbon (AAC) with porous structure derived from renewable camellia pollen grains is synthesized by a facile approach of homogeneous KOH activation followed by carbonization. AAC exhibits the enhancement in specific surface area (660.04 m2 g−1), the enrichment in micro/mesopores (0.29 cm3 g−1 for accumulated pore volume) and the enlargement in interlayer distance (from 0.352 to 0.409 nm) as well. Furthermore, the modified microstructure can facilitate the infiltration for electrolyte into the interior of AAC electrode. On the other hand, it accesses the efficient transportation of Li+/Na+ ions by lowering the charge transfer resistance and shorten the diffusion paths. It is noted that the intercalation and adsorption of Li+/Na+ ion are coincidentally promoted, leading to the outstanding electrochemical performances. As a result, AAC anode surprisingly achieves a reversible charge capacity of 691.7 mAh g−1 after 1200 cycles at 2.0 A g−1 and 213.5 mAh g−1 at ultrahigh rate of 5.0 A g−1 for LIBs. Moreover, AAC delivers reversible capacities of 324.8, 321.6, 312.1, 298.9, 282.3, and 272.4 mAh g−1 at 0.1, 0.2, 0.5, 1.0, 2.0, and 5.0 A g−1 in gradient rate tests for SIBs. It has been discussed in plenty of articles that, concentration of activating reagent (KOH in this article) and activation duration are mutually responsible for the eventualities, which means that with rational adjustments can we delicately tailor the microstructure of pollen-grain-derived amorphous carbon. Alone with the particularity of pollen grains as a renewable and cost-effective biomass, the activated amorphous carbon with such intriguing structure is likely to become the future promising anode materials for lithium and sodium ion batteries.

Author contributions

KX and JX conducted the experiments, CY is the supervisor of this research work. YL and XO helped writing. KX, YL, and XO performed the characterization and data analysis. All authors involved the analysis of experimental data and manuscript preparation.

Conflict of interest statement

KX, YL, WS, and GZ were employed by company Electric Power Research Institute of Guangdong Power Grid Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.