Molecular Dynamics Study on Transmission Mechanism of Torsional Deformation in Cellulose Nanofibers with Hierarchical Structure

Cellulose nanofiber (CNF) is a fibrous and nano-sized substance produced by decomposition of bulk-type cellulose which is a main component of plants. It has high strength comparable to steel, and it shows low environmental load during a cycle of production and disposal. Besides it has many excellent properties and functions such as high rigidity, light-weight, flexibility and shape memory effect, so it is expected as a next-generation new material. Usually it is composed of many cellulose micro fibrils (CMFs) in which molecular chains of cellulose are aggregated in a crystal structure, the knowledge of mechanical properties for each CMF unit is important. Since actual fibrils are complicatedly intertwined, it is also crucial to elucidate the transmission mechanism of force and deformation not only in one fibril but also in between fibrils. How the dynamic and hierarchical structure composed of CMFs responds to bending or torsion is an interesting issue. However, little is known on torsional characteristics (shear modulus, torsional rigidity, etc.) concerning CMF. In general, in a wire-like structure, it is difficult to enhance torsional rigidity and strength, compared with tensile ones. Therefore, in this study, we try to build a hierarchical model of CNF by multiplying CMF fibers and to conduct molecular dynamics simulation for torsional deformation, by using hybrid model between all-atom and united-atoms model. First, shear modulus was estimated for one CMF fibril and it showed a value close to the experimental values. Also, we assume a state in which two CMFs are ideally arranged in parallel, and create a hierarchical structure. We evaluate the dependence on the temperature for the bond strength and toughness in the hieHow to cite this paper: Takada, K., Saitoh, K.-I., Sato, T., Takuma, M. and Takahashi, Y. (2019) Molecular Dynamics Study on Transmission Mechanism of Torsional Deformation in Cellulose Nanofibers with Hierarchical Structure. Soft Nanoscience Letters, 9, 45-57. https://doi.org/10.4236/snl.2019.94004 Received: September 5, 2019 Accepted: October 21, 2019 Published: October 24, 2019 Copyright © 2019 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access


Introduction
Currently, cellulose exists in the largest amount of resources on Earth. In recent years, a technology for manufacturing this cellulose from wood to fiber has been developed, and it is expected that cellulose can be applied to various fields as a kind of sustainable materials [1]. This type is called a cellulose nanofiber (CNF) and has a diameter of nanometer order [2]. CNF has the highly hierarchical structure, and the minimum component of CNF is called cellulose microfibril (CMF). Multiple cellulose molecules are polymerized to form a linear molecular chain, and 30 to 40 chains are regularly formed into a bundle. By intermolecular hydrogen bonds between chains, a molecular sheet is formed in the direction perpendicular to the axis of the chain. CNF is considered to be used as a structural member for automobiles and is expected to construct a fiber-reinforced composite material [3]. It will be an alternative material for carbon or glass fiber, having the advantage concerning resource amount and low environmental load.
The report on the evaluation of experimental mechanical properties of CNF has been increasing trend in recent years [4]. However, the mechanical properties of CMF are more required because CMF is dispersed in a matrix of composite and is nanoscale, and its experiment is not easy [5]. Therefore, CMF has been researched by a numerical simulation. For example, there is a report that estimates Young's modulus and tensile strength by tensile simulation in the axial or molecular sheet direction [6]. It is also pointed out that twist around the axis exists in CMF [7]- [12], and it is reported that the cause of twisting is discussed in terms of hydrogen bonding [13]. In considering that CNF will be integrated in actual structural and mechanical equipment, it will be also important to know the material's behavior in which it is subject to various deformations such as bending and twisting as well as tensile loading. It is reported that CNF's experimental shear modulus is between 1.8 and 3.8 GPa [14], and the simulated one by all-atom model is 1.6 GPa [15]. Although knowledge on tensile properties of CMF is increasing, the material's behavior in torsion or bending deformation is quite different [13] [16]. In addition, the actual CMF is tangled complicatedly to form a hierarchical structure. In natural material, through a lot of hierarchical structures inside one material, torsional force or deformation (i.e. chirality) is successfully transmitted [17]. Thus, it is worthwhile investigating hierarchical transmission mechanism in cellulose materials. Therefore, it is important not

Computational Methods and Models
For the conformation of cellulose nanofibers, we use the potential functions (force field) including stretch, bending, torsion, van der Waals interaction (Buckingham) and electrostatic (Coulomb) interaction [18] [19] [20]. cule is modelled as shown in Figure 1, and its parameters are shown in Table 1.
The cellulose I β type, which is a main component of plant-derived natural cellulose, is the present objective. The crystal structure of cellulose I β is monoclinic.
Therefore, 10 or 15 molecules, using the repeat unit model of Figure 1, are polymerized in the y axis direction to prepare a molecular chain. Subsequently, molecular chains are arranged in a crystal structure around the central molecular chain in the xz plane, and a CMF of five layers is made as shown in Figure   2. The conditions of those CMFs are shown in Table 2. CMFs with degrees of polymerization of 10 and 15 are named Model 1 and Model 2, respectively.
Generally, the actual degree of polymerization of CMF is 200 -300. But, in this study, we assume simple twisting of linear-shaped material with almost circular cross-section. The shorter model is preferable to exclude any complicated mode of deformation other than simple twist and it will make the discussion simple.
We can assume, when simple twist is applied, the cross-sectional shape is still axisymmetric around the rotation axis, so that all points at the same radius are   0.873 The number of atoms 28 The number of united atoms for CH 10 The number of united atoms for CH 2 2 The number of O atoms 10 The number of H atoms 6

Torsional Simulation of Single CMF
At first, structural relaxation to obtain thermal equilibration is performed on the The temperature condition is 10, 100, 200 and 300 K. In this research, to evaluate the shear modulus, it is necessary to calculate the atomic shear stress occurring in the rotation direction of each atom in the CMF. Therefore, we assume that the CMF is placed in a system of cylindrical coordinate. Each atomic stress tensor components in the Cartesian coordinates should be converted into those in the cylindrical coordinates and the average value calculated is obtained for all the atoms for each step.

Hierarchical Structure Simulation Using Two CMFs
Hierarchical structure simulation with multiple CMFs is performed using the Configuration of the torsion simulation model for the hierarchical structure is shown in Figure 5. The CMF on the left is named CMF1 and that on the right CMF2. Both CMF1 and CMF2 are fixed in all the directions along x, y, and z axes at one end (yellow-colored atoms in the figure), while only CMF 1 is given a constant angular velocity around the y axis at the other end (blue-colored atoms). That is, the other end of CMF2 is free of boundary constraint. The given constant angular velocity is ω = +20 rad/ns (it is done only in anticlockwise rotation, viewing from +y), and the calculation time is for 0.03 ns. The temperature condition is 10, 100, 200 and 300 K.
In addition, in this research of the hierarchical structure, "transmission ratio" of stress (force) or strain (deformation) between components (i.e. between CMF1 and CMF2) will be discussed. Here, the transmission ratio is defined as the ratio of the value of CMF2 to the value of CMF1 (since CMF1 is driven first). Suppose, for example, that the averaged shear stress of CMF1 and CMF2 at the final step of the loading are τ 1 and τ 2 , then the transmission ratio of the stress between them R is expressed by R = τ 2 /τ 1 .

Torsional Simulation of Single CMF
Shear stress-shear strain diagrams under positive and negative angular velocities are shown in Figure 6 and Figure 7. Table 3 shows the shear moduli calculated within the strain range from 0 to 0.05. Even when the strain is zero, the stress is not zero. This is because twist angle already occurs during structural relaxation anticlockwisely in viewing from +y axis and the residual stress exists inside. It is observed that, regardless of the twisting direction, the higher the temperature is, the lower stress value is obtained at the same strain. This indicates that, when the temperature rises, the CMF becomes soft and easily deformed. For any angular velocity condition, the shear modulus at the temperature less than 300 K is always smaller than that at 10 K. The value of shear modulus obtained here is about 1.4 GPa at 10 K, and about 0.8 -1.0 GPa between 100 and 300 K. They are somewhat smaller than the experimental value, 1.8 -3.8 GPa [15]. The calculated value reported for CMF using MD with the all atom model is 1.6 GPa [16]. This slight discrepancy is caused by inhomogeneous deformation of the whole system due to the influence of the fixed and velocity-constrained ends.

Tensile Simulation
The stress-strain diagram obtained by tensile simulation is shown in Figure 8. The maximum stress obtained for each temperature is shown in Table 4. The relationship between toughness and temperature is shown in Figure 9. Here, the toughness value is defined as the strain energy until the tensile strain reaches 0.3. The toughness in Figure 8 is shown by the strain energy per the initial volume of CMFs.
With respect to Figure 8, the stress increases with increasing strain up to the maximum stress, but then declines. Even when tensile deformation is applied, the original crystal shape of individual CMF is almost preserved. From this fact,      the hydrogen bond between CMF surfaces tends to be locally stretched and broken there, more possibly than bonds inside the CMF crystal. As a result, a large reduction in stress occurs. Under the temperature between 100 and 400 K, there is no big difference in the maximum stress values. The reason why the maximum stress at 10 K is particularly low is that the intermolecular hydrogen bond between the CMF surfaces is weakly formed there, compared to other temperatures. It means that the deformation resistance between the CMFs has weakened during combining process of the hierarchical structure. At 500 K, no stress reduction is observed after reaching the maximum values. In this temperature, the molecules in the CMF are close to the thermal decomposition, so that the entire hierarchical structure uniformly deforms, without concentration on the hydrogen bond between CMFs. In Table 4, the maximum stress value at 300 K is 0.68 GPa. Using the all-atom model, the tensile strength in the direction of the molecular sheet inside the CMF crystal was estimated from 0.40 to 0.90 GPa [6]. It is recognized that the hybrid model used in this research maintains important hydrogen bonds and accurately reflects the limit value of hydrogen bonding in the cellulose.
Since the manner of degradation of stress after reaching the maximum value is quite different between temperatures as shown in Figure 9, we will mention about the temperature dependence of the toughness value. Generally speaking, in crystals of small molecules, a solid state in low temperature changes to a liquid state in high temperature by increasing diffusion of the molecules inside the structure. On the other hand, in polymer crystals, atomic vibration in a molecule becomes more active as temperature increase, and so each single molecular chain moves in more flexible manner. However, the whole molecular chain is not able to diffuse or move in large distance due to entanglement and so on.
In CMF as well, the toughness also increases as the temperature rises, by virtue of flexibility of molecules. However, when the temperature is 400 K or more, the toughness value saturates at the same limit because the thermal decomposition of CMF into separated molecules has been reached in such high temperature.

Torsion Simulation
Transmission ratio of shear stress and shear strain from CMF1 to CMF2 at temperature from 10 to 300 K are shown in Figure 10. At 10 K, the shear strain is not transmitted at all, and the shear stress induces the counterpart to produce a stress in the opposite rotation direction (from rotational and fibrillar axis, it will be like a pair of mechanical gears transmitting torque). In this case, for very large deformation, the bond between surfaces of two CMFs is broken due to large separation and the CMF2 return to the original shape quickly. The transmission ratio of shear strain increases as the temperature rises. On the other hand, the transmission ratio of the shear stress decreases over 200 K. This indicates that the shear stress is not greatly transmitted at high temperature. Therefore, it is inferred that in a high temperature shear strain or deformation can be transmitted via small stress. Figure 11 shows the atomic arrangement for each time at 300 K and a schematic of torsion deformation in hierarchical structure.
As can be seen in Figure 11, it is confirmed that CMF1 shows torsional deformation, while CMF2 behaves like bending so that deflection occurs in the negative z direction. It shows possibility that the rotation of one CMF (torsional moment) is directly converted to the bending moment of adjacent CMF. Figure 10. Difference in transmission ratio for each temperature. Figure 11. Snapshots of hierarchical structure in torsion simulation at 300 K and a schematic of torsion deformation in hierarchical structure.

Conclusions
In this study, the mechanical properties of cellulose microfibrils (CMF), as the basic component of cellulose nanofibers (CNF), are investigated. In particular, torsional deformation and its transmission mechanism between components in hierarchical structure are discussed based on the results obtained by MD simulations.
• The shear modulus is about 0.7 -1.0 GPa and there is no clear temperature dependence between 100 and 300 K; • The strength of the hierachical structure in the sheet direction does not much depend on the temperature, but its toughness estimated from strain energy until it is broken shows strong temperature dependency; • Transmission of torsional deformation between components in the hierarchical structure is well observed but is dependent on temperature. The transmission ratio of mechanical values (ratio between averaged stress or strain occurring in two adjacent CMFs) increases as temperature rises, but that of shear stress tends to decrease above a certain temperature (in this model, above 200 K).
In our opinion, it is possible to design a new composite material by using the shear modulus of CMF which is difficult to evaluate by experiments but is able to be estimated by MD as done in our study. This study will lead us to a better understanding of the dynamic behavior of CMF in hierarchical structure.