Cellulose has applications in several areas due to its unique properties that are compatible with the growing search for materials from renewable sources. It is a biopolymer that is the main constituent of the cell wall of plants, and its extraction process follows the top-down strategy and can be mechanically or chemically extracted. Cellulose nanocrystals (NCC) are one of the materials obtained from this polymer, which can be used for several studies due to its versatility and characteristics. So, many studies report the obtaining cellulose nanocrystals (NCC) by hydrolysis with strong mineral acids as the most commonly used methodology, such as formic acid, acetic acid, sulfuric acid, hydrochloric acid, phosphoric acid or acid mixtures  - ; in 2012, Espinosa  investigated the influence of the three later acids behaved on the cellulose hydrolysis, as well as the influence on the dispersion of NCCs in different solvents and the differences in thermal stability. It was observed that NCCs obtained from hydrochloric acid treatment, presented higher thermal stability when compared to the ones obtained from the other acids, which was attributed to the absence of residual sulfates and phosphate groups.
Some issues are attributed to acid hydrolysis process, such as the recovery of the used acid, large amounts of base employed for neutralization, large amount of water, equipment corrosion and environmental pollution  . Alternatively, other reagents such as phosphotungstic acid   and dicarboxylic acids are used, which are more easily recoverable by the cold crystallization process and still have enough solubility at elevated temperatures, being able to hydrolyze the amorphous part of the cellulose and thus, be able to produce NCCs . Another major advantage related to the use of these dicarboxylic acids is the possibility of cellulose surface modification through the Fischer-Speier Esterification reaction, as reported by     .
Many works in the literature describe that the hydrolysis of cellulose with sulfuric acid can be found; where concentration, time and temperature of the reaction are cited as fundamental factors for the successful of the process. The break of glycosidic bonds in cellulose occurs preferentially in the amorphous region, while the crystalline part is less affected . Wang and co-workers  studied the kinetics of acid hydrolysis and classified it in three stages: the first one occurred in the amorphous region, which promoted depolymerization of cellulose fiber; in the second one, they observed the glucose     formation from cellulose, needing more severe hydrolysis conditions; the third one directly converts cellulose to monosaccharides.
Sulfuric acid concentration and temperature were crucial in reducing the degree of polymerization (DP) of cellulose and yield of NCC, since with 16 wt% acid, the DP decreases from 700 to 200 along temperature increase from 45˚C to 85˚C, maintaining the yield; by increasing the acid concentration to 64% the DP was not altered, but the yield decreased due to competition between depolymerization and degradation reactions  . Wang’s group  optimized this process by balancing the yield of NCC and the percentage of cellulosic solid residue (CSR) generated, and this solid was transformed into NCC by the addition of one additional reaction step: the mechanical defibration stage. Although 60% acid concentration provided a higher percentage of NCC  , Chen’s group  found that 58 wt% acid concentration was the transition point for the highest yield of NCC, despite of the higher lengths achieved.
NCCs and cellulose nanofibrils (CNFs) have interesting properties beyond their mechanical reinforcement to nanocomposites; those properties are related to the numerous possibilities of chemical modifications over hydroxyls located at cellulose surface, which can improve the polymer-cellulose interaction, decreasing the tendency for agglomeration and can add new functionality to cellulose for different applications such as: biomedical engineering, electronics and energy and water treatment sectors  . Among the many modifications, the anchoring of inorganic particles in cellulose has been investigated to confer antimicrobial, anti-UV and flame retardant properties, as in textiles, and to extend the thermal stability range, which allows to process those materials at higher temperatures   . Due to its good thermal stability and biocidal property, zinc oxide is a promising material, besides having a wide absorption range of UV radiation in the 260 - 370 nm range   it possesses good photo stability, low toxicity and biocompatibility, which makes it an interesting material for packaging and medical area  . Among the methodologies for obtaining inorganic nanoparticles, the sol-gel methodology is advantageous because it is a low-cost process, with good control over the parameters and reaction occurs at lower temperatures when compared to other methodologies. Therefore, zinc oxide can be obtained in situ from a chemical precursor undergoing hydrolysis and condensation reaction. These reactions are controlled by the pH and temperature of the heat treatment performed at the end of the reaction    .
Much research is oriented to increase the yield of NCCs, reduce cellulose decomposition by monitoring the production of NCC and cellulosic solid residues (CSR) obtained from acid hydrolysis processes  ; other efforts are focused in reusing cellulosic tailings from the final bleached kraft pulp process . The CSR is the precipitate obtained after the centrifugation at the end of the hydrolysis reaction, which being submitted to another mechanical step can render NCC, reducing the amount of waste generated, according to Wang . Normally the isolation process could be improved by ultrasonication , ultra-turrax , the mechanical process by ball mill equipped  as well as combine ultrasonication and ultra-turrax showed better results when used simultaneous . The purpose of this work was to perform the acid hydrolysis of microcrystalline cellulose in a single step, in order to achieve that goal, the mechanical stirrer (ultra-turrax) was used before and after the reaction. The chemical modification with zinc oxide is expected to improve the thermal stability of the cellulose and to confer photocatalytic activity to the composite.
2.1. Materials and Reagents
Microcrystalline cellulose (MCC) was purchased from Synth (density 0.26 - 0.32 g/cm3, polymerization grade 350), nanocrystalline cellulose was purchased from Celluforce (density 0.7 g/cm3) (NCC_C), sulfuric acid, sodium hydroxide, ethyl alcohol and di-hydrated zinc acetate were obtained from Sigma-Aldrich. The reagents were used as received.
2.2. Microcrystalline Cellulose Acid Hydrolysis
Cellulose nanocrystals were obtained from acid hydrolisis of microcrystalline cellulose. First, a suspension of MCC (10%w/v) with 32%v/v sulfuric acid was produced and homogenized on a mechanical stirrer (IKA RW20) for 30 min at 17,000 rpm. Then the acid concentration was increased up to 64%v/v, and hydrolysis continued for a further 2 h under reflux at 50˚C. During the addition of sulfuric acid, the system was cooled with a cold ice bath.
By the end of the reaction time, 300 ml of cold water was added to stop the hydrolysis reaction, after which, the suspension (identified as NCC_H) was taken to the mechanical stirrer for a further 20 min at 17,000 rpm. Finally, dialysis (molecular weight cut-off 14,000) was performed to neutralize the acid suspension, and the suspension stored in a refrigerator for approximately 24 h.
2.3. Sol-Gel Reaction at Cellulose Surface
The previously obtained cellulose suspension (200 ml) was submitted to reaction with the same volume of ethanol, adjusting the pH to 8 with sodium hydroxide solution (5 M).
The metal precursor used was zinc acetate dihydrate, which, being previously dissolved in ethanol for 24 h, reacted with cellulose in a 1:1 ratio. The solution containing the precursor was dropped into the cellulose suspension, while pH was monitored and corrected to pH 8. After the addition was complete, the reaction was heated at 80˚C under reflux for 2 h. Then, part of the solvent was removed by evaporation, the precipitate washed with ethanol and centrifuged (4x at 3500 rpm for 10 min). At each washing step the cellulose was re-suspended with the aid of tip ultrasonic homogenizer (Ultronique QR500, frequency of 20 KHz). The sample was oven dried at 120˚C for 1 h .
2.4.1. Zeta Potential
The Zeta potential of the samples was measured on the NICOMP DLS/ZLS, model Nano Z3000, by the Electrophoretic light scattering technique. Concentration suspensions of 0.001 wt% were prepared, which before being analyzed were dispersed for 1 min in Ultrasonic QR500 ultrasound, then tripled analyzed at room temperature.
2.4.2. Scanning Helium Ion Microscope
Helium beam assisted electron microscopy images were acquired on the Zeiss Orion nanofab model equipped with tungsten filament, operating at a constant voltage of 30 kV. The sample preparation involved suspension of 0.5 mg cellulose in 1 ml acetone, which was then diluted 50× and dispersed in ultrasonic bath for 15 min. Then a drop was deposited on the silicon substrate, and the samples analyzed. As these are insulating samples, floodgun was used to neutralize the surface charge.
2.4.3. Atomic Force Microscope
Cellulose samples were prepared in aqueous suspensions (1%w/v), these suspensions were diluted 100× and sonicated for 15 min. A drop of the diluted suspension was deposited on the surface of the freshly cleaved mica, then dried under N2 flux. Images were acquired on a Dimension Icon® (Bruker, Santa Barbara-CA) AFM instrument using Peak Force Tapping QNM scan mode, cantilever model HQ-300-AU (Oxford Instruments®, UK), spring constant 40 N/m, average radius of 10 nm and 300 kHz frequency.
2.4.4. X-Ray Diffraction
A Rigaku Ultima IV X-ray diffractometer, operating with voltage and current of 40 Kv and 30 mA, respectively, was used for X-ray diffraction characterization. X-rays were obtained from a CuKα source (ʎ = 1.5418 Å). The analysis conditions for the MCC, NCC_C and NCC_H samples were made at 0.01˚/s, 2θ between 10˚ and 40˚, with 2θ for the NCC-Zn sample between 10˚ and 80˚. The analyses were carried out at room temperature.
The size of the crystallite, Ʈ, perpendicular to the crystalline (2 0 0) plane, peak 2θ at 22.5˚, was calculated by the Scherrer equation:
where K is the correction value 0.9, λ is the wavelength of the X-rays (1.54 Å), β is the FWHM of the diffractogram peak (radians), and θ is the peak diffraction angle.
The crystallinity index, CrI, was calculated by the Segal Method whose equation is:
where It is the total intensity of peak 200 (2θ at 22.5˚), Ia is the intensity of the amorphous contribution at 2θ = 18.4˚.
2.4.5. Fourier Transform Infrared Spectroscopy (FTIR)
Analyses were performed in a Nicolet 6700-FT-IR equipment, using KBr pellet, 4 cm−1 spectral resolution and 16. Sample was dried for 48 h before analysis.
2.4.6. Thermogravimetric Analysis (TGA)
The thermal stability of the materials was analyzed on TA Instruments model TGA Q-50, under nitrogen atmosphere, heating rate of 10˚C/min in the temperature range of 30˚C - 700˚C.
2.4.7. Nuclear Magnetic Resonance Time Domain (TD-NMR)
Relaxation time measurements were performed on a MARAN Ultra NMR (Oxford Instruments) NMR spectrometer at 30˚C ± 2˚C with 0.54 T (23 MHz for 1H). The samples were dried for 48 h in a vacuum oven at 60˚C and placed directly into the 18 mm diameter glass tube for analysis.
The Magic Sandwich Echo pulse sequence was used to acquire the fully refocused Free Induction Decay (MSE-FID). All relaxation curves were acquired 512 points spaced by 0.5 microseconds, with a repeating time of 1 second, 128 scans and 34 dB receiver gain.
The obtained signals were adjusted according to the following equation:
where FR, and FM are the fractions of the rigid (crystalline) and mobile (amorphous) regions, respectively. T2* is the transverse relaxation time of one of these fractions, ν is a constant component of sine wave oscillation based on Van Vleck’s second and fourth moments. βM corresponds to the Weibullian function exponents ranging from close to zero to 1. μa is centroid of the Abragamian and k is the offset or baseline of the relaxation signal that compensates for the influence of background signal during nonlinear adjustment.
Considering a = 1/T2*FR and b = 2πν, the mean value of the second dipolar coupling moment can be obtained through the equation:
3. Results and Discussion
3.1. Zeta Potential Measurement
Despite being an important parameter to evaluate the stability of suspensions, this technique was employed in order to follow the influence of the surface modification influenced the cellulose surface charge. NCC_H presented a zeta potential around −44 mV and reached a value close to 0 mV when modified with zinc. The value for NCC_H is consistent with the literature, since there is a direct relationship of this value with hydrolysis times . The negative value for NCC_H is the result of the repulsion forces between the sulphate groups, promoting the good dispersion of these nanoparticles in aqueous suspension    . The sol-gel reaction was performed with cellulose in cellulosic suspension, avoiding agglomerations resulting from the lyophilization process   .
3.2. Morphological Studies
The morphology of the materials was examined by helium ion microscopy (HIM) and atomic force microscopy (AFM). Microcrystalline cellulose (Figure 1(a)) had a rod-like shape, rough surface and irregularities, such as the presence of fibrils on smaller scales that pierced the cellulosic fibers (see Figure 1(a)), which may facilitate the accessibility of acid during hydrolysis . After the MCC hydrolysis reaction (Figure 1(d)), cellulose tended to agglomerate in a round-shaped cluster with a thin layer around it. AFM micrograph (Figure 1(e)) isolated two cellulose nanocrystals of 100 nm length and 3 nm diameter, revealing that the experimental hydrolysis conditions were able to produce cellulose nanocrystals, and its morphology, size and diameter were similar to those found in the literature  . Comparing the images of Figure 1(d) and Figure 1(b)
Figure 1. HIM images of MCC (a), NCC_C (b), NCC_H (d), and NCC_ZnO (f); topographic AFM image of sample NCC_C (c), NCC_H (e) and NCC_ZnO (g).
(commercially obtained cellulose nanocrystals), both samples showed the same agglomeration tendency in round-shaped clusters and the presence of a thin layer around the cluster. The topography image for sample NCC_C (Figure 1(c)), revealed that this clustering was formed from the network of randomly interconnected cellulose nanoparticles. This comparison between NCC_H and NCC_C permits to state that the methodology promotes the formation of cellulose nanocrystals, although agglomerated in micrometrical clusters.
HIM image Figure 1(f) confirms the presence of zinc nano oxide on the cellulose surface and in its non-associated form into reaction medium. It was observed that the zinc structure formed is porous and polydisperse (image not shown here), which may be related to the lack of control in the dropping of the metal precursor solution, pH fluctuations during the reaction and the absence of capping agent    . Fibril-like structures were identified by microscopy, which correspond to the type of nanofibrillated cellulose. These structures are probably the result of partially hydrolyzed structures present in the CSR, which were defragmented in fibrils with higher aspect ratio by mechanical stirring of Ultra-turrax and the negative surface charges, which agrees with previous observations by  .
The AFM image (Figure 1(g)) complements the information in Figure 1(f), since the roughness on the cellulose surface as well as distinct morphologies are visible, which were attributed to the zinc oxide (round shape) anchored on the surface of the CNF.
3.3. X-Ray Diffraction
The X-ray diffraction pattern shown in Figure 2 corresponds to type I cellulose; all the samples have peaks at 2θ = 15˚, 16˚, 21˚ and 22.8˚, such peaks correspond to the planes (110), (110), (012), and (200), respectively. The MCC and NCC_H samples exhibit a small 12.4˚ shoulder that can be attributed to type II cellulose. Sèbe  confirmed the polymorphism of in NCCs (cellulose type I and II) obtained from hydrolysis with sulfuric acid, which is related to the concentration and time of acid addition. Jin’s research group  also confirmed the coexistence of the two types of cellulose (I and II) by attributing a peak at 2θ equal to 12˚ to type II pulp.
Gaussian deconvolution of the diffractograms allowed to better evaluate the profile of the crystalline planes before and after hydrolysis, making the changes of the crystallographic patterns more noticeable. There is a greater crystalline contribution to the MCC when compared to the NCC_H, if we consider the peak intensity at (110) and (200). The degree of crystallinity was calculated by the Segal method, which, as described by  is an efficient and widespread method in the literature. The data in Table 1 corroborate the abovementioned hypothesis, showing that the crystallinity is slightly lower for NCC_H. Therefore, although the mechanical effect benefited the disintegration of MCC the crystal structure of NCC_H has been preserved, acid hydrolysis under these
Figure 2. Diffractograms of MCC (a), NCC_C (b), NCC_H (c) and NCC_ZnO (d) samples. The deconvolution of the samples (a), (b) and (c) was performed, which evidenced the crystalline and amorphous contributions in the samples.
Table 1. Crystallinity index (CrI) and crystallite size (Ʈ200) data obtained from XRD.
conditions (64%, 50˚C, 120 min) favored degradation reactions, with cellulose being reduced to sugars or soluble by-products     . The divergence in the value of CrI and Ʈ200 for NCC_H and NCC_C is at the origin of the fiber, since the DP varies from source to source, finally reflecting the length of the obtained cellulose crystal, as observed by Jia’s group . The values found for Ʈ200 (Table 1), agree with report by Hamad , despite the fact of MCC and NCC_H values were identical.
The well-defined diffraction peaks at 2θ = 31.9˚ (100), 34.7˚ (002), 36.4˚ (101), 47.7˚ (102), 56.7˚ (110), 63.1˚ (103), 68.2˚ (200) and 69.3˚ (201) correspond to the hexagonal structure of the wurtzite, as shown in Figure 2(d)     . Still in the same Figure 2(d), the less intense diffraction peaks of cellulose are encountered, and it is possible to state that there is no presence of intermediate structures (ε-Zn (OH) 2), as the ones specified by  .
No adjustments could be made for the NCC_ZnO sample due to the high zinc oxide diffraction signal.
3.4. Fourier Transform Infrared Spectroscopy (FTIR)
Changes in the chemical environment and the appearance of new absorption bands after chemical modifications were observed in the FTIR spectrum (Figure 3). The wide absorptions in the range of 3800 - 3000 cm−1 and in the region of 1638 cm−1 were attributed to the stretching and bending vibration of the OH groups present in the cellulose, even though it was dried   . The bands found in the range between 1420 - 897 cm−1 correspond to the chemical structure of cellulose, such as symmetric CH2 bonding, C-H and C-O vibration, C-O-C bridges, C-O stretching, C-H rocking and C-O-C β-glycosidic bonding   . The emergence of the 540 - 412 cm−1 Zn-O bond broadband confirms that there are cellulose-ZnO structures, as suggested in Figure 1(e), and the characteristic cellulose bands are still present      .
Although negative zeta potential value was observed for NCC_H, cellulose-O-SO3 binding was not detected in the IR spectrum, therefore, the degree of sulfonation is below the FTIR detection limit .
3.5. Thermogravimetric Analysis
The degradation of MCC and NCC_C presents a single degradation pattern, as
Figure 3. FTIR spectrum of MCC, NCC_C, NCC_H and NCC_ZnO samples.
seen in Figure 4, corresponding to the steps of depolymerization, dehydration, decomposition of glycosidic units and ash formation, and such reactions occur concurrently  . However, this profile is altered with the presence of sulphate groups or different surface functionalization, as described by  . Obtaining cellulose nanocrystals by sulfuric acid hydrolysis, NCC_H, reduced the temperature at which 10% of the mass was degraded (T10%, shown in Table 2) by up to 30% when compared to MCC. As observed in this paper, it is described in the literature that the presence of sulphate clusters on the cellulose surface reduces thermal stability, which agrees with Kargarzadesh group , which associated low thermal stability with longer hydrolysis time, making thus evident that there was a change in the degradation profile after the neutralization process, which possibly removed sulfate residues . The NCC_C degraded in a single event, because as reported in the technical sheet, the neutralization of the groups with sodium hydroxide was performed.
In the DTG curve, it is possible to see the changes in the cellulose degradation behavior. The degradation of the NCC_H sample occurs over a wider temperature range, and events that occur at lower temperatures are related to the most accessible sulfated region, which, according to Kargarzadesh , corresponds
Figure 4. Mass loss curve and temperature mass loss derivative for MCC, NCC_C, NCC_H and NCC_ZnO cellulose samples.
Table 2. Data extracted from TGA curves, temperature relative to 10% mass loss (T10%), peak temperature obtained from DTG, and % residue of samples after analysis.
to the amorphous region; also, it was attributed to the increase in sulfated chain terminals, which begin to degrade faster   . The increase in residue percentage after hydrolysis with sulfuric acid compared to MCC is related to the flame-retardant property of sulfate groups   .
Zinc oxide surface modification (NCC_ZnO) improved cellulose thermal stability compared to NCC_H, evidenced by a Tonset increase of 54C, as shown in Table 2. Increase on thermal stability of sulfuric acid hydrolysis produced cellulose was reported by Wang , who neutralized the acidic suspension with sodium hydroxide. The observed Tonset increase of NCC_ZnO is justified by the removal of these acid groups and the increase of zinc nano oxide on the surface. Comparing the DTG curves of NCC_H and NCC_ZnO samples, the peak width was decreased after modification, like MCC. Giraldi  attributed this narrowing of the DTG peak to the conversion of the zinc metal precursor to zinc oxide, where the acetate groups degraded under 200˚C. All samples showed a mass loss of up to 5% at 100˚C and this loss can be attributed to humidity .
The NCC_ZnO sample exhibited 5% mass loss at 100˚C, followed by a single degradation event with 30% mass loss. The residue percentage corresponds to the sulfur residue, derived from the NCC_H sulfate ester groups, and to the modification with zinc oxide. Azizi’s group  tested different proportions of precursor and cellulose, they found that the higher the amount of precursor used, the higher the percentage of waste generated, which agrees with . Arputharaj  observed that the ZnO produced on the surface of the cotton cloth was able to increase its flame-retardant characteristics. Zinc oxide does not lose mass during the temperature range studied.
3.6. Nuclear Magnetic Resonance Time Domain (TD-NMR)
The nuclear magnetic relaxation technique allows the morphological investigation of nanostructured systems such as cellulose, through the relaxation times of protons present in each environment or domain, where the polymer segments assume a specific conformation/arrangement. From the measurements obtained by MSE-FID (Figure 5) and the adjustment of the points based on Equation (3), it is observed that all samples have two well defined domains, one of them rigid (Fr) associated with the crystalline fraction of matrix, and another of higher molecular mobility (Fm) related to the amorphous fraction of cellulose. The crystallinity degree values obtained from Fr are like those found in CrI, calculated by XRD analysis. In Table 3, MCC and NCC_C have similar values of relaxation times, which may be associated with uniformity in particle size, as seen in Figure 1 and discussed above. For NCC_H these relaxation times are shorter since not only cellulose nanocrystals are present but also cellulose nanofibrils which have a highly oriented morphology. For the NCC_ZnO sample, it is observed that the presence of oxide in the cellulose structure reduces its crystalline fraction to 86% .
The mean value of the second dipole coupling moment
Figure 5. MSE-refocused FID curves with nonlinear fit based on Equation (3).
the distance between cellulose chain segments in the crystalline fraction, being proportional to the coupling force between the hydrogen nuclei in the crystals. The value of 9262 rad·ms−2, shown in Table 3, found for NCC_H indicates that while having a slightly lower crystallinity than the MCC and NCC_C samples, the crystals in this matrix have a larger organization, with better packed segments and high coupling  . However, for the NCC_ZnO sample, the value of
Nanocrystals and cellulose nanofibrils were produced in a single step combining acid hydrolysis mechanical stirring with ultra-turrax, which was responsible for the improvement of MCC dispersibility in acid medium, the high shear promoted the mass transfer and helped the disintegrations and isolation to NCC
Table 3. Transversal relaxation time for rigid and mobile fraction.
and NFC. Next, the nanoparticles were superficially modified with nano oxide by sol-gel methodology under low temperature, as found by HIM and AFM, FTIR, XRD, TGA and NMR analyses. The NMR technique clarified the changes in the constituent fractions of cellulose, rigid and mobile fraction, validating the observations made through the XRD technique. Furthermore, the crystalline fraction of NCC_ZnO was detected by this technique. The properties presented by nanocellulose and zinc oxide were combined into a single particle, and the improvement in thermal stability enables the preparation of nanocomposite via melting mixture. Therefore, this nanoparticle is expected to play the reinforcing role and barrier property, characteristic of nanocellulose, and antimicrobial and anti-UV activity, property of ZnO, when added to a polymer. According to all properties studied about this particle we can say that it could be used in packing.
This work was financially supported by CAPES code 001 and CNPQ. The authors thank for the Biological Physical Laboratory, located in Carlos Chagas Institute of Biophysics, special for the Professors Dr. Gilberto Weissmüller and Dr. Paulo Bish. We also thank for Inorganic Chemistry Department, Institute of Chemistry, mainly MSc. Glaucia Martins for all the support and contributions. Both Institutes are in the Federal University of Rio de Janeiro.
 Kargarzadeh, H., Mariano, M., Huang, J., Lin, N., Ahmad, I., Dufresne, A. and Thomas, S. (2017) Recent Developments on Nanocellulose Reinforced Polymer Nanocomposites: A Review. Polymer, 132, 368-393.
 Wang, N., Ding, E. and Cheng, R. (2007) Thermal Degradation Behaviors of Spherical Cellulose Nanocrystals with Sulfate Groups. Polymer, 48, 3486-3493.
 Bondeson, D., Mathew, A. and Oksman, K. (2006) Optimization of the Isolation of Nanocrystals from Microcrystalline Cellulose by Acid Hydrolysis. Cellulose, 13, 171.
 Li, B., Xu, W., Kronlund, D., Maattanen, A., Liu, J., Smatt, J.-H., Peltonen, J., Willfor, S., Mu, X. and Xu, C. (2015) Cellulose Nanocrystals Prepared via Formic Acid Hydrolysis Followed by TEMPO-Mediated Oxidation. Carbohydrate Polymers, 133, 605-612.
 Melo, A.R.A., da Silva, E.O., Menezes, L.R. and Tavares, M.I.B. (2018) The Effect of Modified Cellulose Particles on Morphology and Properties Ethylene Vinyl Acetate Copolymer. Polymer Testing, 68, 333-339.
 Camarero Espinosa, S., Kuhnt, T., Foster, E.J. and Weder, C. (2013) Isolation of Thermally Stable Cellulose Nanocrystals by Phosphoric Acid Hydrolysis. Biomacromolecules, 14, 1223-1230.
 Chen, L., Zhu, J., Baez, C., Kitin, P. and Elder, T. (2016) Highly Thermal-Stable and Functional Cellulose Nanocrystals and Nanofibrils Produced Using Fully Recyclable Organic Acids. Green Chemistry, 18, 3835-3843.
 Xu, W.-J., Zhai, J.-W., Cui, Q., Liu, J.-Z., Luo, M., Fu, Y.-J. and Zu, Y.-G. (2016) Ultra-Turrax Based Ultrasound-Assisted Extraction of Five Organic Acids from Honeysuckle (Lonicera japonica Thunb.) and Optimization of Extraction Process. Separation and Purification Technology, 166, 73-82.
 Liu, Y., Wang, H., Yu, G., Yu, Q., Li, B. and Mu, X. (2014) A Novel Approach for the Preparation of Nanocrystalline Cellulose by Using Phosphotungstic Acid. Carbohydrate Polymers, 110, 415-422.
 Lu, Q., Cai, Z., Lin, F., Tang, L., Wang, S. and Huang, B. (2016) Extraction of Cellulose Nanocrystals with a High Yield of 88% by Simultaneous Mechanochemical Activation and Phosphotungstic Acid Hydrolysis. ACS Sustainable Chemistry & Engineering, 4, 2165-2172.
 Boujemaoui, A., Mongkhontreerat, S., Malmstrom, E. and Carlmark, A. (2015) Preparation and Characterization of Functionalized Cellulose Nanocrystals. Carbohydrate Polymers, 115, 457-464.
 Jia, C., Chen, L., Shao, Z., Agarwal, U.P., Hu, L. and Zhu, J.Y. (2017) Using a Fully Recyclable Dicarboxylic Acid for Producing Dispersible and Thermally Stable Cellulose Nanomaterials from Different Cellulosic Sources. Cellulose, 24, 2483-2498.
 Li, D., Henschen, J. and Ek, M. (2017) Esterification and Hydrolysis of Cellulose Using Oxalic Acid Dihydrate in a Solvent-Free Reaction Suitable for Preparation of Surface-Functionalised Cellulose Nanocrystals with High Yield. Green Chemistry, 19, 5564-5567.
 Tang, L., Huang, B., Lu, Q., Wang, S., Ou, W., Lin, W. and Chen, X. (2013) Ultrasonication-Assisted Manufacture of Cellulose Nanocrystals Esterified with Acetic Acid. Bioresource Technology, 127, 100-105.
 Kargarzadeh, H., Ahmad, I., Abdullah, I., Dufresne, A., Zainudin, S.Y. and Sheltami, R.M. (2012) Effects of Hydrolysis Conditions on the Morphology, Crystallinity and Thermal Stability of Cellulose Nanocrystals Extracted from Kenaf Bast Fibers. Cellulose, 19, 855-866.
 Wang, Q.Q., Zhao, X. and Zhu, J.Y. (2014) Kinetics of Strong Acid Hydrolysis of a Bleached Kraft Pulp for Producing Cellulose Nanocrystals (CNCs). Industrial & Engineering Chemistry Research, 53, 11007-11014.
 Hamad, W.Y. and Hu, T.Q. (2010) Structure-Process-Yield Interrelations in Nanocrystalline Cellulose Extraction. The Canadian Journal of Chemical Engineering, 88, 392-402.
 Wang, Q.Q., Zhu, J.Y., Reiner, R.S., Verrill, S.P., Baxa, U. and McNeil, S.E. (2012) Cellulose Loss in Cellulose Nanocrystal (CNC) Production: Recovery and Characterization of Cellulosic Solid Residues (CSR) and CNCApproaching Zero. Cellulose, 19, 2033-2047.
 Aguayo, M., Fernández Pérez, A., Reyes, G., Oviedo, C., Gacitúa, W., Gonzalez, R. and Uyarte, O. (2018) Isolation and Characterization of Cellulose Nanocrystals from Rejected Fibers Originated in the Kraft Pulping Process. Polymers, 10, 1145.
 Chen, L., Wang, Q., Hirth, K., Baez, C., Agarwal, U.P. and Zhu, J. (2015) Tailoring the Yield and Characteristics of Wood Cellulose Nanocrystals (CNC) Using Concentrated Acid Hydrolysis. Cellulose, 22, 1753-1762.
 Grishkewich, N., Mohammed, N., Tang, J. and Tam, K.C. (2017) Recent Advances in the Application of Cellulose Nanocrystals. Current Opinion in Colloid & Interface Science, 29, 32-45.
 Arputharaj, A., Nadanathangam, V. and Shukla, S.R. (2017) A Simple and Efficient Protocol to Develop Durable Multifunctional Property to Cellulosic Materials Using in Situ Generated Nano-ZnO. Cellulose, 24, 3399-3410.
 Jia, M., Zhang, X., Weng, J., Zhang, J. and Zhang, M. (2019) Protective Coating of Paper Works: ZnO/Cellulose Nanocrystal Composites and Analytical Characterization. Journal of Cultural Heritage, 38, 64-74.
 Khosravian, S., Montazer, M., Malek, R.M. and Harifi, T. (2015) In Situ Synthesis of Nano ZnO on Starch Sized Cotton Introducing Nano Photo Active Fabric Optimized with Response Surface Methodology. Carbohydrate Polymers, 132, 126-133.
 Azizi, S., Ahmad, M., Mahdavi, M. and Abdolmohammadi, S. (2013) Preparation, Characterization and Antimicrobial Activities of ZnO Nanoparticles/Cellulose Nanocrystal Nanocomposites. BioResources, 8, 1841-1851.
 Azizi, S., Ahmad, M.B., Ibrahim, N.A., Hussein, M.Z. and Namvar, F. (2014) Cellulose Nanocrystals/ZnO as a Bifunctional Reinforcing Nanocomposite for Poly (Vinyl Alcohol)/Chitosan Blend Films: Fabrication, Characterization and Properties. International Journal of Molecular Sciences, 15, 11040-11053.
 Anvzlovar, A., Orel, Z.C., Kogej, K. and Zigon, M. (2012) Polyol-Mediated Synthesis of Zinc Oxide Nanorods and Nanocomposites with Poly (Methyl Methacrylate). Journal of Nanomaterials, 2012, Article ID: 760872.
 Jia, W., Dang, S., Liu, H., Zhang, Z., Yu, C., Liu, X. and Xu, B. (2012) Evidence of the Formation Mechanism of ZnO in Aqueous Solution. Materials Letters, 82, 99-101.
 Li, X., Zhang, X., Li, L., Huang, L., Zhang, W. and Ye, J. (2016) Preparation of Nano-ZnO/Regenerated Cellulose Composite Particles via Co-Gelation and Low-Temperature Hydrothermal Synthesis. Materials Letters, 175, 122-125.
 Eyholzer, C., Bordeanu, N., Lopez-Suevos, F., Rentsch, D., Zimmermann, T. and Oksman, K. (2010) Preparation and Characterization of Water-Redispersible Nanofibrillated Cellulose in Powder Form. Cellulose, 17, 19-30.
 Jin, E., Guo, J., Yang, F., Zhu, Y., Song, J., Jin, Y. and Rojas, O. J. (2016) On the Polymorphic and Morphological Changes of Cellulose Nanocrystals (CNC-I) Upon Mercerization and Conversion to CNC-II. Carbohydrate Polymers, 143, 327-335.
 Khoshkava, V. and Kamal, M. (2014) Effect of Drying Conditions on Cellulose Nanocrystal (CNC) Agglomerate Porosity and Dispersibility in Polymer Nanocomposites. Powder Technology, 261, 288-298.
 Sèbe, G., Ham-Pichavant, F., Ibarboure, E., Koffi, A.L.C. and Tingaut, P. (2012) Supramolecular Structure Characterization of Cellulose II Nanowhiskers Produced by Acid Hydrolysis of Cellulose I Substrates. BioMacromolecules, 13, 570-578.
 Jayachandiran, J., Vajravijayan, S., Nandhagopal, N., Gunasekaran, K. and Nedumaran, D. (2019) Fabrication and Characterization of ZnO Incorporated Cellulose Microfiber Film: Structural, Morphological and Functional Investigations. Journal of Materials Science: Materials in Electronics, 30, 6037-6049.
 Singh, A., Viswanath, V. and Janu, V. (2009) Synthesis, Effect of Capping Agents, Structural, Optical and Photoluminescence Properties of ZnO Nanoparticles. Journal of Luminescence, 129, 874-878.
 Nam, S., French, A.D., Condon, B.D. and Concha, M. (2016) Segal Crystallinity Index Revisited by the Simulation of X-Ray Diffraction Patterns of Cotton Cellulose Iβ and Cellulose II. Carbohydrate Polymers, 135, 1-9.
 Ramírez, J.A.á., Suriano, C.J., Cerrutti, P. and Foresti, M.L. (2014) Surface Esterification of Cellulose Nanofibers by a Simple Organocatalytic Methodology. Carbohydrate Polymers, 114, 416-423.
 Giraldi, T.R., Santos, G.V., Mendonca, V.R., Ribeiro, C. and Weber, I.T. (2011) Annealing Effects on the Photocatalytic Activity of ZnO Nanoparticles. Journal of Nanoscience and Nanotechnology, 11, 3635-3640.
 Grunin, L.Y., Grunin, Y.B., Nikolskaya, E.A., Sheveleva, N.N. and Nikolaev, I.A. (2017) An NMR Relaxation and Spin Diffusion Study of Cellulose Structure during Water Adsorption. Biophysics (Russian Federation), 62, 198-206.
 Mackay, A.L., Tepfer, M., Taylor, I.E.P. and Volke, F. (1985) Proton Nuclear Magnetic Resonance Moment and Relaxation Study of Cellulose Morphology. Macromolecules, 18, 1124-1129.
 Papon, A., Saalwachter, K., Schaler, K., Guy, L., Lequeux, F. and Montes, H. (2011) Low-Field NMR Investigations of Nanocomposites: Polymer Dynamics and Network Effects. Macromolecules, 44, 913-922.