Structure, Reactivity, and Mechanical Properties of Sustainable Geopolymer Material: A Reactive Molecular Dynamics Study

Sodium aluminosilicate hydrate (NASH) gel, the primary binding phase in geopolymer, determines the mechanical properties and durability of environment-friendly construction materials. In this work, the models of NASH gel were obtained through a two-step procedure: the temperature quenching method and Grand Canonical Monte Carlo water adsorption. The reactive force field (ReaxFF) molecular dynamics were utilized to investigate the structure, reactivity, and mechanical performance of the NASH gel with Na/Al ratio ranging from 1 to 3. Q species, the connectivity factor, shows that the increase of sodium content in NASH gel leads to depolymerization of the aluminosilicate network and more non-bridging oxygen (NBO) atoms. The adsorbed water molecules dissociate near the NBO with high reactivity in defective aluminosilicate structure. The newly produced hydroxyls associate with the aluminate species, contributing to the formation of the pentahedron local structure. The sodium ions distributed in the cavity of the aluminosilicate skeleton have around 4 7 nearest neighbors. Furthermore, with an increase in sodium, the molecular structure of the aluminosilicate skeleton is transformed from an integrity network to partially destroyed branch structures, which gradually decrease the stiffness and cohesive force of NASH gel, characterized by the uniaxial tensile testing. During the large tensile deformation process, the ReaxFF MD correlates the mechanical response with the chemical reaction pathway. The aluminosilicate skeleton is stretched broken to resist the tensile loading and the hydrolytic reaction of water molecules near the stretched Si-O and Al-O bond further accelerates the degradation of NASH gel. Hopefully, this work can shed light on the material design for a high performance of sustainable geopolymer at the nanoscale.


INTRODUCTION
Ordinary Portland cement (OPC) is ubiquitously utilized as the essential construction and building material worldwide (Li, 2011). Cement demand is estimated to increase from around 3.5 Gt in 2015 to nearly 4.4 Gt per year by 2050 (Luukkonen et al., 2018). The manufacturing of OPC is accompanied with CO 2 emissions such as the calcination of limestones at 1,450 • C and the energy consumed by the cement plant itself, which is supplied by coal combustion (Damtoft et al., 2008). In 2016, the production of OPC contributed to around 1.45 Gt CO 2 , occupying approximately 8% of the global CO 2 release (Andrew, 2018). Even though great efforts have made to improve energy efficiency, the clean production of cement remains a difficult issue to solve, especially considering the increasing cement demand and the inevitable CO 2 emission reactions from cement production. Consequently, to reduce the carbon footprint, it is necessary to develop environment-friendly construction material to serve as an alternative and supplementary binder for OPC. Geopolymer, one type of alkali activated aluminosilicate cement, is used primarily as an environmentally beneficial alternative to OPC (Davidovits, 1991;Duxson et al., 2007a). Geopolymers use industrial byproducts as precursors, and therefore result in dramatically less CO 2 emissions per ton of concrete produced (Duxson et al., 2007a). The industrial byproducts utilized in geopolymer synthesis include fly ash and slag, which when combined with alkaline activators react to form a hardened binder possessing performance characteristics comparable to traditional Portland cement (Davidovits, 1982;Palomo et al., 1999;Van Jaarsveld et al., 2002;Bakharev, 2005). It is synthesized by dissolution of Al and Si in alkali medium, transportation, and polycondensation, forming a three-dimensional network structure. Due to their inorganic three-dimensional network structure, geopolymers show high efficiency of fireproofing, excellent thermal stability, and superior mechanical properties as compared with traditional OPC (Palomo and Glasser, 1992;Xu and Van Deventer, 2000; Barbosa and MacKenzie, 2003;Kriven et al., 2003). Geopolymers have therefore received great growing in recent years in the field of cleaner production of construction materials.
Sodium aluminosilicate hydrate (NASH) gel is a primary binding phase in geopolymers. According to most of the papers published so far (White et al., 2012;Parthasarathy et al., 2017), the mechanical properties, mass transport, ion exchange, and other physicochemical properties of NASH gels are controlled by the chemical composition and nanostructure, which have drawn attention from researchers in related fields. In order to determine the coordination state of aluminum and oxygen, and to explain the abnormal changes of physical properties from the perspective of structure, many experiments and simulations have been carried out (McKeown et al., 1984;Leonelli et al., 2001;Okuno et al., 2005;Sadat et al., 2016). Earlier spectroscopic studies on geopolymer paste show that the structure of NASH gel primarily consists of an alkali aluminosilicate hydrate gel framework (the alkali used commonly is Na) formed by disordered interlinked aluminosilicate tetrahedrons. Al and Si are both present in tetrahedral coordination connected by bridging oxygen (BO), and the negative charge associated with Al substitution for Si is balanced by the alkali cations (Barbosa et al., 2000;Duxson et al., 2005Duxson et al., , 2007aSingh et al., 2005). Compared with the bonding between Na and NBO, the interaction between charge compensating alkali cations (Na + ) and [AlO 4 ] − units are more ionic and weaker (Uchino et al., 1993). The Na/Al ratio in compositions indicates the change in internal structure. In the geopolymers with low sodium content, Na + is bound superficially to the [AlO 4 ] − . When the Na/Al ratio is larger than 1, Na + may break the Si-O-Si linkage, associating it with part of the silicon in the form of SiO-(Na + )-Si − (Uchino et al., 1993). Using molecular dynamics, Xiang et al. (2013) studied the changes in the structure and mechanical properties of sodium aluminosilicate glass as a function of Al/Na ratios, and the results were in good agreement with experimental data. Zhang M. et al. (2018) studied the local structure and dynamics of sodium in NASH gel, demonstrated the loose connection between sodium and aluminate tetrahedron, and modeled the dissociation process of sodium, which is so-called leaching. Meanwhile, the ion immobilizing ability of NASH gel was studied by molecular dynamics, providing a fundamental understanding of the immobilization mechanism of geopolymer materials (Hou et al., 2019b). An empirical force field molecular dynamics method was also employed to investigate the structure, dynamics, energetics, and mechanical properties of calcium silicate hydrate (C-S-H) (Hou et al., 2019a(Hou et al., , 2015bWang et al., 2019). Furthermore, reactive molecular dynamics were used to simulate the polymerization process and the molecular structure of geopolymer gels, and evaluated the influence of the simulation temperature and Si/Al ratio on the geopolymer (Zhang M. et al., 2018).
Typically, the properties of the materials subjected to loading will be greatly reduced largely due to the degradation of the chemical and physical structure. The aluminosilicate skeleton significantly influences the adsorption capability, reactivity, and hydrogen bond of the structural water molecules. The way of binding water in NASH gel is of particular importance, as it is closely related to strength and stiffness. In some cases, the internal damage of NASH gel is caused by the hydrolysis reaction of the silicate network, attacked by dissociated water molecules (Hou et al., 2014a(Hou et al., , 2016a. Sadat et al. (2016) reported that exposure to water promotes the diffusion and dissociation of structural sodium, leading to destabilization of aluminate tetrahedrons and underling mechanical failure of the local structure. However, minute quantities of water could stabilize the sodium in NASH gel, which was attributed to the hydrolysis associated with aluminate tetrahedron and the formation of Na-OH (Zhang Y. et al., 2018). The attack of water molecules on chemical bonds promoted bond breaking in the loading test (Hou et al., 2014a;Zhang Y. et al., 2018). In addition, some studies have focused on the effects of chemical compositions (Hou et al., 2015a(Hou et al., , 2016b(Hou et al., , 2015c. To the best of our knowledge, the mechanical properties and weakening mechanisms have not been fully understood and the case of coupling with hydrolysis and depolymerization, in particular, has not been studied yet. In this study, NASH gel with different compositions of the Na/Al ratio were characterized by molecular dynamics using a reactive force field (ReaxFF), which can perform a physical and chemical process in hydration, aiming to understand the structural role of sodium as a function of the chemical composition. The structural, mechanical, and chemical properties of NASH gel with different Na/Al ratios have been investigated. The influences of an existent form of water within NASH gel were also analyzed to obtain a better understanding of their behavior. The obtained results highlight that NASH gel is weakened by the breakage of aluminosilicate chains and the penetration of water molecules, and show the mechanism of influence on its properties.

Reactive Force Field
The ReaxFF, developed by van Duin et al. (2001), was used to simulate the chemical reaction between sodium aluminosilicate glass and water molecules for both atomic structure construction and uniaxial tensile testing. The ReaxFF provides an advanced description of the interaction between water and NASH surface. In addition, the ReaxFF has been widely used in silica glass (Hou et al., 2014a), C-S-H gel (Hou et al., 2015c), ordered crystal, and disordered glass (Hou et al., 2016b), discussing their structure evolution and tensile behavior under water penetration. The short-range interactions for the ReaxFF are determined by bond length-bond order scheme so that the bonds can be broken and formed with the potential energy transforming into a smooth state (Brenner et al., 2002). On the contrary, the long-range Coulombic interactions are determined by a seventh-order taper function, with an outer cutoff radius of 10 Å. The parameters of the force field for Na, Si, Al, O, and H can be obtained directly from previously published reference data (Cygan et al., 2004).

NASH Models With Different Na/Al Ratios
The LAMMPS package was used to construct three NASH gel models with Na/Al ratios of 1, 2, and 3 via two steps, respectively. First, a temperature quenching process was performed to build the NAS glass model. The water molecules were then inserted into the NAS glass model to obtain the NASH gel model. In order to obtain initial configuration with different Na/Al ratios, four types of atoms (Na, Si, Al, and O) were added in a cubic box with a side length of 20 Å to generate random coordinates. The density was adjusted to a typical value (see Table 1; Sadat et al., 2016). A time step of 0.1 fs was used for the simulation. The system was operated by a heating/cooling cycle process from 300 to 4,000 K under NPT ensemble at a rate of 10 K/ps, then the system was cooled from 4,000 to 300 K at the same rate. Finally, the system was further equilibrated at NVT ensemble with the temperature maintained at 300 K for 300 ps. The simulated equilibrium structure of the aluminosilicate chain with different Na/Al ratios are shown in Figure 1.
The Grand Canonical Monte Carlo (GCMC) method was utilized to adsorb water molecules into the NAS structure to obtain a NASH gel model (Bonnaud et al., 2012). The simulation process is analogous to the water adsorption in a microscopic porous structure such as C-S-H and silica glass (Hou et al., 2014a(Hou et al., , 2015c. A two million step simulation was performed for the system equilibrium, followed by one million steps for data analysis. Water molecules can be inserted, deleted, displaced, and rotated in the constant volume system. After GCMC simulation, the NASH gel model saturated with water molecules will be obtained. The aluminosilicate skeletons with a Na/Al ratio of 1, 2, and 3 are shown in Figure 1. The composition of the NASH gel models is shown in Table 1.

Uniaxial Tension Test
A uniaxial tension test was utilized to detect the mechanical behavior of amorphous NASH gel. In order to reduce error and to ensure stability of the simulation system during stretching, the initial NASH models (20 × 20 × 20 Å) were expanded periodically (40 × 40 × 40 Å) and then used in a uniaxial tension test in the X direction. Moreover, there are 5,000 to 6,000 atoms in every supercell with a specific Na/Al ratio. It should be noted that a large number of atoms in simulation models can provide stable statistical results, especially for the reliable failure modes. To explore the failure mechanism of the NASH gel, the stress-strain relation and the change of molecular structure during loading were investigated.
To obtain the stress-strain relations, the supercell structures were subjected to uniaxial tensile loadings through gradual elongation at a strain rate of 0.08/ps. In the uniaxial tensile simulation process, NPT ensembles were defined for the system. When stretched along the X direction, the supercells were relaxed, coupled with zero external pressure in the x, y, and z dimensions at 300 K for 400 ps. Until the pressures reached equilibrium in three directions, the NASH gel would be elongated in the X direction. Furthermore, the pressure was zero in the Y and Z directions, respectively. Pressure evolution in the X direction was taken as the internal stress σ xx . By setting the pressure to 0, which is perpendicular to the stretch direction, the normal direction can be relaxed toward the anisotropy without any restriction. The Poisson's ratio was considered in the setting, eliminating the artificial constraints of deformation.

Structure Properties
After water adsorption, the molecular structures of the NASH gel were obtained. Water molecules penetrated into the cavity region of NAS glass and randomly distributed in the NASH gel. As shown in Figure 2, for gel with a Na/Al ratio of 1, two forms of the water molecules influencing the NAS glass can be found: chemical decomposition and physical adsorption. The water molecules adsorb and form physical bonds with NASH at first. Very small percentages of water molecules are decomposed to hydrogen and hydroxyl. The hydrogen forms a chemical bond with the oxygen of silicon-oxygen tetrahedron, and the hydroxyl forms a chemical bond with the aluminum. For chemical decomposition (see region I), a small number of water molecules (less than 7.2%) break the aluminosilicate skeleton by occupying the BO site, and construct aluminum hydroxyl or silicate hydroxyl. In the case of physical adsorption water (see region II), the hydroxyl in water molecules point to the aluminosilicate skeleton to form a hydrogen bond with its oxygen atoms, and these water molecules, called physical adsorption water, makes up more than 92.8% of the total. From the molecular structure of the static state, physical adsorption is the main existing form for water molecules in NASH gel, and the proportion of decomposition is very low. This conclusion agrees with previously published works (Hou et al., 2014a;Zhang Y. et al., 2018). The three-dimensional network structure of the NASH gel with a Na/Al ratio of 1 is shown in Figure 2, which is the most representative feature. In the structure, silicon and aluminum atoms are interconnected by BO, and nearby sodium ions play a role in charge balance (Barbosa et al., 2000;Schmücker and MacKenzie, 2005). The coordination number (CN) of NASH gel, obtained using dynamic simulations, shows the local structure. In thermodynamic equilibrium state, as shown in Table 2, the CN of all the silicon in the gel with different Na/Al ratios is fourfold coordination, indicating that silicon in NASH gel exists in the form of tetrahedron. The aluminum coordination in the network structure is also predominantly tetrahedral (∼95.57, 94.14, and 96.51%), and some aluminum (∼2.16, 5.84, and 3.49%) with fivefold coordination, which is related to formation of Al-O w H w . Sodium diffuses in the cavity of the Si-Al framework, and exists in the form of Na(H 2 O) n + rather than Na + , implying a loose interaction between sodium and the neighboring Si-O-Al bond. Sodium is located at a vacant region among the skeleton, showing multiple coordination forms. The results of MD simulation match well with previous experimental results ( Barbosa et al., 2000;Duxson et al., 2005;De Silva et al., 2007;Walkley et al., 2018).
In NASH gel, Si and Al are predominantly tetrahedral, which can also be revealed by  respectively. The simulation results are in agreement with the experimental results given by nuclear magnetic resonance, and X-ray diffractions are between 142 • and 151 • (Mauri et al., 2000). More significantly, the Si-O-Si angle peak value decreased slightly with the increase of the Na/Al ratio due to the change of sodium ions. However, Figure 3C shows that the angle distributions of O-Al-O are far broader than that of O-Si-O, and there are more obvious peaks located at different angles. It can be attributed to two possible factors: aluminum exists in multiple coordination configurations in the NASH gel, and the Al-O bond is relatively unstable in the aluminosilicate skeleton compared with the Si-O bond.

Stress-Strain Relation for NASH Gel
Given the constitutive stress-strain relation, the stress-strain curve can reflect the mechanical properties of NASH gel in the process of uniaxial tension testing. Due to the isotropic structure of NASH, in the directions X, Y, and Z, they have the same stressstrain curves. Figure 4 shows the constitutive relation between stress and strain with a strain rate of 0.08/ps along the X direction.
Although the stress-strain curves have similar evolutionary trends, the intrinsic strength and failure strain of the specimen with a different Na/Al ratio is obviously different from others. The tensile strain, when the stress reaches the maximum value, is gradually decreased from 0.45 Å/Å to 0.3 Å/Å. The NASH gel with a lower Na/Al ratio is more easily stretched broken, and the yield strength of the NASH gel with a Na/Al ratio of 1, 2, and 3 is 17.7, 16.2, and 14.3 GPa, respectively. From the perspective of composition, the stress maximum value decreased as the Na/Al ratio increased. The number of sodium ions beyond that required, compensated for the negatively charged [AlO4] − units. Excess sodium ions with a positive charge can then attack the Si-O-Si and Si-O-Al bond (Irwin et al., 1988;Barbosa et al., 2000;Schmücker and MacKenzie, 2005), leading to aluminosilicate skeleton instability. It means that the increase of sodium content reduces the mechanical strength of the system. This result is consistent with the experimental result that the compressive strengths are decreased with increasing the Na/Al ratios when the Na/Al ratio is greater than 1 (Rowles and O'connor, 2003). The stress-strain relations for different Na/Al ratios show the detailed mechanical performance during the process of tensile loading. The whole process can be divided into several stages. Taking the NASH gel with a Na/Al ratio of 1 as an example, the stress initially increases linearly in the elastic stage and subsequently slowly increases to a maximum value at a strain around 0.30 Å/Å. After the stress varies very slowly with strain from 0.3 to 0.45 Å/Å, there is a region with a dramatic drop in stress as the strain increases. Finally, the NASH gel was completely fractured at a strain of 0.8 Å/Å, corresponding to zero stress. It is obvious that the NASH gel is more prone to yield deformation as the Na/Al ratio increases. As the Na/Al ratio increase from 1 to 3 in the NASH gel, the tensile strength of the gel greatly decreases from 17.69 to 14.30 GPa accompanied by Young's modulus deduction from 71.84 to 47.46 GPa (see Table 3), and results in an obvious decrease of  ultimate tensile strength of the NASH gel. For the specimens with different Na/Al ratios, the difference in the mechanical properties can be explained by different molecular structures and reactivity mechanisms.

Molecular Local Structure Evolution
The evolution of the structure during the tensile process from a strain of 0.2 Å/Å to 0.6 Å/Å is shown in Figure 5, a preliminary qualitative illustration of the damage process in the NASH gel. The aluminosilicate skeleton is stretched broken to resist the tensile loading. In the elastic region, Si-O and Al-O bonds in the aluminosilicate skeleton, subjected to tensile loading, are elongated. The Si-O-Si and Si-O-Al angles are stretched to bear the strain in the NASH gel system. In the subsequent yield region, the Si-O and Al-O bonds break down gradually, causing morphology transformation. As shown in Figure 5, the shape of the aluminosilicate skeleton changes at the strain of 0.4 Å/Å due to the bond breakage and reconnection. Small cracks have been created during the tensile process, but the local structure rearrangement slows down its propagation (Hou et al., 2014b(Hou et al., , 2015c. In the failure stage of strain states 0.4 Å/Å to 0.6 Å/Å, the cracks created before continuing to grow and coalesce fast as strain increases in the area with the defective aluminosilicate skeleton, leading to the stress dropping rapidly (see Figures 4, 5).  It should be noted that the increase of sodium ions at the end of the crack tips prevents the Si-O-Si and Si-O-Al bonds' reconnection, and further decreases the tensile properties of the gel structure. According to Table 2, the higher the Na/Al ratio, the more Na-O bonds and the less Al-O bonds form. The Na and Al are competitive, therefore, we can conclude that the aluminosilicate skeleton is unstable after an increase in the ratio of Na/Al. For the purpose of investigating the morphology evolution quantitatively, the RDF, angle distribution, and Q species changes have been calculated to depict the structural variation of the aluminosilicate skeleton during the tensile process. The RDFs of Si-O and Al-O as a function of strain for the NASH gel with a Na/Al ratio of 1 are shown in Figures 6A,B. Based on Figures 6A,B | Evolution of Q n species (n = 0, 1, 2, 3, and 4) for the NASH gel with a Na/Al ratio of (A) 1, (B) 2, and (C) 3 during the tensile process at a strain rate of 0.08 ps.
of the structure during the tensile process ( Figure 5). In stage 2 (the failure region), the Si-O distance observed in RDF, recovered close to that of the unstrained state (1.62 Å). The trends of RDFs for Al-O are the same, and further verifies the evolution of structure during the tensile process. A similar evolution trend can be obtained in the angle distribution variation as a function of the strain. As shown in Figure 6C, the average value of the Si-O-Si angle is stretched from 139.2 • to 150.7 • in the elastic region. However, in the yield region, the angle remains basically unchanged and finally returns to the unstrained state. These are consistent with the trends of RDFs.
To depict the structural variation of the aluminosilicate skeleton, the change of Q n percentages has been calculated. Q species play an important role in the quantitative analysis of the morphological evolution. A connectivity factor, Q n (n = 0, 1, 2, 3, and 4), is a parameter that estimates the molecular skeleton connection, where n is the number of connected BO of central atoms (Si or Al). Q 0 is the monomer; Q 1 represents the dimmer structure (two connected BO of central atoms); Q 2 is the long chain; Q 3 is the branch structure; and Q 4 is the network structure (Feuston and Garofalini, 1990). When the system reaches the thermodynamic equilibrium state, the basic parameters of the Q species, BO, and NBO in different Na/Al ratio samples are shown in Table 4. It is basically consistent with the experimental results (Sadat et al., 2016). With a Na/Al ratio increase, the percentage of Q 4 species (from 81.22 to 51.25%) and BO are both decreased, implying that the network structure (aluminosilicate skeleton) depolymerizes, further resulting in a decrease of the mechanical properties. Furthermore, the aluminosilicate skeleton depolymerizes to branch structure, long chain, short chains, or monomers, so that the percentage of other phases, including Q 3 , Q 2 , and Q1 obviously increases as the sodium content increases.
During the tensile process, the Q species can depict the structural variation of the aluminosilicate skeleton quantitatively. At the beginning of the elastic stage, the aluminosilicate skeleton takes up the strain by changing the Si-O and Al-O bonds, and the Q 4 percentage remains unchanged at low strain levels for the NASH gel with all Na/Al ratios (Figure 7), implying that the structure only undergoes elastic deformation with no structural damage. As the strain increases, the Q 4 species starts to reduce at a strain of 0.06 Å/Å and continues to decrease from 81 to 39%. This means that some aluminosilicate skeletons were stretched broken, resulting in the morphological transformation from network structure to branch structure, long chain, short chains, or monomers. At the failure stage, the Q 4 species percentage decreases slightly and finally tends to become stable. Furthermore, the increase of the Q 3 and Q 2 species suggests that the network transforms to branch structures and long aluminosilicate chains. The Q 3 and Q 2 species continued to depolymerize to shorter chains and monomers, and the proportion of other Q species (Q 1 , Q 0 ) also increased later. The tensile damage of the NASH gel structures is therefore mainly attributed to the depolymerization of the aluminosilicate network. The Q 4 species percentage increased partially in the NASH gel with a Na/Al ratio of 3 at the failure stage, suggesting that the increase of sodium content rearranges the structure by electrostatic attraction of Coulomb force.

Hydrolytic Reaction Under Tensile Loading
According to previous studies, water molecules play an important role in the depolymerization of the aluminosilicate network during the tensile process, and the chemical dissociations can be observed (Hou et al., 2014a(Hou et al., , 2015c. Water molecules attack the Si-O-Si and Si-O-Al bonds form Si-OH and Al-OH groups, then further weakens the aluminosilicate skeleton during the tensile process. Therefore, the number of Si-OH, Al-OH groups, and water molecules were recorded as functions of strain to analyze the dissociation process quantitatively. The primary factor of the tensile process is the chemical effect. For the sample with a Na/Al ratio of 1 (Figure 8A), the number of water molecules decreases as the strain increases until reaching a strain of 0.5 Å/Å, and then it remains unchanged, implying that a large number of water molecules dissociate with the strain increase. The number of Si-OH and Al-OH groups increases co-instantaneous, and is close to a stable state at the strain of 0.5 Å/Å. This indicates that the hydrolysis reaction has gradually occurred in large quantities, especially after the strain is greater than 0.05 Å/Å. It is consistent with the theory that increasing the stress in the structure can reduce the energy barrier for activating the hydrolytic reaction (Zhu et al., 2005). For samples with different Na/Al ratios, the strain at which the number of water molecules starts to reduce with a Na/Al ratio increase (Figures 8B,C). In addition, at the same strain rate in the failure region, the number of water molecules of NASH gels with higher Na/Al ratios are significantly smaller than that of the lower sample. These are good validations of the inference that increasing the amount of sodium content accelerates the hydrolytic reaction and the silicate structure transformation, which can be interpreted by the attack on Si-O and Al-O bonds from sodium ions. Furthermore, the number of Al-OH groups increases significantly more than that of Si-OH as the Na/Al ratio increases. This is related to the poorer stability of Al-O bonds than Si-O bonds. The hydrolytic reaction greatly affects the silicate structure transformation during the uniaxial tension process. The reaction mechanisms of water dissociation have been investigated, and there are two reaction mechanisms: the first mechanism of reaction started in the yield stage and is listed in the following sequence ( Figure 9A): water diffuses and water molecule is adsorbed close to the neighboring aluminum atom; a bond forms between the aluminum atom and the oxygen atom of the water molecule; a five coordination aluminum structure is formed; the Al-O bond is broken and the aluminum structure restores to four coordination; the water molecule dissociation and proton transfers. In the dissociation process, water molecules lead to the formation of a five coordination aluminum structure, which is disadvantageous in terms of energy. Then the unstable structure separates from the BO. The stretched Si-O-Al bond was attacked by water molecules, and the separation of aluminosilicate network was accelerated. The other mechanism of water dissociation is illustrated in Figure 9B, the aluminosilicate skeleton was stretched broken by increasing the tensile strain. The sequence is displayed as follows: the water molecule diffuses close to the neighboring Si-O-Al bond or Si-O-Si bond; a bond is formed between the BO atom and the hydrogen atom of water molecules; the Al-O bond or Si-O bond is broken; water dissociation and proton transfers. However, it should be noted that Al-O bonds are more unstable and easily broken than Si-O bonds.
On the other hand, the dissociation process of water corresponds to the evolution of BO ( Figure 10A) and NBO ( Figure 10B) as functions of strain. During the tensile process, the percentage of BO decreases, and that of NBO increases as the Na/Al ratio increases from 1 to 3. It means that the aluminosilicate networks with a lower Na/Al ratios have a higher degree of polymerization, confirming the previous conclusion that an increase of sodium content leads to a decrease in the degree of polymerization and in the mechanical strength of system.
With an increase in the strain, the percentage change of BO and NBO shows nearly the opposite tendencies. At the beginning of the elastic region, the percentage of BO and NBO is relatively stable, that is the process of crossing the energy barrier for activating the hydrolytic reaction (Hou et al., 2019c). And then the percentage of BO drops sharply and the percentage of NBO increases up to the tensile strain of around 0.4 Å/Å. With an increase of tensile strain before failure, the percentage of BO reduces slowly and fluctuates on a small-scale, which can also be found in the opposite trend of change for NBO. This is consistent with the analysis of the hydrolysis mechanism above. Thus, the hydrolytic reaction of water molecules near the stretched Si-O and Al-O bond further accelerates the NASH gel degradation to resist the tensile loading.

CONCLUSION
Based on the molecular dynamics by a ReaxFF, the structure, reactivity, and mechanical performance of the NASH gel with a Na/Al ratio ranging from 1 to 3 have been investigated. The conclusions can be drawn as follows.
(1) For NASH gels with different Na/Al ratios, Si and Al are predominantly tetrahedral, and the tetrahedral SiO 4 structure and silicate skeleton are not sensitive to the Na/Al ratio. The sodium ions distributed in the NASH gel have around 4 ∼ 7 nearest neighbors including oxygen atoms in the silicate/aluminate tetrahedron and hydroxyl. The adsorbed water molecules dissociate into the hydroxyl group near the defective aluminosilicate structure, contributing to the formation of an aluminate pentahedron.
(2) With the amount of sodium content increasing, the molecular structure of the aluminosilicate skeleton transforms from an integrity network to partially destroyed branch structures, decreasing the strength and cohesive force of the NASH gel, characterized by the uniaxial tensile testing. It means that the connection of the aluminosilicate skeleton becomes unstable and further affects the tensile properties of the structure, due to the attack on the Si-O-Si and Si-O-Al bond by sodium ions.
(3) The aluminosilicate skeleton is stretched broken to resist the tensile loading. A high degree of polymerization of the aluminosilicate skeleton network enhances the loading resistance. As the Na/Al ratio increases, the polymerization degree of the aluminosilicate skeleton gradually decreases, and the defective aluminosilicate chains produce more NBO.
(4) Water molecules play an important role in the depolymerization of the aluminosilicate network during the tensile process, and the chemical dissociations can be observed. The adsorption and dissociation are promoted by hydrophilic NBO, further forming Si-OH and Al-OH with the surrounding aluminosilicate skeleton. And the number of Al-OH changes more than that of Si-OH at high Na/Al ratios due to the increase of the sodium content.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/supplementary material.