Next Article in Journal
Sustainable Plant Growth Promotion and Chemical Composition of Pyroligneous Acid When Applied with Biochar as a Soil Amendment
Next Article in Special Issue
Prediction of pH Value of Aqueous Acidic and Basic Deep Eutectic Solvent Using COSMO-RS σ Profiles’ Molecular Descriptors
Previous Article in Journal
Essential Oil Content and Compositional Variability of Lavandula Species Cultivated in the Mid Hill Conditions of the Western Himalaya
Previous Article in Special Issue
Enhanced Extraction Efficiency of Flavonoids from Pyrus ussuriensis Leaves with Deep Eutectic Solvents
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Neutron Total Scattering Investigation of the Dissolution Mechanism of Trehalose in Alkali/Urea Aqueous Solution

1
Institute of High Energy Physics, Chinese Academy of Sciences (CAS), Beijing 100049, China
2
Spallation Neutron Source Science Center, Dongguan 523803, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
4
STFC ISIS Facility, Rutherford Appleton Laboratory, Didcot OX11 0QX, UK
5
Institute for Advanced Study, Shenzhen University, Shenzhen 508060, China
*
Author to whom correspondence should be addressed.
Molecules 2022, 27(11), 3395; https://doi.org/10.3390/molecules27113395
Submission received: 28 April 2022 / Revised: 20 May 2022 / Accepted: 23 May 2022 / Published: 25 May 2022
(This article belongs to the Special Issue Recent Advances in Green Solvents)

Abstract

:
The atomic picture of cellulose dissolution in alkali/urea aqueous solution is still not clear. To reveal it, we use trehalose as the model molecule and total scattering as the main tool. Three kinds of alkali solution, i.e., LiOH, NaOH and KOH are compared. The most probable all-atom structures of the solution are thus obtained. The hydration shell of trehalose has a layered structure. The smaller alkali ions can penetrate into the glucose rings around oxygen atoms to form the first hydration layer. The larger urea molecules interact with hydroxide groups to form complexations. Then, the electronegative complexation can form the second hydration layer around alkali ions via electrostatic interaction. Therefore, the solubility of alkali aqueous solution for cellulose decreases with the alkali cation radius, i.e., LiOH > NaOH > KOH. Our findings are helpful for designing better green solvents for cellulose.

1. Introduction

Cellulose is the world’s most produced natural polymer, and it is a potential candidate to replace petroleum-based materials. To achieve high performance, cellulose has to be dissolved first [1].
Cellulose is difficult to dissolve. Strong inter- and intra-chain interactions prevent its structure from being deconstructed. In industry, the widely used viscose method produces a large amount of alkaline and acidic waste, carbon disulfide and hydrogen sulfide gases. It pollutes the environment. N-methylmorpholine noxide (NMMO) and ionic liquids (ILS) are environment friendly, but they have high costs [1].
Around 2000, Prof. Zhang proposed the use of precooled NaOH/LiOH urea aqueous solution to dissolve cellulose [1,2]. This solvent has the advantages of fast dissolution, low cost and low pollution and, thus, has good application prospects. Unfortunately, the solubility of cellulose in green solvent is still too low to reach the requirements of industrial production. We need to know the dissolution mechanism first to increase its solubility. Lots of methods, such as NMR, FTIR, DSC, TEM, et al., have been used. These studies qualitatively showed that cationic hydrates are more easily adsorbed around cellulose molecules to form a new, stable hydrogen bond network at low temperatures, and hydrates of urea molecules form a sheath-like inclusion complex (IC) around their periphery [3,4,5]. Bjorn et al. summarized the dissolution mechanism of cellulose and proposed that cellulose is amphiphilic and that hydrophobic interactions are important for its solubility [6,7,8,9]. Wolfgang et al. agreed that hydrophobic interactions are the driving forces in an amorphous system, but all types of cellulose are highly crystalline so hydrogen bonds must be broken before hydrophobic interactions can be effective, and this requires a strong alkaline medium [10]. To reveal the dissolution mechanism of cellulose in alkaline/urea aqueous solution, a detailed atomic picture needs to be presented.
The combination of neutron total scattering and empirical potential structure refinement (EPSR) can observe in situ the most probable all-atom structure of the liquids [11,12]. In a previous study, we used the combined methods to study the atomic structure of trehalose in NaOH/urea aqueous solution [13]. We used trehalose as the model molecule for cellulose, because it has similar glucose rings and is one of the disaccharide molecules that has no reducibility (other binary sugars, such as cellobiose, glucose, maltose and lactose, are oxidized in the alkaline solution during neutron total scattering experiments). We found that NaOH, urea and water work cooperatively to dissolve trehalose. Na+ accumulates around electronegative oxygen atoms in the hydration shell, while urea molecules only participate in the dissolution process via Na+ bridging. Additionally, we predicted that alkali with smaller ions, such as LiOH, have better solubility for cellulose.
To prove this, we further observed the microscopic dissolution pictures in two different alkali/urea aqueous solutions, i.e., LiOH and KOH. Then, we compared them with previous results in NaOH/urea aqueous solution. The hydration shell of trehalose has a layered structure; cations directly interact with the glucose rings to form the first hydration layer. It destroys their intra- and intermolecular hydrogen bonds. Thus, the smaller the radius of the cation, the easier it approaches the inside of the glucose ring. From K+ to Na+ to Li+, its ability to dissolve cellulose gradually increases. Urea does not directly interact with glucose rings, and it forms strong complexations with hydroxide groups. The urea hydration complexation forms the second hydration layer via electric interaction. It prevents it from re-aggregating. Temperature effect was also investigated. The atomic structure of the solution did not change when it was cooled down to −10 °C (or −5 °C). Taking into account the fact that cellulose only dissolves in green solvent at lower temperatures, the dissolution had to be a dynamic process.

2. Theory and Methods

2.1. Neutron Scattering Method and SANDALS

The total neutron scattering experiment is an important method for the study of the atomic structure of liquids. In neutron scattering, the observed neutron structure factor (F(Q)) and the atomic structure of the sample have the relationship:
F ( Q ) = α = 1 N c α b α 2 + α = 1 , β α ( 2 δ α β ) c α c β b α b β { 4 π ρ 0 r 2 ( g α β ( r ) 1 ) sin ( Q r ) Q r d r }  
where Q and r are the momentum transfer and the distance between the two atoms in the sample; ci and bi (i = α, β) are the quantity ratio and neutron scattering length of i species; δαβ is the Kronecker δ function; ρ is the atomic number density of the sample; and gαβ (r) is the radial distribution function (hereinafter referred to RDF), which reflects the microscopic atomic structure of amorphous matter.
g α β ( r ) = n α β 4 π r 2   d r ρ β  
where nαβ is an average number of β atoms around an α atom contained in a spherical shell with r radius and dr thickness; ρβ is the average number density of β atoms in the sample; and gαβ(r) describes the number density change of the β atom as the function of the distance from the α atom, which reflects the interaction between atoms α and β. We used g(1)αβ(r) to mark the first peak of gαβ(r). It represents the closest, most probable distance between atoms β and α. The peak height indicates the ratio of the atomic number density of atom β to its average value at this position.
The neutron total scattering experiments were performed at SANDALS in ISIS (Didcot, UK) and Multi-Physics Instrument (MPI) in CSNS (Dongguan, China). The scattering vector (Q) range was from 0.1 Å−1 to 50 Å−1. This means that the instrument could only measure the micro-structure of samples from ~0.1 Å to ~30 Å. Therefore, we chose trehalose as the model molecule for cellulose. The size of trehalose is ~11.6 Å; we could, thus, use a neutron total scattering instrument to observe its microscopic atomic structure in alkali/urea aqueous solution.

2.2. EPSR Simulations

The empirical potential structure refinement (EPSR) is a program developed by Prof. Soper to explore the most probable, all-atom structure of an experimental sample based on neutron scattering [11]. EPSR is essentially a Monte Carlo simulation program. In EPSR, there are two kinds of potential energy: reference potential and experimental potential. The reference potential is taken from the molecular dynamics simulation force field, which is used to realize the basic structural constraints of the simulation system, such as molecular structure, the minimum distance between atoms, etc.; the experimental potential is the structural data observed in the neutron scattering experiment [11,12]. After the experiment, we used EPSR simulation to reconstruct the atomic structure of the experimental sample.

2.3. Experiment Samples

The trehalose, LiOH, KOH and urea used in the experiment were purchased from Shanghai Aladdin Reagent Company. The purity of trehalose was over 99%, and the purity of alkali and urea was over 99.9% for both. Both deuterated urea and heavy water were purchased from Sigma-Aldrich China. The deuteration rate of heavy water was over 99.9%, and the deuteration rate of urea was over 98%.
In this study, we used SANDALS and MPI to observe 4 kinds of experimental sample with different chemical components. They were aqueous solutions of trehalose/LiOH/urea, trehalose/KOH/urea, trehalose/LiOH and trehalose/KOH. The molar ratio of each component in the sample was consistent with that in Prof. Zhang’s experiment [2]. Each kind of solution included three different deuterated ratio test samples, i.e., full hydrogen, full deuterium and half deuterium. All samples were measured at room temperature (25 °C) and −10 °C. However, the trehalose LiOH/urea aqueous solution froze at −10 °C in the neutron scattering experiment, so we raised the temperature of this sample to −5 °C. Each sample was measured in a flat sample cell of titanium–zirconium alloy with a capacity of 1.3 mL for approximately 6 h. Empty sample cells were measured for approximately 4 h for background subtraction. The sample composition and symbols are shown in Table 1.
For the EPSR simulation and the following discussion of the results, the symbols of the atoms in trehalose were as shown in Figure 1. Oxygen and hydrogen atoms of H2O were labeled as OW and HW, respectively; alkali cations were labeled as Li, Na and K, respectively; the oxygen and hydrogen atom of hydroxide anion were labeled as OOH and HOH, respectively; and the carbon, oxygen, nitrogen and hydrogen atoms in the urea molecule were labeled as CU, OU, NU and HU, respectively. The label of each atom and the parameters of reference potential used in EPSR simulation are listed in Table 2(1,2). Among them, epsilon and sigma are the parameters of the Lennard-Jones reference potential (L-J), and q is the charge of the atoms [13,14,15,16,17]. The information for the EPSR simulation box is shown in Table 3. The data about NaOH aqueous solutions can be found in our previous study [13].

3. Results

The neutron scattering structure factor and EPSR simulation results of different samples are shown in Figure 2. As shown in the figure, the EPSR simulation results and the neutron scattering profiles were in good agreement. Therefore, the structure reconstructed by EPSR represented the most probable atomic structure of the solutions.
EPSR gave us the most probable, all-atom structures in trehalose alkali/urea aqueous solution. The simulation, as shown in Figure 3, qualitatively showed the distribution diagram of Li+, Na+, K+ ions and urea molecules around trehalose. In order to improve the statistics, the ions and atoms in the figure were the number of atoms gathered in the 100-frame EPSR simulation conformation. Three things could be seen directly. The first thing was that those ions and urea molecules were anisotropically distributed around trehalose because of steric repulsion, so we built a coordinate system around its O1 atom (Figure 1). Thus, RDF could be calculated at different space angles. The second thing was that ions concentrated around the hydrophilic oxygen atoms of trehalose, and smaller ions were much more easily penetrated into the glucose ring, from K+ to Na+ to Li+ (Figure 3a). The third thing was that larger urea molecules could not directly interact with trehalose, and they were further away from trehalose than the ions (Figure 3b).
Then, we quantitatively analyzed the distribution of alkali, urea and water, as shown in Figure 3. We paid specific attention to the hydration shell around the oxygen atoms (we call them Os hereafter) of the glucose rings. They included all of the hydrophilic atoms on the glucose rings.
Because Os are electronegative, the first hydration shell had to be electropositive. There are three kinds of electropositive atom in an alkaline system, i.e., hydrogen atoms of urea and water and cations of alkali. Thus, there were 3 × 4 = 12 kinds of g(r) between them and the Os. There was only one g(1)(r) larger than 1, i.e., it was between cation and the Os. From the previous study, we knew that cations play a role in breaking the intra- and inter-molecular hydrogen bonds of cellulose. So, we first examined the RDF distribution between the cations and oxygen atoms of trehalose. Figure 4 compares the RDFs of Li+/Na+/K+ ions and O1~O4. Here, we separated Os into two groups. O1 and O2 were in one group, while O3 and O4 were in the other group. We discussed their RDFs with cations, respectively. In group one, both O1 and O2 were inside the glucose rings. g(1)O1-K(r)~5, centered at r1~2.6 Å. This means that the first K+ shell around O1 was centered about 2.6 Å, where the K+ concentration was five times larger than the average K+ concentration in the solution. g(1)O1-Na(r)~10, located at r1~2.4 Å, and g(1)O1-Li(r)~20 was at r1~2.0 Å. Thus, smaller cations were more close to O1, and their concentrations were 10 and 20 times their averaged concentration in solution. Similar things happened around O2. In group two, both O3 and O4 were in the periphery of the glucose rings. The spaces around them were, thereby, relatively open. g(1)O-K(r)~1.8, centered at r1~2.8 Å; g(1)O-Na(r)~1.9, centered at r1~2.45 Å; and g(1)O3-Li(r)~1.4, centered at r1~2.1 Å. Larger ions more easily concentrated around them.
From these comparisons, it was found that the Li+ ions were closest to the two oxygen atoms inside the glucose ring, i.e., O1 and O2, while larger ions, such as Na+ and K+, were more likely to be enriched near O4. Because the LiOH urea aqueous solution had the most powerful solubility, these comparisons confirmed that the cations in the alkaline solution directly interact with the glucose ring. An interesting result was that K+ ions with a relatively larger ionic radius were more concentrated around oxygen atoms away from the center of the glucose ring. This may explain why the KOH/urea aqueous solution had a better dissolving effect on chitin [18]. We suppose that, in the dissolution process of chitin, breaking the association between the side groups may be crucial, while, in the dissolution process of cellulose, breaking the interaction between the glucose backbones is a prerequisite.
Then, we looked for the components in the second layer of the hydration shell. The first hydration shell was electropositive, so the second shell had to be electronegative. There were four electronegative solvent atoms in the system, i.e., OW, NU, OU and OOH. They themselves or their complexations composed the second hydration layer. We listed all of their g(r) with the Os and compared their structures with gO1-Li(r) in LiOH/urea aqueous solution. Figure 5 assumes that they were isotropically distributed. We added two dash lines to indicate their possible boundary. It is easy to see that all of those g(r), except g(1)OW-O3 /O4(r), were smaller than 1. Therefore, three conclusions can be made, i.e., the second hydration layer almost did not exist; it was anisotropically distributed, or it only existed around O3 and O4 on the periphery of the glucose ring. We checked them one by one.
As shown in Figure 3, the distributions of cations and urea around O1 and O2 were very asymmetric. To quantitatively study the asymmetry of cationic and electronegative atoms around O1, we calculated their g(r) as a function of spatial angle. In the calculation, the setting of the coordinate axis was as shown in Figure 1. We first set a cone with an apex angle of 60° and the Y coordinate axis as the axis of symmetry. Then, we placed the apex of this cone at O1. Finally, we let the cone take O1 as the center of rotation and rotated it around the X-axis and calculated the g(r) changes with the rotation angle. The calculation results showed that g(1)O1-Li(r) reached its maximum value between 90° and 180°, which coincided with the Z-axis. Figure 6 shows gO1-Li(r), gO1-OW(r), gO1-OU(r), gO1-NU(r) and gO1-OOH(r) along the Z-axis. g(1)O1-OW(r), g(1)O1-OOH(r) and g(1)O1-NU(r) were larger than 1.0. Therefore, electronegative OW, OOH, NU or their complexation composed the second hydration layer.
The existence of urea is crucial for dissolving cellulose in green solvent. Previous NMR observations claimed that the amino group of urea forms an electronegative complexation with OH through hydrogen bonding. Such complexations wrap around cations and eventually form cellulose–NaOH–urea–H2O inclusion complexes (ICs) [5,19]. To clarify this, we listed all of the possible complexations between OH and urea (Figure 7a). There were four possibilities, i.e., HU as the proton donor to form a hydrogen bond with OW (HB1) and OOH (HB2), NU as the proton acceptor (HB3) and OU as the proton acceptor (HB4). Figure 7b gives the RDFs between OW and NU, OW and HU and HW and NU. We used 2.9 ± 0.3Å as the donor–acceptor distance constraint and linear bond (±20°) as the angle constraint to define the hydrogen bond. The possibility of forming a hydrogen bond for HB1 and HB3 was 78.5% and 10.9%, respectively. Thus, HU prefers forming a hydrogen bond with OW. Figure 7c shows the RDFs between OOH and NU, OOH and HU and HOH and NU. If we used 2.5 ± 0.4 Å as the donor–acceptor distance constraint and (±20°) as the angle constrain to define the hydrogen bond, the possibility of forming a hydrogen bond for HB2 was 51.9%. Thus, HU can also form a hydrogen bond with OOH. Finally, Figure 7d is the RDF between OW and OU and HW and OU. They could form a hydrogen bond (HB4), but they were away from the second hydration layer, as shown in Figure 6. Therefore, as a proton donor, the amino group of urea can hydrogen bond with hydroxy group to form complexations. The electronegative complexation forms the second hydration layer.
Temperature effect on dissolution was also investigated. Figure 8 shows the neutron scattering profiles of trehalose in LiOH urea aqueous solutions at 25 °C and −5 °C and their differences. It can be seen that there was no obvious difference between the two different temperatures. Therefore, low temperature does not change the atomic structure of the solutions. It only increases the stability of the solution, making it difficult for the dissolved cellulose to re-aggregate.

4. Conclusions

In this research work, we used trehalose as a model molecule and neutron scattering and EPSR simulation as the main tools to study the rapid dissolution mechanism of cellulose in alkali/urea aqueous solution. The three-dimensional atomic structures of trehalose in three different alkali/urea aqueous solutions were thus compared. Alkali, urea and water work cooperatively to dissolve trehalose. A layered hydration shell is crucial. Cations directly interact with the Os of the glucose rings. They first break the inter- and intra-hydrogen bonding. The smaller the ions, the more easily they penetrate into the glucose rings. Urea molecules are too large to approach the glucose ring. As proton donor, their amino group can form a hydrogen bond with the hydroxyl group. The resultant electronegative complexation constitutes the second hydration layer via bridging force. They further stabilize trehalose and prevent it from re-aggregating.
We also found an interesting phenomenon, that is, Li+ ions are more concentrated around O1 and O2, while K+ ions are more concentrated in the vicinity of O3 and O4, especially near O4. Although the low-temperature aqueous solution of KOH urea could not dissolve cellulose, it could dissolve chitin; at the same time, the low-temperature aqueous solution of LiOH/NaOH urea could dissolve cellulose, but could not dissolve chitin [18]. Based on this experimental evidence, we speculate: in the dissolution of cellulose, breaking the hydrogen bond formed by O1 and O2 is the key factor, while, in the dissolution of chitin, breaking the hydrogen bonds formed by its hydroxyl and amide groups is crucial. These findings are expected to be verified in future experiments.

Author Contributions

Conceptualization, H.C. (He Cheng) and C.M.; methodology, H.C. (He Cheng); software, C.M.; validation, T.Z., Z.H. and Y.L.; formal analysis, C.M.; investigation, H.C. (He Cheng); resources, W.Y.; data curation S.G. and H.C. (Huaican Chen); writing—original draft preparation, C.M.; writing—review and editing, H.C (He Cheng). and C.C.H.; visualization, C.M.; supervision, H.C. (He Cheng); project administration, H.C. (He Cheng); funding acquisition, H.C. (He Cheng) All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (NSFC) (grant numbers U1932161, U1830205 and 11805035); the National High Energy Physics Data Center of China (grant number NHEPSDC-OP-2021-001); and the National Key Research and Development Program of China (grant number 2017YFA-0403703).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are openly available in ISIS Neutron and Muon Source Data at http://doi.org/10.5286/ISIS.E.RB1820090 (accessed on 27 November 2018) and http://doi.org/10.5286/ISIS.E.RB1920190 (accessed on 15 December 2019).

Acknowledgments

The experiments were conducted using SANDALS at the ISIS Neutron, Muon Source and NOVA of J-PARC and the Multi-Physics Instrument (MPI) of CSNS. The authors thank the staff of the ISIS Disordered Materials Group, Materials and Life Science Experimental Facility of J-PARC and the instrument scientists of SANDALS and NOVA for the support with experiment design, operation and data analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, S.; Lu, A.; Zhang, L. Recent advances in regenerated cellulose materials. Prog. Polym. Sci. 2016, 53, 169–206. [Google Scholar] [CrossRef]
  2. Cai, J.; Zhang, L. Rapid Dissolution of Cellulose in LiOH/Urea and NaOH/Urea Aqueous Solutions. Macromol. Biosci. 2005, 5, 539–548. [Google Scholar] [CrossRef] [PubMed]
  3. Cai, J.; Zhang, L.; Chang, C.; Cheng, G.; Chen, X.; Chu, B. Hydrogen-Bond-Induced Inclusion Complex in Aqueous Cellulose/Lioh/Urea Solution at Low Temperature. ChemPhysChem 2007, 8, 1572–1579. [Google Scholar] [CrossRef] [PubMed]
  4. 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. [Google Scholar] [CrossRef]
  5. Jiang, Z.; Fang, Y.; Xiang, J.; Ma, Y.; Lu, A.; Kang, H.; Huang, Y.; Guo, H.; Liu, R.; Zhang, L. Intermolecular Interactions and 3d Structure in Cellulose-Naoh-Urea Aqueous System. J. Phys. Chem. B 2014, 118, 10250–10257. [Google Scholar] [CrossRef] [PubMed]
  6. Lindman, B.; Karlström, G.; Stigsson, L. On the mechanism of dissolution of cellulose. J. Mol. Liq. 2010, 156, 76–81. [Google Scholar] [CrossRef]
  7. Alves, L.; Medronho, B.; Filipe, A.; Antunes, F.E.; Lindman, B.; Topgaard, D.; Davidovich, I.; Talmon, Y. New Insights on the Role of Urea on the Dissolution and Thermally-Induced Gelation of Cellulose in Aqueous Alkali. Gels 2018, 4, 87. [Google Scholar] [CrossRef] [Green Version]
  8. Lindman, B.; Medronho, B. The Subtleties of Dissolution and Regeneration of Cellulose: Breaking and Making Hydrogen Bonds. BioResources 2015, 10, 3811–3814. [Google Scholar] [CrossRef]
  9. Medronho, B.; Lindman, B. Brief overview on cellulose dissolution/regeneration interactions and mechanisms. Adv. Colloid Interface Sci. 2015, 222, 502–508. [Google Scholar] [CrossRef] [PubMed]
  10. 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. [Google Scholar] [CrossRef]
  11. Soper, A.K. Empirical Potential Monte Carlo Simulation of Fluid Structure. Chem. Phys. 1996, 202, 12. [Google Scholar] [CrossRef]
  12. Ma, C.L.; Cheng, H.; Zuo, T.S.; Jiao, G.S.; Han, Z.H.; Qin, H. Neudatool: An Open Source Neutron Data Analysis Tools, Supporting Gpu Hardware Acceleration, and across-Computer Cluster Nodes Parallel. Chin. J. Chem. Phys. 2020, 33, 727–732. [Google Scholar] [CrossRef]
  13. Qin, H.; Ma, C.; Gärtner, S.; Headen, T.F.; Zuo, T.; Jiao, G.; Han, Z.; Imberti, S.; Han, C.C.; Cheng, H. Neutron total scattering investigation on the dissolution mechanism of trehalose in NaOH/urea aqueous solution. Struct. Dyn. 2021, 8, 014901. [Google Scholar] [CrossRef] [PubMed]
  14. Soper, A.K.; Castner, E.W., Jr.; Luzar, A. Impact of Urea on Water Structure: A Clue to Its Properties as a Denaturant? Biophys. Chem. 2003, 105, 649–666. [Google Scholar] [CrossRef]
  15. McLain, S.E.; Imberti, S.; Soper, A.K.; Botti, A.; Bruni, F.; Ricci, M.A. Structure of 2 molar NaOH in aqueous solution from neutron diffraction and empirical potential structure refinement. Phys. Rev. B 2006, 74, 094201. [Google Scholar] [CrossRef]
  16. Imberti, S.; Botti, A.; Bruni, F.; Cappa, G.; Ricci, M.A.; Soper, A.K. Ions in water: The microscopic structure of concentrated hydroxide solutions. J. Chem. Phys. 2005, 122, 194509. [Google Scholar] [CrossRef] [PubMed]
  17. Soper, A.K.; Ricci, M.A.; Bruni, F.; Rhys, N.H.; McLain, S.E. Trehalose in Water Revisited. J. Phys. Chem. B 2018, 122, 7365–7374. [Google Scholar] [CrossRef] [PubMed]
  18. Gong, P.; Wang, J.; Liu, B.; Ru, G.; Feng, J. Dissolution of chitin in aqueous KOH. Cellulose 2016, 23, 1705–1711. [Google Scholar] [CrossRef]
  19. Xiong, B.; Zhao, P.; Hu, K.; Zhang, L.; Cheng, G. Dissolution of Cellulose in Aqueous Naoh/Urea Solution: Role of Urea. Cellulose 2014, 21, 1183–1192. [Google Scholar] [CrossRef]
Figure 1. Structure of a trehalose molecule and its atomic labels used in EPSR simulation. All carbon atoms are labeled as C; all hydrogen atoms connected to the oxygen atom are labeled as H; the remaining hydrogen atoms connected to the carbon atom are labeled as M; the oxygen atom linking the two glucose rings is labeled as O1; the oxygen atoms on the glucose ring is labeled as O2; the oxygen atoms on the hydroxyl group connected to the glucose ring are labeled as O3; and the oxygen atom on the methyl group is labeled as O4. In order to describe the distribution of atoms around trehalose with spatial angles, we set up a coordinate system with O1 as the coordinate origin. The blue arrow connecting two C atoms at the glucose rings represents the direction of the X-axis; the red arrow pointing vertically to the blue arrow from O1 is the direction of the Y-axis; the Z-axis is the cross product of the X-axis and the Y-axis, which is approximately perpendicular to the paper.
Figure 1. Structure of a trehalose molecule and its atomic labels used in EPSR simulation. All carbon atoms are labeled as C; all hydrogen atoms connected to the oxygen atom are labeled as H; the remaining hydrogen atoms connected to the carbon atom are labeled as M; the oxygen atom linking the two glucose rings is labeled as O1; the oxygen atoms on the glucose ring is labeled as O2; the oxygen atoms on the hydroxyl group connected to the glucose ring are labeled as O3; and the oxygen atom on the methyl group is labeled as O4. In order to describe the distribution of atoms around trehalose with spatial angles, we set up a coordinate system with O1 as the coordinate origin. The blue arrow connecting two C atoms at the glucose rings represents the direction of the X-axis; the red arrow pointing vertically to the blue arrow from O1 is the direction of the Y-axis; the Z-axis is the cross product of the X-axis and the Y-axis, which is approximately perpendicular to the paper.
Molecules 27 03395 g001
Figure 2. Experimental neutron total scattering profiles (dashed black lines), the EPSR fitted neutron structure factors (red lines) and their differences (blue dots). (a) Trehalose in LiOH/urea aqueous solution, (b) trehalose in LiOH aqueous solution, (c) trehalose in KOH/urea aqueous solution, (d) trehalose in KOH aqueous solution. Data have been offset for clarity.
Figure 2. Experimental neutron total scattering profiles (dashed black lines), the EPSR fitted neutron structure factors (red lines) and their differences (blue dots). (a) Trehalose in LiOH/urea aqueous solution, (b) trehalose in LiOH aqueous solution, (c) trehalose in KOH/urea aqueous solution, (d) trehalose in KOH aqueous solution. Data have been offset for clarity.
Molecules 27 03395 g002
Figure 3. The schematic diagram of the accumulation of alkali metal ions and urea molecules around trehalose. (a) Li+ (blue dots), Na+ (dark blue dots), K+ (red dots) ions are distributed around a trehalose molecule. (b) Distribution of Li+ and the oxygen atoms of urea (red dots) around a trehalose molecule.
Figure 3. The schematic diagram of the accumulation of alkali metal ions and urea molecules around trehalose. (a) Li+ (blue dots), Na+ (dark blue dots), K+ (red dots) ions are distributed around a trehalose molecule. (b) Distribution of Li+ and the oxygen atoms of urea (red dots) around a trehalose molecule.
Molecules 27 03395 g003
Figure 4. RDF distributions between the alkali metal ions and the solute oxygen atoms. g(r)s between Li+, Na+, K+ and O1 (a), O2 (b), O3 (c), O4 (d) of trehalose.
Figure 4. RDF distributions between the alkali metal ions and the solute oxygen atoms. g(r)s between Li+, Na+, K+ and O1 (a), O2 (b), O3 (c), O4 (d) of trehalose.
Molecules 27 03395 g004
Figure 5. RDF distributions of some solvent atoms, i.e., Li, OW, OU, NU and OOH, around Os atoms in trehalose LiOH urea aqueous solution. (ad) is the RDF distribution around O1 to O4, respectively. The dashed lines are used to indicate the boundary of the second hydration shells.
Figure 5. RDF distributions of some solvent atoms, i.e., Li, OW, OU, NU and OOH, around Os atoms in trehalose LiOH urea aqueous solution. (ad) is the RDF distribution around O1 to O4, respectively. The dashed lines are used to indicate the boundary of the second hydration shells.
Molecules 27 03395 g005
Figure 6. RDF distributions of Li, OW, OU, NU and OOH around O1 atom along Z-axis. The dashed lines are used to indicate the boundary of the second hydration layer.
Figure 6. RDF distributions of Li, OW, OU, NU and OOH around O1 atom along Z-axis. The dashed lines are used to indicate the boundary of the second hydration layer.
Molecules 27 03395 g006
Figure 7. Schematic representation of the four important hydrogen bonds that urea forms with water and OH (a). The colors of the four atoms C, H, O and N are cyan, white, red and blue, respectively. RDF distributions between the atoms of urea and H2O/OH in trehalose LiOH urea aqueous solution (bd).
Figure 7. Schematic representation of the four important hydrogen bonds that urea forms with water and OH (a). The colors of the four atoms C, H, O and N are cyan, white, red and blue, respectively. RDF distributions between the atoms of urea and H2O/OH in trehalose LiOH urea aqueous solution (bd).
Molecules 27 03395 g007
Figure 8. Neutron total scattering profiles of trehalose in LiOH/urea aqueous solution at 25 °C and −5 °C and the difference between them (green lines).
Figure 8. Neutron total scattering profiles of trehalose in LiOH/urea aqueous solution at 25 °C and −5 °C and the difference between them (green lines).
Molecules 27 03395 g008
Table 1. Sample labels, chemical components, deuterium ratios and molar ratios of the samples.
Table 1. Sample labels, chemical components, deuterium ratios and molar ratios of the samples.
Sample LabelsLi(K)TrH2O/HDO/D2OLi(K)TrUrH2O/HDO/D2O
Chemical ComponentLi(K)OH Trehalose WaterLi(K)OH Urea Trehalose Water
Deuterium Ratio0.0/0.5/1.00.0/0.5/1.0
Molar RatioLi(K)OH:Urea:Trehalose:Water = 222:254:64:5716
Table 2. Lennard-Jones (L-J) reference potential parameters and charge (q) of atoms used in EPSR.
Table 2. Lennard-Jones (L-J) reference potential parameters and charge (q) of atoms used in EPSR.
(1)
Atom LabelOWHWLiNaKOOHHOHCUOU
e [KJ/mole]0.6500.0000.6900.1250.5000.2510.1840.4390.878
σ [Å]3.1660.0001.5102.5003.0002.7501.4430.3752.960
q [e]−0.8480.4240.6790.6790.679−1.1030.4240.142−0.390
(2)
Atom LabelNUHUCO1O2O3O4HM
e [KJ/mole]0.7110.0000.2760.5860.5860.7110.7110.0500.121
σ [Å]3.2500.0003.5003.1002.9003.1003.1001.7001.700
q [e]−0.5420.3330.258−0.500−0.500−0.500−0.5000.3010.000
Table 3. Cubic simulation box atomic number density and number of molecules used in EPSR.
Table 3. Cubic simulation box atomic number density and number of molecules used in EPSR.
Sample LabelsLi(K)OHUreaTrehaloseWaterDensity (Atoms/Å3)
Li(K)TrH2O22206457160.106531 (Li) 0.103323 (K)
Li(K)TrUrH2O2222546457160.106405 (Li) 0.103437 (K)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, C.; Zuo, T.; Han, Z.; Li, Y.; Gärtner, S.; Chen, H.; Yin, W.; Han, C.C.; Cheng, H. Neutron Total Scattering Investigation of the Dissolution Mechanism of Trehalose in Alkali/Urea Aqueous Solution. Molecules 2022, 27, 3395. https://doi.org/10.3390/molecules27113395

AMA Style

Ma C, Zuo T, Han Z, Li Y, Gärtner S, Chen H, Yin W, Han CC, Cheng H. Neutron Total Scattering Investigation of the Dissolution Mechanism of Trehalose in Alkali/Urea Aqueous Solution. Molecules. 2022; 27(11):3395. https://doi.org/10.3390/molecules27113395

Chicago/Turabian Style

Ma, Changli, Taisen Zuo, Zehua Han, Yuqing Li, Sabrina Gärtner, Huaican Chen, Wen Yin, Charles C. Han, and He Cheng. 2022. "Neutron Total Scattering Investigation of the Dissolution Mechanism of Trehalose in Alkali/Urea Aqueous Solution" Molecules 27, no. 11: 3395. https://doi.org/10.3390/molecules27113395

APA Style

Ma, C., Zuo, T., Han, Z., Li, Y., Gärtner, S., Chen, H., Yin, W., Han, C. C., & Cheng, H. (2022). Neutron Total Scattering Investigation of the Dissolution Mechanism of Trehalose in Alkali/Urea Aqueous Solution. Molecules, 27(11), 3395. https://doi.org/10.3390/molecules27113395

Article Metrics

Back to TopTop