The initial model described above was generated such that water was 3 weight percent of the total simulation system. To generate alternative systems with different hydration levels, the hemicellulose and lignin polymers were moved such that the center of geometry each individual molecule was multiplied by a constant factor, effectively moving each polymer away from cellulose fibril without introducing ring penetrations. Additional water molecules were added into the interstitial spaces after rewrapping the lignocellulose using the solvate plugin to VMD. By empirically varying the constant factor by which the system was expanded, we obtain models that are 0, 3, 5, 10, 15, 20, 25, or 30 percent water by weight (100% x mass water / mass total) (Figure ) and range from approximately 120,000 atoms to 190,000 atoms in size. Converting to equilibrium moisture content, these systems have an equilibrium moisture content of 0, 3. .9%, respectively. This range of hydration spans hygroscopic measurements for wood at across all humidity levels up to the fiber saturation point . To track ion diffusion within the system, ions were added by replacing water molecules within the system. Three ionization models were prepared for each level of hydration, shown in Figure . In the first case, only neutralizing sodium ions were added to each model. To directly compare against experimental findings, 100mM KCl, and CuCl2 salts were added to the hydrated neutralized models using the autoionize plugin to VMD to create two alternative systems.