Molecular simulation and microtextural characterization of quartz dissolution in sodium hydroxide

This study uses empirical experimental evidence and Material Studio simulations to explain the interaction of sodium hydroxide (NaOH) with quartz. Density functional theory (DFT) calculations were carried out using the Cambridge Serial Total Energy Package. In addition, quartz grains subjected to dissolution in NaOH were characterized using scanning electron microscopy. The so-called O-middle termination in the quartz tetrahedron structure, typified by a solitary exposed oxygen atom at the surface, is the most susceptible SiO2 terminations to NaOH attack, as it is associated with the lowest surface energy. The adsorption energy values are − 1.44 kcal/mol and − 5.90 kcal/mol for a single atom layer and five-layered atomic structure, respectively. The DFT calculation reveals intramolecular energy is the dominant adsorption energy, followed by a weak van der Waals energy. The NaOH adsorbed on quartz (001) surface constitutes a lower band gap of 0.138 eV compared to cleaved quartz (001) surface (0.157 eV). In addition, the energy range of NaOH adsorbed on quartz is wider (− 50 to 10 eV), compared to (001) quartz (− 20 to 11 eV). The dissolved quartz showed the precipitation of sorbed silicate phases due to incongruent reactions, which indicates new voids and etch pits can be created through the cleaving of the sodium silicates sorbed into the quartz surface. The adsorption energy for NaOH interactions with reservoir sandstone was significantly higher compared to the solitary crystal grains, which can be attributed to the isotropic deformation of a single crystal, and non-uniform deformations of adjacent grains in granular quartz of sandstone reservoir. It can be inferred that exposure to NaOH will affect the structure and reactivity of quartz. The quartz surface textural study indicates that dissolution of crystalline (granite) and clastic rocks (sandstone) is critical to the development of voids, which will improve permeability by providing channels and routes for the passage of hydrothermal and reservoir fluids.


Introduction
Quartz is an abundant silicate mineral that plays a significant part in geological studies, as oxygen and silicon are the key elements of Earth's crust. Understanding the kinetics of quartz-water reaction is critical in elucidating a number of important topics in geochemistry such as soil formation and degradation, chemical weathering, enhanced porosity in reservoirs via dissolution and CO 2 sorption (White and Brantley 1995;White et al. 2001;Zhang and Liu 2014). The dissolution of quartz mineral controls the formation of voids in aquifers and sandstone reservoirs, transport and sorption of noxious waste, trace metals and CO 2 and ore formation (Brantley 2004). Quartz dissolution is crucial to chemical reactions that control the porosity, permeability and sorption capacity of sandstone reservoirs (Schott and Oelkers 1995;Lasaga and Luttge 2001;Brantley 2008). In geoscience, deciphering the reactions that control the dissolution of quartz is imperative since developing secondary pores in siliciclastics is a constraint in the exploitation of precious metals and hydrocarbon (HC) recovery (Smith and Steele 1984). The precipitates derived from pore waters following dissolution reactions may plug/block in situ pore spaces, thus impeding HC recovery efforts. Accordingly, the porosityincreasing route of dissolving labile framework grains in 1 3 sandstones is dependent on the dissolution of quartz. The structure-property of quartz is also intrinsically relevant for various technological applications such as pressure gauges, oscillators, resonators and wave stabilizers (Bosak et al. 2011).
Quartz is known to dissolve in alkali solutions such as NaOH. Hence, subjecting quartz to alkali dissolution is a viable approach to expanding existing pores and creating new ones. Moreover, the injection of alkali fluids into HC reservoirs improves recovery by lowering the interfacial tension (IFT) of oil via a blend of alkali chemicals and surfactants (Hornof et al. 2000). Alkali reaction with siliciclastics between the injection and production wells can expand the pores and improve pore connectivity along the flood path (Labrid and Duquerroix 1991;Bagci et al. 2000). NaOH reduces the interfacial tension (IFT) considerably by developing in situ surfactants (Wang et al. 2010). However, despite the several experimental studies (Rimstidt 2015; Dove 1999) on quartz dissolution in NaOH, few studies have researched molecular-level elucidation of the dissolution mechanism with regard to singular SiO 2 grains (Wang et al. 2018;Yanina et al. 2006). Deciphering this interaction using density functional theory (DFT) modelling is vital since traditional microscopic and spectroscopic analyses have not been completely successful in providing comprehensive structural properties of the quartz surface at atomic/molecular level due to their restricted resolutions or precisions. Therefore, this study is an attempt to use simulation (DFT calculations) carried out in Material Studio software to explain the interaction of NaOH with single crystalline SiO 2 via state-of-the-art ab initio calculations (Nilges et al. 2009). The surface chemical activities of quartz are closely related to its surface structure and properties (Bosak et al. 2011). Thus, studies on the surface features of cleaved quartz with the use of computational techniques unanimously agree that the (001) surface is the thermodynamically predominant plane (Dove and Rimstidt 1994;Yang and Wang 2006;Chen et al. 2008;Wander and Clark 2008). The quartz crystal will be cleaved along the 001 plane for NaOH adsorption. DFT calculations will then be carried out on the cleaved surface. This study will also provide morphological evidence of NaOH-dissolved quartz surface. In addition to the dissolution of a solitary quartz crystal, this paper provides a realistic mechanism of alkali-induced quartz dissolution, which was achieved by creating a model of reservoir sandstone that integrates quartz, feldspar, clay, oil and brine into a single unit. This approach considers reservoir minerals and fluids that were disregarded in previous simulations and experimental studies of quartz dissolution. Thus, the novelty of this research is based on simulating a polymineralic scenario (reservoir sandstone) for quartz dissolution.

Mechanism of quartz dissolution
Quartz comprises Si-O-Si linkage structures; however, the surface terminations are composed of hydrophilic hydroxyl groups (Nangia and Garrison 2008). The metal oxide is hydrolysed in the presence of H 2 O to yield the surface hydroxide ligands (-Si-OH) (Degen and Kosec 2000;Yang et al. 2004). The surface ligands may then develop into a charged one through reaction with hydroxyl (OH-) or hydronium (H + ) ions due to the amphoteric characteristic of surface reactions. The quartz reacts with water via two contrasting processes: physical adsorption (Eq. 1) and chemical dissolution (Eq. 2). Physiochemical adsorption of H 2 O molecules onto the quartz's surface via hydrogen bonding leads to the formation of stable tessellated arrangements (Goumans et al. 2007;Nangia et al. 2007). Protons are adsorbed on surface of the Si-OH under acidic conditions to form a positively charged surface (Eq. 1), whereas under alkali pH conditions, the surface of Si-OH deprotonates to induce a negatively charged surface (Eq. 2). Furthermore, Si-OH adsorption, in addition to other charges present on the metal oxide surface, produces deprotonated/hydroxylated metal species in the solution: Generally, the hydrolysis of quartz can involve four possible dissolution mechanisms (Nangia et al. 2007;Bickmore et al. 2008;Pelmenschikov et al. 2000): (a) hydrolysis of the Si centre by water (H 2 O) to form a 5-coordinated specie, accompanied by weakening of the contiguous linking/bridging bond (Fig. 1a); (b) the hydrolysis of Si centre by H 2 O and rupture of the neighbouring bridging bond, but here the reaction is catalysed by > SiO − (Fig. 1b); (c) elimination of the H + from >SiOH by OH-, thus forming H 2 O, which moves on to hydrolyse the Si centre in a way similar to the earlier scenario, i.e. mechanism 2 (Fig. 1c); and (d) protonation of a bridging O by H 3 O + , followed by the hydrolysis of a contiguous Si centre by H 2 O, formation of a 5-coordinated specie and weakening of the bridging bond (Fig. 1d). Nonetheless, all the dissolution mechanisms involve the interaction of water with quartz surface, resulting in the cleaving of the siloxane (Si-O-Si) bonds to form hydrolysed species. Therefore, surface hydroxyl is key to this quartz dissolution process (Bickmore et al. 2008;Pelmenschikov et al. 2000;Xiao and Lasaga 1996). The -OH group may be physically sorbed to several sheets of polar H 2 O molecules.
Numerous mechanisms have been developed to explain the influence of basic/alkali ions on quartz dissolution. Rate limiting hydrolysis steps reported in a related study (Bickmore et al. 2008) proposed that OH ions first hydrolyse a (1) Si centre, creating a 5-coordinated structure and weakening the neighbouring/adjacent bridging bond (Fig. 2a). Similarly, in the case of NaOH, the -OH ligand hydrolyses the Si centre after creating a 5-coordinated specie and weakening the bridging bond (Fig. 2b). Moreover, the inclusion of hydrated cations to the quartz surface either enhances the physical access of H 2 O molecules to Si-O bonds through reorientation of hydration waters to more productive sites, or transforms the mobility of connected H 2 O, or increases the incidence of effective Si-O hydrolysis reactions by enabling the discrete formation of OH ions at the interface, thereby increasing the nucleophilic nature of water (Zhang and Liu 2014). Although the development of only outer-sphere surface complexes by monovalent cations (Na + and K + ) wields no effect on the hydrolysis of Si-O-Si linkages (Wallace et al. 2010), the inner-sphere surface complexes have been affirmed to be able to modify the activation energies considerably for such hydrolysis reactions (Dove 1999). On the other hand, secondary surface complexes of alkali silicates are adsorbed on quartz surface (Ali et al. 2017), accompanied by the cleaving of these sorbed cationic secondary phases, which accounts for quartz dissolution in alkali. However, the Si-O-Si hydrolysis-initiated surface complexes/ ligands are inconsistent phenomena in dissolution.

Simulation of NaOH adsorption on quartz (001)
The adsorption of NaOH was carried out using the Accelrys Material Studio software 6.0. Density functional theory (DFT) is a potent, quantum mechanics-based tool that predicts and calculates variations in internal energy of a crystal lattice with compression under varying conditions (Eschrig 2003). The DFT calculations were performed using the Cambridge Serial Total Energy Package (CASTEP) (Clark et al. 2005) contained in Material Studio software. The software implements the density functional theory (DFT) using periodic boundary conditions and a plane-wave basis set  (Bickmore et al. 2008;Pelmenschikov et al. 2000;Xiao and Lasaga 1996) (Vanderbilt 1990). The DFT calculations utilized the Perdew-Burke-Ernzerhof function. The modelling involved several steps. The adsorbate (NaOH) was drawn, cleaned, symbolically labelled and subsequently placed in a cubic cell for geometric optimization using the Forcite module before carrying out the adsorption. Then, unit cells of quartz (SiO 2 ) were introduced separately from material studio library and subsequently cleaved to (001) plane, as illustrated in Fig. 2. Thus, the quartz (001) surface (Fig. 3a) was cleaved (Fig. 3b) to form the basis for subsequent adsorption. The unit cell of cleaved quartz was enclosed in vacuum slabs with precise volumetric proportion, which was followed by Forcite geometry optimization. Afterwards, the optimized cleaved quartz was subjected to the Adsorption Locator module to determine the most stable active adsorption sites on the optimized cleaved quartz's surface.

Surface models and surface energy calculation
The silicon atom contained in SiO 2 is sited at the tetrahedron centre of the structure with each silicon atom synchronized with four atoms of oxygen. For this study, the valence electron configurations of SiO 2 include: Si 3s 2 3p 2 , O 2s 2 2p 4 and H 1s 1 . The convergence tolerances for geometry optimization calculations include maximum displacement, maximum force, maximum energy change, maximum stress and SCF convergence tolerance of 0.002 Å, 0.05 eV Å −1 , 2.0 × 10 −5 eV atom −1 , 0.1 GPa and 2.0 × 10 −6 eV atom −1 , respectively. Prior to adsorption, the NaOH was geometrically optimized. The quartz cell parameters obtained after minimization were a = b = 4.913 Å, c = 14.76 Å. The (001) cleaved quartz plane reveals three different exposed atoms. The cleaved surface results from the breaking of Si-O bonds to produce highly reactive Si and Si-O free radicals or ligands (Fig. 3b). It is commonly acknowledged that the stability of the cleaved surface increases with lower surface energy (Qiu et al. 2017). Compared to the Si termination and O-rich termination, cleaving on the O-middle termination requires less energy, signifying that it is potentially the most stable termination model of the quartz (001) surface. The single point energy type of calculation was employed.
The next stage involved NaOH adsorption on the cleaved quartz. Given the real mineral-water interface is hydrated, H 2 O molecules were first adsorbed on the cleaved quartz surface (Fig. 4a). A supercell surface was built to accommodate the isolated H 2 O. Prior to the NaOH adsorption configuration, the amorphous cell of the hydrated quartz was transformed into a supercell and a vacuum slab was created afterwards. As illustrated in Fig. 4b, c, H 2 O was adsorbed on NaOH. The adsorbate (NaOH·H 2 O) was then imported into the vacuum slab of hydrated quartz, and the target atoms were marked. The adsorption module was then run. This procedure was repeated, but this time atomic layer number of the quartz along with the loadings of H 2 O and NaOH was increased to generate a five-layered structure. The resulting molecular models are presented in Fig. 5. This was done to determine possible variation in adsorption energy with number atomic layers.
Afterwards, DFT was also applied to simulate the adsorption energy between the molecules of the NaOH adsorbate and sandstone reservoir rock surfaces. Sandstone with porosity of 24% and composition of quartz, feldspar and kaolinite was used. The procedural steps for constructing reservoir sandstone in Material Studio are presented in Fig. 6. The sandstone components (quartz, feldspar and kaolinite) were separately constructed and then geometrically optimized (step 1). In step 2, an amorphous cell was built. The sandstone components were subsequently imported into the amorphous cell. For step 3, oil (C 8 H 18 ) and brine (H 2 O + NaCl) were incorporated into the cell to fully replicate the reservoir environment. The adsorption energies of the NaOH on the sandstone in the presence of brine and oil were measured using the adsorption configuration locator.

Experimental quartz dissolution
Laboratory experiments of quartz dissolution in NaOH involved the use of Teflon reaction vessels sheathed in ceramic bombs containing a few quartz grains (≤ 2 mm) and 50 ml of the alkali solution (NaOH). The bombs were then placed in a water bath shaker set at 100 °C for a duration of 20 h and subsequently cooled in cold running water. Afterwards, the quartz grains were removed and placed under a stream of distilled water for approximately a quarter of an hour to help dispel the adsorbed artefacts and products of the reaction. Some of the quartz grains were directly studied after drying to provide a realistic view of alkali (NaOH)-dissolved quartz. The surface microtextural features and sorbed phases of the dissolved quartz substrate were studied with the use of field emission scanning electron microscopy (Carl Zeiss Supra 55VP FESEM; operated at 5-20 kV).

Adsorption locator calculations and CASTEP simulation
CASTEP simulation was used to elucidate the action of NaOH molecules adsorbed on the optimized SiO 2 surfaces, by acquiring data on orbital overlapping, electron density, band structures, band gap, orbital and spin density of states (DOS). From the DFT calculations, the O-middle termination exhibited the lowest surface energy at 1.969 J m −2 compared to the others (surface energies of 2.896 J m −2 and 2.892 J m −2 for Si and O-rich terminations, respectively.

Adsorption of an isolated NaOH on quartz (001) surface
NaOH adsorbed on the cleaved quartz (001) surface was analysed according to the O-middle termination because it is the most stable quartz (001) termination. A supercell surface was constructed to incorporate the discrete adsorption of NaOH on quartz surfaces. The optimized geometry of NaOH adsorbed on quartz is depicted in Fig. 7, and their adsorption energies are presented in Table 1. The adsorption characteristics of NaOH on the hydrated quartz were derived by determining the preferred adsorption geometries and their equivalent adsorption energies calculated from the equation below: where E surface+NaOH is the total energy of the adsorbate-substrate system, E surf is the energy of the adsorbent (hydrated quartz) surface, and E NaOH is the energy of the adsorbate (NaOH·H 2 O). The adsorption energy values are outlined in Table 1. The calculated adsorption energy values were − 1.44 kcal/mol and − 5.90 kcal/mol for a single atom layer and five-layered atomic structure, respectively. A negative value for the absorption energy means that the corresponding adsorbed state is thermodynamically more stable E ads = E surface+NaOH − E surf + E NaOH than the unbound state. Based on the adsorption energies obtained in the DFT calculations, molecular adsorption (intramolecular energy) is favourable for isolated NaOH on the O-middle termination of a quartz (001) surface. DFT calculation shows intramolecular energy is the dominant adsorption energy, followed by weak van der Waals energy (Fig. 8). Thus, the NaOH adsorption transforms the structure, electronic properties of the quartz surface using both molecular and van der Waals energies. Accordingly, NaOH molecule adsorption is dependent on both the structure and electronic properties of the quartz.

Band structure
The locus of the Fermi level (E F ) is 0 eV. The region close to the Fermi level is crucial for different physiochemical properties that include conductor type, surface reaction activity and electrochemical potential of minerals (Escamilla-Roa et al. 2017). Figure 9 shows the reaction mechanism with and without the NaOH adsorption in terms of band structure.

3
The blue lines indicate spin lines. Figure 9a shows band structure of cleaved SiO 2 (001) surface, whereas Fig. 9b shows the band structure for NaOH adsorbed on quartz (001) surface. The cleaved SiO 2 (001) surface constitutes a higher band gap of 0.157 eV compared to 0.138 eV for NaOH adsorbed on quartz (001) surface. In addition, the energy range of NaOH adsorbed on quartz is wider (− 50 to 10 eV) compared to (001) quartz (− 20 to 11 eV). Top of the valence band and bottom of the conductive bands are occupied by Si p states. The narrow band at the bottom of Fig. 8 indicates oxygen's s state combined with Si p state. The adsorption of NaOH molecules on the SiO 2 (001) surface tightens the band gap. The wave functions of the electrons bound to the NaOH atoms create an overlap that forces the energies to  form an energy band rather than a discrete level. However, after NaOH absorption, occupancy of energy levels remains almost same because the 'doping' density of NaOH is not high enough to generate a band. Hydrogen's s state can be clearly seen at the lowest energy level, as shown in Fig. 9b.

Density of states (DOS)
Density of states (DOS) is a vital parameter for defining the states of motion of electrons in solid state physics, defining how many energy states exist at a given energy. The DOS in the conduction band depends on available states for occupation, while the valence band depends on the probability than an electron occupies an available quantum state, taking into account Pauli's exclusion principal. Figure 10 presents the partial orbital of states of sub-orbitals (s, p, d and sum) of the cleaved quartz (001) surface and NaOH adsorbed on quartz (001) surface. From Fig. 10, it is evident that the existing energy orbital states below Fermi level are ascribed to s and p orbitals. It can be deduced that the activity of electrons at p states is significantly larger compared to the s states. The p states occupy the top of the valence band and bottom of the conductance band, particularly in the region of the Fermi (EF) level. In addition, occupied energy levels for both NaOH absorbed on quartz surface and the cleaved SiO 2 are almost same. However, total orbital, as well as p orbital, DOS for NaOH adsorbed on quartz (001) surface is higher than cleaved SiO 2 (001) surface below and above Fermi level. Therefore, considering details at the vicinity of Fermi level, activity of NaOH adsorbed on quartz (001) surface is greater than that of cleaved SiO 2 (001).

Simulated NaOH interaction with reservoir sandstones
The adsorption module was then run to create a combined adsorbent-adsorbate structure (Fig. 11). The Adsorption Locator module performs a simulation of the reservoir environment, which comprises quartz substrates loaded with brine, oil, feldspars and clays of fixed composition. The adsorption energies of the different reservoir components were measured using the adsorption configuration locator. The adsorption energy calculation combines the bonding (stretching, bending, torsion energy and the diagonal and off-diagonal cross-coupling terms) and nonbonding terms. Similar to the cleaved single quartz, the dominant energies involved in NaOH interaction with the reservoir sandstones are intramolecular. However, an energy difference was observed. The intramolecular energy for adsorption increased drastically to 270 kcal/mol. The average total energy also increased significantly from 20 to 210 kcal/mol (Fig. 12). This significantly higher adsorption energy for NaOH interactions with reservoir sandstone compared to the solitary crystal grains can be attributed to the different components incorporated into the system (Aufort et al. 2015;Pertsev and Salje 2000;Zubko et al. 2013).

Stress autocorrelation function (SACF) for shear viscosity behaviour
Stress autocorrelation function (SACF) is pressure function that denotes shearing force, and thus used to explain the stresses the NaOH-dissolved quartz is subjected to due to the diffusion of Na and OH molecules, resulting in intermolecular kinetics. The shearing force refers to molecular stresses caused by attraction and repulsion of Na, OH and Si and O molecules. The SACF values are used to infer viscosity of reservoir fluids. However, as the SACF reaches stability, it begins to congregate to attain a monotonic/harmonic level, which indicates that the analysed viscosity is acceptable. The SACF was performed in the Materials Studio simulation software, at reservoir temperature of 120 °C and pressure of 70 MPa. The SACF plot for reservoir sandstone and with the addition of NaOH is shown in Fig. 13. For the ideal SiO 2 system (SACF2), the SACF remains relatively stable and undulating at zero (0) with time. With the incorporation of NaOH (SACF), the SACF slightly increases and resonates between 0 and 0.02, indicating improved stress in the system. This implies stress is imparted on the viscosity of fluids in the reservoir with the addition of NaOH. The imparted stress reduces viscosity, thereby improving fluid mobility in the reservoir.

Surface microtextural characterization of dissolved quartz
The surface of the NaOH-dissolved quartz was characterized using scanning electron microscopy (SEM). The morphological dissolution features include surface smoothening, etch pits, striations and chatter mark/curved trails. The SEM images of the in situ unwashed samples depict fibrous strands of sodium silicate, which are derived from the reaction of dissolved silica with sodium ions from the hydroxide, enmeshed into the quartz grains (Fig. 14a). These precipitated laths expand the crystal lattice of the quartz. The laths of sodium silicates are also shown as fibrous bundles that are stacked in between adjoining grains (Fig. 14b). This phase of sorbed silicates generally shows another perspective of etch pit development on quartz surface, where a new void is created through the 'cleaving off' of the sodium silicates sorbed into the quartz surface, instead of the traditional concept of dissolution of a pre-existent crystal void.
In the case of washed quartz grains, continuous elongated striations and curved trails rather than distinct etch marks can be observed on the quartz surface (Fig. 14c). The microstructures merge to such an extent that their size cannot be determined precisely. Figure 14d shows etch pits resulting from intense quartz dissolution. The walls of the etch pits will gradually advance into unreactive negative crystal faceted forms. The quartz textural disparities can be ascribed to density variations of pre-existing steps in the original crystal surface or dislocation points (Yanina et al. 2006) as well as the differential sizes of sodium silicates. The precipitation of sorbed phases suggests that the dissolution reaction needs to be incongruent (the forward reaction of dissolution should be faster than the backward reaction, i.e. disequilibria) to ensure the development of etch pits, which will serve as micro-and nanopores (Pollard et al. 1989). This phase or layer of adsorbed silicate (Fig. 14a, b) represents an intermediate species between the quartz solid and dissolved silica. The quartz surface textural study indicates that dissolution of crystalline (granite) and clastic rocks (sandstone) is critical to the development of voids, which will improve permeability by providing channels and routes for the passage of hydrothermal and reservoir fluids.

Discussion
The selectivity of NaOH absorption on quartz surface can be inferred. The NaOH adsorbate is dissociated into Na and OH (Fig. 15a), after being ionized in solution. From the dissociated adsorption model, Na prefers to adsorb on the most stable quartz (001) termination (O-middle termination) via Na-O bond at 3.828 Å, while the OH selects the Si termination via Si-OH bond at 4.605 Å (Fig. 15b). This implies that the surface of quartz is hydroxylated after adsorption of the H 2 O molecule. The bond lengths of the newly formed bonds are highlighted in green (Fig. 15b).
Previous studies reported that the O-middle termination is the most stable model of quartz (001) surface, with H 2 O adsorption on quartz achieved through van der Waals force and weak hydrogen bond (Wang et al. 2018). Although this current study established the contribution of van der Waals force, intramolecular energy was found to be the dominant energy controlling the adsorption process. Most importantly, an energy difference was observed between the adsorption of NaOH on single quartz crystal (as performed in earlier studies) and that of quartz in a polymineralic scenario (reservoir sandstone). The intramolecular energy for adsorption increased drastically to 270 kcal/mol. The average total energy also increased significantly from 20 to 210 kcal/mol. Thus, contrary to the homogenous deformation of a solitary crystal, dissolution reactions in granular quartz are affected by anisotropic deformations of contiguous grains (Aufort et al. 2015). It has been posited that such inhomogeneous deformations in polymineralic sandstone result in phase variations (Pertsev and Salje 2000), where the stepwise transition would be levelled, and the transition temperature would rise. Porosity of sandstone may, in principle, increase the energy distribution. It is likely that the combination of polarization from NaOH adsorption with strain gradient derived from motion of adjacent quartz grains can disrupt the local symmetry (Zubko et al. 2013). High porosity subsequently infers larger surface vibrations and thus higher energy distributions. The quartz grains in the sandstone are arbitrarily positioned, so the bulk material is subject to isotropic expansion that is averagely consistent with that of the singular grain. Zhang and Liu (2014) and Ali et al. (2017) reported that surface complexes and morphological features of dissolved quartz are dependent on pH and electrolyte conditions, while the computational study by Nangia and Garrison (2008) noted that quartz dissolution is controlled by the distribution of protonated, deprotonated or neutral species. The adsorption of alkali metal ions was reported to weaken the Si-O bond by reducing the reaction barrier heights thereby increasing hydrolysis (Ali et al. 2017). Rimstidt (2015) posited the integration of existing surface complexation, surface charge, terrace-ledge-kink and Lewis acid-base models into a reliable model of quartz dissolution reaction. This study takes a different approach to elucidating quartz dissolution. Based on the quartz dissolution review (Nangia et al. 2007;Bickmore et al. 2008;Pelmenschikov et al. 2000;Xiao and Lasaga 1996;Wallace et al. 2010), microtextural analysis of dissolved quartz and DFT simulation data, the mechanism of quartz dissolution is being reassessed to reflect realistic geological scenarios. Figure 16 illustrates the proposed model. First, the siloxane bonds (Si-O-Si) in quartz form hydrolysed silica species (Si-OH) when reacting with water at neutral pH [Eq. (3)]. Subsequently, continuous reaction with water causes the protonation of Si-OH to form SiOH 2 ligands on the quartz surface [Eq. (4)]. At this point, the adsorption energies from intramolecular forces facilitate the bonding of the H + ions with the silica species. As pH increases with the immersion in alkali solution (NaOH), the hydrolysed silica species begins to deprotonate due to interactions with OH − ions [Eq. (5)]. This deprotonation mechanism is enabled by desorption energies derived from endothermic reactions between the monovalent cations, hydroxyl ions and quartz surface ligands, as indicated by the high intermolecular energy of 250 kcal/mol and increase in shearing force and perturbation (i.e. increase in SACF). The deprotonation reaction reduces SiOH 2 to SiO − [Eq. (6)]. The desorbed H ions create vacancies in the initial siloxane and SiOH bonds, which are represented as the defects (etch pits) observed on the quartz surface. However, the band structure shows the occupancy of energy levels remains almost same after NaOH absorption, which can be attributed to the stable Si-O quartz bonds. Nonetheless, the slight drop in band gap from 0.157 to 0.138 eV and expanded energy range from (− 20 to 11 eV) to (− 50 to 10 eV) are suggestive of increasing interactions, density variations, improved conductivity and molecular changes of surface ligands on the quartz. Moreover, the anisotropic deformations derived in polymineralic scenario (sandstone) increase the energy in the system, which points to possible faster quartz dissolution. Thus, despite the strong covalent bonds of quartz, the formation of surface ligands due to adsorption of cations and hydroxyls on quartz and subsequent cleaving of these silica species allow for quartz dissolution. The dominant intramolecular energy indicates chemical bond formation of the hydroxide adsorbate with the oxide substrate, while the van der Waals energy stabilizes the adsorbed NaOH molecules (Dzade et al. 2014).

Conclusion
Quartz dissolution involves the weakening of siloxane (Si-O-Si) bonds by OH-ions to create defects in the quartz tetrahedron structure and form new surface ligands (SiOH, SiO − ). Density function theory (DFT) calculations on NaOH interaction with quartz showed that the O-middle termination of SiO 2 is potentially the most thermodynamically stable termination, since it has the lowest associated surface energy. Intramolecular energy was revealed to be the predominant adsorption energy in NaOH-quartz interaction. Thus, molecular NaOH adsorption accounts for the transformation of the quartz structure. The DFT calculation also indicates a considerable level of variability in intramolecular adsorption energies at different sites. The adsorption energy for NaOH interactions with reservoir sandstone was significantly higher compared to the solitary crystal grains, which can be attributed to the isotropic deformation of a single crystal, and non-uniform deformations of adjacent grains in granular quartz of sandstone. The reduced band gap with the addition of NaOH indicates improved conductivity. The laboratory dissolution experiments showed the precipitations of sorbed phases, indicating that the dissolution reaction needs to be incongruent to ensure the development of etch pits, which will serve as micro-and nanopores. Quartz dissolution is also dependent on in situ mineralogy (i.e. quartz ligands, silica flakes, dissolved silica), in addition to concentration of the alkali. The DFT calculations and experimental data clearly show that exposure to NaOH affects the structure and reactivity of quartz, although selectively.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creat iveco mmons .org/licen ses/by/4.0/.