Crystalline cellulose elastic modulus predicted by atomistic models of uniform deformation and nanoscale indentation

The elastic modulus of cellulose Iβ in the axial and transverse directions was obtained from atomistic simulations using both the standard uniform deformation approach and a complementary approach based on nanoscale indentation. This allowed comparisons between the methods and closer connectivity to experimental measurement techniques. A reactive force field was used that explicitly describes hydrogen bond, coulombic and van der Waals interactions, allowing each contribution to the inter- and intra-molecular forces to be analyzed as a function of crystallographic direction. The uniform deformation studies showed that the forces dominating elastic behavior differed in the axial and transverse directions because of the relationship between the direction of the applied strain and the hydrogen bonding planes. Simulations of nanoscale indentation were then introduced to model the interaction between a hemispherical indenter with the $(1\bar{1}0)$ surface of a cellulose Iβ rod. The role of indenter size, loading force and indentation speed on the transverse elastic modulus was studied and, for optimized parameters, the results found to be in good agreement with experimentally-measured transverse elastic modulus for individual cellulose crystals.     

Effect of hydrogen bonds on the axial stiffnes of crystalline native cellulose

Previous theoretical calculations of elastic constants for cellulose based on force constants for bond stretching and bending of valence angles have yielded axial stiffness values admittedly too low. The present analysis accounts for a hitherto unexamined geometrical effect associated with deformation of interchain hydrogen bonds. To do this, most primary bond deformations are neglected so the resulting calculation gives an upper bound for the axial stiffness. By using two different sets of hydrogen bond force constants, values of 24.6 and 31.9 × 10¹¹ dyne/cm² were obtained for Young's modulus in the chain direction. These values are very much larger than earlier calculations and experimental determinations from cellulosic fibers, indicating both the importance of the effect considered here and the likelihood of an exact analysis yielding an acceptable result.     

Short cellulose nanofibrils as reinforcement in polyvinyl alcohol fiber

Short cellulose nanofibrils (SCNF) were investigated as reinforcement for polyvinyl alcohol (PVA) fibers. SCNF were mechanically isolated from hard wood pulp after enzymatic pretreatment. Various levels of SCNF were added to an aqueous PVA solution, which was gel-spun into continuous fibers. The molecular orientation of PVA was affected by a combination of wet drawing during gel spinning and post-hot-drawing at a high temperature after drying. A maximum total draw ratio of 27 was achieved with various SCNF contents investigated. The PVA crystal orientation increased when small amounts of SCNF were added, but decreased again as the SCNF content was increased above about 2 or 3 %, likely due to SCNF percolation resulting in network formation that inhibited alignment. SCNF fillers were effective in improving PVA fiber tensile properties (i.e., ultimate strength and elastic modulus). For example, the ultimate strength and modulus of PVA/SCNF composite fiber with a SCNF weight ratio of 6 were nearly 60 and 220 % higher than that of neat PVA. Shifts in the Raman peak at ~1,095 cm^?1, which were associated with the C–O–C glycosidic bond of SCNF, indicated good stress-transfer between the SCNF and the PVA matrix due to strong interfacial hydrogen bonding. Cryogenic fractured scanning electron microscopy images of filled and unfilled PVA fibers showed uniform SCNF dispersion.     

Nanofibrillation of dried pulp in NaOH solutions using bead milling

The drying process in typical pulp production generates strong hydrogen bonding between cellulose microfibrils in refined cell walls and increases the difficulty in obtaining uniform cellulose nanofibers. To investigate the efficacy of alkaline treatment for cellulose nanofibrillation, this study applied a bead-milling method in NaOH solutions for the nanofibrillation of dried pulps. NaOH treatments loosened the hydrogen bonding between cellulose microfibrils in dried pulps and allowed preparation of cellulose nanofibers in 8 % NaOH with a width of approximately 12–20 nm and a cellulose I crystal form. Both the nanofiber suspensions prepared in 8 and 16 % (w/w) NaOH were formed into hydrogels by neutralization because of surface entanglement and/or interdigitation between the nanofibers. When the dried pulp was fibrillated in 16 % (w/w) NaOH, the sample after neutralization had a uniquely integrated continuous network. These results can be applied to the preparation of high-strength films and fibers with cellulose I crystal forms without prior dissolution of pulps.     

Young’s modulus calculations for cellulose I<Subscript>β</Subscript> by MM3 and quantum mechanics

Quantum mechanics (QM) and molecular mechanics (MM) calculations were performed to elucidate Young’s moduli for a series of cellulose I_β models. Computations using the second generation empirical force field MM3 with a disaccharide cellulose model, 1,4′-O-dimethyl-β-cellobioside (DMCB), and an analogue, 2,3,6,2′,3′,6′-hexadeoxy-1,4′-O-dimethyl-β-cellobioside (DODMCB), that cannot make hydrogen bonds reveal a considerable contribution of intramolecular hydrogen bonding to the molecular stiffness of cellulose I_β; the moduli for DMCB and DODMCB being 85.2 and 37.6 GPa, respectively. QM calculations confirm this contribution with modulus values of 99.7 GPa for DMCB and 33.0 GPa for DODMCB. However, modulus values for DMCB were considerably lower than values previously reported for cellulose I_β. MM calculations with extended cellulose chains (10–40 glucose units) resulted in modulus values, 126.0–147.5 GPa, more akin to the values reported for cellulose I_β. Comparison of the cellodecaose model, 1,4′-O-dimethyl-β-cellodecaoside (DMCD), modulus with that of its hydrogen bonding-deficient analogue, 2,3,6,2′,3′,6′-hexadeoxy-1,4′-O-dimethyl-β-cellodecaoside (DODMCD), corroborates the observed stiffness conferred by intramolecular hydrogen bonds; the moduli for DMCD and DODMCD being 126.0 and 63.3 GPa, respectively. Additional MM3 determinations revealed that modulus values were not strongly affected by intermolecular hydrogen bonding, with multiple strand models providing values similar to the single strand models; 87.5 GPa for a 7-strand DMCB model and 129.5 GPa for a 7 strand DMCD model.     

Preparation of tough cellulose II nanofibers with high thermal stability from wood

Well-dispersed cellulose II nanofibers with high purity of 92 % and uniform width of 15–40 nm were isolated from wood and compared to cellulose I nanofibers. First, ground wood powder was purified by series of chemical treatments. The resulting purified pulp was treated with 17.5 wt% sodium hydroxide (NaOH) solution to mercerize the cellulose. The mercerized pulp was further mechanically nanofibrillated to isolate the nanofibers. X-ray diffraction patterns revealed that the purified pulp had been transformed into the cellulose II crystal structure after treatment with 17.5 wt% NaOH, and the cellulose II polymorph was retained after nanofibrillation. The cellulose II nanofiber sheet exhibited a decrease in Young’s modulus (8.6 GPa) and an increase in fracture strain (13.6 %) compared to the values for a cellulose I nanofiber sheet (11.8 GPa and 7.5 %, respectively), which translated into improved toughness. The cellulose II nanofiber sheet also showed a very low thermal expansion coefficient of 15.9 ppm/K in the range of 20–150 °C. Thermogravimetric analysis indicated that the cellulose II nanofiber sheet had better thermal stability than the cellulose I nanofiber sheet, which was likely due to the stronger hydrogen bonds in cellulose II crystal structure, as well as the higher purity of the cellulose II nanofibers.     

Comparison of polyethylene glycol adsorption to nanocellulose versus fumed silica in water

The recent intensification of industrially produced cellulose nanocrystals (CNCs) and cellulose nanofibrils has positioned nanocelluloses as promising materials for many water-based products and applications. However, for nanocelluloses to move beyond solely an academic interest, a thorough understanding of their interaction with water-soluble polymers is needed. In this work, we address a conflicting trend in literature that suggests polyethylene glycol (PEG) adsorbs to CNC surfaces by comparing the adsorption behaviour of PEG with CNCs versus fumed silica. While PEG is known to have strong hydrogen bonding tendencies and holds water tightly, it is sometimes (we believe erroneously) presumed that PEG binds to cellulose through hydrogen bonding in aqueous media. To test this assumption, the adsorption of PEG to CNCs and fumed silica (both in the form of particle films and in aqueous dispersions) was examined using quartz crystal microbalance with dissipation, isothermal titration calorimetry, rheology and dynamic light scattering. For all PEG molecular weights (300–10,000 g/mol) and concentrations (100–10,000 ppm) tested, strong rapid adsorption was found with fumed silica, whereas no adsorption to CNCs was observed. We conclude that unlike silanols, the hydroxyl groups on the surface of CNCs do not readily hydrogen bond with the ether oxygen in the PEG backbone. As such, this work along with previous papermaking literature supports the opinion that PEG does not adsorb to cellulose surfaces.     

Structure conversions of cellulose IIII crystal models in solution state: a molecular dynamics study

This paper re-examines our previous molecular dynamics (MD) study on cellulose III_I crystal models with finite dimensions solvated in explicit water molecules. Eight crystal models, differing in a constituent lattice plane and dimensions, were studied. One calculation allowed for O–H and C–H bond stretching, and had a small time step of 0.5 fs. The other calculation adopted non-scaling factors of the 1–4 non-bonded interactions. As in our previous study, in the former MD calculations, six of the eight crystal models exhibited structure conversion with cooperative chain slippages generated by a progressive fiber bend. This converted the initial non-staggered chain packing of cellulose III_I into a near one-quarter staggering and gave the crystal model a triclinic-like configuration. In contrast, in the non-1–4 scaling MD calculations, all of the eight crystal models retained the initial cellulose III_I crystal structure. Another series of non-1–4 scaling MD calculations were performed for the four crystal models containing chains with a degree of polymerization (DP) of 40 at 370 K, which simulated hot water treatment to convert cellulose III_I to Iβ. Some of the hydroxymethyl groups irreversibly rotated from gt into tg conformation. This accompanied exchange of the intrasheet hydrogen bonding scheme along the (1 ?1 0) lattice plane from O2–O6 to O3–O6. The original corrugated (1 ?1 0) chain sheet was partly converted into a cellulose I-like flat chain sheet.     

Application of polarizable ellipsoidal force field model to pnicogen bonds

Noncovalent interactions, such as hydrogen bonds and halogen bonds, are frequently used in drug designing and crystal engineering. Recently, a novel noncovalent pnicogen bonds have been identified as an important driving force in crystal structures with similar bonding mechanisms as hydrogen bond and halogen bond. Although the pnicogen bond is highly anisotropic, the pnicogen bond angles range from 160° to 180° due to the complicated substituent effects. To understand the anisotropic characters of pnicogen bond, a modification of the polarizable ellipsoidal force field (PEff) model previously used to define halogen bonds was proposed in this work. The potential energy surfaces (PESs) of mono‐ and polysubstituted PH₃–NH₃ complexes were calculated at CCSD(T), MP2, and density functional theory levels and were used to examine the modified PEff model. The results indicate that the modified PEff model can precisely characterize pnicogen bond. The root mean squared error of PES obtained with PEff model is less than 0.5 kcal/mol, compared with MP2 results. In addition, the modified PEff model may be applied to other noncovalent bond interactions, which is important to understand the role of intermolecular interactions in the self‐assembly structures. © 2015 Wiley Periodicals, Inc.     

Diversity of potential hydrogen bonds in cellulose I revealed by molecular dynamics simulation

We have performed molecular dynamics calculations using a revised version of the Gromos56A_carbo force field to understand the consequences of the different potential hydrogen bonding patterns on the structural stability and thermal behavior of the I_α and I_β forms of native cellulose. For each allomorph, we considered three patterns of hydrogen bonds: two patterns obtained from neutron diffraction data refinement and a regular mixture of the two. Upon annealing, the hydrogen bonding schemes of cellulose I_β, irrespective of the starting structure, re-arranged into the main hydrogen bond pattern experimentally observed (pattern A). On the other hand, the I_α structures, irrespective of the starting hydrogen bonding pattern, converged to a non-experimental structure where the adjacent chains are shifted along the chain direction by 0.12 nm in the hydrogen-bonded plane, and the hydroxymethyl group conformation alternates between gt and tg along the chain. The exotic structure in I_α might be a consequence of a deficiency in force field parameters and/or potential molecular arrangement in less crystalline cellulose.     

Neutron crystallographic and molecular dynamics studies of the structure of ammonia-cellulose I: rearrangement of hydrogen bonding during the treatment of cellulose with ammonia

The hydrogen bond arrangement in a complex of cellulose with ammonia has been studied using neutron crystallography in combination with molecular dynamics simulations. The O6 atom of the hydroxymethyl group is donor in a highly occupied hydrogen bond to an ammonia molecule. This rotating ammonia molecule is donor in partially occupied and transient hydrogen bonds to the O2, O3 and O6 atoms of the hydroxyl groups of other chains. The hydrogen atom bound to the O3 atom is disordered but it is almost always involved in some type of hydrogen bonding. It is donated in a hydrogen bond most of the time to the O5 atom on the same chain. However, it also rotates away from this O5 atom to be donated to an ammonia molecule part of the time. On the other hand the hydrogen atom bound to the O2 atom is free from hydrogen bonding most of the time. It is donated in a hydrogen bond to the O6 atom on a neighboring chain only with a relatively small probability. These results provide new insights into how hydrogen bonds are rearranged during the conversion of cellulose I to cellulose III_I by ammonia treatment.     

An empirical potential for the O-H ⋯ O hydrogen bond: Part 2. Comparison with observed hydrogen bond geometries

The room temperature distribution of O-H ? O hydrogen bond geometries has been predicted by a Monte Carlo calculation, with an empirical potential energy function for the hydrogen bond. The results are compared with a recent survey of hydrogen bonds in carbohydrate crystal structures. The calculated and observed distributions of the O-H ? O angle have mean values of 165.5° and 167.1° respectively. Both the theoretical and experimental results suggest that short O ? H hydrogen bonds tend to be more linear than long O ? H bonds. The distribution of hydrogen bonding within the lone pair plane of the acceptor oxygen atom is predicted to be broader than the distribution perpendicular to this plane, in agreement with the experimental data. The empirical hydrogen bond function, in conjunction with the molecular mechanics program MMI, has also been used to predict the geometries of inter-residue hydrogen bonds in five disaccharides. The O ? H distances and O-H ? O angles are reproduced with r.m.s. deviations of 0.06 Å and 9° respectively.     

Gelation behavior of cellulose in NaOH/urea aqueous system via cross-linking

Sol–gel transition of cellulose solution in NaOH/urea aqueous solution with the addition of epichlorohydrin (ECH) was investigated by rheological means. The gelation was controlled by a synergy of chemical and physical cross-linking processes, namely, the etherification reaction between cellulose and ECH as well as the self-association and entanglement of cellulose chains via hydrogen bonding re-construction in NaOH/urea. The results revealed that the cross-linker concentration, cellulose concentration and temperature played important roles in the gelation behavior. The gel time decreased with increasing either ECH or cellulose concentration, and the gel temperature dropped from 38 to 28 °C with an increase of cellulose concentration from 4 to 6 wt%, i. e. easier gelation was reached with higher cross-linker concentration, cellulose concentration or temperature, since higher cross-linker or cellulose concentration led to more network junctions via chemical or physical cross-linking, while higher temperature was favorable to both the etherification reaction and re-construction of cellulose hydrogen bonds. The compressive modulus of cellulose/ECH hydrogels was improved a lot by increasing either cellulose or ECH concentration, indicating the chemical cross-linking obviously improved the mechanical property, on the other hand, the swelling property could be tunable by changing the gelation parameter. This work supplied useful information to the control and optimization of the structure and properties of cellulose based hydrogels.     

2‐(2‐Oxazolin‐2‐yl)benzene‐1,4‐diol: X‐ray and density functional theory studies

In the crystal structure of the title compound, C₉H₉NO₃, there are strong intra­molecular O—H⋯N and inter­molecular O—H⋯O hydrogen bonds which, together with weak inter­molecular C—H⋯O hydrogen bonds, lead to the formation of infinite chains of mol­ecules. The calculated inter­molecular hydrogen‐bond energies are −11.3 and −2.7 kJ mol⁻¹, respectively, showing the dominant role of the O—H⋯O hydrogen bonding. A natural bond orbital analysis revealed the electron contribution of the lone pairs of the oxazoline N and O atoms, and of the two hydr­oxy O atoms, to the order of the relevant bonds.     

Hydrogen bonding interactions between adsorbed polymer molecules and crystal surface of acetaminophen

The objective of this work was to investigate whether or not the hydrogen bonding interaction between polymer and crystal surface can be detected by the etching pattern changes in the presence of polymers. The (010) face of acetaminophen single crystal was used as a model solid surface. The etching patterns on the (010) face of acetaminophen crystal by water are in the directions of a- and c-axes, which are the same as the directions of the dominant attachment energies on the (010) face. In the presence of polymer, the hydrogen bonding interactions between adsorbed polymer and crystal surface can affect surface diffusion of acetaminophen molecules and change the etching patterns in the direction of a-axis, i.e., the direction of one hydrogen bond chain. Studies with 2-hydroxyethyl cellulose (HEC), hydroxypropyl cellulose (HPC), hydroxypropyl methylcellulose (HPMC) and poly(vinyl alcohol) (PVA) showed that polymers, which can form hydrogen bonds with acetaminophen crystal surface, can change etching patterns in the direction of a-axis. Study with Dextran suggested that if a polymer cannot form hydrogen bonds with crystal surface due to steric repulsion, it will not change the etching pattern in the direction of a-axis. Studies with poly(ethylene glycol) (PEG) and poly(propylene glycol) (PPG) further confirmed that only if a polymer can form hydrogen bonds with acetaminophen on crystal surface, the etching patterns in the direction of a-axis will be affected. The study results suggest that in the presence of polymers, the etching pattern change in the direction of hydrogen bond chain, the a-axis of acetaminophen crystals, can be used to indicate the existence of the hydrogen bonding interactions between adsorbed polymers and acetaminophen crystal surface.     

Electron (charge) density studies of cellulose models

Introductory material first describes electron density approaches and demonstrates visualization of electron lone pairs and bonding as concentrations of electron density. Then it focuses on the application of Bader’s Quantum Theory of Atoms-in-Molecules (AIM) to cellulose models. The purpose of the work is to identify the various interactions that stabilize cellulose structure. AIM analysis aids study of non-covalent interactions, especially those for which geometric criteria are not well established. The models were in the form of pairs of cellotriose molecules, methylated at the O1 and O4 ends. Based on the unit cell of cellulose Iβ, there were corner–corner, and center–center pairs that correspond to (200) sheets, and corner–center pairings that corresponded to (1–10) and (110) stacks. AIM analysis (or charge-density topology analysis) was applied before and after minimization in vacuum and in continuum solvation. Besides the conventional O–H···O hydrogen bonds, all of which were known from geometric criteria, C–H···O hydrogen bonds (some previously reported), and some O···O and H···H interactions were found. Non-covalent bonds in the (200) sheets were maintained in all calculations with the exception of a weak, bifurcated O6–H···O2′′ bond that was not found in the corner–corner pair model and did not survive minimization. Nor did the O6···O4 interactions on the reducing ends of the triosides. Pairs of molecules along the (110) plane had an equal number (12) of non-covalent bonds compared to the pairs along the (1–10) plane, but the AIM parameters indicated the bonds between the pairs in the (110) plane were weaker. Intra-molecular O–H···O hydrogen bonds survived in these minimized pairs, but the relative chain alignments usually did not.     

Crystal Structure and Conformation of a 10-Thiabilirubin

Summary.  A crystal structure determination of a bilirubin analog with a sulfur instead of a C(10)–CH₂ linking the two dipyrrinones is reported. Conformation-determining torsion angles and key hydrogen bond distances and angles are compared to those obtained from molecular dynamics calculations as well as to the corresponding data from X-ray determinations and molecular dynamics calculations of bilirubin. Like other bilirubins, the component dipyrrinones of the analog are present in the bis-lactam form with (Z)-configurated double bonds at C(4) and C(15). Despite the large differences in bond lengths and angles at –S–vs.–CH₂–, the crystal structure shows considerable similarity to bilirubin: both pigments adopt a folded, intramolecularly hydrogen-bonded ridge-tile conformation stabilized by six hydrogen bonds – although the interplanar angle of the ridge-tile conformation of the title compound is smaller (∼ 86°) than that of bilirubin (∼ 98°). The collective data indicate that even with long C–S bond lengths and a smaller C–S–C bond angle at the pivot point on the ridge-tile seam, intramolecular hydrogen bonding persists. Received August 16, 2001. Accepted September 12, 2001     

Concentration effects on the isolation and dynamic rheological behavior of cellulose nanofibers via ultrasonic processing

Chemical pretreatment combined with high-intensity ultrasonication was performed to disintegrate cellulose nanofibers from poplar wood powders. The cellulose content in each suspension was treated as the control variable because the suspension concentration significantly influences the properties of the resultant cellulose nanofibers via ultrasonic processing. The as-obtained cellulose nanofibers were characterized by fiber diameter distribution, crystal structure, and rheological analysis. An increase of not more than 1.2 % of the cellulose content resulted in finer nanofibers. Both storage modulus and loss modulus of cellulose nanofiber suspensions rapidly increased with increasing concentration because of the gradual formation of a stronger network structure. In addition, the dynamic mechanical behavior of suspensions with fiber contents lower than 0.8 % was affected by the frequency and temperature alteration in contrast with the suspension with higher fiber contents. The sol–gel transformation and the visco-elastic transition depend on the hydroxyl bonding and the cross-linking extent of cellulose nanofibers in various concentration environments.     

Understanding the mechanism of cellulose dissolution in 1-butyl-3-methylimidazolium chloride ionic liquid via quantum chemistry calculations and molecular dynamics simulations

While N,N′-dialkylimidazolium ionic liquids (ILs) have been well-established as effective solvents for dissolution and processing of cellulose, the detailed mechanism at the molecular level still remains unclear. In this work, we present a combined quantum chemistry and molecular dynamics simulation study on how the ILs dissolve cellulose. On the basis of calculations on 1-butyl-3-methylimidazolium chloride, one of the most effective ILs dissolving cellulose, we further studied the molecular behavior of cellulose models (i.e. cellulose oligomers with degrees of polymerization n = 2, 4, and 6) in the IL, including the structural features and hydrogen bonding patterns. The collected data indicate that both chloride anions and imidazolium cations of the IL interact with the oligomer via hydrogen bonds. However, the anions occupy the first coordination shell of the oligomer, and the strength and number of hydrogen bonds and the interaction energy between anions and the oligomer are much larger than those between cations and the oligomer. It is observed that the intramolecular hydrogen bond in the oligomer is broken under the combined effect of anions and cations. The present results emphasize that the chloride anions play a critically important role and the imidazolium cations also present a remarkable contribution in the cellulose dissolution. This point of view is different from previous one that only underlines the importance of the chloride anions in the cellulose dissolution. The present results improve our understanding for the cellulose dissolution in imidazolium chloride ILs.     

AIM and NBO analyses on hydrogen bonds formation in sugar-based surfactants (α/β-d-mannose and n-octyl-α/β-d-mannopyranoside): a density functional theory study

Density functional theory calculations on α/β-d-mannose (α/β-d-Man) and the corresponding glycosides of n-octyl-α/β-d-mannopyranoside (C₈O-α/β-d-Man) were carried out for geometrical optimisation and stability predictions at the B3LYP/6-31G level of theory. These compounds are related anomerically, since they differ by only the orientation of the hydroxyl group at the C1 position. The aim of this study is to investigate the effect of the hydroxyl group's orientations (axial vs. equatorial) at the C1 position on the intra-molecular interactions and the conformational stability of these isomers. The structural parameters of X-H???Y intra-molecular hydrogen bonds were analysed, while the nature of these bonds was considered using the atoms-in-molecules (AIM) approach. Natural bond orbital (NBO) analysis was used to determine bond orders and the effective non-bonding interactions. We have also reported thermodynamic properties and the electronic properties, such as the highest occupied molecular orbital, lowest unoccupied molecular orbital, ionisation energy, electron affinity, electronic chemical potential, chemical hardness, softness and electrophilicity index in the gas phase for all compounds. These results showed that while α-anomers possess only one intra-molecular hydrogen bond, β-anomers possess two intra-molecular hydrogen bonds, which further confirms the anomalous stability of the latter in the self-assembly phenomena.