Molecular Dynamics Simulation of the Effect of Low Temperature on the Properties of Lignocellulosic Amorphous Region
Abstract
:1. Introduction
2. Materials and Methods
2.1. Model Establishment
2.2. Dynamic Simulation
- Create a molecular model structure of the studied system, including the simulation box size, the number of model molecules in the simulation box, and the specific atom types, as in Section 2.1.
- Perform initial geometry optimization of the initial molecular simulation model by the computational options in the Forcite module. The algorithm is chosen as the Smart algorithm [35], and the total number of iterations was 5000 steps. This step is to balance the free motion state of the molecules in the model so as to achieve energy minimization.
- The kinetic relaxation of the geometrically optimized model by the computational option in the Forcite module so that the model is in a lower energy stable state and the initial internal stress is reduced. The target temperature of the system is set to room temperature (25 °C), and the initial velocity is random. A total length of 1 ns is simulated under the macroscopic rule (NVT) system and the time step is set to 1 fs.
- After the system reached equilibrium and then entered into formal kinetic simulations for six sets (20 °C, 0 °C, −30 °C, −70 °C, −110 °C, and −150 °C). All six sets of simulations were performed at atmospheric pressure (0.1 MPa) and 1 ns under isothermal isobaric system synthesis (NPT). The pressure of the simulated experiments was controlled using the Berendsen method, chosen to be suitable for calculating the PCFF force field of natural polymer materials [30]. The electron summation was controlled by the Ewald method, and the van der Waals force calculation was controlled by the Atom-based method [37,38,39,40].
3. Results and Discussion
3.1. The Energy
3.1.1. System Energy
3.1.2. Interaction Energy
3.2. Model Volume and Density
3.3. Hydrogen Bonding
3.4. Mechanical Properties
4. Conclusions
- The absolute value of the interaction energy between water molecules and cellulose chains increases as the temperature decreases. It indicates that the interaction between water molecules and cellulose amorphous region is stronger by decreasing the temperature. The interaction energies of the six groups of models at different temperatures are all negative, indicating that water molecules and cellulose chains are attracted to each other, and the system can exist stably.
- From 20 °C to −150 °C, the model density increased from 1.413 g/cm3 to 1.459 g/cm3, an increase of 3.26%. The volume decreases from 9098.866 Å3 to 8803.711 Å3, a decrease of 3.24%. The internal structure of the model is influenced by the temperature, and the changes in the model parameters can be verified with the interaction energy with each other. As the temperature decreases, the model interaction energy increases, the volume decreases, and the density increases.
- The decrease in temperature makes cellulose molecular activity weaker. The total number of hydrogen bonds and the number of hydrogen bonds between water molecule–cellulose chains for each model increased with decreasing temperature. The new interchain hydrogen bonds enhanced the restraining effect on the arrangement of cellulose molecular chains. The values of G, E, and K increased with decreasing temperature, and K/G decreased with decreasing temperature. The low-temperature treatment increased the stiffness and reduced the toughness of the wood. The changes in mechanical properties were characterized by the number of hydrogen bonds to each other, and the increase in the number of hydrogen bonds will inevitably increase the mechanical properties and chemical stability of the cellulose material.
Author Contributions
Funding
Data Availability Statement
Conflicts of Interest
References
- Esteves, B.; Pereira, H. Wood modification by heat treatment: A review. BioResources 2009, 4, 370–404. [Google Scholar] [CrossRef]
- Ma, X.; Zhang, F.; Wei, L. Effect of wood charcoal contents on the adsorption property, structure, and morphology of mesoporous activated carbon fibers derived from wood liquefaction process. J. Mater. Sci. 2014, 50, 1908–1914. [Google Scholar] [CrossRef]
- Wei, L.; Liang, S.; McDonald, A.G. Thermophysical properties and biodegradation behavior of green composites made from polyhydroxybutyrate and potato peel waste fermentation residue. Ind. Crops Prod. 2015, 69, 91–103. [Google Scholar] [CrossRef]
- Wei, L.; McDonald, A.G.; Stark, N.M. Grafting of bacterial polyhydroxybutyrate (PHB) onto cellulose via in-situ reactive extrusion with dicumyl peroxide. Biomacromolecules 2015, 16, 1040–1049. [Google Scholar] [CrossRef]
- Carlmark, A.; Larsson, E.; Malmström, E. Grafting of cellulose by ring-opening polymerization—A review. Eur. Polym. J. 2012, 48, 1646–1659. [Google Scholar] [CrossRef] [Green Version]
- Green, D.W.; Evans, J.W.; Logan, J.D.; Nelson, W.J. Adjusting modulus of elasticity of lumber for changes in temperature. For. Prod. J. 1999, 49, 82–94. [Google Scholar]
- Lv, J.X.; Lin, Z.Y.; Jiang, J.L. Influence of different drying methods on the infiltration properties of wood from fir plantations. For. Sci. 2006, 4, 85–90. [Google Scholar]
- Schaffer, E.L. Effect of pyrolytic temperatures on longitudinal strength of dry Douglas-fir. J. Test. Eval. 1973, 1, 319–329. [Google Scholar]
- Moraes, P.D.; Rogaume, Y.; Bocquet, J.F. Influence of temperature on the embedding strength. Holz Roh-Werkst 2005, 63, 297–302. [Google Scholar] [CrossRef]
- Adewopo, J.B.; Patterson, D.W. Effects of heat treatment on the mechanical properties of Loblolly Pine, Sweetgum, and Red Oak. For. Prod. J. 2011, 61, 526–535. [Google Scholar] [CrossRef]
- Zhou, J.H.; Hu, C.S.; Hu, S.F. Effects of temperature on the bending performance of wood-based panels. BioResources 2012, 7, 3597–3606. [Google Scholar]
- Kubojima, Y.; Okano, T.; Ohta, M. Bending strength and toughness of heat-treated wood. J. Wood Sci. 2000, 46, 8–15. [Google Scholar] [CrossRef]
- Moraes, P.; Rogaume, Y.; Triboulot, P. Influence of temperature on the modulus of elasticity (MOE) of Pinus sylvestris L. Holzforschung 2004, 58, 143–147. [Google Scholar] [CrossRef]
- Windeisen, E.; Bchle, H.; Zimmer, B. Relations between chemical changes and mechanical properties of thermally treated wood. Holzforschung 2009, 63, 773–778. [Google Scholar] [CrossRef]
- Awoyemi, L.; Jones, I.P. Anatomical explanations for the changes in properties of western red cedar (Thujaplicata) wood during heat treatment. Wood Sci. Technol. 2011, 45, 261–267. [Google Scholar] [CrossRef]
- Kamody, D.J. Using deep cryogenic to advantage. Adv. Mater. Process. 1998, 154, 215–218. [Google Scholar]
- Yi, X.; Chen, D. Research progress on deep cooling treatment of composite materials. Met. Heat Treat. 2012, 37, 73–76. [Google Scholar]
- Owaku, S. Cryogenic treatment. Heat. Treat. 1981, 21, 44–48. [Google Scholar]
- Kalia, S. Cryogenic processing: A study of materials at low temperatures. J. Low Temp. Phys. 2010, 158, 934–945. [Google Scholar] [CrossRef]
- Yamada, T. Temperature Dependency of Physical Properties of Wood at Low Temperature; Kyoto University Research Information Repository: Kyoto, Japan, 1971. [Google Scholar]
- Berner, D.; Gerwick, B.C.; Polivka, M. Static and cyclic behavior of structural lightweight concrete at cryogenic temperatures. In Temperature Effects on Concrete ASTM International; Springer: Tokyo, Japan, 1985. [Google Scholar]
- DeGeer, D.; Bach, L. Machine Stress Grading of Lumber at Low Temperatures; Alberta Research Council: Edmonton, AB, Canada, 1995. [Google Scholar]
- Xu, H.D.; Wang, Y.T.; Wang, L.H.; Wang, X.P. A review on ice formation and propagation in wood cells at subzero temperatures. J. Nanjing For. Univ. 2017, 41, 169–174. [Google Scholar]
- Li, Z.; Jiang, J.L.; Lu, J.X. Moisture-dependent orthotropic viscoelastic properties of Chinese fir wood in low temperature environment. J. Wood Sci. 2018, 64, 515–525. [Google Scholar] [CrossRef] [Green Version]
- Ozkan, O.E. Effect of Freezing Temperature on Impact Bending Strength and Shore-D Hardness of Some Wood Species. BioResources 2022, 17, 6123–6130. [Google Scholar] [CrossRef]
- Paschalina, T.; Vasiliki, K. Chemical characterization of Wood and Bark biomass of the invasive species of Tree-of-heaven (Ailanthus altissima (Mill.) Swingle), focusing on its chemical composition horizontal variability assessment. Wood Mater. Sci. Eng. 2022, 17, 469–477. [Google Scholar]
- Jiang, X.W.; Wang, W.; Guo, Y.Y.; Dai, M. Effect of pressurized hydrothermal treatment on the properties of cellulose amorphous region based on molecular dynamics simulation. Forestry 2023, 14, 314. [Google Scholar] [CrossRef]
- Khazraji, A.C.; Robert, S. Interaction effects between cellulose and water in nanocrystalline and amorphous regions: A novel approach using molecular modeling. J. Nanomater. 2013, 2013, 44. [Google Scholar] [CrossRef] [Green Version]
- Hou, T.J.; Zhang, W.; Xu, X.J. Binding affinities for a series of selective inhibitors of gelatinase-A using molecular dynamics with a linear interaction energy approach. J. Phys. Chem. B 2005, 105, 5304–5315. [Google Scholar] [CrossRef]
- Wang, W.; Ma, W.; Wu, M.; Sun, L. Effect of Water Molecules at Different Temperatures on Properties of Cellulose Based on Molecular Dynamics Simulation. BioResources 2020, 17, 269–280. [Google Scholar] [CrossRef]
- Karplus, M.; McCammon, J.A. Molecular dynamics simulations of biomolecules. Nat. Struct. Biol. 2002, 9, 646–652. [Google Scholar] [CrossRef]
- Steinhauser, M.O.; Hiermaier, S. A review of computational methods in materials science: Examples from shock-wave and polymer physics. Int. J. Mol. Sci. 2009, 10, 5135–5216. [Google Scholar] [CrossRef] [Green Version]
- Hollingsworth, S.A.; Dror, R.O. Molecular dynamics simulation for all. Neuron 2018, 99, 1129–1143. [Google Scholar] [CrossRef] [Green Version]
- Liu, C.B.; Sun, Z. Molecular dynamics simulation of phase transition process in nematic liquid crystal materials. Mech. Sci. Technol. 2023, 33, 1–9. [Google Scholar]
- Mazeau, K.; Heux, L. Molecular dynamics simulations of bulk native crystalline and amorphous structures of cellulose. J. Phys. Chem. B 2003, 107, 2394–2403. [Google Scholar] [CrossRef]
- Allen, M.P.; Tildesley, D.J. Computer Simulation of Liquids; Clarendon Press: Oxford, UK, 2017. [Google Scholar]
- Ewald, P. Evaluation of optical and electrostatic lattice potentials. Ann. Phys. 1921, 64, 253–287. [Google Scholar] [CrossRef] [Green Version]
- Andersen, H.C. Molecular dynamics simulations at constant pressure and/or temperature. J. Chem. Phys. 1980, 72, 2384–2393. [Google Scholar] [CrossRef] [Green Version]
- Andrea, T.A.; Swope, W.C.; Andersen, H.C. The role of long ranged forces in determining the structure and properties of liquid water. J. Chem. Phys. 1983, 79, 4576–4584. [Google Scholar] [CrossRef]
- Berendsen, H.J.; Postma, J.V.; Van Gunsteren, W.F.; DiNola, A.; Haak, J.R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684–3690. [Google Scholar] [CrossRef] [Green Version]
- Ouyang, D.; Zhang, L.; Mao, R.; Qin, X.; Pang, W. Application Process of Coating Agent and the Coating Effect Evaluation Based on Molecular Dynamics. Langmuir 2023, 39, 3411–3419. [Google Scholar] [CrossRef]
- Hanus, J.; Mazeau, K. The xyloglucan-cellulose assembly at the atomic scale. Biopolymers 2006, 82, 59–73. [Google Scholar] [CrossRef]
- Bodig, J.; Jayne, B.A. Mechanic of Wood and Wood Composites; Van Nostrand Reinhold: New York, NY, USA, 1982. [Google Scholar]
- Zhang, Y.; Duan, Y.; Liu, J. The effect of intermolecular hydrogen bonding on the polyaniline water complex. J. Cluster Sci. 2017, 28, 1071. [Google Scholar] [CrossRef] [Green Version]
- Zheng, Y.; Deng, G.; Guo, R. A DFT-based study of the hydrogen-bonding interactions between myricetin and ethanol/water. J. Mol. Model. 2019, 25, 67. [Google Scholar] [CrossRef]
- Mishra, S.; Suryaprakash, N. Intramolecular hydrogen bonding involving organic fluorine: NMR investigations corroborated by DFT-based theoretical calculations. Molecules 2017, 22, 423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Juan, Y.; An, Y.L. Theoretical study of hydrogen bonding excited states of fluorenone with formaldehyde. Comput. Theor. Chem. 2016, 1101, 62. [Google Scholar]
- Yang, J.; Li, A.Y. Hydrogen bond strengthening between o-nitroaniline and formaldehyde in electronic excited states:A theoretical study. Spectrochim. Acta A 2018, 199, 194. [Google Scholar] [CrossRef]
- Du, D.Y.; Tang, C.; Zhang, J.W. Effects of hydrogen sulfide on the mechanical and thermal properties of cellulose insulating paper: A molecular dynamics simulation. Mater. Chem. Phys. 2020, 240, 122153. [Google Scholar] [CrossRef]
- Ozer, I. Some Mechanical Properties of Frozen Juvenile and Mature Wood of Pinus Nigra Var pallasiana Aral. and Pinus brutia Ten; Institute of Natural Sciences, Istanbul University: Istanbul, Turkey, 2002; Volume 75. [Google Scholar]
- Schmidt, R.A.; Pomeroy, J.W. Bending of a conifer branch at subfreezing temperatures: Implications for snow interception. Can. J. For. Res. 1990, 20, 1250–1253. [Google Scholar] [CrossRef]
Temperature (°C) | 20 | 0 | −30 | −70 | −110 | −150 |
---|---|---|---|---|---|---|
1476 | 1417 | 1336 | 1158 | 1078 | 1002 | |
2170 | 2138 | 2075 | 1944 | 1893 | 1839 | |
−93 | −88 | −89 | −110 | −136 | −148 | |
−601 | −633 | −650 | −676 | −679 | −689 |
Temperature (°C) | Cell Parameters (Å) | ||
---|---|---|---|
The Length | The Width | The Height | |
20 | 20.88 | 20.88 | 20.88 |
0 | 20.82 | 20.82 | 20.82 |
−30 | 20.89 | 20.89 | 20.89 |
−70 | 20.78 | 20.78 | 20.78 |
−110 | 20.68 | 20.68 | 20.68 |
−150 | 20.65 | 20.65 | 20.65 |
Temperature (°C) | 20 | 0 | −30 | −70 | −110 | −150 |
---|---|---|---|---|---|---|
229.20 | 232.01 | 235.22 | 244.80 | 244.64 | 251.59 | |
87.00 | 83.85 | 86.37 | 89.21 | 89.60 | 85.89 | |
32.69 | 30.23 | 28.62 | 32.96 | 38.66 | 34.99 | |
109.51 | 117.94 | 120.23 | 122.63 | 116.38 | 130.71 |
Temperature (°C) | 20 | 0 | −30 | −70 | −110 | −150 |
---|---|---|---|---|---|---|
8.29 | 8.01 | 10.31 | 10.30 | 9.73 | 11.00 | |
8.38 | 7.83 | 9.82 | 10.54 | 12.58 | 11.18 | |
7.49 | 8.02 | 9.72 | 11.71 | 11.64 | 11.62 | |
4.07 | 5.11 | 8.42 | 5.28 | 6.17 | 9.06 | |
1.97 | 4.29 | 7.15 | 6.32 | 6.24 | 7.26 | |
5.61 | 3.38 | 8.22 | 6.79 | 7.60 | 7.32 | |
8.05 | 7.95 | 9.95 | 10.85 | 11.32 | 11.27 | |
3.88 | 4.26 | 7.93 | 6.13 | 6.67 | 7.88 | |
10.39 | 11.30 | 20.27 | 16.18 | 17.54 | 20.40 | |
8.23 | 8.11 | 10.04 | 10.96 | 11.42 | 11.35 | |
2.12 | 1.90 | 1.27 | 1.79 | 1.71 | 1.44 |
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Jiang, X.; Wang, W.; Guo, Y.; Dai, M. Molecular Dynamics Simulation of the Effect of Low Temperature on the Properties of Lignocellulosic Amorphous Region. Forests 2023, 14, 1208. https://doi.org/10.3390/f14061208
Jiang X, Wang W, Guo Y, Dai M. Molecular Dynamics Simulation of the Effect of Low Temperature on the Properties of Lignocellulosic Amorphous Region. Forests. 2023; 14(6):1208. https://doi.org/10.3390/f14061208
Chicago/Turabian StyleJiang, Xuewei, Wei Wang, Yuanyuan Guo, and Min Dai. 2023. "Molecular Dynamics Simulation of the Effect of Low Temperature on the Properties of Lignocellulosic Amorphous Region" Forests 14, no. 6: 1208. https://doi.org/10.3390/f14061208
APA StyleJiang, X., Wang, W., Guo, Y., & Dai, M. (2023). Molecular Dynamics Simulation of the Effect of Low Temperature on the Properties of Lignocellulosic Amorphous Region. Forests, 14(6), 1208. https://doi.org/10.3390/f14061208