Extraction of Cellulose Nano-Whiskers Using Ionic Liquid-Assisted Ultra-Sonication: Optimization and Mathematical Modelling Using Box–Behnken Design

This study focuses on the extraction of cellulose nano-whiskers (CNWs) from the leaves of Adansonia kilima (AK), usually known as African baobab, using a combination of a microwave-assisted alkali (KOH) pre-treatment with subsequent bleaching process prior to ultra-sonication. Ultra-sonication was carried out using the ionic liquid (IL) 1-butyl-3-methylimidazolium hydrogen sulfate (Bmim-HSO4). Process parameters for ultra-sonication were optimized using a two-level factorial Box–Behnken design (BBD). Process variables such as ultra-sonication power (x1), hydrolysing time (x2) and temperature (x3) were varied. Responses selected were percentage crystallinity index, CrI% (y1) and yield% (y1) for the finally procured CNWs sample. Regression analysis was carried out to develop quadratic model to analyze the effect of process variables on IL-assisted ultra-sonication process. Analysis of variance (ANOVA) showed that ultra-sonication power was the most influential aspect for hydrolyzing the amorphous segments of crude cellulose extracted from baobab leaves. A relative study of the physio-chemical properties of the starting lignocellulosic substrate (AK), KOH pre-treated, bleached and IL-assisted ultra-sonicated CNWs was conducted. The synthesized samples were characterized using Fourier transform infrared spectroscopy, Scanning electron microscopy, atomic force microscopy, high resolution transmission electron microscopy, X-ray diffraction and thermo-gravimetric and zeta potential analysis. Under optimum condition, the extracted CNWs showed an average width of 15–20 nm; with high crystallinity index of 86.46%. This research provides an insight about the delignification of Adansonia kilima (AK) leaves and its effective conversion to CNWs having high crystallinity.


Introduction
Application of renewable lignocellulosic feedstock to produce nanocellulose has gained extensive interest over the last few decades.Cellulose nano-whiskers (CNWs) can exhibit some unique features of having low density, high aspect ratio with extended surface area, tunable surface chemistry containing negative hydroxyl groups (-OH) with biodegradability and superior thermal and mechanical properties [1][2][3][4].They are considered as a promising candidate for versatile application in opto-electronic devices, smart composite materials for environmental as well as biomedical application [5][6][7].Several preparation techniques are available to obtain CNWs from ligno-cellulosic residues.Various mechanical or chemical approaches can be undertaken to obtain the nanocellulose having various morphological and chemical properties [8][9][10][11][12][13].Usually preliminary treatment of raw biomass substrates is carried out using alkali treatment and bleaching first to separate the lignin and hemicellulose fraction to obtain crude cellulose.After that, subsequent traditional acid hydrolysis with concentrated Lewis acid (H 2 SO 4 , HNO 3 , HCl) is carried out for selective hydrolysis of amorphous region of cellulose to yield nanocellulose.The cellulosic matrix in ligno-cellulosic biomass contain lignin polymer having complex structure in the primary cell wall which causes recalcitrance.This results in difficulties to yield nanocellulose [14,15].Eliminating lignin using energy efficient, green technology to intensify the degree of nano-fibrillation is often critical.Commonly used chemical reagents, especially corrosive acids, sulphite treatment and TEMPO mediated oxidation have been widely used to delignify the lignocellulosic substrates.Although successful extraction has been carried out by following those methods, some practical obstacles remain for their large-scale production, such as prolonged reaction time with lower yield percentages, use of concentrated corrosive mineral acids and environmental pollution.In this regard, some novel protocols consisting of microwave-assisted heating in presence of alkali (NaOH, KOH, Ca(OH) 2 ) and peroxide bleaching process with subsequent hydrolysis using ionic liquids (ILs)-as green solvent during the ultrasonication have been recognized as eco-friendly approaches [16][17][18].Application of ionic liquids (ILs) can prevent the use of corrosive mineral acids and reduce environmental pollution [8,19,20].
Ultrasonication can be carried out in different mediums to disperse the nanocellulose crystals in homogeneous suspension with a desired aspect ratio.The ultrasonication process also improves the mass transfer efficiency of nanocellulose extraction process in the presence of suitable solvent under optimized condition.Ultrasonic waves can cause mechanical oscillation inside the reaction medium to generate a cavitation effect.This will initiate formation, growth and collapse of microscopic gas bubbles.It will disintegrate the intra-and intermolecular hydrogen bonding in cellulosic fiber.At the secondary stage, the cellulosic microfiber will further disintegrate in the presence of suitable medium under optimum condition to yield cellulose nano-whiskers (CNWs).Furthermore, it does not alter the fibers' properties significantly [21][22][23][24].Application of acid hydrolysis reduces the yield of CNWs to a greater extent making its large-scale commercial production difficult.Thus, overall the process becomes time consuming which is not economically feasible [14,15].In this context, production efficiency of CNWs should be upgraded by combination of physiochemical treatment as well as avoiding the use of hypochlorite bleaching and application of harsh, corrosive Lewis acids like H 2 SO 4 , H 3 PO 4 , HCl, etc.In this research, the acidic IL of [Bmim]HSO 4 was chosen to extract CNWs during ultrasonication process.[Bmim]HSO 4 is classified as Bronsted acidic ionic liquid having pKa values of 2 in presence of water [16][17][18].It was expected that [Bmim]HSO 4 will play dual role for swelling and selective hydrolysis of crude cellulose obtained from pretreated and bleached AK leaves.The presence of shorter alkyl chain (Butyl) with alkylimidazolium cations of [Bmim]HSO 4 will have improved solvation capability with higher anion concentration.The presence of negative anions of [HSO 4 ]-will make it acidic and enhance the interlayer distance between the cellulosic chain by swelling and will cause destruction of strong hydrogen bond [16][17][18].
Compared to conventional heating, microwave-assisted heating process can be much more efficient based on starting lignocellulosic residues.Microwave-assisted pretreatment will initiate the breakdown of the cell walls of lignocellulosic residues which involves some biochemical conversion process.This will lead to disruption as well as disintegration of the recalcitrance structures and open some channels for subsequent hydrolysis processes.The cavitation effect can significantly improve the diffusion of heat through convection process to intensify the delignification process.The electromagnetic energy of microwave can transfer the heat inside the interior region of lignocellulosic substrate to separate the lignin from the cellulosic matrix within a shorter reaction time.Thus, the process also becomes energy efficient.It will liberate the cellulose fiber from lignin matrix, enhance the porosity of cellulose and increase the crystallinity index.Until recently, response surface methodology (RSM) has been widely used by the researchers to optimize the complex process of synthesis, modeling, and simulation with statistical analysis.Relatively lesser number of experiments are needed to determine the optimum condition for the reaction, thereby reducing the cost and time of the synthesis process effectively.It also can provide the insight about the relationship between the reaction variables (inputs) with responses (desired output) [25,26].
To the best of our knowledge, until now no research has been conducted to extract CNWs with high crystallinity index from Adansonia kilima (AK) leaves.In this research, dried leaves of Adansonia kilima (AK) were solvent-extracted to reduce pectin, wax, pigments, etc. Subsequently after that, microwave-assisted KOH pre-treatment and peroxide bleaching were carried out to extract CNWs in presence of Bmim-HSO 4 IL using ultrasonication.The ultrasonication process was optimized and relevant mathematical polynomial models were developed for crystallinity index, CrI% (y 1 ) and yield% (y 2 ) of CNWs.Statistical analysis with model validation were carried out using RSM technique based on Box-Behnken (BBD) experimental design.Finally extracted CNWs samples were characterized in terms of FESEM, HRTEM, AFM, XRD, TGA, FTIR and Zeta potential analysis.The output of this research clearly shows the role of Bmim-HSO 4 IL in extracting highly purified crystals of CNWs having high crystallinity index (CrI%).

Method
The dried leaves (AKL-1) were crushed and sieved to 200 µm.They were extracted with a 1:2 solution of ethyl alcohol and toluene for 8 h to remove pectin, wax, pigments, etc. [27,28].The sample thus obtained was washed with hot deionized water (DI) for several times and dried using a vacuum oven overnight at 60 • C. The sample was labeled as AKL-2 and stored in an airtight container to prevent fungal development.It was sent for subsequent characterizations.

Microwave Pre-Treatment Using KOH
Microwave heating of AKL-2 sample was carried out at constant temperature of 90 • C and the power was fixed to 300 watts to avoid char formation.The required amount (gm) of AKL-2 was added with KOH pellets and the ratio between dewaxed residues (AKL-2) with KOH was maintained at 1:0.25.The mixture was transferred to microwave glass tube with addition of 100 mL deionized water (DI).Microwave heating was carried out for 30 min under nitrogen flow.After 30 min, the sample was cooled to room temperature and filtered to separate the liquid fraction from the solid one.The solid fraction was the partially delignified sample (AKL-3) which was washed repeatedly with DI water and ethanol.The sample thus obtained was dried overnight in vacuum oven at 55 • C and labeled as AKL-3 for further characterizations.

Peroxide Bleaching
AKL-3 was refluxed with 100 mL of 30% H 2 O 2 for 6 h at 60 • C to bleach the sample completely.The sample was filtered and washed with DI water several times until the pH became 6.5-7.5.It was dried overnight in in vacuum oven at 55 • C and labeled as AKL-4 for subsequent characterizations.

Extraction of CNWs Using Ultrasonication
After alkaline peroxide treatment, the sample AKL-4 was ultra-sonicated in presence of ionic liquid, maintaining the reaction condition suggested by the design of experiment based on BBD design.The reaction conditions along with the design matrix obtained for ultrasonication are summarized in Tables 1 and 2. A total of 10 g of crude cellulosic sample (AKL-4) was added with 100 g of ILs solvent.Ultra-sonication power (x 1 ), time (x 2 ) and temperature (x 3 ) were varied accordingly (Table 2).After the predetermined period (hydrolysis time) of ultrasonication at different conditions (Table 2), 30 mL cold deionized waster (DI) was added with the mixture and was stirred until the temperature reached room temperature.The milky white, colloidal suspension of cellulose thus obtained was centrifuged at 6500 rpm for 30 min.The white cellulosic precipitate obtained was separated from the mixture and washed with hot DI water several times until the pH became 5.5-6.The sample thus obtained was freeze dried at −4 • C for two days and labelled as AKL-5 which was the finally extracted CNWs sample.In this research, ultrasonication experiments were conducted using response surface methodology using BBD design where the reaction variables (Power = x 1 , time = x 2 and temperature = x 3 ) were varied to observe its impact on preselected responses of crystallinity index, CrI% (y 1 ) and yield% (y 2 ) of finally extracted cellulose nano-whiskers (CNWs).The experimental data were analyzed using Design the Experiment (Design-Expert V9.0).A relevant mathematical model was developed and the accuracy of the model was validated by analysis of variance (ANOVA) test.The process optimization was carried out to get maximum crystallinity index, CrI% and yield% of CNWs (AKL-5).

Characterizations
The surface morphological features of all the samples synthesized here (AKL-2, AKL-3, AKL-4 and AKL-5) including the starting lignocellulosic residues of dried leaves (AKL-1) were observed using a field emission scanning electron microscope (FESEM) (Leo-Supra 55, Carl Zeiss, London, UK).The sample was placed over the conductive carbon tape and it was coated with gold having a thickness of 5 nm before FESEM analysis.From FESEM images, a software (Olympus, Tokyo, Japan) was used to measure the length of the extracted sample (AKL-5).After ultrasonication, HRTEM analysis was carried out for optimum sample (AKL-5) (JEOL-2100 F, Tokyo, Japan).The samples were mixed with C 2 H 5 OH and sonicated for 10 min and then a drop of sample was placed over the copper grid and allowed to dry before taking the HRTEM images.Atomic force microscopic (AFM) analysis was carried out for finally extracted CNWs sample (AKL-5) to observe the shape and surface topography (Multimode Nano-Scope-IIIA, Bruker, Billerica, USA).The samples were sonicated for 15 min to prevent aggregation and a drop was placed over the glass slides.After it was dried at ambient temperature, AFM analysis was carried out.Thermogravimetric analysis with DTG max (TGA-Q0500, Shimadzu, New York, USA) was carried out to observe the thermal degradation profile of all the extracted samples where 5 mg of all the sample was heated from 30 • C to 1000 • C under the nitrogen flow (150 mL/min) at a heating rate of 5 • C per minute.The crystalline phase, crystallinity index (CrI%) including the d-spacing of the samples, was analyzed by X-ray diffraction technique using CuKα radiation at a voltage around 2.7 Kw (Shimadzu-XRD-6000, Kyoto, Japan).The crystallinity index was calculated using Segal's method using Equation (1): CrI% = (I 002 −I am )/I 002 × 100%. (1) Here, I 002 represents the intensity of both crystalline and amorphous regions of the sample at 2θ = 22-24 • and I am represents the intensity of amorphous region of the sample [29].The samples were ground to powder and mixed with dried KBr and pressed to ultrathin films (approximately 4 mm) and sent for FTIR analysis (Perkin Elmer, Tokyo, Japan).

Mathematical Modeling and Statistical Analysis
For crystallinity index, CrI% (y 1 ) and yield% (y 2 ) of CNWs, two polynomial equations were suggested after regression analysis based on the data obtained from the basic design matrix (Table 2).The design matrix provided the reactions conditions for each experimental run with the output responses of crystallinity index, CrI% (y 1 ) and yield% (y 2 ).The models were chosen based on their polynomial order which should be highest here [28].Preselected input variables such as power (x 1 ), time (x 2 ) and temperature, (x 3 ) were significant for developing the models.All these input variables had positive impact on crystallinity index, CrI% (y 1 ) and yield% (y 2 ) up to certain extent.After that range, enhancing power (x 1 ), time (x 2 ) and temperature (x 3 ) would reduce the crystallinity index, CrI% (y 1 ) and yield% (y 2 ) both due to extensive hydrolysis of the amorphous and crystalline domain of cellulose to form other organic liquid fractions.Based on the magnitude of the sequential model sum of the squares, the models were not aliased [1,29].
For crystallinity index, CrI% (y 1 ) and yield% (y 2 ), following models were proposed, and the empirical equations obtained for that is given bellow: The numerical coefficients before the linear terms depicting power (x 1 ), time (x 2 ) and temperature (x 3 ) illustrated the impact of these specific parameters for ultra-sonication of crude cellulose (AKL-4).However, the terms obtained after multiplication of two different variables such as x 1 x 2 , x 2 x 3 and x 3 x 1 represent the interaction effects on the desired output/responses.The terms x 2 1 , x 2 2 and x 2 3 denote the quadratic effects.A positive sign in front of each term specifies the synergistic effect while the negative symbol shows the antagonistic effect [1,29].
Figure 1a,b displays the linear regression plots for predicted versus actual/experimentally analyzed results for crystallinity index, CrI% (y 1 ) and CNWs yield% (y 2 ), respectively.The plots revealed that the experimental data obtained here for both the responses were closer to the predicted values exhibiting R 2 values for the Equations ( 2) and (3) around 0.987 and 0.986 for crystallinity index, CrI% (y 1 ) (Figure 1a) and CNWs yield%, (y 2 ) (Figure 1b); correspondingly.This evidences the outstanding adjustment of the proposed models with the experimental observation.
The numerical coefficients before the linear terms depicting power (  ), time (  ) and temperature ( ) illustrated the impact of these specific parameters for ultra-sonication of crude cellulose (AKL-4).However, the terms obtained after multiplication of two different variables such as   ,   and   represent the interaction effects on the desired output/responses.The terms  ,  and  denote the quadratic effects.A positive sign in front of each term specifies the synergistic effect while the negative symbol shows the antagonistic effect [1,29].
Figure 1a,b displays the linear regression plots for predicted versus actual/experimentally analyzed results for crystallinity index, CrI% (y1) and CNWs yield% (y2), respectively.The plots revealed that the experimental data obtained here for both the responses were closer to the predicted values exhibiting R 2 values for the Equations ( 3) and ( 4) around 0.987 and 0.986 for crystallinity index, CrI% (y1) (Figure 1a) and CNWs yield%, (y2) (Figure 1b); correspondingly.This evidences the outstanding adjustment of the proposed models with the experimental observation.The studentized residuals obtained here for CrI% (y1) and yield% (y2) of CNWs (AKL-5) throughout the 17 experimental runs (Table 2) are shown by Figure 2a,b.Residual error is the difference between experimental data and predicted data for modelling.Studentized residual is obtained by dividing the residual error with standard deviation.Figure 2a,b showed that the points are arbitrarily scattered but were within the range of ±3.00.Thus, transformation of responses is not required, and the models used here are suitable to analyze the overall extraction process using IL-assisted ultra-sonication [1,29].The studentized residuals obtained here for CrI% (y 1 ) and yield% (y 2 ) of CNWs (AKL-5) throughout the 17 experimental runs (Table 2) are shown by Figure 2a,b.Residual error is the difference between experimental data and predicted data for modelling.Studentized residual is obtained by dividing the residual error with standard deviation.Figure 2a,b showed that the points are arbitrarily scattered but were within the range of ±3.00.Thus, transformation of responses is not required, and the models used here are suitable to analyze the overall extraction process using IL-assisted ultra-sonication [1,29].

Analysis of Variance (ANOVA) Test and Statistical Analysis
An analysis of variance (ANOVA) test was done to observe the competency of the developed model Equations ( 3) and ( 4).The statistical parameters with relevant correlation coefficient (actual and adjusted) after the linear regression analysis are listed in Table 3.The appropriateness of the developed models can be additionally confirmed by observing the values of standard deviation, coefficient of variation (CV) and adequate precision.The experimental R 2 values were in close agreement with the adjusted R 2 values.The magnitudes of standard deviation observed here for both the models were relatively smaller.Furthermore, the values for coefficient of variation (CV) were also only 1.22 and 1.06 for the empirical models Equations ( 3) and (4).This represents the reproducibility of the experimental data with the developed model.Adequate precision is defined as the ratio between signal to noise and for effective development of the polynomial model, the magnitude must be greater than 4.Here for both the responses, it was around 29.82 and 29.20, respectively, reflecting suitable navigation of the design [29].
Experimental data obtained here were analyzed using F-test and p-test to observe the consequence of the linear, interaction and quadratic model terms over the responses preselected here (CrI% and yield%).Tables 4 and 5 summarize the ANOVA results obtained for CrI%, (y1) and yield%, (y2); respectively.It was observed based on the magnitude of the linear term that the

Analysis of Variance (ANOVA) Test and Statistical Analysis
An analysis of variance (ANOVA) test was done to observe the competency of the developed model Equations ( 2) and ( 3).The statistical parameters with relevant correlation coefficient (actual and adjusted) after the linear regression analysis are listed in Table 3.The appropriateness of the developed models can be additionally confirmed by observing the values of standard deviation, coefficient of variation (CV) and adequate precision.The experimental R 2 values were in close agreement with the adjusted R 2 values.The magnitudes of standard deviation observed here for both the models were relatively smaller.Furthermore, the values for coefficient of variation (CV) were also only 1.22 and 1.06 for the empirical models Equations ( 2) and (3).This represents the reproducibility of the experimental data with the developed model.Adequate precision is defined as the ratio between signal to noise and for effective development of the polynomial model, the magnitude must be greater than 4.Here for both the responses, it was around 29.82 and 29.20, respectively, reflecting suitable navigation of the design [29].
Experimental data obtained here were analyzed using F-test and p-test to observe the consequence of the linear, interaction and quadratic model terms over the responses preselected here (CrI% and yield%).Tables 4 and 5 summarize the ANOVA results obtained for CrI%, (y 1 ) and yield%, (y 2 ); respectively.It was observed based on the magnitude of the linear term that the ultrasonication power (x 1 ) had a greater impact on CrI% rather than the time (x 2 ) and temperature (x 3 ) (Table 4).However, the yield% (y 2 ) was mostly affected by temperature (x 3 ) (Table 5).Based on the values of F-test as illustrated by Tables 4 and 5; both the regression models (Equations ( 2) and ( 3)) developed here were significant by considering a 5% confidence level.The probability p-values were less than 0.0001 which also ensures that the developed models were significant and confirms an acceptable fitness of the polynomial regression models as illustrated by Equations ( 2) and (3) with the actual experimental observation for each run (Table 2).
The analysis of variance test (ANOVA) observed from Table 4, shows that ultrasonication power (x 1 ), time (x 2 ), temperature (x 3 ), and their interaction terms x 1 x 2, x 2 x 3 and x 1 x 3 together with the quadratic terms of (x 2 1 ) and (x 2 3 ) were significant model terms.Ultrasonication power (x 1 ) had played most influential role over the CrI% (y 1 ) by showing the highest F-value of 444.77.However, time (x 2 ) had least impact on CrI% (y 1 ) compared to the other two factors of time (x 2 ) and temperature (x 3 ), as observed from F-values of Table 4.As illustrated by Table 5, ultrasonication power (x 1 ), time (x 2 ), temperature (x 3 ), and their interaction terms x 2 x 3 together with the quadratic terms of (x 2  3 ) are significant model terms for CNWs yield% (y 2 ).The F-values for temperature (x 3 ) was 218.09 reflecting its prominent effect on yield% (y 2 ) (Table 5).However, power (x 1 ) and time (x 2 ) had relatively less remarkable impact on yield% (y 2 ) (Table 5) as observed from the magnitude of F-test.The interaction effect of power and time (x 1 x 2 = 0.43) had negligible impact whereas power and temperature (x 1 x 3 = 1.81) had reasonable impact on yield% of CNWs (y 2 ).Time and temperature combined (x 2 x 3 ) played a vital role for the yield% (y 2 ) of CNWs.

Process Variables Optimization
For successful commercialization of finally extracted CNWs, both yield% and crystallinity index, CrI% should be maximum.Thus, the numerical optimization test was conducted for both the responses.The optimum condition thus obtained is listed in Table 6.The input variables of power (x 1 ), time (x 2 ) and temperature (x 3 ) were kept within the range selected earlier (Table 1).The experiment was conducted under optimal condition.The analytical data obtained from experiments under optimum condition was compared with the predicted results as suggested by the software (Table 6).The percentage error between the predicted and actual/experimental results was determined.As observed from Table 6, the percentage error for both the responses obtained were negligible supporting suitability of the polynomial models developed earlier (Equations ( 2) and ( 3)) additionally.The optimization ramp is shown by Figure 3.The desirability value obtained was 0.921 reflecting reproducibility of the developed model (Equations ( 2) and ( 3)).

Process Variables Optimization
For successful commercialization of finally extracted CNWs, both yield% and crystallinity index, CrI% should be maximum.Thus, the numerical optimization test was conducted for both the responses.The optimum condition thus obtained is listed in Table 6.The input variables of power (x1), time (x2) and temperature (x3) were kept within the range selected earlier (Table 1).The experiment was conducted under optimal condition.The analytical data obtained from experiments under optimum condition was compared with the predicted results as suggested by the software (Table 6).The percentage error between the predicted and actual/experimental results was determined.As observed from Table 6, the percentage error for both the responses obtained were negligible supporting suitability of the polynomial models developed earlier (Equations ( 3) and ( 4)) additionally.The optimization ramp is shown by Figure 3.The desirability value obtained was 0.921 reflecting reproducibility of the developed model (Equations ( 3) and ( 4)).

Effect of Process Variables for Synthesis of Cellulose Nano-Whiskers (CNWs)
The cumulative effects of ultrasonication power (x1) and time (x2) on CrI% (y1) of finally extracted CNWs is illustrated by a 3D surface mesh with contour plots where the temperature (x3) was kept at center point, in other words constant at 105 °C (Figure 4a).Similarly, Figure 4b was constructed to observe the combined effect of power (x1) and temperature over (x3) the CrI% (y1) where the ultrasonication time (x2) was fixed at 30 min.Overall, the three process parameters used here for ultrasonication have substantial effect on CrI% (y1) as observed from the basic design matrix of Table 2. CrI% of extracted CNWs increased with successive increase of ultrasonication power (x1) and time (x2) in presence of ILs as hydrolyzing solvent (Figure 2a).This was anticipated as the enhancement of power (x1) with hydrolyzing time (x2) would increase the kinetic force of ILs and

Effect of Process Variables for Synthesis of Cellulose Nano-Whiskers (CNWs)
The cumulative effects of ultrasonication power (x 1 ) and time (x 2 ) on CrI% (y 1 ) of finally extracted CNWs is illustrated by a 3D surface mesh with contour plots where the temperature (x 3 ) was kept at center point, in other words constant at 105 • C (Figure 4a).Similarly, Figure 4b was constructed to observe the combined effect of power (x 1 ) and temperature over (x 3 ) the CrI% (y 1 ) where the ultrasonication time (x 2 ) was fixed at 30 min.Overall, the three process parameters used here for ultrasonication have substantial effect on CrI% (y 1 ) as observed from the basic design matrix of Table 2. CrI% of extracted CNWs increased with successive increase of ultrasonication power (x 1 ) and time (x 2 ) in presence of ILs as hydrolyzing solvent (Figure 2a).This was anticipated as the enhancement of power (x 1 ) with hydrolyzing time (x 2 ) would increase the kinetic force of ILs and improve its diffusion rate inside the amorphous segment of cellulose.Enhancing power (x 1 ) during ultrasonication would produce cavitation bubbles which would travel fast towards the interior region of cellulose and collapse there suddenly to initiate the disruption of intra-and intermolecular hydrogen bonds of cellulose.Extended ultrasonication time (x 2 ) would facilitate the physical swelling of crude cellulose (AKL-4) Symmetry 2019, 11, 1148 10 of 20 obtained from the previous pretreatment process.Thus, more surface area would be exposed for ILs to enter and hydrolyze the amorphous domain resulting in higher crystallinity index (CrI%) of the finally extracted CNWs sample (AKL-5).Increasing the temperature up to a certain level increased the CrI% of the CNWs.In that case, viscosity of the ILs will be reduced and its solvation properties will be improved [16,17].Due to reduction of viscosity, significant amount of IL molecules will be dissociated to give [Bmim] + cations which can form complex with negative oxygen atom of -OH groups and form intermediate complexes [16].Relatively higher temperature would reduce the viscosity of ILs which would subsequently increase the diffusion rate of them inside the compact cellulosic matrix [18].The negative part of -HSO 4 would interact with positive H + ions.This would disintegrate the glycosidic bonds between the monosaccharide units of cellulose resulting in nano-dimensional cellulose having higher CrI% [16][17][18].However, whenever the temperature was increased up to the maximum point of 120 • C at different power (x 1 ) and time (x 2 ); the crystallinity index dropped significantly.This happened due to excessive dissolution of cellulose and char formation [17,18].improve its diffusion rate inside the amorphous segment of cellulose.Enhancing power (x1) during ultrasonication would produce cavitation bubbles which would travel fast towards the interior region of cellulose and collapse there suddenly to initiate the disruption of intra-and intermolecular hydrogen bonds of cellulose.Extended ultrasonication time (x2) would facilitate the physical swelling of crude cellulose (AKL-4) obtained from the previous pretreatment process.Thus, more surface area would be exposed for ILs to enter and hydrolyze the amorphous domain resulting in higher crystallinity index (CrI%) of the finally extracted CNWs sample (AKL-5).Increasing the temperature up to a certain level increased the CrI% of the CNWs.In that case, viscosity of the ILs will be reduced and its solvation properties will be improved [16,17].Due to reduction of viscosity, significant amount of IL molecules will be dissociated to give [Bmim] + cations which can form complex with negative oxygen atom of -OH groups and form intermediate complexes [16].
Relatively higher temperature would reduce the viscosity of ILs which would subsequently increase the diffusion rate of them inside the compact cellulosic matrix [18].The negative part of -HSO4 would interact with positive H + ions.This would disintegrate the glycosidic bonds between the monosaccharide units of cellulose resulting in nano-dimensional cellulose having higher CrI% [16][17][18].However, whenever the temperature was increased up to the maximum point of 120 °C at different power (x1) and time (x2); the crystallinity index dropped significantly.This happened due to excessive dissolution of cellulose and char formation [17,18].
(a) (b) Figures 5a,b demonstrate the 3D response surface mesh with contour plots to exhibit the effect of process variables on CNWs yield% (y2).Figure 5a illustrates the collective effect of ultrasonication power (x1) and time (x2) on the yield% (y2) of the finally extracted CNWs where temperature (x3) was kept at fixed level of center point (105 °C). Figure 5b exhibits the effect of ultrasonication power (x1) and temperature (x3) over the yield% (y2), where ultrasonication time (x2) was fixed at center level (30 min).However, yield% (y2) of CNWs showed a decreasing trend with increasing power (x1), time (x2) and temperature (x3).Both the plots exhibited that (Figure 3a,b), temperature (x3) played most vital role for yield% (y2) whereas the other two variables of power (x1) and time (x2) had a relatively moderate impact on it.The yield% (y2) was lowest when power was maintained at 350 watts for time period of 45 min in presence of ILs at a maximum temperature of 120 °C (65.69%,Sample 7) as summarized by Table 2.This was evident due to dissolution of cellulose to other liquid fractions rather than forming solid fraction containing nanocellulose [1,29].Figure 5a,b demonstrate the 3D response surface mesh with contour plots to exhibit the effect of process variables on CNWs yield% (y 2 ). Figure 5a illustrates the collective effect of ultrasonication power (x 1 ) and time (x 2 ) on the yield% (y 2 ) of the finally extracted CNWs where temperature (x 3 ) was kept at fixed level of center point (105 • C). Figure 5b exhibits the effect of ultrasonication power (x 1 ) and temperature (x 3 ) over the yield% (y 2 ), where ultrasonication time (x 2 ) was fixed at center level (30 min).However, yield% (y 2 ) of CNWs showed a decreasing trend with increasing power (x 1 ), time (x 2 ) and temperature (x 3 ).Both the plots exhibited that (Figure 3a,b), temperature (x 3 ) played most vital role for yield% (y 2 ) whereas the other two variables of power (x 1 ) and time (x 2 ) had a relatively moderate impact on it.The yield% (y 2 ) was lowest when power was maintained at 350 watts for time period of 45 min in presence of ILs at a maximum temperature of 120 • C (65.69%, Sample 7) as summarized by Table 2.This was evident due to dissolution of cellulose to other liquid fractions rather than forming solid fraction containing nanocellulose [1,29]

Surface Morphological Studies
The surface morphological features of AKL-1, AKL-2, AKL-3, AKL-4 and AKL-5 samples are illustrated by Figure 6.The surface of AKL-1 was relatively even and smooth as there was a coating of pectin, wax, pigments, lignin, and hemicellulose over its surface (Figure 6a).After solvent extraction process, the smooth coating of pectin and wax was removed to obtain AKL-2 [28].The surface of AKL-2 became rough and the sample contained some minor amounts of flakes or depositions (Figure 6b).Microwave heating in presence of KOH caused partial delignification and removed the residual lignin and hemicellulose fraction from AKL-2.Thus, the sample AKL-3 became uneven and to some extent, microfibrils were released from the compact binding of lignin (Figure 6c).Bleaching process had further removed the lignin and hemicellulose.Thus, long cellulosic strands were separated from each other, defragmented to smaller sizes and are visible in Figure 6d.At this stage, more micro-fibrillated cellulose was being released (AKL-4) [29,30].

Surface Morphological Studies
The surface morphological features of AKL-1, AKL-2, AKL-3, AKL-4 and AKL-5 samples are illustrated by Figure 6.The surface of AKL-1 was relatively even and smooth as there was a coating of pectin, wax, pigments, lignin, and hemicellulose over its surface (Figure 6a).After solvent extraction process, the smooth coating of pectin and wax was removed to obtain AKL-2 [28].The surface of AKL-2 became rough and the sample contained some minor amounts of flakes or depositions (Figure 6b).Microwave heating in presence of KOH caused partial delignification and removed the residual lignin and hemicellulose fraction from AKL-2.Thus, the sample AKL-3 became uneven and to some extent, microfibrils were released from the compact binding of lignin (Figure 6c).Bleaching process had further removed the lignin and hemicellulose.Thus, long cellulosic strands were separated from each other, defragmented to smaller sizes and are visible in Figure 6d.At this stage, more micro-fibrillated cellulose was being released (AKL-4) [29,30].
After ultrasonication, the cavitation effect caused defibrillation as well as defragmentation via ruptures of intra-molecular and inter-molecular hydrogen bonding to obtain CNWs (AKL-5).However, the samples were aggregated after freeze drying and can be observed from FESEM images of AKL-5 (Figure 6e).The surfaces of AKL-4 and AKL-5 were eroded due to hydrolysis.The mechanical force of ultrasonication caused micro bubbles to form which moved with high speed and aided in defragmentation of the micro-fibers to nano size.The nano-dimensional sample thus obtained (CNWs) contained lot of -OH groups (AKL-5) and was found to be agglomerated after the drying process [29][30][31][32][33].After ultrasonication, the cavitation effect caused defibrillation as well as defragmentation via ruptures of intra-molecular and inter-molecular hydrogen bonding to obtain CNWs (AKL-5).However, the samples were aggregated after freeze drying and can be observed from FESEM images of AKL-5 (Figure 6e).The surfaces of AKL-4 and AKL-5 were eroded due to hydrolysis.The mechanical force of ultrasonication caused micro bubbles to form which moved with high speed and The HRTEM and AFM images after the ultrasonication process in presence of ILs were taken and shown by Figure 7.The CNWs obtained were needle-shaped and overlapping each other to form larger aggregates (Figure 7a).Microwave-assisted heating with KOH and subsequent peroxide bleaching process yielded crude cellulose which contained amorphous and crystalline domain.The anions and Symmetry 2019, 11, 1148 13 of 20 cations present in hydrolyzing medium have strong influence for dissolution of amorphous region of cellulose to obtain purified CNW crystals having high crystallinity index.During ultrasonication BmimHSO 4 would form [Bmim] + cations which would interact with the negative -OH groups of cellulose.The negative anions of [HSO 4 ] -would interact β-1-4 glycosidic bonds.The positive cations of the ILs have electron rich π system which could attack the negative oxygen atoms of β-1-4 glycosidic bonds [34].This would cause disintegration of intra-molecular hydrogen bonding between two cellulosic chains.It would further initiate the rupture of β-1-4 glycosidic linkages inside the single, long chain of cellulose resulting defragmented, smaller CNW crystals [35,36].Most of the CNW crystals finally obtained had the length of 80-110 nm and width of 15-20 nm (Figure 7c).The AFM images (Figure 7b) of AKL-5 clearly showed the presence of needle-shaped cellulosic nano crystal (CNWs) samples finally after ultrasonication at optimal condition [37].
The HRTEM and AFM images after the ultrasonication process in presence of ILs were taken and shown by Figure 7.The CNWs obtained were needle-shaped and overlapping each other to form larger aggregates (Figure 7a).Microwave-assisted heating with KOH and subsequent peroxide bleaching process yielded crude cellulose which contained amorphous and crystalline domain.The anions and cations present in hydrolyzing medium have strong influence for dissolution of amorphous region of cellulose to obtain purified CNW crystals having high crystallinity index.During ultrasonication BmimHSO4 would form [Bmim] + cations which would interact with the negative -OH groups of cellulose.The negative anions of [HSO4] -would interact β-1-4 glycosidic bonds.The positive cations of the ILs have electron rich ᴫ system which could attack the negative oxygen atoms of β-1-4 glycosidic bonds [34].This would cause disintegration of intra-molecular hydrogen bonding between two cellulosic chains.It would further initiate the rupture of β-1-4 glycosidic linkages inside the single, long chain of cellulose resulting defragmented, smaller CNW crystals [35,36].Most of the CNW crystals finally obtained had the length of 80-110 nm and width of 15-20 nm (Figure 7c).The AFM images (Figure 7b) of AKL-5 clearly showed the presence of needle-shaped cellulosic nano crystal (CNWs) samples finally after ultrasonication at optimal condition [37].
(c) The mechanical process of ultra-sonication leads to solvo-dynamic shear force, that is, the cavitation effect inside the hydrolyzing medium of ILs used here.Thus, lots of micro-bubbles are formed in presence of heat during ultrasonication.The growth of micro-bubbles increases with lapse of time and after a certain time, it collapses causing shock waves inside the hydrolyzing medium.This initiates the erosion as well as disintegration of micro-sized fiber to nano-size during the ultrasonication process [35,36,38].

XRD Analysis
Figure 8 exhibits the XRD pattern of starting lignocellulosic residues (AK), solvent extracted, microwave-assisted KOH pretreated, bleached and ultrasonicated samples in presence of Bmim-HSO 4 ILs under optimum condition.The crystallinity index along with the crystallite size of the extracted sample after each treatment step was calculated and the results obtained are summarized in Table 7.The mechanical process of ultra-sonication leads to solvo-dynamic shear force, that is, the cavitation effect inside the hydrolyzing medium of ILs used here.Thus, lots of micro-bubbles are formed in presence of heat during ultrasonication.The growth of micro-bubbles increases with lapse of time and after a certain time, it collapses causing shock waves inside the hydrolyzing medium.This initiates the erosion as well as disintegration of micro-sized fiber to nano-size during the ultrasonication process [35,36,38].

Sample Crystallinity Index (%) Crystallites Sizes (nm)
AKL It was observed that after each treatment the crystallinity index was increased subsequently but the crystallites sizes were reduced.Solvent extraction process had removed pectin, wax, pigments, etc. from untreated sample of AKL-1 to yield AKL-2 with negligible change in crystallinity index and crystallite sizes.KOH pretreatment using microwave heating, lignin, hemicellulose and amorphous regions of cellulose were partially degraded (AKL-3).Peroxide bleaching of AKL-3 sample further reduced the amount of lignin to yield AKL-4 where the amorphous region of cellulose was more exposed for selective hydrolysis during ultrasonication.After ultrasonication, in presence of Bmim-HSO 4 , ILs, the sample AKL-5 (CNWs) showed comparatively high crystallinity index with smaller crystallite size than the other samples.The cavitation effect of ultrasonication had initiated the removal of the intra-and intermolecular hydrogen bonds, releasing the micro fibrillated cellulosic fiber resulting in a higher crystallinity index than the untreated (AKL-1), solvent extracted (AKL-2), KOH treated (AKL-3) and bleached (AKL-4) samples.
All the XRD diffraction patterns showed two major peaks approximately at 2θ = 14 • to 16 • (110) and 22 • to 24 • (200) representing typical cellulose I structure [39,40].The results demonstrated the presence of monoclinic crystalline cellulose having Iβ structure.It was reported earlier that cellulose I structure contains around 85 repetitive units β-D-glucopyranose molecules linked together by strong networking of intermolecular hydrogen bonds [41].However, after ultrasonication in presence of ILs, the intra-and intermolecular hydrogen bonds will be broken, and the amorphous domain will be hydrolyzed to finally yield the cellulose nano-whiskers (CNWs) with a higher crystallinity index.Thus, a sharper peak at 2θ = 22-24 • (200) was observed for the samples after being ultrasonicated in presence of ILs.

Thermogravimetric Analysis
Thermogravimetric analyses of all the samples (AKL-1, AKL-2, AKL-3, AKL-4 and AKL-5) were carried out and are shown by Figure 9.The thermal stability of each sample changed after the consecutive steps of treatment.Thermal stability of the synthesized sample basically depends on the extent of lignin proportion and intra-and intermolecular hydrogen bonding present inside the lignocellulosic matrix of the biomass residues.
the crystallites sizes were reduced.Solvent extraction process had removed pectin, wax, pigments, etc. from untreated sample of AKL-1 to yield AKL-2 with negligible change in crystallinity index and crystallite sizes.KOH pretreatment using microwave heating, lignin, hemicellulose and amorphous regions of cellulose were partially degraded (AKL-3).Peroxide bleaching of AKL-3 sample further reduced the amount of lignin to yield AKL-4 where the amorphous region of cellulose was more exposed for selective hydrolysis during ultrasonication.After ultrasonication, in presence of Bmim-HSO4, ILs, the sample AKL-5 (CNWs) showed comparatively high crystallinity index with smaller crystallite size than the other samples.The cavitation effect of ultrasonication had initiated the removal of the intra-and intermolecular hydrogen bonds, releasing the micro fibrillated cellulosic fiber resulting in a higher crystallinity index than the untreated (AKL-1), solvent extracted (AKL-2), KOH treated (AKL-3) and bleached (AKL-4) samples.
All the XRD diffraction patterns showed two major peaks approximately at 2θ = 14° to 16° (110) and 22° to 24° (200) representing typical cellulose I structure [39,40].The results demonstrated the presence of monoclinic crystalline cellulose having Iβ structure.It was reported earlier that cellulose I structure contains around 85 repetitive units β-D-glucopyranose molecules linked together by strong networking of intermolecular hydrogen bonds [41].However, after ultrasonication in presence of ILs, the intra-and intermolecular hydrogen bonds will be broken, and the amorphous domain will be hydrolyzed to finally yield the cellulose nano-whiskers (CNWs) with a higher crystallinity index.Thus, a sharper peak at 2θ = 22-24° (200) was observed for the samples after being ultrasonicated in presence of ILs.

Thermogravimetric Analysis
Thermogravimetric analyses of all the samples (AKL-1, AKL-2, AKL-3, AKL-4 and AKL-5) were carried out and are shown by Figure9.The thermal stability of each sample changed after the consecutive steps of treatment.Thermal stability of the synthesized sample basically depends on the extent of lignin proportion and intra-and intermolecular hydrogen bonding present inside the lignocellulosic matrix of the biomass residues.Initially all the samples lost some weight at temperature 100-120 °C which was due to the evaporation of adsorbed moisture.The first steps of weight loss observed for all the samples (AKL-1, AKL-2, AKL-3, AKL-4 and AKL-5) were 3.5%, 3.1%, 2.8%, 2.3% and 1.9%, respectively.During the second degradation step, around 76.74%, 80.68%, 81.31%, 83.26% and 85.02% of weight loss, respectively, were observed for all the samples (AKL-1, AKL-2, AKL-3, AKL-4 and AKL-5).After that, decomposition of cellulose, hemicellulose and residual lignin took place around the Initially all the samples lost some weight at temperature 100-120 • C which was due to the evaporation of adsorbed moisture.The first steps of weight loss observed for all the samples (AKL-1, AKL-2, AKL-3, AKL-4 and AKL-5) were 3.5%, 3.1%, 2.8%, 2.3% and 1.9%, respectively.During the second degradation step, around 76.74%, 80.68%, 81.31%, 83.26% and 85.02% of weight loss, respectively, were observed for all the samples (AKL-1, AKL-2, AKL-3, AKL-4 and AKL-5).After that, decomposition of cellulose, hemicellulose and residual lignin took place around the temperature range 230-370 • C. At a higher temperature range, around 300-400 • C, pyrolysis of cellulose mainly takes place [39].The third degradation step is observed after that, and was recorded until 1000 • C. Lignin decomposition is relatively difficult, and it starts from 200 • C until 1000 • C.This is due to presence of high temperature resistant phenyl groups inside the lignin molecules.
Table 8 summarizes the major decomposition peak (DTG max ), crystallinity index, amount of char and weight loss for all the samples obtained here.Table 8 illustrates that the thermal decomposition temperature of the synthesized cellulosic substrates was reduced after each treatment.This is obviously due to elimination of lignin as well as dissolution of hemicellulose from the lignocellulosic samples [1,29].However, the IL-assisted ultrasonicated sample (AKL-5) showed lowest decomposition temperature compared to the other samples of AKL-1, AKL-2, AKL-3 and AKL-4.The mechanical force of ultrasonication with the application of ILs was able to defibrillate the crude cellulose sample obtained after the bleaching (AKL-4) process.Presence of acidic ILs during the ultrasonication process would further enhance the breaking of intra-molecular hydrogen bond to release the microfiber and convert it to nano size by disrupting the intermolecular hydrogen bonds as well.Thus, the finally extracted CNWs (AKL-5) sample would have a larger number of free ends of chain, resulting lower thermal stability [40][41][42].The amount of char was maximum for untreated AKL-1 sample reflecting the presence of hemicellulose and lignin inside the cellulosic fiber.However, after the successive treatment of solvent extraction would reduce the amount of pectin and wax from AKL-1.Microwave heating in presence of KOH and bleaching would decompose the hemicellulose and lignin fraction from the starting sample of AKL-1 to yield AKL-3 and AKL-4, respectively.Thus AKL-3 and AKL-4 would have lower char residues.The mechanical force of ultrasonication would further reduce the trace amount of lignin present in AKL-4 and would partially hydrolyze the amorphous segment of cellulosic chain to obtain the sample of AKL-5 having cellulose nano-whiskers (CNWs) with high crystallinity index.The combined force of mechanical treatment (ultrasonication) and chemical (ionic liquid, IL) hydrolyzing medium would enhance defibrillation and defragmentation of long cellulosic chains by disrupting the intraand intermolecular hydrogen bonds, respectively.The sample AKL-5 thus obtained will have lower amount of char residues.Similar observation was reported earlier for extracting nano-cellulosic fiber from sisal fiber and pulp derived from sugar beet [30,43].

Surface Functional Groups Analysis
FTIR spectroscopic analysis was carried out to observe the surface functional groups of starting lignocellulosic residues of AKL-1 and other samples of AKL-2, AKL-3, AKL-4 and AKL-5 after each treatment.All the samples used here had some common peaks which were almost similar [28,34,44].The absorbance peak around 3400-3300 cm −1 represented the stretching of -OH groups [9].The presences of -OH groups in cellulosic materials enable them to form different kinds of intra-and intermolecular hydrogen bonding between the different cellulosic strands and each cellulosic chain, respectively [45].The stretching of -CH functional groups was visible in all the samples around 2900-2800 cm −1 [46].Peaks around 1620-1650 cm −1 were related to O-H bending vibrations which showed the presence of absorbed moisture in all the samples [47].The minor peaks obtained around 1420-1450 cm −1 showed the intermolecular hydrogen bonding in C-6 groups of cellulosic substrates [44].Peaks around 1510 to 1530 cm −1 exhibited the presence of aromatic ring vibration in AKL-1, AKL-2 and AKL-3.In AKL-4 and AKL-5, the peak disappeared showing complete elimination of lignin from the synthesized sample.The bending vibration of -CH groups as well as presence of -CO groups were visible for all the samples in the range 1320-1380 cm −1 [44].All the samples contained the peaks around 1045 cm −1 reflecting the presence of C-O-C stretching vibration of pyranose ring and glycosides linkages between two monosaccharaides units of cellulose [48].The peaks around 800-900 cm −1 typically represented the cellulosic chain which can be ascribed to β-glycosidic linkages between glucose units in cellulosic chain [44].During microwave-induced heating, the alkali acts as microwave energy absorber and initiates easy penetration of solvent inside the cellulosic substrates.Thus, partial dissolution of hemicellulose and lignin takes place due to uniform heating process of microwaves.This makes the surface area of AKL-3 more accessible for the peroxide bleaching process which initiates the defibrillation process.Therefore, in AKL-3 and AKL-4, more micro-and nano-fibrillated fiber were released, which was supported earlier by our FESEM images.Table 9 summarizes the major peaks with their frequency regions usually observed for ligno-cellulosic biomass and the micro-and nano-dimensional cellulose extracted from them.
Table 9. List of FTIR peaks with their frequency region for synthesized samples (AKL-1, AKL-2, AKL-3, AKl-4 and AKL-5).Zeta potential analysis was carried out for the samples (AKL-3, AKL-4, AKL-5).It was observed that the surface charge of the synthesized samples was increasing after each consecutive treatment.Initially the ligno-cellulosic residue (AKL-1) was covered with pectin and wax and the cellulosic matrix was embedded inside the lignin molecules.Thus, for sample AKL-1 and AKL-2, zeta potential values were not measured.However, after KOH pre-treatment, delignification was carried out partially and a fractional amount of cellulosic chain was released containing negative -OH groups over its surface (AKL-3) exhibiting zeta potential values of −19.56 mV.After bleaching, more micro-fibrillated fiber was further released from compact binding of lignin.Thus, more negative-OH groups were exposed resulting in higher zeta potential values (AKL-4) of −23.54 mV.Finally extracted CNWs sample (AKL-5) showed the highest zeta potential value of -29.9 mV.This is expected as ultrasonication in presence of ILs caused disruption of intra-and intermolecular hydrogen bonding of long cellulosic chains, leading to the formation of milky, white, colloidal suspension of CNWs (AKL-5) containing defragmented smaller versions of cellulosic chain having lots of -OH groups [42].The strong mechanical force of ultrasonication and acidic hydrolyzing medium of [Bmim]HSO 4 ILs had initiated the defibrillation as well as defragmentation of long cellulosic chain resulting CNWs [28,37].Thus, more negative -OH groups were exposed resulting higher zeta potential values [49,50].

Conclusions
Ultrasonication in presence of ILs can effectively extract CNWs with high crystallinity index after the microwave-assisted delignification and bleaching of the dried leaves of Adansonia kilima (AK).Presence of ILs under optimum condition of ultrasonication can selectively degrade the amorphous domain of crude cellulose obtained from AK. Application of ILs can enhance the crystallinity index of CNWs up to a satisfactory level.X-ray crystallographic analysis demonstrated that the CrI% (y 1 ) was increasing with successive decrease in particle size of the synthesized sample.HRTEM and AFM images showed development of needle like CNWs.Highly negative zeta potential values of finally extracted sample (AKL-5) clearly reflected higher density of -OH groups over the defragmented nano-sized crystals (CNWs).The experimental design matrix was made based on Box-Behnken (BBD) design and ultrasonication variables were optimized.CrI% (y 1 ) was influenced mostly by ultrasonication power (x 1 ) and temperature (x 3 ), whereas ultrasonication time (x 2 ) had a moderate impact on it.On the other hand, yield% (y 2 ) was affected more prominently by temperature (x 3 ) and time (x 2 ) rather than the power (x 1 ) itself.Model validation with analysis of variance (ANOVA) test was carried out to check the adequacy of the proposed polynomial Equations ( 2) and (3).Theoretical optimization to obtain maximum CrI% (y 1 ) and yield% (y 2 ) was done.The suggested values by the software for CrI% (y 1 ) and yield% (y 1 ) were very close with the actual ones illustrating negligible error percentages of 1.61% and 1.30% for both the responses, respectively.The findings of this research revealed that the BBD design using mathematical modeling and statistical analysis concurrently is appropriate for detecting and optimizing the input variables influencing the IL-assisted ultrasonic extraction process of CNWs from the delignified, dried leaves of Adansonia kilima (AK) plants.

Figure 4 .
Figure 4. 3D surface mesh and contour plots.(a) Combined effects of power (x1) and time (x2); (b) combined effect of power (x1) and temperature (x3) on percentage crystallinity, CrI% (y1) of CNWs when the other two variables were at center level.

Figure 4 .
Figure 4. 3D surface mesh and contour plots.(a) Combined effects of power (x 1 ) and time (x 2 ); (b) combined effect of power (x 1 ) and temperature (x 3 ) on percentage crystallinity, CrI% (y 1 ) of CNWs when the other two variables were at center level.
Figure5a,b demonstrate the 3D response surface mesh with contour plots to exhibit the effect of process variables on CNWs yield% (y 2 ).Figure5aillustrates the collective effect of ultrasonication power (x 1 ) and time (x 2 ) on the yield% (y 2 ) of the finally extracted CNWs where temperature (x 3 ) was kept at fixed level of center point (105 • C).Figure5bexhibits the effect of ultrasonication power (x 1 ) and temperature (x 3 ) over the yield% (y 2 ), where ultrasonication time (x 2 ) was fixed at center level (30 min).However, yield% (y 2 ) of CNWs showed a decreasing trend with increasing power (x 1 ), time (x 2 ) and temperature (x 3 ).Both the plots exhibited that (Figure3a,b), temperature (x 3 ) played most vital role for yield% (y 2 ) whereas the other two variables of power (x 1 ) and time (x 2 ) had a relatively moderate impact on it.The yield% (y 2 ) was lowest when power was maintained at 350 watts for time period of 45 min in presence of ILs at a maximum temperature of 120 • C (65.69%, Sample 7) as summarized by Table2.This was evident due to dissolution of cellulose to other liquid fractions rather than forming solid fraction containing nanocellulose[1,29].

Figure 5 .
Figure 5. 3D surface mesh and contour plots.(a) Combined effects of power (x1) and time (x2); (b) combined effect of power (x1) and temperature (x3) on yield% (y2) of CNWs when the other two variables were at center level.

Figure 5 .
Figure 5. 3D surface mesh and contour plots.(a) Combined effects of power (x 1 ) and time (x 2 ); (b) combined effect of power (x 1 ) and temperature (x 3 ) on yield% (y 2 ) of CNWs when the other two variables were at center level.
Author Contributions: Z.Z.C. conceived and designed the experiments.A.J., K.K., O.A. have performed the experiment.The data was analyzed by Z.Z.C., M.M.R. and M.A.N.A.H., Y.A.W. and M.R.B.J. has provided all reagents and chemicals support.I.A.B., T.M.Y.K., G.A.A. and S.K. had contributed for manuscript editing, publication charges and characterizations.Z.Z.C., R.R.R.C., M.A. and A.J. had compiled the results and written the paper.R.R.R.C and G.A.A. had calculated particle size distribution and provided the software to analyze the data.Funding: The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work through a research groups program under grant number R.G.P. 2/11/39.The research was also conducted under 17AET-RP44C under Dr. Zaira Zaman Chowdhury as Principal Investigator from University of Malaya, Malaysia.

Table 1 .
Input and output parameters with their levels using BBD design for extraction of cellulose nano-whiskers (CNWs) (AKL-5).

Table 3 .
Statistical accuracy test for proposed polynomial models (crystallinity index% and CNWs yield%).

Table 3 .
Statistical accuracy test for proposed polynomial models (crystallinity index% and CNWs yield%).

Table 4 .
Analysis of variance test (ANOVA) test for regression polynomial equation for crystallinity index, CrI% (y 1 ) of CNWs.

Table 5 .
Analysis of variance test (ANOVA) test for regression polynomial equation for yield% (y 2 ) of CNWs.

Table 6 .
Process parameter optimization for ultrasonication in presence of ILs.

Table 6 .
Process parameter optimization for ultrasonication in presence of ILs.