Quantum mechanical modeling of the structures, energetics and spectral properties of Iα and Iβ cellulose

Periodic planewave and molecular cluster density functional theory (DFT) calculations were performed on Iα and Iβ cellulose in four different conformations each. The results are consistent with the previous interpretation of experimental X-ray and neutron diffraction data that both Iα and Iβ cellulose are dominantly found in the tg conformation of the hydroxymethyl group with a H-bonding conformation termed “Network A”. Structural and energetic results of the periodic DFT calculations with dispersion corrections (DFT-D2) are consistent with observation suggesting that this methodology is accurate to within a few percent for modeling cellulose. The structural and energetic results were confirmed by comparison of calculated vibrational frequencies against observed infrared and Raman frequencies of Iα and Iβ cellulose. Structures extracted from the periodic DFT-D2 energy minimizations were used to calculate the ¹³C nuclear magnetic resonance chemical shifts (δ¹³C), and the tg/Network A conformations of both Iα and Iβ cellulose produced excellent correlations with observed δ¹³C values.     

Influence of drying method on the material properties of nanocellulose I: thermostability and crystallinity

The effect of drying method on selected material properties of nanocellulose was investigated. Samples of nanofibrillated cellulose (NFC) and cellulose nanocrystals (CNC) were each subjected to four separate drying methods: air-drying, freeze-drying, spray-drying, and supercritical-drying. The thermal stability and crystallinity of the dried nanocellulose were evaluated using thermogravimetric analysis (TGA) and X-ray diffraction. Supercritical-drying produced NFCs with the least thermal stability and the lowest crystallinity index. Air-drying or spray-drying produced NFCs which were more thermally stable compared with freeze-dried NFCs. The CNCs dried by the three methods (air-drying, freeze-drying, and spray-drying) have similar onset temperature of thermal degradation. The different drying methods resulted in various char weight percentages at 600 °C for the dried NFCs or CNCs from TGA measurements. The dried NFCs are pure cellulose I while the dried CNCs consist of cellulose I and II. The calculated crystallinity indices differ with each drying method. The cellulose II content in CNCs changes as a function of drying method. For the application of nanocellulose in non polar thermoplastics, spray-dried products are recommended according to their higher thermal stability and higher crystallinity index.     

Fluorescence quenching of three molecular weight fractions of a soil fulvic acid by UO(2)(II)

A soil fulvic acid isolated from a Korean forest was divided into three different molecular weight fractions (F1: less than 220 Da; F2: 220-1000 Da; and F3: 1000-4000 Da) by gel filtration chromatography and the fractions were studied by synchronous fluorescence (SyF) spectroscopy. Analysis of the SyF spectra for the fulvic acid fractions showed that the fractions with molecules of larger sizes have a higher content of condensed aromatic compounds. The information about their interaction with UO(2)(II) ions in an aqueous solution (100 mg l(-1) of fulvic acid, in 0.1 M NaClO(4) at pH 3.5) was obtained from the measurement of SyF spectra at increasing concentrations of metal ions. Self-modeling mixture analysis of the quenching spectra gives two distinct peak components having a maximum peak position of 386 (type I) and 498 nm (type II) for all the size-fractionated fulvic acids. From the analysis of the quenching profiles of the peaks, using a non-linear method, the concentration of binding sites (C(L)), and the corresponding stability constants (logK) were calculated. The stability constants of the UO(2)(II)-fulvate complexes ranged from 4.10 to 5.33, and increased with higher molecular weight fractions, which indicates a stronger affinity for UO(2)(II) in the fraction with molecules of larger size.     

Theoretical investigations on the structure and physical properties of cellulose

The structure of cellulose is investigated using a recently derived force field. Published experimental data is taken only as a starting point for purely theoretical investigations. The reliability of the method is validated by calculating physical properties of the obtained geometries. In the course of the investigations, the geometries of cellulose Iα, cellulose Iβ and cellulose II are derived. For these geometries the Young's-modulus is calculated. The structure of cellulose in aqueous solution is investigated, using cellohexaose (a hexamer of β-D -glucose) as a fragment of a cellulose chain. Here, the diffusion coefficient is calculated.     

Mercerized cellulose biocomposites: a study of influence of mercerization on cellulose supramolecular structure, water retention value and tensile properties

In this study the effect of the mercerization degree on the water retention value (WRV) and tensile properties of compression molded sulphite dissolving pulp was evaluated. The pulp was treated with 9, 10, or 11 % aqueous NaOH solution for 1 h before compression molding. To study the time dependence of mercerization the pulp was treated with 12 wt% aqueous NaOH for 1, 6 or 48 h. The cellulose I and II contents of the biocomposites were determined by solid state cross polarization/magic angle spinning carbon 13 nuclear magnetic resonance (CP/MAS ¹³C NMR) spectroscopy. By spectral fitting of the C6 and C1 region the cellulose I and II content, respectively, could be determined. Mercerization decreased the total crystallinity (sum of cellulose I and cellulose II content) and it was not possible to convert all cellulose I to cellulose II in the NaOH range investigated. Neither increased the conversion significantly with 12 wt% NaOH at longer treatment times. The slowdown of the cellulose I conversion was suggested as being the result from the formation of cellulose II as a consequence of coalescence of anti-parallel surfaces of neighboring fibrils (Blackwell et al. in Tappi 61:71–72, 1978; Revol and Goring in J Appl Polym Sci 26:1275–1282, 1981; Okano and Sarko in J Appl Polym Sci 30:325–332, 1985). Compression molding of the partially mercerized dissolving pulps yielded biocomposites with tensile properties that could be correlated to the decrease in cellulose I content in the pulps. Mercerization introduces cellulose II and disordered cellulose and lowered the total crystallinity reflected as higher water sensitivity (higher WRV values) and poorer stiffness of the mercerized biocomposites.     

CP/MAS ~(13)C NMR技术对木浆纤维微观结构的研究

利用交叉极化结合魔角旋转技术~(13)C核磁共振法(CP/MAS ~(13)C NMR)对桉木浆纤维的微观结构进行研究,为进一步研究木质纤维素材料开发过程中反应障碍特征奠定基础.通过对NMR光谱C1区(δ 102～108)进行洛仑兹拟合,得到桉木浆纤维中纤维素Iα的相对含量为26.92%,纤维素Iβ的相对含量为52.04%,主要以纤维素Iβ晶体形式为主.通过计算纤维素C4结晶区(δ 86～92)和非结晶区(δ 80～86)的相对含量得到桉木浆的纤维素结晶度为47%.通过洛仑兹和高斯函数的混合模型对NMR光谱C4区(δ 80～92)进行拟合得到基原纤尺寸和微原纤横向尺寸分别为4.0与17.9 nm,并通过计算不同形态的结晶纤维素的相对含量得到纤维素结晶度为51%,证实了在微原纤内部次晶纤维素的存在.     

Theoretical investigation of azo dyes adsorbed on cellulose fibers: 2. Spectroscopic study

Electronic absorption spectra and the frontier orbitals of 1-arylazo-2-naphtol dyes are computed and analyzed in four models, namely in the gas phase (model I), in a solvent (model I + CPCM), adsorbed on the cellulose surface (model II), and model II in the presence of solvent (model II + CPCM) via time-dependent density functional theory (TD-DFT) and conductor-like polarizable continuum model (CPCM) at the B3LYP/6-31G** level of theory. A bathochromic shift is observed for the λ_max peak due to both short-range and long-range interactions of the non-ionic dyes with cellulose, while the ionic dyes exhibit hypsochromic shift in their λ_max peak. The results predict that the studied dyes should be nearly yellow after being adsorbed on cellulose with excellent color strength. Furthermore, the ionic dyes are suitable for the dyeing of cellulose fibers. The nuclear magnetic resonance (NMR) chemical shieldings calculated for the azo dyes in the gas phase and adsorbed states and for their tautomeric equilibrium mixtures show that the NMR technique can be used successfully to follow the dyeing process.     

Two-dimensional Rietveld analysis of celluloses from higher plants

The Rietveld method is a versatile tool to parameterize the fine structure of crystallites analyzed by diffraction. The method relies on a crystallographic model representing what is known a priori, and free coefficients determined from fits to experimental data. This article provides an introduction to Rietveld analysis of celluloses from higher plants that are adequately described by the cellulose Iβ crystal structure. Possibilities of Rietveld analysis have been recently enhanced by a tailored crystallographic model and computational algorithm, named Cellulose Rietveld Analysis for Fine Structure (CRAFS). From each two-dimensional diffraction pattern, CRAFS automated analysis outputs unit cell parameters, crystallite sizes, peak profile functions, integrated crystalline intensity (proportional to cellulose degree of crystallinity), and crystallite orientation distribution function. Two of the major hurdles for analysis of plant cellulose—overlapping of diffraction peaks and preferred crystallite orientation—are consistently treated by the two-dimensional Rietveld analysis. Hence, the method is a unique tool to explore cellulose fine structural variability, with differences arising from specimen conditioning, processing, and biological origins.     

X-ray diffraction from faulted cellulose I constructed with mixed Iα–Iβ stacking

Cellulose from higher plants is usually thought to be a composite of the Iα and Iβ allomorphs, with predominance of the latter. Instead of the pure allomorphs, this article proposes that Iα and Iβ stacking patterns coexist within each crystallite, forming a type of crystallographic defect known as stacking fault. Models of faulted crystallites are constructed with mixed Iα–Iβ stacking and their X-ray diffraction intensities are calculated using the Diffracted Intensities From Faulted Xtals (DIFFaX) computer program. Simulated powder diffractograms from faulted crystallites compare favorably with experimental data, modifying diffractogram regions that have been misfit by models based on the Iβ crystal structure. Calculations also reveal that stacking faults generate a signature in the (hkl) dependence of diffraction line broadening, guiding further experimental verification and eventual quantification of stacking faults. Our findings bring an alternative view of native cellulose polymorphism and suggest that the proposed stacking faults are ubiquitous crystallographic defects in cellulose from higher plants.     

Partial crystalline transformation of solvated cellulose IIII crystals, reproduced by theoretical calculations

In hot-water molecular dynamics simulation at 370 K, four cellulose III_I crystal models, with different lattice planes and dimensions, exhibited partial crystalline transformations of (1 ?1 0) chain sheets, in which hydroxymethyl groups were irreversibly rotated from gt into tg conformations, accompanied by hydrogen-bond exchange from the original O3–O6 to cellulose-I-like O2–O6 bonds. The final hydrogen-bond exchange ratio was about 95 % for some of the crystal models after 50 ns simulation. The corrugated (1 ?1 0) chain sheet was converted to a cellulose-I-like flat chain sheet with a slightly right-handed twist. The 3D structures of the three types of isolated chain sheet models were optimized using density functional theory calculations to compare their stabilities without crystal packing forces. The cellulose Iβ (1 0 0) models were more stable than the cellulose III_I (1 ?1 0) models. The optimized structure of cellulose III_I (1 0 0) models deviated largely from the initial sheet form. It was proposed to the crystalline transformation from cellulose III_I to Iβ that conversion of the chain sheet structure first take place, followed by sliding of the chain sheet along the fiber axis.     

Simulating infrared spectra and hydrogen bonding in cellulose Iβ at elevated temperatures

We have modeled the transformation of cellulose Iβ to a high temperature (550 K) structure, which is considered to be the first step in cellulose pyrolysis. We have performed molecular dynamics simulations at constant pressure using the GROMOS 45a4 united atom forcefield. To test the forcefield, we computed the density, thermal expansion coefficient, total dipole moment, and dielectric constant of cellulose Iβ, finding broad agreement with experimental results. We computed infrared (IR) spectra of cellulose Iβ over the range 300-550 K as a probe of hydrogen bonding. Computed IR spectra were found to agree semi-quantitatively with experiment, especially in the O-H stretching region. We assigned O-H stretches using a novel synthesis of normal mode analysis and power spectrum methods. Simulated IR spectra at elevated temperatures suggest a structural transformation above 450 K, a result in agreement with experimental IR results. The low-temperature (300-400 K) structure of cellulose Iβ is dominated by intrachain hydrogen bonds, whereas in the high-temperature structure (450-550 K), many of these transform to longer, weaker interchain hydrogen bonds. A three-dimensional hydrogen bonding network emerges at high temperatures due to formation of new interchain hydrogen bonds, which may explain the stability of the cellulose structure at such high temperatures.     

The swelling and dissolution of cellulose crystallites in subcritical and supercritical water

The swelling and dissolution phenomena of microcrystalline cellulose (MCC) were investigated in subcritical and supercritical water. Commercial MCC was treated in water at temperatures of 250–380 °C and a pressure of 250 bar for 0.25–0.75 s. As reaction products, undissolved but depolymerised cellulose residue, short-chain cellulose precipitate, water-soluble cello-oligosaccharides and monosaccharides, as well as their degradation products, were detected. The highest yield of the cellulose II precipitate was obtained after a reaction time of 0.25 s at 360 °C. Our hypothesis was that if the crystallites were swollen, the depolymerization pattern would be that of homogeneous reaction and the cellulose Iβ to cellulose II transformation would be observed. The changes in the structure of the undissolved cellulose residue were characterised by size exclusion chromatography, wide-angle X-ray scattering and ¹³C solid-state NMR techniques. In many cases, the cellulose residue samples contained cellulose II; however, due to experimental limitations, it remains unclear whether it was formed through the swelling of crystallites or the partial readsorption of the dissolved cellulose fraction. The molar mass distributions of untreated MCC and after low intensity treatments showed a bimodal shape. After high intensity treatments the high molar mass chains disappeared which indicated a complete swelling or dissolution of the crystallites.     

Simulation of X-ray diffractograms relevant to the purported polymorphs cellulose IV<Subscript>I</Subscript> and IV<Subscript>II</Subscript>

Diffractograms were simulated for model nanocrystals of cellulose I_β, using numerical summation of radiation scattered from all carbon and oxygen atoms in the nanocrystal. Diffractogram peaks were sometimes displaced by a few degrees from positions calculated by the Bragg equation, as predicted in a published study based on a different mathematical approach. Simulated diffractograms showed 2 or 3 peaks, depending on the cross-sectional size and shape of the model nanocrystal. Some of the 2-peak diffractograms resembled published results for the purported polymorph cellulose IV_I, or for cell-wall cellulose, supporting suggestions that cellulose IV_I is simply cellulose I fragmented into nanocrystals with relatively small cross-sectional dimensions. A published diffractogram for cellulose IV_II could not be simulated with acceptable precision, suggesting that this polymorph might have a crystal structure distinctly different from that of cellulose I_β.     

Hydrolysis of steam-pretreated lignocellulose

The mechanism of hydrolysis of cellulose is important for improving the enzymatic conversion in bioprocesses based on lignocellulose. Adsorption and hydrolysis experiments were performed with cellobiohydrolase I (CBH I) and endoglucanase II (EG II) from Trichoderma reesei on a realistic lignocellulose substrates: steam-pretreated willow. The enzymes were studied both alone and in equimolar mixtures. Adsorption isotherms were determined at 4 and 40 degrees C during 90-min reaction times. Both CBH I and EG II adsorbed stronger at 40 than at 4 degrees C. The time course of adsorption and hydrolysis, 3 min to 48 h, was studied at 40 degrees C. About 90% of the cellulases were adsorbed within 2 h. The hydrolysis rate was high in the beginning but decreased during the time course. Based on adsorption data, the hydrolysis and synergism were analyzed as function of adsorbed enzyme. CBH I showed a linear correlation between hydrolysis and adsorbed enzyme, whereas for EG II the corresponding curve leveled off at both 4 and 40 degrees C. At low conversion, below 1%, EG II produced as much soluble sugars as CBH I. At higher conversion, CBH I was more efficient than EG II. The synergism as function of adsorbed enzyme increased with bound enzyme before reaching a stable value of about 2. The effect of varying the ratio of CBH I:EG II was studied at fixed total enzyme loading and by changing the ratio between the enzymes. Only a small addition (5%) of EG II to a CBH I solution was shown to be sufficient for nearly maximal synergism. The ratio between EG II and CBH I was not critical. The ratio 40% EG II:60% CBH I showed similar conversion to 5% EG II:95% CBH I. Modifications of the conventional endo-exo synergism model are proposed.     

The Structure of Cellulose by Conformational Analysis. Part 5. The Cellulose Ii Water Sorption IsothermXS

Exact values of the sorption energies of single molecules of water on all available sorption sites of crystalline cellulose II have been obtained by conformational analysis. The sorption energies are equated to the total energy (E_tot ) of interaction between the water molecule and all the suitable atomic groups of the cellulose. E_tot is composed of van der Waals, H-bond, and electrostatic energies. The interferences of water molecules on vicinal sorption sites were obtained. Sites in which such interference can occur were identified for crystalline cellulose II. Sorption energy in crystalline cellulose II appears to depend only on the interaction of water with surface sorption sites of the crystal. There appears to be favorable sorption on 1) sites exerting high attractive forces, and 2) sites which are exposed and protrude from the crystal surface. Sites recessed from the crystal surface are generally repulsive due to strong interactions with neighboring groups. All the sorption energies of the “monolayer” were calculated. Very strong sorption sites cannot always form a second layer because of strong steric hindrance from vicinal groups. Sorption capacities of crystalline cellulose II were calculated, and the isotherm of the schematic five chain crystallite used was constructed by theoretical means. The results obtained were briefly compared with those for cellulose I crystallites and amorphous cellulose. The inflection points of the isotherm and the variability of Dent's k ₁ constant for the water monolayer with relative humidity for the cellulose I and II isotherms were also calculated by theoretical means.     

Characterization of low-frequency modes in aqueous peptides using far-infrared spectroscopy and molecular dynamics simulation

Far-infrared spectroscopy was used to study the dynamics of three aqueous peptides having varied helicity. Experimental data were compared to the molecular dynamics simulated far-infrared absorbance spectrum derived from the dipole time correlation function. Vibrational density of state (VDOS) simulation was then used to analyze the contribution of different structural elements to the bands. Frozen aqueous peptide samples were studied in the frequency range between 325 and 540 cm(-1) where the ice absorbance is low. Three resonances were identified; band I centered at approximately 333 cm(-1), band II centered at approximately 380 cm(-1), and band III comprising two constituent bands at approximately 519 and 528 cm(-1). The peak height and frequency of the maximum absorbance of bands I and II varied depending on the helicity of the peptide. VDOS of the far-infrared absorbance spectrum confirmed that bands I and II were associated with the peptide backbone and that band III had both potential backbone and side chain components.     

Enzymatic hydrolysis of cellulose I is greatly accelerated via its conversion to the cellulose II hydrate form

Cellulose II hydrate was prepared from microcrystalline cellulose (cellulose I) via its mercerization with 5 N NaOH solution over 1 h at room temperature followed by washing with water. The structure of cellulose II hydrate changed to that of cellulose II after drying. Compared with cellulose II, cellulose II hydrate exhibited a slightly (8.5%) expanded structure only along the direction. The hydrophobic stacking sheets of the cellulose II were conserved in the cellulose II hydrate, and water molecules could be incorporated in the inflated two-chain unit cell of cellulose II hydrate. Enzymatic hydrolysis of cellulose I, cellulose II hydrate, and cellulose II was carried out at 37 °C using solutions comprising a mixture of cellulase and β-glucosidase. The hydrolysis of cellulose II hydrate proceeded much faster than the hydrolysis of the other two substrates, while the saccharification ratio of cellulose II was only slightly higher than that of cellulose I. The alkaline mercerization treatment was also applied to sugarcane bagasse. After its direct mercerization, the cellulose in bagasse was converted from cellulose I to cellulose II hydrate, and then to cellulose II after drying. Similar to in the case of microcrystalline cellulose, the rate of the enzymatic hydrolysis of the mercerized bagasse without drying (cellulose II hydrate) was much faster than the enzymatic hydrolysis of the other two substrates. Thus, the wet forms of cellulose and cellulosic biomass after mercerization, and after hydrolysis with cellulolytic enzymes, afforded superior products with extremely high degradability.     

Preparation and characterization of cellulose nanofibers from partly mercerized cotton by mixed acid hydrolysis

Cellulose nanofibers with a diameter of 70 nm and lengths of approximately 400 nm were fabricated from partly mercerized cotton fibers by acid hydrolysis. Morphological evolution of the hydrolyzed cotton fibers was investigated by powder X-ray diffraction, Fourier transform infrared analysis and field emission scanning electron microscopy. The XRD results show that the cellulose I was partially transformed into cellulose II by treatment with 15 % NaOH at 150° for 3 h. The crystallinity of this partially mercerized sample was lower than the samples that were converted completely to cellulose II by higher concentrations of NaOH. The intensities of all of the diffraction peaks were noticeably increased with increased hydrolysis time. Fourier transform infrared results revealed that the chemical composition of the remaining nanofibers of cellulose I and II had no observable change after acidic hydrolysis, and there was no difference between the hydrolysis rates for cellulose I or II. The formation of cellulose nanofibers involves three stages: net-like microfibril formation, then short microfibrils and finally nanofibers.     

Resolution method for overlapping peaks based on the fractional-order differential

Equations between the differential order and the maximum of the fractional-order differential for the specified peak signals are developed based on the variation of the maximum of the specified peak signals at different orders. Also, equations between the differential order and the zero-crossing of the fractional-order differential of the specified peak signals are proposed according to the variation of the zero-crossing of the specified peak signals at different orders. Characteristic paramters of the Gaus- sian peak, Lorentzian peak, and Tsallis peak can be estimated using estimator I and estimator II which are obtained by the equations above. As a result, a new method is presented to resolve the overlapped peaks signal. Firstly, a fractional-order differential of the specified peak signals is obtained with the fractional-order differentiation filter. Then, characteristic paramters of the specified peak signals can be extracted using estimator I and estimator II. Finally, the Tsallis peak is used as a model to assign the overlapping peak signals correctly. Experimental results show that the proposed method is efficient and effective for the simulated overlapping peaks and detected overlapping voltammetric peak signals.