Intermolecular Interactions and 3D Structure in Cellulose–NaOH–Urea

Aug 11, 2014 - This material is available free of charge via the Internet at ..... and take-up speed, and dope characteristics, e.g. cellulose DP and ...
0 downloads 3 Views 2MB Size
Article pubs.acs.org/JPCB

Intermolecular Interactions and 3D Structure in Cellulose−NaOH− Urea Aqueous System Zhiwei Jiang,†,‡ Yan Fang,‡ Junfeng Xiang,† Yanping Ma,† Ang Lu,‡ Hongliang Kang,† Yong Huang,† Hongxia Guo,† Ruigang Liu,*,† and Lina Zhang*,‡ †

State Key Laboratory of Polymer Physics and Chemistry, Beijing National Laboratory of Molecular Sciences, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China ‡ Department of Chemistry, Wuhan University, Wuhan 430072, China S Supporting Information *

ABSTRACT: The dissolution of cellulose in NaOH/urea aqueous solution at low temperature is a key finding in cellulose science and technology. In this paper, 15N and 23Na NMR experiments were carried out to clarify the intermolecular interactions in cellulose/NaOH/urea aqueous solution. It was found that there are direct interactions between OH− anions and amino groups of urea through hydrogen bonds and no direct interaction between urea and cellulose. Moreover, Na+ ions can interact with both cellulose and urea in an aqueous system. These interactions lead to the formation of cellulose−NaOH−urea−H2O inclusion complexes (ICs). 23Na relaxation results confirmed that the formation of urea−OH− clusters can effectively enhance the stability of Na+ ions that attracted to cellulose chains. Low temperature can enhance the hydrogen bonding interaction between OH− ions and urea and improve the binding ability of the NaOH/urea/H2O clusters that attached to cellulose chains. Cryo-TEM observation confirmed the formation of cellulose− NaOH−urea−H2O ICs, which is in extended conformation with mean diameter of about 3.6 nm and mean length of about 300 nm. Possible 3D structure of the ICs was proposed by the M06-2X/6-31+G(d) theoretical calculation, revealing the O3H···O5 intramolecular hydrogen bonds could remain in the ICs. This work clarified the interactions in cellulose/NaOH/urea aqueous solution and the 3D structure of the cellulose chain in dilute cellulose/NaOH/urea aqueous solution.



INTRODUCTION Cellulose is the most abundant renewable biopolymer on the earth and the potential resource of biofuels1 and new functional materials.2−4 The stereoregularity of the linear poly(1 → 4)-βD-glucan chain and the abundant hydroxyl groups of cellulose result in the strong intra- and intermolecular hydrogen-bonding network in cellulose (Scheme 1). Therefore, native cellulose is resistant to chemical accessibility and enzymatic degradation.5,6 On the other hand, due to the strong hydrogen-bonding network, cellulose cannot be melted or be dissolved in common

solvents. The traditional approach for producing regenerated cellulose products is mainly based on viscose process, which is still dominated nowadays. In the viscose process, the usage of CS2 and the byproducts cause serious environmental pollution. Therefore, the finding of new environmental friendly cellulose solvents is always the key issue in both producing regenerated cellulose materials and synthesis new cellulose derivatives. Several cellulose solvent systems, such as N-methylmorpholineN-oxide/water (NMMO/H2O),7,8 ionic liquids (ILs),9−11 lithium chloride/N,N-dimethylacetamide (LiCl/DMAc),12−14 alkali aqueous solutions,15 and alkali/urea aqueous solution,16 have been discovered for dissolving cellulose, in which the NMMO/H2O system has been commercialized for producing regenerated cellulose fibers with the category name of Lyocell.7 NaOH/urea aqueous solution is a cheap, nonvolatile, nontoxic, and environmentally friendly cellulose solvent with quick dissolving power and has been acknowledged to be a milestone in the history of cellulose processing technology.17 Based on this solvent system, regenerated cellulose materials including fiber,18 film,19 particle,20 aerogel,21 and hydrogel22,23 have been prepared. However, the origin of the rapid

Scheme 1. Schematic Representation of Cellulose and Its Hydrogen Bond Network

Received: February 9, 2014 Revised: August 8, 2014 Published: August 11, 2014 © 2014 American Chemical Society

10250

dx.doi.org/10.1021/jp501408e | J. Phys. Chem. B 2014, 118, 10250−10257

The Journal of Physical Chemistry B

Article

TEM Observation. The cryo-TEM samples were prepared in a controlled environment vitrification system (CEVS) at 277 K. The cellulose/NaOH/urea aqueous solution with the cellulose mass fraction of 1 wt % was first filtered through a 0.22 μm filter. A micropipet was utilized to load 2.8 μL of the resulting solution onto a carbon-coated copper grid, which was blotted with two pieces of filter paper, resulting in the formation of thin films suspended on the mesh holes. After waiting for about 20 s, the samples were quickly plunged into a reservoir of liquid ethane (cooled by nitrogen) at 108 K. The vitrified samples were stored in liquid nitrogen until they were transferred to a cryogenic sample holder (Gatan 626) and then were observed with a FEI Tecnai 20 TEM (200 kV) at about 99 K. The images were recorded with a Gatan Ultrascan 894 CCD. TEM images were obtained by using a JEM-2010 (HT) transmission electron microscope (JEOL TEM, Japan). The samples were prepared by casting the diluted cellulose solution with the cellulose mass fraction of 0.10% onto a holey carbon film supported on copper grid. To form extremely thin liquid layer for TEM observation, the specimen was tilted during dried at room temperature.

dissolution of cellulose in NaOH/urea aqueous solution at low temperature and the intermolecular interactions in the cellulose/NaOH/urea aqueous solution are still unclear. Hydrogen bonding is a critical factor in condensed structures and dissolution of cellulose. Among the methods for the investigation of the interactions, nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for the investigation of hydrogen bonding interactions, which has been widely applied to investigate structures as well as the dissolution behaviors and mechanism of cellulose,24−29 while the intermolecular interactions between cellulose and solvents were somewhat insufficiently discussed. In present work, NMR experiments, including temperature-dependent 1H NMR, 13C NMR, 23Na, 15 N NMR, and 23Na relaxation, and transmission electron microscope (TEM) were used to study the intermolecular interactions in the cellulose−NaOH−urea aqueous system. The effect of temperature on the intermolecular interactions was also discussed. Moreover, the intermolecular interactions and 3D structure of cellulose in the solutions were further clarified by density functional theory (DFT), a powerful method for simulating the intermolecular interactions and the structure,30,31 based on the experimental data. Our findings would reveal the nature of cellulose dissolution in NaOH/urea solvent at low temperature and provide significant information for the chain structure of complicated macromolecule in solution.



RESULTS AND DISCUSSION Interactions between Urea and NaOH Revealed by 15 N NMR. Cellulose has a low solubility in NaOH aqueous solution without urea, and the resultant solutions are less stable.16 However, cellulose can be readily dissolved in NaOH/ urea aqueous solution at relatively high solubility, and the cellulose solutions are more stable than those in NaOH aqueous system. To clarify the role of urea on the dissolution of cellulose in NaOH/urea aqueous solution, 15N NMR spectra of the samples Ur12, Na7Ur12, and Na7Ur12-Cell were measured at the temperature of 261, 273, and 298 K (Figure 1). Urea and H2O were replaced by 15N-urea and D2O, respectively, in all the samples. Table 1 summarizes the chemical shifts of 15N NMR spectra of all samples at different temperatures. In the 15N NMR spectra of Ur12, the signal of urea-15N is a multiplet at 298 and 273 K, which is due to the 15N−2D (quintet) and 15 N−1H (triplet) spin−spin coupling. It is difficult to obtain the exact H/D ratio. The 15N NMR signal in the 15N NMR spectra of urea-15N in H2O is a triplet (Figure S1), which confirmed the rapid exchange of the N−H (N−D) proton in the present of base. At 261 K, the NMR signal of 15N becomes a singlet (Figure 1a), which is due to that the solution becoming frozen at this temperature and the resolution of the NMR signal of the frozen solution would be decreased.32 Interestingly, the NMR signal of 15N is a singlet in the 15N NMR spectra of all the samples in the presence of NaOH at all experimental temperatures (Figure 1b,c). This could be explained that the exchange rate of N−H (or N−D) proton in the presence of base is typically rapid.33 A similar phenomenon was observed in the 1H NMR spectra of urea/D2O and NaOH/urea/D2O (Figure S2). Moreover, the 15N resonance of urea with NaOH shifted to low field by 0.49 ppm at 298 K and 0.60 ppm at 273 K compared to that in urea-15N/D2O solution. The low-field shift of the 15N NMR signal is attributed to the formation of hydrogen bonds between amino groups of urea and OH− ions of NaOH in the system. As a result, the C−N bonds are in a more double bond character as shown in Scheme 2.34 Because of the cellulose cannot dissolve in urea aqueous solution, cellobiose, a constitutional unit of cellulose, was used as a model compound of cellulose to investigate the



EXPERIMENTAL SECTION Materials. Microcrystalline cellulose was purchased from Alfa Aesar. NaOH, NaCl, and urea were of analytical reagent (Shanghai Chemical Reagent Co. Ltd., China) and were used without further purification. The cellulose (cotton linter pulp) with Mw = 9.6 × 104 g/mol was provided by Hubei Chemical Fiber Group Ltd. (Xiangfan, China) and was only used in cryoTEM and TEM experiments. D2O of 99.9% purity for NMR was purchased from Cambridge Isotope Laboratories, Inc. Urea-15N with 98 atom % 15N for NMR was purchased from Sigma-Aldrich. Deionized water was used in all experiments. Sample Preparation. NaOH and urea were dissolved in D2O to obtain 12 wt % urea, 7 wt % NaOH, and 7 wt % NaOH/12 wt % urea D2O solutions, denoted as Ur12, Na7, and Na7Ur12 for short, respectively. Cellulose was added to Na7 and Na7U12 solutions, and the mixtures were stirred for 2 h at room temperature. The mixtures were then transferred into a freezer to become solid for 30 min. The resultant samples were placed at room temperature and thus obtained clear cellulose solutions, denoted as Na7-Cell and Na7Ur12-Cell, respectively. The content of cellulose was kept at 4 wt % for NMR experiments. NMR Measurement. 15N NMR spectra were recorded on a Bruker AVANCE 600 NMR spectrometer with a TBO probe at a settled temperature range from 261 to 298 K (−12 to 25 °C). 23 Na and 13C NMR measurements were carried out on a Bruker AVIII 500WB spectrometer with a TBO probe at certain temperature. 1H NMR measurements were carried out on a Bruker AVANCE 400 at 298 K. Deuterium oxide (D2O) was used as lock signal. Longitudinal relaxation times (T1) were measured by using the inversion recovery sequence. 13C NMR chemical shifts were externally referenced to CDCl3 at 75.38 ppm. Nitromethane (CH3NO2) was used as the external chemical-shift reference at 382 ppm for 15N NMR experiments. 23 Na NMR chemical shifts were externally referenced to methanol solution of 0.7 M NaBPh4 at 25 °C at −3.087 ppm. 10251

dx.doi.org/10.1021/jp501408e | J. Phys. Chem. B 2014, 118, 10250−10257

The Journal of Physical Chemistry B

Article

Scheme 2. Structure of OH−−Urea Complex

urea-15N/D2O solutions are shown in Figure S3. The results show that the triplet and quintet NMR signals of 15N still exist in 15N NMR spectra of cellobiose/urea-15N/D2O solution, although the peaks are not clearly separated from each other. Furthermore, there is small variation in chemical shift when cellobiose was added to urea/D2O solution. The results confirmed that there are weak interactions between cellulose and urea. Generally, electrons in p-orbitals of the double bond can give rise to a very large paramagnetic effect, resulting in a downfield shift. Increasing in hydrogen bonding interaction can cause the nitrogen lone pairs to a more delocalized state and makes the nitrogen atoms more positive. However, positive nitrogen atoms reduce the second-order paramagnetic effect at the 15N nucleus and lead to an upfield shift of the NMR signal. With decreasing temperature, the nitrogen resonance of urea aqueous solutions with and without NaOH shifted to upfield, which indicates the enhancement of the intermolecular hydrogen bonding interactions for urea alone and those associated with OH− ions of NaOH. Our results are similar to the association interaction between urea and acetate at different temperatures.34 Moreover, the 15N chemical shift in the NMR spectrum of Na7Ur12 is similar to that of Na7Ur12-Cell at the temperature investigated. In addition, the 13C NMR chemical shifts of cellulose in Na7Ur12-Cell remained fundamentally unchanged compared with those in cellulose/NaOH/D2O solution (Figure S4 and Table S1), which also confirmed that there is no strong interaction between urea and cellulose. Combined with the above results, the role of urea in cellulose−NaOH−urea aqueous system can be summarized as follows. The amino groups of urea bound to OH− ions of NaOH through hydrogen bonds to form complex in aqueous solution. There is almost no direct interaction between urea and cellulose. Low temperature can enhance the hydrogen bonding interactions between OH− ions of NaOH and urea. Cluster Formation Revealed by 23Na NMR. For the clarification of the role of Na+ ions on the dissolution of cellulose, 23Na NMR measurements were carried out for Na7, Na7Ur12, Na7-Cell, and Na7Ur12-Cell solutions, and the results are shown in Figure 2. Only the singlet peak, assigned to 23Na characteristic peak, can be observed on the 23Na NMR spectra of all the samples at 298 K (Figure 2). Adding urea or cellulose in NaOH aqueous solution caused the upfield shift of 23Na NMR signal. Moreover, the presence of both urea and cellulose in NaOH/D2O further resulted in the upfield shift of 23Na resonance signal. The results indicate that Na+ ion can interact with both urea and cellulose in the system. In general, the position of the NMR signal is determined by the total shielding that the nucleus under investigation receives from various sources.35 The shielding can be expressed by the screening constant, σ. On a practical basis, the screening constant can be divided into three independent parts by36 σ = σd + σo + σp (1)

Figure 1. 15N NMR spectra of (a) Ur12, (b) Na7Ur12, and (c) Na7Ur12Cell at the temperature of 261, 273, and 298 K. 15N-urea and D2O were used instead of urea and H2O in all the samples.

Table 1. Chemical Shift of 15N NMR Spectra of Several Solutions chemical shift (ppm) sample

261 K

273 K

298 K

NaOH/urea-15N cellulose/NaOH/urea-15N urea-15N

75.90 75.92

75.99 76.02 75.42

76.12 76.12 75.63

where σd and σo correspond to the contribution from local diamagnetic and other atoms, respectively, which are generally

interactions between urea and cellulose in aqueous solution. The 15N NMR spectra of urea-15N/D2O and cellobiose/ 10252

dx.doi.org/10.1021/jp501408e | J. Phys. Chem. B 2014, 118, 10250−10257

The Journal of Physical Chemistry B

Article

shift decreases linearly with the increase in temperature, which suggested that low temperature could enhance the interactions between Na+ ions and surround solvent or ions (water, urea, or OH− ions). The slopes of the lines obtained by linear fitting experimental data are denoted in Figure 3. The slope for Na7 is nearly equal to that for Na7Ur12, which confirmed that the Na+−urea interaction is independent of the temperature. It has been reported that the interaction of Na+−Z− ion pairs (Z− represents the anions such as OH−, Cl−, NO3−, etc.) has no effect on 23Na shift at various temperature.37 Therefore, the shift of 23Na NMR signal as a function of temperature is mainly caused by the Na+−H2O interaction. Combined with 23Na and 15N NMR results, it can be concluded that NaOH can directly interact with urea and cellulose and there was no direct association interaction between urea and cellulose. Based on the above results and discussion, the interactions in the cellulose−urea−NaOH− H2O system can be depicted as follows. Inclusion complexes (ICs) associated with NaOH, urea, cellulose chains, and water clusters were formed in the system and urea hydrates located at shell of the ICs, by which to prevent the self-association of the separated cellulose chains. Thereafter, cellulose molecules can be dispersed at the molecular level to form a stable solution in the NaOH/urea aqueous system. Low temperature can promote the hydrogen-bonding interactions in the system and mainly enhance the Na+−H2O interactions as well as the hydrogen-bonding interactions between OH− ions and urea. Intermolecular Interaction Revealed by 23Na Relaxation. The 23Na nucleus has a spin quantum number I = 3/2 and an electric quadrupole moment Q. The interaction between the Q and the electric field gradient (dE/dr) at the nucleus will give rise to the nuclear magnetic relaxation mechanism.38 Sodium is monatomic, and the electric field gradient originates from intermolecular interactions. Therefore, the relaxation rate of 23Na was very sensitive to its environment. To further investigate the role of Na+ ion in different cellulose solutions and solvent systems, the longitudinal relaxation time T1 of 23Na in Na7, Na7Ur12, Na7-Cell, and Na7Ur12-Cell was determined at different temperature (263−298 K). For all experiments, the 23Na resonance longitudinal relaxation times were monoexponential. The temperature dependence of the relaxation rate (1/T1) of 23Na for Na7, Na7Ur12, Na7-Cell, and Na7Ur12-Cell is shown in Figure 4, and the corresponding data are listed in Table S2. The results indicated that the addition of

Figure 2. 23Na NMR spectra of Na7, Na7Ur12, Na7-Cell, and Na7Ur12Cell at 298 K.

minor on nuclei other than protons. This is a result of the anisotropy in the magnetic susceptibilities of the distant atoms from the resonant nucleus. σp is the paramagnetic term and is usually the dominant shielding term. This contribution is due to the net orbital angular momentum of electrons in the valence orbitals of the nucleus under investigation and makes a negative contribution to the screening. The most important factor affecting these chemical shifts is the orbital angular momentum of electrons occupying the n = 3 shell of sodium and its effect upon σp. This arises from the overlapping of the outer p orbitals of the central sodium ion with the outer s and p orbitals of the neighboring solvent molecules or ions.35 The 23 Na shift values increase with increasing the σp term. As a result, the presence of urea or cellulose in NaOH aqueous solution led to the upfield shift of 23Na resonance signal, indicating that the overlap between the orbitals of Na+ and H2O is larger than that between Na+ and urea or cellulose. This also meant that the Na+−H2O interaction is stronger than the Na+−cellulose or Na+−urea interaction.37 Because urea can dramatically improve the dissolution of cellulose in NaOH aqueous solution, especially at low temperature,16 the temperature-dependent 23Na NMR for Na7 and Na7Ur12 solutions were also used to explore the effect of temperature on the Na+−urea interaction. The 23Na chemical shift for Na7 and Na7Ur12 as a function of temperature is shown in Figure 3. The results show that the 23Na chemical

Figure 3. Temperature dependences of the 23Na NMR chemical shift for Na7 and Na7Ur12.

Figure 4. Temperature dependences of the Na7, Na7Ur12, Na7-Cell, and Na7Ur12-Cell. 10253

23

Na relaxation rates for

dx.doi.org/10.1021/jp501408e | J. Phys. Chem. B 2014, 118, 10250−10257

The Journal of Physical Chemistry B

Article

urea and/or cellulose to the NaOH aqueous solutions or the decrease in the temperature of the samples caused a significant increase in 1/T1. Generally, a large 1/T1 corresponds to a low mobility. Therefore, the complex structure formed by attaching the NaOH−urea−H2O cluster on cellulose chains is relatively more stable than that cellulose−NaOH−H2O. These results confirmed that urea can effectively enhance the stability of NaOH attached to the cellulose chains and thus favor to the breakdown of intermolecular hydrogen bonds among cellulose chains. Low temperature favors for the stabilization of the complex structure in the system. Arrhenius plots for the dependence of the relaxation rate (1/T1) of Na7, Na7Ur12, Na7Cell, and Na7Ur12-Cell on temperature are shown in Figure 5.

Figure 5. Arrhenius plot for the dependence of the relaxation rate (1/ T1) of Na7, Na7Ur12, Na7-Cell, and Na7Ur12-Cell on temperature. The number represents the slope of the line.

Figure 6. Typical images of cryo-TEM (a) and HR TEM (b) of cellulose in NaOH/urea aqueous solution and the model of the cellulose ICs (c).

Based on the slopes of the relevant fitting lines, the relative activation energies of the Na7, Na7Ur12, Na7-Cell, and Na7Ur12Cell were calculated to obtain 21.0, 22.3, 20.1, and 20.5 kJ/mol, respectively. The results further confirmed that Na+ ions in Na7, Na7Ur12, Na7-Cell, and Na7Ur12-Cell possess different interactions with the surrounding species in solutions. Intermolecular Hydrogen Bonds in Cellulose/NaOH/ Urea Aqueous Solution. It is known that the partial dissociation of the C6-OH is the crucial factor for cellulose dissolved in NaOH aqueous solution.39 The binding of NaOH to hydroxyl groups of cellulose breaks the intermolecular hydrogen bonds in cellulose and leading to the dissolution of cellulose. In previous work, cellulose chains were proposed in the extended wormlike conformation in NaOH/urea aqueous solution.16 In order to clarify the morphology of cellulose chains in cellulose/NaOH/urea aqueous solutions, cryo-TEM observation was carried out. The samples for cryo-TEM observation were prepared by a rapid deep freezing procedure in liquid nitrogen at 99 K, by which the cellulose−NaOH− urea−H2O ICs in the solution can be successfully maintained in their original morphology and the formation of salt crystals can be prevented. Figure 6 shows the typical images of cryo-TEM and TEM of cellulose/NaOH/urea aqueous solution. As expected, extended wormlike structure was observed in both cryo-TEM and TEM images, which proved that ICs of cellulose−NaOH−urea exist in the solution. The extended wormlike ICs are about 250−350 nm in length and about 3.6 ± 0.4 nm in diameter evaluated by cryo-TEM. The length extended wormlike ICs is comparable with the contour length

(Lc = 300 nm) of cellulose calculated from the length of a projection of the repeating unit (Lunit = 0.5 nm) and the molecular weight of the repeating unit (Munit = 162 g/mol). The diameter of the ICs estimated by TEM is about 4.7 ± 0.6 nm, which is larger than that of cryo-TEM observation due to the collapse of the ICs during drying procedure for preparing TEM samples. Moreover, the cellulose ICs are well dispersed in the solution, indicating that the ICs are very stable in the solution system. It can be explained that urea was as shell of the ICs to prevent the cellulose chains to aggregate each other. On the basis of the above results and discussion, the steric structure of the ICs could be estimated by the M06-2X/631+G(d) level of theory, which is due to that M06-2X has been proven to have better performance than B3LYP for main-group noncovalent interactions.40−42 The M06-2X/6-31+G(d) method is implemented in the Gaussian 03 package.43 Previous work proved that the Na+ or OH− ions are generally in their hydrated forms, and the OH− ions could directly bind to the amino groups of urea molecules as confirmed by 15N NMR and 23Na NMR results. To further reveal the microscopic interaction among the solvent small molecules in NaOH/urea aqueous solutions, several kinds of clusters, such as NaOH(H2O)5·urea, NaOH(H2O)6·urea, NaOH(H2O)7·urea, NaOH(H2O)8·urea, and NaOH(H2O)9·urea, were designed and optimized (Figure S5). The successive binding energies were calculated to be 49.5, 64.0, 41.3, and 30.1 kJ/mol for NaOH(H2O)6·urea, NaOH(H2O)7·urea, NaOH(H2O)8·urea, and NaOH(H2O)9·urea, respectively. The results indicated that the optimized structure 10254

dx.doi.org/10.1021/jp501408e | J. Phys. Chem. B 2014, 118, 10250−10257

The Journal of Physical Chemistry B

Article

of NaOH(H2O)7·urea cluster appears as a perfect cage and is more stable than the others, which is due to that the NaOH(H2O)7·urea cluster possesses the biggest successive binding energy (Table S3). In this cage structure, the proton on the amino group of urea is easily transferred to the OH− ion, nearly forming a water molecule, and the C−N bond turns into partly double bond character, which is consistent with the 15N NMR results in this work. To explore the interaction between the hydroxyl groups in cellulose and the NaOH/urea/water clusters, cellobiose with tg conformation (shown in Figure S4) was used as the candidate of cellulose to suppose the possible structures. As for the conformation of cellobiose, there have been a great number of reports to investigate. French et al. have extensively investigated the structures and energies of cellobiose.44 However, in the case of conformations of cellulose in solvated environment, it is complex to confirm all the possible compound structures. In this work, the tg conformation was chose referred to the structure proposed by Nishiyama et al.5 Many researchers from various aspects have indicated that the hydroxyl groups in cellulose tend to interact with OH− ions in the NaOH solutions.39 On the basis of this point, we constructed two possible structures for NaOH/urea/H2O interacting with the hydroxyl groups in C-6 and C2, C3 positions of cellobiose. The optimized structures are shown in Figure 7. As depicted in Figure 7a, the distance between the

Cluster-O2,3 is weaker than that with water in NaOH(H2O)7· urea, which is consistent with the 23Na NMR results in this work. In Cell-Cluster-O2,3 cluster, the intramolecular hydrogen bond between H atom of hydroxyl groups in the C3 position and O5 atom of adjacent AHG unit is about 1.909 Å that are not disrupted when interacting with NaOH/urea/H2O cluster. It should be noted that the conformations shown in Figure 7 are just the typical representatives and not the only way by which the NaOH/urea/H2O cluster interacts with cellulose. This is due to that the minimization in the gas phase may not reflect accurately the complicated conformations of the cellulose in aqueous conditions.45,46 However, due to the complex situation in cellulose/NaOH/urea aqueous solutions and the considerable consuming of the computational time, our theoretical motivation in this work is only primarily for proposing possible interaction models for NaOH/urea/H2O with the hydroxyl groups at C6 and C2, C3 positions by using cellobiose in the tg conformation proposed by Nishiyama et al.5 The primary results are good agreement with the experimental data in this work. Therefore, no more attempts have been tried for determining the global minimized structures of cellulose/ NaOH/urea aqueous solution at the present state. Furthermore, the binding energy of NaOH(H2O)6·urea cluster with cellobiose to generate Cell-Cluster-O6 in Figure 7b is 120.1 kJ/ mol, which is higher than those of cellobiose with cellobiose (109.2 kJ/mol, in Figure S7) and NaOH(H2O)6·urea with water (81.4 kJ/mol, in Figure S5c). This result suggested that the NaOH/urea/water clusters prefer to attach to cellulose chains. Moreover, the diameter of the ICs estimated by quantum chemistry calculation is about 2.8 nm (Figure S8), which is comparable with the result obtained by the cryo-TEM. The strong hydrogen bonding in cellulose was mainly referred for the condensed structure and the dissolution of cellulose in the literature and so did in this work. However, on the basis of reviewing basic physicochemical aspects, Lindman et al. hypothesized that cellulose is significantly amphiphilic and that the low aqueous solubility must have a marked contribution from hydrophobic interactions,47 which was debated by cellulose scientists with a wide range of experiences representing a variety of scientific disciplines.48 Lindman’s hypothesis provided the general perspective of cellulose as a polymer in which intermolecular stress transfer involves more than hydrogen bonds. Hydrophobic and amphiphilic behaviors have been acknowledged for some time but may have been underestimated in conventional considerations of structure, solubility, etc.48 Lindman’s hypothesis and the following debates provided a new viewpoint for clarifying the fundamental problems of cellulose. The effects of the interactions besides the hydrogen in cellulose−NaOH−urea aqueous system will be investigated and discussed based on experimental results in the future.

Figure 7. Optimized structures of (a) Cell-Cluster-O2,3 and (b) CellCluster-O6 by the M06-2X/6-31+G(d) theoretical method.

hydroxyl groups at C2 and C3 is too close to bind two NaOH/ urea/H2O clusters simultaneously. Therefore, two H2O molecules in NaOH(H2O)7·urea cluster were replaced by the hydroxyl groups at C2 and C3 to form the complex structure, denoted as Cell-Cluster-O2,3. For the hydroxyl group at C-6 of cellulose, one H2O molecule in NaOH(H2O)7·urea cluster was replaced by the hydroxyl group at C-6 to result the complex structure, defined as Cell-Cluster-O6 (Figure 7b). In both cases, the cage structure of the primary NaOH(H2O)7·urea cluster are preserved. It is noticed that the conformation of cellulose moiety in the structure of Cell-Cluster-O6 still has the tg conformation with the χ (O6−C6−C5−O5 torsion angle) value of 167.61°. In addition, the distance of Na···O6 in CellCluster-O6 is about 3.929 Å, which is larger than that (2.378 Å) between Na+ and O atom of the corresponding H2O in NaOH(H2O)7·urea. In Cell-Cluster-O2,3, the distances of Na··· O2 and Na···O3 increase to 2.448 and 3.446 Å from 2.378 and 3.227 Å, respectively. These results indicate that the interaction of Na+ ions with hydroxyl groups in Cell-Cluster-O6 and Cell-



CONCLUSIONS The intermolecular interactions among NaOH, urea, and cellulose in aqueous system were investigated in detail and discussed. Temperature-dependent 23Na and 15N NMR results revealed that OH− ions of NaOH interacted directly with amino groups of urea through hydrogen bonds to form a cluster, and there was no direct interaction between urea and cellulose. Adding urea can effectively enhance the stability of the NaOH−cellulose complex as a result of the strong interaction between NaOH and urea, leading to low mobility of the NaOH−urea−cellulose complex. Low temperature can 10255

dx.doi.org/10.1021/jp501408e | J. Phys. Chem. B 2014, 118, 10250−10257

The Journal of Physical Chemistry B

Article

(2) Klemm, D.; Heublein, B.; Fink, H. P.; Bohn, A. Cellulose: Fascinating Biopolymer and Sustainable Raw Material. Angew. Chem., Int. Ed. 2005, 44, 3358−3393. (3) Klemm, D.; Kramer, F.; Moritz, S.; Lindström, T.; Ankerfors, M.; Gray, D.; Dorris, A. Nanocelluloses: A New Family of Nature-Based Materials. Angew. Chem., Int. Ed. 2011, 50, 5438−5466. (4) Kang, H. L.; Liu, R. G.; Huang, Y. Cellulose Derivatives and Graft Copolymers as Blocks for Functional Materials. Polym. Int. 2013, 62, 338−344. (5) Nishiyama, Y.; Langan, P.; Chanzy, H. Crystal Structure and Hydrogen-Bonding System in Cellulose Iβ from Synchrotron X-Ray and Neutron Fiber Diffraction. J. Am. Chem. Soc. 2002, 124, 9074− 9082. (6) Nishiyama, Y.; Sugiyama, J.; Chanzy, H.; Langan, P. Crystal Structure and Hydrogen Bonding System in Cellulose Iα from Synchrotron X-Ray and Neutron Fiber Diffraction. J. Am. Chem. Soc. 2003, 125, 14300−14306. (7) Fink, H.-P.; Weigel, P.; Purz, H.; Ganster, J. Structure Formation of Regenerated Cellulose Materials from NMMO-Solutions. Prog. Polym. Sci. 2001, 26, 1473−1524. (8) Liu, R. G.; Shen, Y. Y.; Shao, H. L.; Wu, C. X.; Hu, X. C. An Analysis of Lyocell Fiber Formation as a Melt-Spinning Process. Cellulose 2001, 8, 13−21. (9) Wang, H.; Gurau, G.; Rogers, R. D. Ionic Liquid Processing of Cellulose. Chem. Soc. Rev. 2012, 41, 1519−1537. (10) Zhang, H.; Wu, J.; Zhang, J.; He, J. 1-Allyl-3-Methylimidazolium Chloride Room Temperature Ionic Liquid: A New and Powerful Nonderivatizing Solvent for Cellulose. Macromolecules 2005, 38, 8272−8277. (11) Swatloski, R. P.; Spear, S. K.; Holbrey, J. D.; Rogers, R. D. Dissolution of Cellose with Ionic Liquids. J. Am. Chem. Soc. 2002, 124, 4974−4975. (12) McCormick, C. L.; Callais, P. A.; Hutchinson, B. H., Jr. Solution Studies of Cellulose in Lithium Chloride and N,N-Dimethylacetamide. Macromolecules 1985, 18, 2394−2401. (13) Morgenstern, B.; Kammer, H. W.; Berger, W.; Skrabal, P. 7Li NMR Study on Cellulose LiCl N,N-Dimethylacetamide Solutions. Acta Polym. 1992, 43, 356−357. (14) Zhang, C.; Liu, R. G.; Xiang, J. F.; Kang, H. L.; Liu, Z. J.; Huang, Y. Dissolution Mechanism of Cellulose in N,N-Dimethylacetamide/ Lithium Chloride: Revisiting through Molecular Interactions. J. Phys. Chem. B 2014, 118, 9507−9514. (15) Isogai, A.; Atalla, R. Dissolution of Cellulose in Aqueous NaOH Solutions. Cellulose 1998, 5, 309−319. (16) Cai, J.; Zhang, L.; Liu, S.; Liu, Y.; Xu, X.; Chen, X.; Chu, B.; Guo, X.; Xu, J.; Cheng, H. Dynamic Self-Assembly Induced Rapid Dissolution of Cellulose at Low Temperatures. Macromolecules 2008, 41, 9345−9351. (17) http://cell.sites.acs.org/. (18) Cai, J.; Zhang, L. N.; Zhou, J. P.; Qi, H. S.; Chen, H.; Kondo, T.; Chen, X. M.; Chu, B. Multifilament Fibers Based on Dissolution of Cellulose in NaOH/Urea Aqueous Solution: Structure and Properties. Adv. Mater. 2007, 19, 821−825. (19) Qi, H.; Chang, C.; Zhang, L. Properties and Applications of Biodegradable Transparent and Photoluminescent Cellulose Films Prepared via a Green Process. Green Chem. 2009, 11, 177−184. (20) Luo, X.; Liu, S.; Zhou, J.; Zhang, L. In Situ Synthesis of Fe3O4/ Cellulose Microspheres with Magnetic-Induced Protein Delivery. J. Mater. Chem. 2009, 19, 3538−3545. (21) Cai, J.; Liu, S.; Feng, J.; Kimura, S.; Wada, M.; Kuga, S.; Zhang, L. Cellulose−Silica Nanocomposite Aerogels by in Situ Formation of Silica in Cellulose Gel. Angew. Chem., Int. Ed. 2012, 51, 2076−2079. (22) Chang, C.; Peng, J.; Zhang, L.; Pang, D.-W. Strongly Fluorescent Hydrogels with Quantum Dots Embedded in Cellulose Matrices. J. Mater. Chem. 2009, 19, 7771−7776. (23) He, M.; Zhao, Y.; Duan, J.; Wang, Z.; Chen, Y.; Zhang, L. Fast Contact of Solid−Liquid Interface Created High Strength MultiLayered Cellulose Hydrogels with Controllable Size. ACS Appl. Mater. Interfaces 2014, 6, 1872−1878.

accelerate the breakdown of the intermolecular hydrogen bonds among cellulose and thus prevent the aggregation of the separated cellulose ICs. The stability of the ICs at the lower temperature is a key to the cellulose dissolution. TEM observation provided a direct evidence of the stiffness of cellulose chains and the formation of ICs in the solution. The ICs of cellulose in the solution have an extended wormlike morphology with a diameter of 3.6 ± 0.4 nm and a length about 300 nm estimated by cryo-TEM observation. The model theoretic calculation confirmed that the structure of the NaOH−urea−water cluster maintained when cellulose interacted with the cluster. Moreover, the O3H···O5 intramolecular hydrogen bonds of cellulose remained in the cellulose ICs. The diameter of the ICs estimated theoretically is about 2.8 nm, which was comparable with the data obtained by cryo-TEM. The combination of NMR and quantum chemistry calculation revealed that the key point of cellulose dissolution was cleavage of the intermolecular hydrogen bonds of cellulose as well as that the stabilizing of the NaOH−cellulose complex by the surrounding urea. Low temperature improved the formation and stability of small molecular cluster as well as complex associated with small molecules and cellulose chains through the hydrogen bonding. This work clarified the intermolecular interactions in cellulose/NaOH/urea aqueous solution, which is particularly important for promoting the applications of this green cellulose solvent system in commercial producing regenerated cellulose products.



ASSOCIATED CONTENT

* Supporting Information S

NMR spectra and data of cellulose solution, the detailed Cartesian coordinates of the optimized structures of NaOH/ urea/H2O clusters and the cellulose IC, total binding energies and successive binding energies of NaOH/urea/H2O clusters, and the electronic energies of the optimized geometries. This material is available free of charge via the Internet at http:// pubs.acs.org.



AUTHOR INFORMATION

Corresponding Authors

*E-mail [email protected] (R.L.). *E-mail [email protected] (L.Z.). Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by National Basic Research Program of China (973 Program, 2010CB32203), the Major Program of National Natural Science Foundation of China (21334005), the National Natural Science Foundation of China (21274154, 51203122, and 51373191), and the Fundamental Research Fund for the Central Universities (201220302020205). We also thank Dr. Xiaojun Huang and Dr. Gang Ji from Institute of Biophysics, Chinese Academy of Science, for technical support in cryo-TEM samples preparation and data collection.



REFERENCES

(1) Himmel, M. E.; Ding, S.-Y.; Johnson, D. K.; Adney, W. S.; Nimlos, M. R.; Brady, J. W.; Foust, T. D. Biomass Recalcitrance: Engineering Plants and Enzymes for Biofuels Production. Science 2007, 315, 804−807. 10256

dx.doi.org/10.1021/jp501408e | J. Phys. Chem. B 2014, 118, 10250−10257

The Journal of Physical Chemistry B

Article

(24) Jiang, Z.; Lu, A.; Zhou, J.; Zhang, L. Interaction between −OH Groups of Methylcellulose and Solvent in NaOH/Urea Aqueous System at Low Temperature. Cellulose 2012, 19, 671−678. (25) Horii, F.; Yamamoto, H.; Kitamaru, R.; Tanahashi, M.; Higuchi, T. Transformation of Native Cellulose Crystals Induced by Saturated Steam at High Temperatures. Macromolecules 1987, 20, 2946−2949. (26) Xiong, B.; Zhao, P. P.; Cai, P.; Zhang, L. N.; Hu, K.; Cheng, G. Z. NMR Spectroscopic Studies on the Mechanism of Cellulose Dissolution in Alkali Solutions. Cellulose 2013, 20, 613−621. (27) Yamashiki, T.; Kamide, K.; Okajima, K.; Kowsaka, K.; Matsui, T.; Fukase, H. Some Characteristic Features of Dilute Aqueous Alkali Solutions of Specific Alkali Concentration (2.5 mol l−1) Which Possess Maximum Solubility Power against Cellulose. Polym. J. 1988, 20, 447− 457. (28) Remsing, R. C.; Swatloski, R. P.; Rogers, R. D.; Moyna, G. Mechanism of Cellulose Dissolution in the Ionic Liquid 1-N-Butyl-3methylimidazolium Chloride: A 13C and 35/37Cl NMR Relaxation Study on Model Systems. Chem. Commun. 2006, 1271−1273. (29) Yamamoto, H.; Horii, F. CP/MAS 13C NMR Analysis of the Crystal Transformation Induced for Valonia Cellulose by Annealing at High Temperatures. Macromolecules 1993, 26, 1313−1317. (30) Coomber, A. T.; Beljonne, D.; Friend, R. H.; Bredas, J. L.; Charlton, A.; Robertson, N.; Underbill, A. E.; Kurmoo, M.; Day, P. Intermolecular Interactions in the Molecular Ferromagnetic NH4Ni(mnt)2·H2O. Nature 1996, 380, 144−146. (31) Vaidhyanathan, R.; Iremonger, S. S.; Shimizu, G. K. H.; Boyd, P. G.; Alavi, S.; Woo, T. K. Direct Observation and Quantification of CO2 Binding within an Amine-Functionalized Nanoporous Solid. Science 2010, 330, 650−653. (32) Cho, H.; Shepson, P. B.; Barrie, L. A.; Cowin, J. P.; Zaveri, R. NMR Investigation of the Quasi-Brine Layer in Ice/Brine Mixtures. J. Phys. Chem. B 2002, 106, 11226−11232. (33) Yavari, I.; Roberts, J. D. Nitrogen-15 Nuclear Magnetic Resonance Spectroscopy. A Simple Procedure for Determination of the Rates of N-H Proton Exchange of cis- and trans-1-Aza-2cyclononanone. J. Am. Chem. Soc. 1978, 100, 5217−5218. (34) Haushalter, K. A.; Lau, J.; Roberts, J. D. An NMR Investigation of the Effect of Hydrogen Bonding on the Rates of Rotation About the CN Bonds in Urea and Thiourea. J. Am. Chem. Soc. 1996, 118, 8891− 8896. (35) Bloor, E.; Kidd, R. Solvation of Sodium Ions Studied by 23Na Nuclear Magnetic Resonance. Can. J. Chem. 1968, 46, 3425−3430. (36) Saika, A.; Slichter, C. P. A Note on the Fluorine Resonance Shifts. J. Chem. Phys. 1954, 22, 26−28. (37) Templeman, G. J.; Van Geet, A. L. Sodium Magnetic Resonance of Aqueous Salt Solutions. J. Am. Chem. Soc. 1972, 94, 5578−5582. (38) Mason, J. Multinuclear NMR; Plenum Press: New York, 1987. (39) Isogai, A. NMR Analysis of Cellulose Dissolved in Aqueous NaOH Solutions. Cellulose 1997, 4, 99−107. (40) Zhao, Y.; Truhlar, D. G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215−241. (41) Zhao, Y.; Truhlar, D. G. Density Functionals with Broad Applicability in Chemistry. Acc. Chem. Res. 2008, 41, 157−167. (42) Parthasarathi, R.; Bellesia, G.; Chundawat, S. P. S.; Dale, B. E.; Langan, P.; Gnanakaran, S. Insights into Hydrogen Bonding and Stacking Interactions in Cellulose. J. Phys. Chem. A 2011, 115, 14191− 14202. (43) Frisch, M. J.; Trucks, G. W.; et al. Gaussian 03, RevisionD.02; Gaussian, Inc.: Pittsburgh, PA, 2003. (44) French, A. D.; Johnson, G. P. Roles of Starting Geometries in Quantum Mechanics Studies of Cellobiose. Mol. Simul. 2008, 34, 365−372. (45) French, A. D.; Csonka, G. I. Hydroxyl Orientations in Cellobiose and Other Polyhydroxyl Compounds: Modeling Versus Experiment. Cellulose 2011, 18, 897−909.

(46) French, A. D.; Johnson, G. P.; Cramer, C. J.; Csonka, G. I. Conformational Analysis of Cellobiose by Electronic Structure Theories. Carbohydr. Res. 2012, 350, 68−76. (47) Medronho, B.; Romano, A.; Miguel, M. G.; Stigsson, L.; Lindman, B. Rationalizing Cellulose (in)Solubility: Reviewing Basic Physicochemical Aspects and Role of Hydrophobic Interactions. Cellulose 2012, 19, 581−587. (48) Glasser, W. G.; Atalla, R. H.; Blackwell, J.; Brown, R. M.; Burchard, W.; French, A. D.; Klemm, D. O.; Nishiyama, Y. About the Structure of Cellulose: Debating the Lindman Hypothesis. Cellulose 2012, 19, 589−598.

10257

dx.doi.org/10.1021/jp501408e | J. Phys. Chem. B 2014, 118, 10250−10257