## 2 Model and Methods
### Methods for Calculating the Elastic Constants  
Acting on external forces, the hydrated MMT system will be in a state of stress. If the system is in equilibrium, the external forces must be exactly balanced by internal stress. In general, stress is a second rank tensor with nine components as in Equation (1):  
\[\overset{\rightarrow}{\sigma}=\left[\sigma_{ij}\right]=\left[ \begin{array}{ ccc}\sigma_{xx}&\sigma_{xy}&\sigma_{xz}\\ \sigma_{xy}&\sigma_{yy}&\sigma_{yz}\\ \sigma_{xz}&\sigma_{yz}&\sigma_{zz}\end{array} \right] \tag{1}\]  
In an atomistic calculation, the internal stress tensor can be obtained using the so-called virial expression as in Equation (2):  
\[\overset{\rightarrow}{\sigma}=-\frac{1}{V_{0}}\left(\sum_{i=1}^{N}m_{i}( \overset{\rightarrow}{v}_{i}\overset{\rightarrow}{v}_{i})+\sum_{i<j} \overset{\rightarrow}{r}_{ij}\overset{\rightarrow}{F}_{ij}^{T}\right) \tag{2}\]  
where index \(i\) runs over all particles 1 through \(N\); \(m_{i}\), \(\overset{\rightarrow}{v}_{i}\), and \(\overset{\rightarrow}{F}_{i}\) denote the mass, velocity, and force acting on particle \(i\), respectively; \(V_{0}\) denotes the initial system volume. The application of stress to the hydrated MMT system results in a change in the relative positions of particles within the system expressed quantitatively via the strain tensor as in Equation (3):  
\[\overset{\rightarrow}{\varepsilon}=\left[\varepsilon_{ij}\right]=\left[ \begin{array}{ccc}\varepsilon_{xx}&\varepsilon_{xy}&\varepsilon_{xz}\\ \varepsilon_{xy}&\varepsilon_{yy}&\varepsilon_{yz}\\ \varepsilon_{xz}&\varepsilon_{yz}&\varepsilon_{zz}\end{array} \right] \tag{3}\]For a parallelepiped (for example, a periodic simulation cell) characterized by the three column vectors \([\stackrel{{\rightharpoonup}}{{a}}_{0}\stackrel{{ \rightharpoonup}}{{b}}_{0}\stackrel{{\rightharpoonup}}{{c}}_{0}]\) in some reference state, and by the vectors \([\stackrel{{\rightharpoonup}}{{a}}\stackrel{{ \rightharpoonup}}{{b}}\stackrel{{\rightharpoonup}}{{c}}]\) in the deformed state, the strain tensor is given by Equation (4):  
\[\stackrel{{\rightharpoonup}}{{\varepsilon}}=\frac{1}{2}\left[ \left(\stackrel{{\rightharpoonup}}{{h}}_{0}^{T}\right)^{-1} \stackrel{{\rightharpoonup}}{{G}}\stackrel{{ \rightharpoonup}}{{h}}_{0}^{-1}-1\right] \tag{4}\]  
where \(\stackrel{{\rightharpoonup}}{{h}}_{0}\) denotes the matrix formed from the three column vectors \([\stackrel{{\rightharpoonup}}{{a}}_{0}\stackrel{{ \rightharpoonup}}{{b}}_{0}\stackrel{{\rightharpoonup}}{{c}}_{0}]\), \(\stackrel{{\rightharpoonup}}{{h}}\) denotes the corresponding matrix forms from \([\stackrel{{\rightharpoonup}}{{a}}\stackrel{{ \rightharpoonup}}{{b}}\stackrel{{\rightharpoonup}}{{c}}]\), \(T\) denotes the matrix transpose, and \(\stackrel{{\rightharpoonup}}{{G}}\) denotes the metric tensor \(\stackrel{{\rightharpoonup}}{{h}}\stackrel{{ \rightharpoonup}}{{h}}\). Therefore, the elastic stiffness coefficients, relating the various components of stress and strain, are defined by Equation (5):  
\[\stackrel{{\rightharpoonup}}{{C}}=\left[C_{ijkl}\right]=\left. \frac{\partial\sigma_{ij}}{\partial\varepsilon_{kl}}\right|_{T,\varepsilon_{ kl}}=\left.\frac{1}{V_{0}}\frac{\partial^{2}A}{\partial\varepsilon_{ij}\partial \varepsilon_{kl}}\right|_{T,\varepsilon_{ij},\varepsilon_{kl}} \tag{5}\]  
where \(A\) denotes the Helmholtz free energy, and \(T\) denotes temperature. For small deformations, the relationship between the stresses and strains may be expressed in terms of a generalized [PERSON]'s law:  
\[\sigma_{ij}=C_{ijkl}\varepsilon_{kl} \tag{6}\]  
and, alternatively:  
\[\varepsilon_{ij}=S_{ijkl}\sigma_{kl} \tag{7}\]  
where \(C_{ijkl}\) and \(S_{ijkl}\) denote stiffness tensor and compliance tensor, respectively. Since both the stress and strain tensors are symmetric, it is often convenient to simplify these expressions by making use of Voigt vector notation. The stress is represented as in Equation (8):  
\[\left[