Nanosponges by the oxo-Michael polyaddition of cyclodextrins as sorbents of water pollutants: the o-toluidine case

Hydrophilic cyclodextrin nanosponges were prepared by the oxo-Michael polyaddition in an aqueous solution at pH > 10 of α-, β-, and γ-cyclodextrin with 1,4-bisacryloylpiperazine or 2,2-bisacrylamidoacetic acid. These nanosponges and, for comparison purposes, their precursor cyclodextrins were tested as sorbents of o-toluidine, a carcinogenic wastewater contaminant, by monitoring the depletion of o-toluidine from a 10−4 M (10 ppm) aqueous solutions. To this aim, an innovative analytical procedure was used: The voltammetric peak currents of o-toluidine in linear sweep voltammetry experiments were registered using multi-walled carbon nanotubes-modified glassy carbon electrodes. The experimental sorption curves fitted a mono-exponential kinetic model, and the residual o-toluidine was 0.16 ppm, one order of magnitude lower than those of all other sorbents reported so far. The sorption capacities ranged from 88 to 199 µmol g−1 (10–21.3 mg g−1), equal to or higher than those of the parent cyclodextrins. All nanosponges were completely regenerated by extracting with methanol. After regeneration, the sorption capacity slightly improved, suggesting a rearrangement of the nanosponge network. Overall, it may be reasonably concluded that the cyclodextrin nanosponges reported in this paper warrant potential as o-toluidine exhaustive sorbents.


Introduction
o-Toluidine is an important organic intermediate used, inter alia, in the production of dyes, synthetic rubber, chemicals, and, as a curing agent, epoxy resins. In addition, it is a major component of tobacco smoke (OECD 2004) and, consequently, ubiquitous in the human environment. o-Toluidine manufacturing and processing have produced increasing amounts of toluidine-polluted wastewaters, eventually merging into water bodies. o-Toluidine is classified by the International Agency for Research on Cancer (IARC) as "carcinogenic to humans" (group 1) (Baan et al. 2008); therefore, its efficient removal from water is required. Different technologies are currently available to remove o-toluidine from wastewater, including photodegradation Pifferi et al. 2015;San and Çjnar 2016), electrochemical treatment (López-Grimau et al. 2013), oxidation (Anotai et al. 2011), microbial digestion (Liu et al. 2002), and adsorption. Among them, adsorption is reckoned as one of the most attractive and effective purification techniques in wastewater treatment. Various sorbents have been reported in the literature (Elgarahy et al. 2021;El-Sayed et al. 2019) and in particular for anilines, including o-toluidine, such as organoclay (Bishop et al. 2002;Zhang and Sparks 1993), montmorillonite (Essington 1994), silica (Voumard et al. 1995), rubber tire (Gupta et al. 2012), zeolites (Titus et al. 2002), activated carbon Responsible Editor: George Z. Kyzas in powder, granular, or fiber form (Duman and Ayranci 2005) and polymers (Jianguo et al. 2005).
Cyclodextrins (CDs) are cyclic oligosaccharides formed by 6 (α-CD), 7 (β-CD), or 8 (γ-CD) glucopyranose units. They are well known to possess a hydrophobic cavity forming inclusion complexes with organic molecules through host-guest interactions, with dimensional selectivity due to the different sizes of their inner cavity (Dodziuk 2006;Szejtli 1998). CD-containing polymers, either linear  or immobilized on solid supports (Allabashi et al. 2007;Fan et al. 2003;Faraji et al. 2011) or crosslinked CD resins García-Zubiri et al. 2007;Romo et al. 2008;Wilson et al. 2011;Yamasaki et al. 2006;Yu et al. 2003) were extensively studied as sorbents of organic pollutants with removal efficacy ranging within ample limits (Crini 2005, Yadav et al. 2022. Since the late nineties, hyper-crosslinked CDs have commonly been referred to as nanosponges (Ahmed et al. 2013;Li and Ma 1999). Nanosponges are supramolecular cage-like architectures in which the CD units are connected by many short moieties forming nanochannels (Bilensoy 2011;Krabicov et al. 2020;Mhlanga et al. 2007;Petitjean et al. 2021;Rizzi et al. 2021;Sherje et al. 2017;Taka et al. 2017;Trotta et al. 2012;Varana et al. 2020). Nanosponges are nanoporous by definition and show high inclusion constants with many organic pollutants, including aromatic and chlorinated compounds (Berto et al. 2007a and2007b;Gupta et al. 2012;Sun et al. 2010;Rizzi et al. 2021;Yadav et al. 2022). Most nanosponges are based on β-CD, the cheapest CD on the market, and only a few examples of α-and γ-CD nanosponges have been reported . Nanosponges were obtained by interconnecting the CD molecules with different crosslinking agents, such as N,N'-carbonyldiimidazole (Ansari et al. 2011), triphosgene (Trotta 2011), diphenyl carbonate (Swaminathan et al. 2013, or organic dianhydrides (Mognetti et al. 2021). For all of them, the coupling reactions demanded organic solvents. Subsequently, β-CD-based nanosponges were synthesized in water at pH > 12 by the Michael oxo-polyaddition of β-CD, or β-CD/2-methylpiperazine mixtures, with 2,2-bisacrylamidoacetic acid (Swaminathan et al. 2010). Under these conditions, β -CD acted as a multifunctional monomer and gave rise to a previously undescribed type of nanosponges. In the present work, this preparation method was extended to α-and γ-CD, and a small library of CD nanosponges was prepared by polyaddition of all three CD's with either 1,4-bisacryloylpiperazine or 2,2-bisacrylamidoacetic acid. These nanosponges were then tested as sorbents to remove o-toluidine from water and compared in this regard with their CD precursors. o-Toluidine sorption was monitored by using a previously optimized electroanalytical method based on linear sweep voltammetry (LSV) at multi-walled carbon nanotubes-modified glassy carbon electrode (Mardegan et al. 2014;Pifferi et al. 2014). The accuracy, sensitivity, and robustness of this optimized method were already demonstrated in previous papers on the degradation of o-toluidine Pifferi et al. 2015).

Synthesis of BAC-based CD nanosponges
In a typical procedure, a solution of α-CD (3.125 g, 2.9 mmol) and LiOH•H 2 O (0.357 g, 8.4 mmol) in water (1.9 mL) was mixed with an aqueous solution of BAC (2.135 g, 10.6 mmol) and LiOH•H 2 O (0.445 g, 10.6 mmol) in 2.2 mL water. The mixture was allowed to react for 24 h at 25 ℃. The final product appeared as a homogeneous, transparent, and soft nanosponge, which was first soaked for 2 h in 50 mL deionized water, allowed for swelling for a further 2 h, brought to pH 5 with 37% HCl, and soaked for 2 h in ethanol (50 mL). This extracting procedure was repeated three times. The extracted product was finally dried to constant weight. Yield: 5.26 g (69%).
BAC-β-CD and BAC-γ-CD were prepared similarly. The reagents used in the preparation of all BAC nanosponges and their amounts were as follows.

Synthesis of BP-based CD nanosponges
BP-α-CD and BP-β-CD were prepared following the same procedure described for BAC-α-CD by substituting BP for equimolar amounts of BAC, apart from the fact that BP was dissolved in plain water before mixing with the CD/LiOH solution. The same reaction performed with γ-CD failed to produce insoluble crosslinked nanosponges. In detail, the reagents and amounts used were as follows: α-CD (3.125 g, 2.9 mmol), LiOH•H 2 O (0.235 g, 5.5 mmol), H 2 O (1.7 mL), BP (2.072 g, 10.7 mmol), H 2 O (1.9 mL). Yield: 52% (5.20 g).

Elemental analysis
Elemental analyses were performed by the Analytical Laboratory of the Polytechnic University of Milano, Italy, using a Mettler Toledo Elemental Analyzer.

FTIR analysis
The FTIR spectra were collected on solid samples dispersed in KBr pellets using a Perkin Elmer 100 spectrometer, as the average of 16 individual scans at 2 cm −1 resolution in the 4000-600 cm −1 interval and with corrections for atmospheric water and carbon dioxide.

Swelling tests
Dry finely ground samples (30 mg) were suspended in ethanol (7 mL) and allowed to settle. The supernatant liquid was then carefully removed and the solid suspended three times in either water or 0.01 M PBS pH 7.4, gently shaken for 15 min, and then allowed to settle. Finally, the swelling degree (SD) was evaluated at equilibrium using Eq. (1): where V t is the volume of the swollen nanosponge at equilibrium and V t0 is the volume of the dry sample.

Thermal analyses
Differential scanning calorimetry (DSC) analyses were performed with a heat flux Mettler Toledo DSC823 calorimeter on 10 mg samples sealed in standard aluminum pans under a 20 mL min −1 nitrogen flow and with a 20 ℃ min −1 heating/ cooling rate. Empty pans were used as references. The analyses consisted of first heating run from 25 to 200 ℃, followed by a cooling run from 200 to 25 ℃ and second heating run from 25 to 300 ℃. Thermogravimetric (TG) analyses were performed with a Perkin Elmer TGA analyzer 4000 on 10 mg samples, in the 30-600 ℃ range, under a 50 mL min −1 nitrogen flow and with a 20 ℃ min −1 heating rate.

Assessment of the o-toluidine sorption ability of CD nanosponges o-Toluidine detection
All experiments were carried out using an Autolab PG-Stat 12 (Ecochemie, the Netherlands) potentiostat/galvanostat, according to a previously optimized method ). The voltammetric cell was conical (5 mL operating volume), including an Ag|AgCl|saturated KCl electrode as a reference, a Pt wire as a counter, and a multi-walled carbon nanotubes (MWCNT)-modified glassy carbon (GC) electrode as working electrodes, respectively. The GC electrode surface was initially cleaned with synthetic diamond powder (Aldrich, diameter 1 µm) on a Struers DP Nap wet cloth. Subsequently, 20 μL of a suspension of MWCNT in dimethylformamide (0.5 mg mL −1 ) was deposited on a glassy carbon electrode and dried at 25 ℃ until complete evaporation of the casting solvent. The working cell was thoroughly cleaned with nitric acid and ultrapure water (MilliQ ® Millipore) before each experiment. The sorption kinetics was assessed by monitoring the corresponding voltammetric peak currents during linear sweep voltammetry (LSV) experiments. To each sample of sorption solution (0.6 mL), 0.6 mL of 0.2 M HCl (ultrapure Fluka) was added as a supporting electrolyte before measurement. Voltammograms were recorded from the lowest concentration to the highest one. Peak heights were obtained after background subtraction.

Assessment of o-toluidine sorption ability
The ability to sorb o-toluidine from the CD nanosponges and the parent soluble α-, β-, and γ-CDs was determined by monitoring o-toluidine depletion from an aqueous solution of 10 −4 M (10 ppm) in batch reactors under slight magnetic stirring. The selected starting concentration of o-toluidine was one order of magnitude higher than the UK maximum exposure limit (Gregg et al. 1998).
As regards α-, β-, and γ-CDs, the sorbent (30 mg) was dissolved in 100 mL 10 −4 M o-toluidine solution. As regards nanosponges, the sorbent (30 mg) was put in a small punctured paper bag (tea bag) and subsequently inserted in 100 mL 10 −4 M o-toluidine solution. This procedure allowed a facile sorbent separation from the decontaminated water at the end of the experiment.

Nanosponge regeneration after o-toluidine sorption
The o-toluidine-loaded samples (30 mg) were soaked in 70:30 methanol/water (20 mL) and shaken for 20 min. The sorbent was then recovered, thoroughly rinsed with water to completely remove methanol, and finally re-used for o-toluidine sorption. The sorption-desorption cycles were repeated twice.

Synthesis of CD nanosponges
The synthetic procedure adopted for preparing CD nanosponges involved the activation of the CD hydroxyl groups by deprotonation in water at pH ≥12 (Scheme 1). Under these conditions, α-and β-CD's behaved as multifunctional monomers giving crosslinked resins (i.e., nanosponges) by one-pot polyaddition with bisacrylamides. No additional solvents or catalysts were added. A neutral bisacrylamide (1,4-bisacryloylpiperazine, BP) and a strongly acidic carboxylated bisacrylamide (2,2-bisacrylamidoacetic acid, BAC) were employed. The latter was chosen to ascertain whether its acidic nature could favor the sorption of o-toluidine by establishing ionic interaction. With γ-CD, the same reaction succeeded only in the case of BAC, while BP failed to produce insoluble resins. Therefore, BAC-γ-CD nanosponge was the only γ-CD-based one prepared.
For the sake of simplicity, Scheme 1 assumes that the polyaddition involved only the primary CD hydroxyl groups, neglecting the secondary ones. This, however, was not proven, and no reasonable predictions in this regard can be made. In fact, whereas the primary hydroxyl groups are less sterically hindered, the secondary ones are stronger acids, hence preferentially ionized (Loftsson and Brewster 1996). The bisacrylamide/CD molar ratio in the feed (Table 1) was adjusted to prepare nanosponges with maximum CD content, hence, presumably, maximum sorption capacities. As stated above, only BAC gave nanosponges with all three CDs, while BP failed in giving nanosponges with γ-CD. All nanosponges were purified by repeatedly soaking with water, then with ethanol, and finally dried to constant weight at 50 ℃ and 0.2 tor. Their Scheme 1 Synthesis of BAC-CD and BP-CD nanosponges. As an example, the structure of the reaction product of β-CD with BAC is shown in detail compositions, as determined by elemental analysis, are reported in Table 1. On the whole, the CD content of the resultant resins ranged between 0.60 mmol CD g −1 and 0.69 mmol CD g −1 nanosponge (67.0 wt% and 78.0 wt%) and, for all nanosponges, was higher than in the reactant mixtures, probably due to the loss of residual unreacted bisacrylamides or soluble CD/acrylamide oligomers. This loss was apparently slightly more pronounced for α-CDbased nanosponges. The bisacrylamide moiety/CD molar ratio, calculated from elemental analyses data, ranged between 1.80 and 2.54. The lowest ratio was observed for BAC-γ-CD, the only γ-CD resins obtained, while in all other cases, it was greater than 2.21. This result is in line with the higher water sorption of BAC-γ-CD (see below), due to a lower crosslink extent of this nanosponge.
All nanosponges were characterized by FTIR spectroscopy. The spectra (Fig. 1) showed intense diagnostic bands typical of parent CDs and their addition products with either BAC or BP.
In particular, the characteristic -OH, -CH 2 , and C-O stretching bands of the CD patterns were observed at 3283 cm −1 , 2924 cm −1 , and 1024 cm −1 , respectively. Meanwhile, in the BAC-based nanosponges, the diagnostic amide C=O stretching and NH bending bands of BAC addition residues were observed at 1616 cm −1 and 1515 cm −1 , respectively. The same bands for BP addition residues, with the exception of the NH bending bands of BAC, were observed at 1616 cm −1 and 1436 cm −1 . The -CH 2 -groups, present in all amide-deriving structural units, gave rise to the C-H stretching and bending bands observed at 2924 cm −1 and 1380 cm −1 , respectively. The broad band typical of the N-H amide stretching, expected at 3283 cm −1 for BAC-deriving nanosponges, overlapped with the OH stretching band of CD.
All nanosponges were characterized by differential scanning calorimetry (DSC) and thermogravimetric analysis (TGA). The DSC traces of the second heating cycle following a cooling cycle carried out with the same rate (20 ℃ min −1 ) (Fig. 2) showed substantially flat curves in the 20-150 ℃ range, with no evidence of transition temperatures. The thermal decomposition onset was observed at temperatures higher than 200 ℃, in line with previously  The thermal stability of BAC-CD and BP-CD nanosponges was compared with that of parent CDs, as shown by the TG traces displayed in Fig. 3 and by the collected data reported in Table 2. In line with previous observations on the thermal behavior of CD nanosponges obtained using different crosslinkers (Kumar et al. 2018;Salgin et al. 2017), the decomposition patterns of BAC-CD and BP-CD appeared to be far more complex than those of plain cyclodextrins.
In agreement with literature reports (Trotta et al. 2000), the TG curves of cyclodextrins had a similar three-stage profile, including (i) a first stage of dehydration, up to 120 ℃, in which sorbed and crystallization water was removed; (ii) a second stage, from 250 to 350 ℃, associated with a T onset10% around 300 ℃, a maximum weight loss of 75-80% at T max1 332-333 ℃, and the formation of a thermally stable residue, called char; and (iii) a third stage at temperatures above 400 ℃ in which char formed in the previous stage progressively degraded down to a residual mass fraction at 600 ℃, RMF 600 , of 11.0-12.5%. The main decomposition phase centered at T max1 was attributed to the opening of the cyclodextrin ring, followed by structural transformations similar to those of cellulose, which include the loss of glycosidic and hydroxyl groups and the formation of double bonds, carbonyl groups, and aromatic structures.
The TG patterns of BAC-CD and BP-CD nanosponges exhibited T onset10% values invariably lower ( Table 2) than those of cyclodextrins. Moreover, the weight loss patterns proceeded through four decomposition steps (Fig. 3) down to a significantly higher RMF 600 of 23.0-33.0%. In detail, apart from the dehydration stage, up to 120 ℃, a second weight loss of approximately 5-15% occurred in the 150-300 ℃ range, ascribed to the cleavage of the C-C bond between CD and the linker. The third stage of degradation, ascribed to CD decomposition, was found at slightly different temperatures depending on the nature of the linker. All BAC-CD nanosponges had similar T max2 , around 345 ℃, with a weight loss of around 45% in the range of 300-400 ℃. Regarding BP-CD nanosponges, BP-α-CD exhibited a T max2 value of 445 ℃ and BP-β-CD of 342 ℃. The fourth step (400-600 ℃ and 500-600 ℃ for BAC-and BP-based cyclodextrins, respectively) was ascribed to the relatively slow thermal degradation of the char formed in the previous step.
The swelling ability of nanosponges was determined both in water and 0.01 M PBS pH 7.4. The results are shown in Fig. 4. The swelling in water was higher than in buffer and, for nanosponges with the same crosslinking arm, increased in the order α-CD < β-CD < γ-CD. The swelling degree is also influenced by the structure of the crosslinking arm. Under the same conditions, BAC nanosponges constantly swelled more than those based on the same CD, but deriving from BP, possibly due to the ionic nature of BAC.

o-Toluidine sorption
The o-toluidine sorption ability of CD nanosponges was assessed by monitoring the o-toluidine depletion from 10 −4 M (about 10 ppm) aqueous solutions in linear sweep voltammetry experiments, making use of MWCNT-modified glassy carbon (GC) electrodes, a robust and sensitive analytical method that has proven to detect o-toluidine down to 0.16 ppm with an apparent recovery factor of about 100% ). The sorption study was made at the natural pH in which o-toluidine is present in water-polluted samples. Moreover, the range of pollutant concentration investigated in this work was really low, thus not influencing the solution pH during the sorption.
The observed sorption capacities were compared with those of native soluble α-, β-, and γ-CDs. In all cases, the sorption kinetics, expressed both referring to the sorbent mass unit (Fig. 5a) and the molar CD content of the sorbent (Fig. 5b), were found to fit a mono-exponential kinetic model derived from a more conventional doubleexponential model in which the second term prevails over the second one, indicating a good homogeneity of the sorbent materials (Ferruti et al. 2006), and the absence of adsorptive phenomena in competition with sorption. The mono-exponential model implies the following equations (Eqs. (2) and (3)): where m sorb is the weight of the sample of the sorbing material, μmol CD is the quantity of CD present in the sorbing material, k D is the relevant sorption rate constant, q e and q e ' are the equilibrium o-toluidine sorption capacities, and D and D' are further parameters modulating the sorption rate. The relevant kinetic parameters obtained for all tested samples are reported in Table 3.   The sorption efficiency (see q e and D values) of β-CD nanosponges is higher than that of α-CD-and γ-CD-nanosponges in line with the higher efficiency of β-CD compared with that of both α-CD and γ-CD, due to the specific capacity of the internal cavity size of β-CD to create stable inclusion complexes with monocyclic aromatic compounds (Szejtli 1998).
The equilibrium sorption capacities of both α-CD and β-CD nanosponges, expressed as mol o-toluidine mol −1 CD (q e ' factors, accounting for the number of sorption sites), were remarkably higher than those of parent cyclodextrins, particularly for α-CD nanosponges. This was ascribed to two concurrent factors: the cooperation of correctly oriented cyclodextrin groups in forming inclusion complexes with o-toluidine and the additional uptake of o-toluidine

Sample
Sorption performance per g of sorbent Sorption performance per µmol of CD in the sorbent q e (μmol g −1 ) D (μmol L −1 ) q e ' (μmol μmol −1 ) D' × 10 −3 (μmol L −1 ) k D × 10 −3 (min −1 ) α-CD 47 ± 6 1.6 ± 0.2 0.05 ± 0.01 0.0031 ± 0.0004 3.6 ± 0.9 BAC-α-CD 104 ± 7 3.5 ± 0. by the interstitial channels created by the interconnected arms. On the opposite, the sorption capacity of BAC-γ-CD was approximately the same as that of γ-CD. This could be attributed to the fact that due to the larger size of the internal cavity size of γ-CD and to the lower crosslinking degree that probably leads to wider interstitial channels in γ-CD nanosponges, o-toluidine complexation does not benefit from the establishment of cooperative effects between adjacent CD macrocycles. It may also be observed that the BP-based nanosponges performed better than their BAC counterparts, more significantly in the case of BP-β-CD. This behavior was ascribed to the higher affinity of the BP-deriving interconnecting arms for o-toluidine. The above-reported assumption that owing to its acidic nature, the BAC interconnecting arms could favor the sorption of o-toluidine by establishing ionic interactions did not materialize.
As regards o-toluidine removal from water, it is worth mentioning that the best results obtained in this work are significantly superior to those of other o-toluidine sorbents. In particular, the o-toluidine abatement dropped to 0.16 ppm, which represents the detection limit of the adopted analytical technique, lower by one order of magnitude than the figures reported so far for all other sorbents (Sun et al. 2010). In addition, all nanosponges, irrespective of the size of the CD inner cavity, proved capable of removing o-toluidine from dilute aqueous solutions under conditions similar to those usually found on the field in most o-toluidine-polluted aquifers, abating it down to at least 0.16 ppm. Remarkably, the fact that all three CDs lead to effective nanosponges allows us to conclude that high-performance o-toluidine sorbents can be prepared from commercially available crude CD blends, which are very cheap compared to pure samples of α, β, and γ-CD used in this pioneering work. This is of paramount importance in view of the mass production of CDbisacrylamide nanosponges for water purification.

Nanosponge regeneration
Regeneration studies were performed for BP-β-CD, the sorbent with the highest o-toluidine equilibrium sorption. The nanosponge was simply regenerated by methanol extraction of the sorbed material. Three o-toluidine sorption-desorption cycles were performed. The sorption capacity after regeneration was monitored each time and compared with that of the pristine nanosponge. Interestingly, the sorption capacity ( Fig. 6) slightly improved after the first regeneration treatment, possibly due to a rearrangement of the nanosponges network that facilitated the o-toluidine inclusion in a sort of imprinting of the polymer. More investigations should be done to prove this preliminary idea.

Conclusions
Tailor-made α-, β-, and γ-CD-based nanosponges characterized by high CD content, in which the CD moieties were covalently connected by either ionizable or neutral bisacrylamide-deriving arms, were synthesized from cyclodextrins and bisacrylamides by base-catalyzed stepwise oxo-Michael polyaddition in aqueous solution. These nanosponges exhibited high water sorption ability and were thermally stable at least up to 250 ℃. Their TG patterns significantly differed from those of the parent CD and were in line with those of previously reported CD nanosponges. These CD nanosponges removed o-toluidine from 10 −4 M (~10 ppm) aqueous solutions (the reference UK-approved maximum exposure limit) down to the detection limit of the adopted analytical technique (0.16 ppm).
The nanosponges' sorption capacities were in the range of 88-199 µmol g −1 (10-21.3 mg g −1 ), one order of magnitude higher with respect to other CD-containing sorbents reported in the literature (Sun et al. 2010). Interestingly, the sorption capacity of nanosponges toward o-toluidine was even higher than that of the corresponding cyclodextrins. This can be attributed to the participation in the sorption process of the interstitial channels between the CD units.
The sorbents were simply regenerated by extracting with methanol and recycled with no loss of performance. On the opposite, their sorption capacity slightly increased, possibly due to some molecular imprinting during the first sorption/ desorption cycle. Also considering the simplicity, potential scalability (the preparation process can be extended to any compound containing two activated acrylamide double bonds (Ferruti et al. 2019)), environmental sustainability of their preparation, and no phyto-toxicity (Alongi et al. 2022), these cyclodextrin-based materials warrant potential as o-toluidine sorbents competitive with most materials presented so far to the same purpose and deserve further investigations aimed at widening their scope as sorbents of organic water pollutants other than o-toluidine.
Funding Open access funding provided by Università degli Studi di Milano within the CRUI-CARE Agreement.

Data availability
The main data generated within this research are included in the paper. Some further data that support the findings are available from the corresponding author upon reasonable request.

Declarations
Ethics approval This work does not contain any studies with human participants or animals.
Consent to participate All authors provided informed consent to participate in this research study.

Consent for publication
The authors transfer to the journal the publication rights and guarantee that this contribution is original. All authors mutually agree on its submission.

Competing interests The authors declare no competing interests.
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 original 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:// creat iveco mmons. org/ licen ses/ by/4. 0/.