We have studied the dehydroxylation reaction in a smectite model, with the purpose of determining the influence of the octahedral layer composition and the possible role of the interlayer cations. Employing ab initio molecular dynamics and with the help of metadynamics for accelerating the reaction, we study the two possible mechanisms for the formation of the first water molecule. In the first mechanism, the migration of the proton across the octahedral vacancy takes place by means of a stable intermediate in which the proton is coordinated with an apical oxygen. The formation of this intermediate is not possible without a local deformation of the tetrahedral sheet, in which a Si−O basal bond is broken so that the Si coordinates with the O that releases the H, stabilizing its residual charge. In turn, the basal oxygen that loses the Si is stabilized by the nearest interlayer Na atoms. This tetrahedral deformation is permanent at this step of the reaction. In the second mechanism the proton migrates to the neighboring OH group in the same octahedral units. This takes place without reaction intermediates and it does not imply any permanent distortion in the structure, though it is observed that the environment of a tetrahedral Si spontaneously fluctuates between a tetrahedral and a ditrigonal bipyramidal arrangement.
■ INTRODUCTIONThe high temperature transformations of phyllosilicates play an important role in ceramic technology, in the design of new materials, in chemical and nuclear waste management, as well as in many geological processes such as earthquakes and vulcanism. Of special relevance is dehydroxylation, which has been linked to the phenomenology of diverse sediments overpressurization, to petroleum diffusion, to the smectite-illite transition, and to the rheology of sediments in the subduction zones. 1,2 Phyllosilicates are mainly composed of two structural units: SiO 4 tetrahedra and MO 6 octahedra, where M = Al 3+ , Fe 3+ , Mg 2+ , and Fe 2+ . The SiO 4 tetrahedra are linked and arranged in hexagons, forming the so-called tetrahedral or T-sheet. The MO 6 octahedra are edge linked, forming the octahedral or Osheet, and one distinguishes between the trioctahedral and the dioctahedral series. In the former all octahedral sites are occupied by divalent cations (mainly Mg 2+ or/and Fe 2+ ), whereas in the latter one out of three octahedra is left empty while the other two are occupied by trivalent cations (mainly Al 3+ and Fe 3+ ). T-and O-sheets are linked by sharing a plane of oxygens, called apical oxygens (O ap ). Nonapical oxygens of the O-sheet form hydroxyl groups, which can be located at the opposite sides of the octahedral vacancy or in the same side (trans-vacant and cis-vacant configurations, respectively). The T-and O-sheets can be associated in TO or TOT units, which are stacked along the c crystallographic axis leaving an interlayer space. TO and TOT phyllosilicates are also referred to as 1:1 and 2:1 phyllosilicates, respectively. Isomorphic cation substitutions are possible: for instance...