Dissolution of cellulose in ionic liquids: an ab initio molecular dynamics simulation study

Optimization of Dry-Jet Wet Spinning of Regenerated Cellulose Fibers Using [mTBDH][OAc] as a Solvent

Superbase-based ionic liquids (ILs) have demonstrated excellent dissolution capability for cellulose, and employing the dry-jet wet spinning process, high-tenacity regenerated textile fibers have been made. Among a range of superbase-based ILs, [mTBDH][OAc] exhibited not only good spinnability but also exceptional recyclability, making it highly suitable for a closed-loop production of regenerated cellulose fibers. To further optimize the spinning process, we investigated the influence of the cellulosic raw materials and the IL with residual water on spinnability and fiber properties. In addition, single-filament spinning and multifilament spinning using spinnerets with different hole densities were investigated to reveal the upscaling challenges of the dry-jet wet spinning process. The air gap conditions, for example, temperature and moisture concentration were simulated using COMSOL multiphysics. The results indicate that the presence of a small amount of water (3 wt%) in the IL has a positive effect on spinnability, while the mechanical properties of the fibers remain unchanged.

Understanding X-ray Photoelectron Spectra of Ionic Liquids: Experiments and Simulations of 1-Butyl-3-methylimidazolium Thiocyanate

We demonstrate a combined experimental and computational approach to probe the electronic structure and atomic environment of an ionic liquid, based on core level binding energies. The 1-butyl-3-methylimidazolium thiocyanate [C₄C₁Im][SCN] ionic liquid was studied using ab initio molecular dynamics, and results were compared against previously published and new experimental X-ray photoelectron spectroscopy (XPS) data. The long-held assumption that initial-state effects in XPS dominate the measured binding energies is proven correct, which validates the established premise that the ground state electronic structure of the ionic liquid can be inferred directly from XPS measurements. A regression model based upon site electrostatic potentials and intramolecular bond lengths is shown to account accurately for variations in core-level binding energies within the ionic liquid, demonstrating the important effect of long-range interactions on the core levels and throwing into question the validity of traditional single ion pair ionic liquid calculations for interpreting XPS data.

Molecular simulations explain the exceptional thermal stability, solvent tolerance and solubility of protein-polymer surfactant bioconjugates in ionic liquids

Proteins complexed electrostatically with polymer surfactants constitute a viscous liquid by themselves, called the solvent-free protein liquid (SFPL). A solution of SFPL in a room temperature ionic liquid (PS-IL) offers the protein hyperthermal stability, higher solubility and greater IL tolerance. A generic understanding of these protein-polymer systems is obtained herein through extensive atomistic molecular dynamics simulations of three different enzymes (lipase A, lysozyme and myoglobin) under various conditions. Along with increased intra-protein hydrogen bonding, the surfactant coating around the proteins imparts greater thermal stability, and also aids in screening protein-IL interactions, endowing them IL tolerance. The reduced surface polarity of the protein-polymer bioconjugate and hydrogen bonding between the ethylene glycol groups of the surfactant and the IL cation contribute to the facile solvation of the protein in its PS-IL form. The results presented here rationalize several experimental observations and will aid in the improved design of such hybrid materials for sustainable catalysis.

Influence of Calcium Silicate and Hydrophobic Agent Coatings on Thermal, Water Barrier, Mechanical and Biodegradation Properties of Cellulose

Thin films of cellulose and cellulose–CaSiO₃ composites were prepared using 1-ethyl-3-methylimidazolium acetate (EMIMAc) as the dissolution medium and the composites were regenerated from an anti-solvent. The surface hydrophilicity of the resultant cellulose composites was lowered by coating them with three different hydrophobizing agents, specifically, trichloro(octadecyl)silane (TOS), ethyl 2-cyanoacrylate (E2CA) and octadecylphosphonic acid (ODPA), using a simple dip-coating technique. The prepared materials were subjected to flame retardancy, water barrier, thermal, mechanical and biodegradation properties analyses. The addition of CaSiO₃ into the cellulose increased the degradation temperature and flame retardant properties of the cellulose. The water barrier property of cellulose–CaSiO₃ composites under long term water exposure completely depends on the nature of the hydrophobic agents used for the surface modification process. All of the cellulose composites behaved mechanically as a pure elastic material with a glassy state from room temperature to 250 °C, and from 20% to 70% relative humidity (RH). The presence of the CaSiO₃ filler had no effect on the elastic modulus, but it seemed to increase after the TOS surface treatment. Biodegradability of the cellulose was evaluated by enzyme treatments and the influence of CaSiO₃ and hydrophobic agents was also derived.

Glucose in dry and moist ionic liquid: vibrational circular dichroism, IR, and possible mechanisms

Ionic liquids and their mixtures with water show remarkable features in cellulose processing. For this reason, understanding the behavior of carbohydrates in ionic liquids is important. In the present study, we investigated three d-glucose isomers (α, β and open-chain) in 1-ethyl-3-methylimidazolium acetate in the presence and absence of water, through ab initio molecular dynamics simulations. In the complex hydrogen bonding network of these mixtures, the most interesting observation is that upon water addition every hydrogen bond elongates, except the glucose-glucose hydrogen bond for the open-chain and the α-form which shortens, clearly showing the beginning of the crystallization process. The ring glucose rearranges from on-top to in-plane and the open form changes from a coiled to a more linear arrangement when adding water which explains the contradiction that the center of mass distances of the glucose molecules with other glucose molecules grow while the hydrogen bonds shorten. The appearance of coiled open forms indicates that the previously suggested isomerization between these forms is possible and might play a role in the solubility of the related carbohydrates. The calculated IR and VCD spectra reveal insight into the intermolecular interactions, with good to excellent agreements with experimental spectra. Investigating the role of the cation, distances between the acidic carbon atom of the cation and the glucose carbon atom where ring closure and opening occurs are found, which are way shorter than dispersion-like interactions between aliphatic hydrocarbons.

Combining random walk and regression models to understand solvation in multi-component solvent systems

Polysaccharides, such as cellulose, are often processed by dissolution in solvent mixtures, e.g. an ionic liquid (IL) combined with a dipolar aprotic co-solvent (CS) that the polymer does not dissolve in. A multi-walker, discrete-time, discrete-space 1-dimensional random walk can be applied to model solvation of a polymer in a multi-component solvent mixture. The number of IL pairs in a solvent mixture and the number of solvent shells formable, x, is associated with n, the model time-step, and N, the number of random walkers. The mean number of distinct sites visited is proportional to the amount of polymer soluble in a solution. By also fitting a polynomial regression model to the data, we can associate the random walk terms with chemical interactions between components and probe where the system deviates from a 1-D random walk. The 'frustration' between solvents shells is given as ln x in the random walk model and as a negative IL:IL interaction term in the regression model. This frustration appears in regime II of the random walk model (high volume fractions of IL) where walkers interfere with each other, and the system tends to its limiting behaviour. In the low concentration regime, (regime I) the solvent shells do not interact, and the system depends only on IL and CS terms. In both models (and both regimes), the system is almost entirely controlled by the volume available to solvation shells, and thus is a counting/space-filling problem, where the molar volume of the CS is important. Small deviations are observed when there is an IL-CS interaction. The use of two models, built on separate approaches, confirm these findings, demonstrating that this is a real effect and offering a route to identifying such systems. Specifically, the majority of CSs - such as dimethylformide - follow the random walk model, whilst 1-methylimidazole, dimethyl sulfoxide, 1,3-dimethyl-2-imidazolidinone and tetramethylurea offer a CS-mediated improvement and propylene carbonate results in a CS-mediated hindrance. It is shown here that systems, which are very complex at a molecular level, may, nonetheless, be effectively modelled as a simple random walk in phase-space. The 1-D random walk model allows prediction of the ability of solvent mixtures to dissolve cellulose based on only two dissolution measurements (one in neat IL) and molar volume.

Accurate prediction of energetic properties of ionic liquid clusters using a fragment-based quantum mechanical method

Accurate prediction of physicochemical properties of ionic liquids (ILs) is of great significance to understand and design novel ILs with unique properties.

Ion-water wires in imidazolium-based ionic liquid/water solutions induce unique trends in density

Ionic liquid/water binary mixtures are rapidly gaining popularity as solvents for dissolution of cellulose, nucleobases, and other poorly water-soluble biomolecules. Hence, several studies have focused on measuring the thermophysical properties of these versatile mixtures. Among these, 1-ethyl-3-methylimidazolium ([emim]) cation-based ILs containing different anions exhibit unique density behaviours upon addition of water. While [emim][acetate]/water binary mixtures display an unusual rise in density with the addition of low-to-moderate amounts of water, those containing the [trifluoroacetate] ([Tfa]) anion display a sluggish decrease in density. The density of [emim][tetrafluoroborate] ([emim][BF4])/water mixtures, on the other hand, declines rapidly in close accordance with the experimental reports. Here, we unravel the structural basis underlying this unique density behavior of [emim]-based IL/water mixtures using all-atom molecular dynamics (MD) simulations. The results revealed that the distinct nature of anion-water hydrogen bonded networks in the three systems was a key in modulating the observed unique density behaviour. Vast expanses of uninterrupted anion-water-anion H-bonded stretches, denoted here as anion-water wires, induced significant structuring in [emim][Ac]/water mixtures that resulted in the density rise. Conversely, the presence of intermittent large water clusters disintegrated the anion-water wires in [emim][Tfa]/water and [emim][BF4]/water mixtures to cause a monotonic density decrease. The differential nanostructuring affected the dynamics of the solutions proportionately, with the H-bond making and breaking dynamics found to be greatly retarded in [emim][Ac]/water mixtures, while it exhibited a faster relaxation in the other two binary solutions.

The effect of pressure on cation-cellulose interactions in cellulose/ionic liquid mixtures

Cation-cellulose interactions in binary mixtures of [EMIM][OAc] and cellulose have been investigated using high-pressure infrared spectroscopy. At low concentrations of cellulose, almost no changes were observed in the imidazolium C(2)-H frequency; on the other hand, at high concentrations of cellulose, increases in the C(2)-H vibration frequency were observed under ambient pressure. As the pressure was elevated, the imidazolium C(2)-H absorption of the [EMIM][OAc]/cellulose mixtures underwent band-narrowing and blue-shifts in the frequency. These observations suggest that high pressures may strengthen the hydrogen bonds formed between C(2)-H and cellulose, possibly forcing the cellulose to dissociate clusters of ionic liquid through enhanced cation-cellulose interactions. In contrast to the cation-cellulose interaction results, the COO(-) absorption of the anion does not show dramatic changes under high pressures. Our results indicate the possibility of enhanced cation-cellulose interactions through pressure elevation, demonstrating that high pressures may have the potential to tune the relative contributions of cation-cellulose and anion-cellulose interactions in cellulose/ionic liquid mixtures.

Chemocatalytic hydrolysis of cellulose at 37 °C, 1 atm

The metal salt – Brønsted acidic ionic liquid system composed of ZnCl₂·1.74H₂O-1-(1-propylsulfonic)-3-methylimidazolium chloride can directly hydrolyze untreated cellulose in 78% total reducing sugar and 19% glucose yield at 37 °C, 1 atm in 4.0 days.

Solution Properties of Hemicellulose Polysaccharides with Four Common Carbohydrate Force Fields

Hemicellulose polysaccharides play an important role in the swelling behavior of the primary plant cell wall, and molecular dynamics simulations provide the means of gaining a concise understanding of the interactions of hemicellulose polysaccharides with water. Here, we compare four of the main polysaccharide force fields (CHARMM36 TIP3P, GROMOS56A6(CARBO) SPC, GLYCAM06h TIP3P, and GLYCAM06h TIP5P) for the most abundant hemicellulose backbone components. In particular, we compare aggregation, diffusion coefficients, system density, and investigate the free energy of hydration of saccharides in water. We find that the saccharides show nonphysical aggregation at low concentrations with the GLYCAM06h TIP3P force field, which can be rectified by the use of the TIP5P water model. As a result of the aggregation, GLYCAM06h TIP3P does not lead to reasonable diffusion coefficients whereas the diffusion coefficients, as well as the system density, agrees best with experimental data for the GLYCAM06h TIP5P force field. Overall, GLYCAM06h TIP5P gives good agreement with experimental free energy of hydration data for small saccharides. In addition, the free energy of hydration for short polysaccharides calculated with the GLYCAM06h TIP5P force field is consistent with the radial distribution functions between the polysaccharides and water, the hydration number of the polysaccharides, and the hydrogen bonds formed in the system.

Dissolution of cellulose in room temperature ionic liquids: anion dependence

The dissolution of cellulosic biomass in room temperature ionic liquids (RTILs) is studied through free energy calculations of its monomer, viz., cellobiose, within a molecular dynamics simulation approach. The solvation free energy (SFE) of cellobiose in ionic liquids containing any of seven different anions has been calculated. The ranking of these liquids based on SFE compares well with experimental data on the solubility of cellulose. The dissolution is shown to be enthalpically dominated, which is correlated with the strength of intermolecular hydrogen bonding between cellobiose and the anions of the IL. Large entropic changes upon solvation in [CF3SO3](-) and [OAc](-) based ionic liquids have been explained in terms of the solvent-aided conformational flexibility of cellobiose.

Mechanisms of hydrogen bond formation between ionic liquids and cellulose and the influence of water content

We explore the complex network of transitions occurring between different hydrogen bonding states within ionic liquids and cellulose.

Multiscale modeling of the trihexyltetradecylphosphonium chloride ionic liquid

Hierarchical trihexyltetradecylphosphonium cationic and chloride anionic models.