6H2O (M2+= Mg, Zn and A+= Na - ACS Publications

Subscriber access provided by UNIV OF LOUISIANA

C: Surfaces, Interfaces, Porous Materials, and Catalysis

DTF Study of Layered Double Hydroxides with Cation Exchange Capacity:(A(HO))[M +AlOH) (SO)]6HO (M = Mg, Zn and A= Na, K) +

2

6

26

3

18

4

2

2

2+

+

Pedro Ivo Rodrigues Moraes, Fernando Wypych, and Alexandre Amaral Leitão J. Phys. Chem. C, Just Accepted Manuscript • Publication Date (Web): 27 Mar 2019 Downloaded from http://pubs.acs.org on March 27, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

DTF Study of Layered Double Hydroxides with Cation Exchange Capacity: 2+ (A+(H2O)6)[M2+ 6 Al3 (OH)18 (SO4 )2 ]·6H2 O (M

= Mg, Zn and A+ = Na, K) †

Pedro Ivo R. Moraes,

Fernando Wypych,

†Departamento

‡

and Alexandre A. Leitão

∗,†

de Química

Universidade Federal de Juiz de Fora Juiz de Fora, MG, CEP-36036-330, Brazil

‡Departamento

de Química

Universidade Federal do Paraná P.O. Box 19032, Curitiba, PR, CEP-81531-980, Brazil

E-mail: *[email protected] Phone: *55-32-21023310

Abstract Layered double hydroxides (LDH) equivalent to the minerals Motukoreaite and Natroglaucocerinite structures were investigated by ab initio calculations as well the capacity to exchange the cations without removing the intercalated sulphate. Due to the M2+ :M3+ molar ratio in all minerals structures, the layer domains contains three posi3+ 3+ tive charges ([M2+ 6 M3 (OH)18 ] ) which after the intercalation with two sulfate anions 3+ − build a negative charged compound ([M2+ 6 M3 (OH)18 ](SO4 )2 ) and these compounds

1

ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 22

are intercalated with hydrated sodium cations. Structures with dierent layered compositions (Mg/Al and Zn/Al) and with dierent intercalated cations (Na+ and K+ ) were constructed. In all cases the 3R polytype was the most stable. The thermodynamic of the cation exchange reaction was performed and the sodium cations were replaced by potassium, maintainig the intercalated sulfate. LDH are traditionally anion exchangers but the studied formulations can also exchange cations, which open new perspectives of applications.

Introduction Layered double hydroxides (LDH) or hydrotalcite-like compounds 1 are natural or synthetic layered materials whose structure is similar to Brucite (Mg(OH)2 ) and are very well known as 3+ anion exchangers. LDHs can be described by a generic chemical composition [M2+ 1−x Mx (OH)2 ]

[(An− )x/n ·yH2 O], where M2+ is a divalent cation that are partially replaced by M3+ trivalent cation, resulting an excess of positive charges, which need to be balanced by the intercalation of hydrated anions (An− )x/n ·yH2 O. The most common example of LDH minerals are Hydrotalcite (Mg6 Al2 (OH)16 CO3 ·4H2 O) and Pyroaurite (Mg6 Fe2 (OH)16 CO3 ·4.5H2 O), where in both cases the intercalated anion is carbonate, due to the solubility of carbon dioxide in water and formation of carbonate under strongly alkaline conditions, and the molar ratio between M2+ /M3+ is of 3:1. Other LDHs minerals with a M2+ :M3+ 2:1 molar ratio were also reported in the literature, like in Motukoreaite ((Na(H2 O)6 )[Mg6 Al3 (OH)18 (SO4 )2 ]·6H2 O), 2 Natroglaucocerinite ((Na(H2 O)6 )[Zn6 Al3 (OH)18 (SO4 )2 ]·6H2 O), 3 Shigaite ((Na(H2 O)6 )[Mn6 Al3 (OH)18 (SO4 )2 ]·6H2 O) 4 5 and Nikisherite (Na(H2 O)6 )[Fe2+ 6 Al3 (OH)18 (SO4 )2 ].6H2 O). In the last cases, instead of car-

bonate, the counter cation allocated in the interlayer domain is sulfate. All these compounds 2+ have a generic formulation (Na(H2 O)6 )[M2+ 6 Al3 (OH)18 (SO4 )2 ]·6H2 O; M = Mg, Zn, Mn and

Fe) or in a reduced way ([Na0.111 (H2 O)0.666 ]M+2 0.666 Al0.333 (OH)2 [SO4 ]0.222 ·0.666H2 O). The Shigaite structure, 4 the interlayer domain is composed of a trivalent anion [Na(SO4 )2 (H2 O)12 ]3− 2


Page 3 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


+ that compensate three edge-sharing octahedral [AlMn2+ 2 ](OH)6 ] units, where an octahedral

arrangement is formed due to the sodium cations coordinated by six water molecules. In the layer top each sulfate anions is surrounded by three sodium cations, each cation is surrounded by three sulfate anions and each sulfate is linked to one aluminum cation and the aluminum cation is coordinated octahedrally by six hydroxide anions. +3 3+ In the generic formulation [M2+ 6 M3 (OH)18 ] , the layer domains contains three positive

charges and in the case of synthetic phases, only sulfate is make available as counter ion, normally supplied by a soluble sulfate, typically Na2 SO4 . The problem is the location of the three positive charges of the layer domain by a double negative charged sulfate and 1.5 sulfate anions would be necessary. Due to the ordered positioning of the cations in the layers and consequently the local positive charge distribution, sometimes the charge of the layers cannot be directly compensated by the intercalated anions, especially with double or triple charged anions. This mismatch of the charges spatial positioning results in a net negative charge, which needs to be compensated by the intercalation of cations. To solve the problem of the layer domain neutralization, one sulfate anion and one single charged species, namely [NaSO4 ]− ·yH2 O are used and minerals like Motukoreaite, Natroglaucocerinite, Shigaite and Nikisherite are obtained. Other intercalated cations together with sulfate anions are also reported in the literature like in the Wermlandite structure ((Ca(H2 O)6 )[Mg7 Al2 (OH)18 (SO4 )2 ]·6H2 O, 6,7 but in this case the M2+ /M3+ molar ratio is of 3.5 and the intercalated cation is hydrated calcium. According to the current classication proposed by S. J. Mills and co-workers, 8 all these compounds belongs to the Wermlandite group were the basal spacing is close to 11 Å and hydrated cations occur in the presence of anions (chiey sulfate and carbonate), which are allocated between the Brucite-like layers. Another more complex structures were also reported like in Karchevskyite (Mg18 Al9 (OH)54 Sr2 (CO3 ,PO4 )9 ·6H2 O·5H2 O) 9 or synthetic mixed-valence Fe2+ ,Fe3+ -layered double hydroxide, known as green rusts (GR), a family of Fe2+ /Fe3+ LDH intercalated with hydrated anions 3



2− such as Cl− , CO2− 3 and SO4 . When GR are synthesized in the presence of NaOH, sodium

cations are intercalated between the layers to partially compensate the intercalated sulfate 3+ anion and compounds with a generic formulation (Na(H2 O)6 )[Fe2+ 6 Fe3 (SO4 )2 (OH)18 ]·6H2 O

are obtained, 10 which presents also the 3a × 3a superstructure (a=b= 9.528 Å) and consequently belongs to the Motukoreaite, Natroglaucocerinite, Shigaite and Nikisherite family. The GR are dicult to be investigated because they decompose through the contact with air into iron oxides such as magnetite (Fe3 O4 ), goethite and lepidocrocite, two FeO(OH) polymorphs. 11 In general, in the LDH having 2:1 M2+ /M3+ molar ration observed in nature, there is a √ √ ordering of the M2+ and M3+ in the structure, originating superstructure of a 3 × a 3 and/or ordering of the sodium cation in a 3a × 3a structure, which stabilizes the structure, 2,46 while in the LDH with the 3:1 M2+ /M3+ , the superstructure is of 2a × 2a type. All these mineral structures have the superstructure 3a × 3a, attributed to the ordering of sodium cations between the layers, which are bonded to sulfate anions (Shigaite: a=b= 9.512 Å; Nikisherite: a=b= 9.347 Å; Motukoreaite: a=b= 9.172 Å; Wermlandite: a=b= 9.303 Å) √ √ but in synthetic phases with similar compositions, the a 3 × a 3. superstructures were observed, attributed to the ordering of the M2+ and M3+ metals in the layers. 12 Density Functional Theory (DFT) calculations are a exceptional tool to provide reliable electronic properties, assist experimental characterization and predict properties. DFT calculations have already shown great eciency in the simulation of layered compounds. The exchange-correlation functionals together with van der Walls corrections and strongly correlated electrons corrections are archiving substantial results and agreement with experimental data. 1318 The present work aims, by means of ab initio calculations, investigate the electronic structure of the Motukoreaite and Natroglaucocerinite. The simulated 3R and 1H polytypes have the follow chemical compositions [M6 Al3 OH)18 ][Na(H2 O)6 (SO4 )2 ]·6H2 O and [M6 Al3 (OH)18 ] [K(H2 O)6 (SO4 )2 ]·6H2 O where M is Mg and Zn. For the sake of simplicity these models will 4


Page 4 of 22

Page 5 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


be called as: M-Al-(NaSO4 ) and M-Al-(KSO4 ). It is also a goal to evaluate thermodynamically a possible cation exchange. The experimental results involving the cation exchange reactions of compounds having alkaline metals and sulfate intercalated, including the phases investigated in the present manuscript, were recently investigated 12

Computational Details All ab initio calculations were performed using the codes available in the Quantum Espresso package, 19 which implements the DFT under periodic boundary conditions with plane wave functions as a basis set. 20,21 The geometry optimizations of the structures were performed with the generalized gradient approximation (GGA/PW91). 22 The ion cores were described using ultrasoft pseudopotential, 23 and the Kohn-Sham oneelectron states were expanded in a planewave basis set with a kinetic cuto energy of 60 Ry (480 Ry for the density). A Monkhorst-Pack mesh of 2x2x2 k-point sampling was used. 24 The equilibrium cell parameters and the atomic positions of all systems were found by minimizing the total energy gradient. For each set of cell parameters, the relative ion positions were relaxed until all of the force components were lower than 0.001 Ry/Bohr. All of the molecular graphics were generated by the XCRYSDEN graphical package. 25 The Motukoreaite structures were built in this work correspond to the 3x3 3R polytype supercell and Natroglaucocerinite 1H polytype as reported by Costa et al. 13 Motukoreaite and Layered Double Hydroxides (LDHs) are compounds belonging to the same family, their crystal system, class and space group are the same. The chemical formula of the Motukoreaite and Natroglaucocerinite structures were described as [M6 Al3 (OH)18 ][Na(H2 O)6 (SO4 )2 ]·6H2 O and [M6 Al3 (OH)18 ][K(H2 O)6 (SO4 )2 ]·6H2 O where M is Mg and Zn. The calculation of the thermodynamic properties like Gibbs free energy, entropy and the enthalpy of the cationic exchange reaction and Gibbs free energy increments for the studied compounds was carried out with the vibrational modes obtained from phonon calculations

5



Page 6 of 22

at the Γ-point using the harmonic approximation by the density functional perturbation theory (DFPT). 26,27 The Gibbs free energy (G) can be calculated using the following equation: G = H - TS. The enthalpy and the entropy were calculated for the Motukoreaite and Natroglaucocerinite models in the solid state by the following approximations, similarly as in previous work for other materials. 2832

Hsolid (T ) = Eele + EZP E + Evib (T )

(1)

Ssolid (T ) = Svib (T )

(2)

where Eele , EZP E , Evib (T) and Svib (T) are the total electronic energy at 0 K, the zero point energy and the vibrational contributions for the enthalpy and the entropy, respectively. The contribution of the PV term for the enthalpy of the solid in the Eq. 1 was neglected due to its small contribution. The enthalpy and the entropy contributions for the ions in the gas state were calculated by

Hgas (T ) = Eele + ZP E + Evib (T ) + Etrans (T ) + Erot (T ) + RT

(3)

Sgas (p, T ) = Svib (T ) + Strans (p, T ) + Srot (T )

(4)

where Etrans (T) and Erot (T) are the translational and rotational contributions of the enthalpy. RT is equivalent to the term PV, which is necessary to obtain the enthalpy of the gases. Strans (p, T) and Srot (T) are the translational and rotational contributions of the entropy. The EZP E , Evib (T), Erot (T), Svib (T) and Srot (T) contributions for a monoatomic ion are zero. The cations were calculated in a cubic cell 20 Å. The long-range Coulomb interactions in the periodic charged species were corrected by the Makov-Payne 33 correction. (aq)

(aq)

(aq)

The contributions of H, S, and G of the aqueous cation (HC , SC , GC , respectively) were computed by the following equations:

6


Page 7 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


C + (g) → C + (aq)

◦(hyd)

(5)

∆EA

Thus, we considered:

= (HC + ∆HC

(aq)

= (SC + ∆SC

(aq)

= (GC + ∆GC

∆SC

∆GC (g)

(g)

(g)

◦(hyd)

(aq)

∆HC

(g)

◦(hyd)

(g)

◦(hyd)

(6)

)

(7)

)

(8)

)

(g)

where HC , SC , and GC are the calculated enthalpy, the entropy, and the Gibbs free energy ◦(hyd)

of the cations in the gas phase. ∆HC

◦(hyd)

, ∆SC

◦(hyd)

, and ∆GC

are the standard enthalpy,

entropy, and Gibbs free energy of the formation of the hydrated ion C+ from its gaseous state, respectively, which were taken from experimental measurements. 34 The cationic exchange reactions considered are the following where M is Mg and Zn.

+ [M6 Al3 (OH)18 ][N a(H2 O)6 (SO4 )2 ] · 6H2 O(s) + K(aq) →

[M6 Al3 (OH)18 ][K(H2 O)6 (SO4 )2 ] · 6H2 O(s) + N a+ (aq) (R1)

Results and discussion Structural analysis The simulated 3R structures cell parameters are shown in Table 1 and the simulated and experimental Mg-Al-(NaSO4) Motukoreaite X-ray diraction patterns are shown in Figure 1,

7



Page 8 of 22

for 1H structures the cell parameters are shown in Table 2. As can be seen a good agreement between the simulated and experimental values of the cell parameters and XRD patterns was obtained. Figure 2 show the optimized structures of 3R Mg-Al-(NaSO4 ), Mg-Al-(KSO4 ) and 1H Zn-Al-(NaSO4 ), Zn-Al-(KSO4 ) structures. The sulfates anion and the cation were positioned linked to the octahedral Al cation and after the structure optimization they remained in these positions.

Table 1: Cell parameters of the simulated 3R polytype of the M-Al-(NaSO4 ) and M-Al(KSO4 ) structures where M = Mg and Zn with their respective relative errors. Parameters a (Å) b (Å) c (Å)

a Mineral

data 35

Mg-Al-(NaSO4 )a 9.17 9.17 33.51

b Synthesized

materials data 12

Mg-Al-(NaSO4 )b 33.03

Mg-Al-(KSO4 )b 32.61

Zn-Al-(NaSO4 )b 33.42

Zn-Al-(KSO4 )b 34.20

Mg-Al-(NaSO4 ) 9.31 / 1.52% 9.29 / 1.31% 33.42 / 0.27-1.18%

Mg-Al-(KSO4 ) 9.29 9.31 34.11 / 4.60%

Zn-Al-(NaSO4 ) 9.33 9.36 33.27 / 0.45%

Zn-Al-(KSO4 ) 9.34 9.38 33.96 / 0.70%

Table 2: Cell parameters of the simulated 1H polytype of the M-Al-(NaSO4 ) and M-Al(KSO4 ) structures where M = Mg and Zn with their respective relative errors Parameters a (Å) b (Å) c (Å)

a Mineral

data 3

Mg-Al-(NaSO4 )b 11.01

b Synthesized

materials data 12

Mg-Al-(KSO4 )b 10.87


Zn-Al-(NaSO4 )a 9.24 9.24 11.12


Mg-Al-(NaSO4 ) 9.30 9.30 11.20 / 1.72%

Mg-Al-(KSO4 ) 9.30 9.33 11.43 / 5.15%

Zn-Al-(NaSO4 ) 9.36 / 1.30% 9.36 / 1.30% 11.16 / 0.18-0.36%

Figure 1: XRD patterns of the Motukoreaite [Mg6 Al3 (OH)18 ][Na(H2 O)6 (SO4 )2 ]·6H2 O. The experimental XRD was simulated from mineral data taken from ref. 35 The diraction patterns show good agreement between the experimental and simulated 8


Zn-Al-(KSO4 ) 9.36 9.39 11.39 / 0.09%

Page 9 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


structures and the polytype of the compound was maintained. The higher basal spacing for K+ ion than Na+ ion would be expected because the cationic radius and the alkaline metal is surrounded octahedrally with water molecules and also with octahedral sulfate groups but experimentally this was conrmed only for Natroglaucocerinite (Na+ = 11.14 Åand K+ = 11.40 Å) but not for Motukoreaite (Na+ = 11.01 Åand K+ = 10.87 Å), in spite of having diculties to determine the values in the last case due to the diraction peaks broadening. 12 The structures in which the layer are composed by Zn (Natroglaucocerinite) showed a higher a and b parameters than the Mg (Motukoreaite) compounds, but the structures with Mg shows higher c parameter for both polytypes. The electronic total energy shows that the most stable polytype is the 3R for both compositions. The total energy dierence between 3R and 1H polytype is 16.08 and 17.41 kJ mol−1 for Mg-Al-(NaSO4 ) and Mg-Al-(KSO4 ) respectively. For the Zn structures the dierences are 17.55 and 18.84 kJ mol−1 for Zn-Al-(NaSO4 ) and Zn-Al-(KSO4 ) respectively.

Figure 2: Optimized structure of the simulated layered double hydroxydes (a) 3R Mg-Al(NaSO4 ), (b) 3R Mg-Al-(KSO4 ), (c) 1H Zn-Al-(NaSO4 ) and (d) 1H Zn-Al-(KSO4 ). Figure 3 shows the increments of ∆G for structures. As can be noticed, the 3R polytype structures are the most stable structures in all the studied temperature range, demonstrat9



ing that the vibrational contribution to the thermodynamical potentials does not change the stability. The variation of the Gibbs free energy for each temperature is approximately constant, thus the enthalpy variation is also almost constant. The increments of ∆G show that both polytypes could be accessed experimentally, since the average dierences between the 1H and 3R polytypes curves are 11.15, 18.22, 13.11 and 16.21 kJ mol−1 for Mg-Al-(NaSO4 ), Mg-Al-(KSO4 ), Zn-Al-(NaSO4 ) and Zn-Al-(KSO4 ) respectively. It is also evidenced that the 1H polytype structures reach ∆G < 0 above 15 ◦ C while the 3R polytype Gibbs free energy is less than zero throughout studied range. Although the structure of Natroglaucocerinite was indexed as 1H polytype, 3 its polytype is still doubtful in the LDH minerals community and needs to be revaluated. 8

Figure 3: Gibbs free energy increments for the M-Al-(NaSO4 ) and M-Al-(KSO4 ) structures where M = Mg and Zn.

Thermodynamic study of the cation exchange reaction The main thermodynamic potentials for the cationic exchange reaction described by the reaction R1 are shown in Table 3. The electronic total energy rules the proposed exchange cation reaction, in other words, the higher contribution for the thermodynamic partial functions is 10


Page 10 of 22

Page 11 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


the electronic total energy. Table 3: Thermodynamic Potentials for the Cation Exchange Reactions (M-Al-(NaSO4 ) → M-Al-(KSO4 )) of Motukoreaite and Natroglaucocerinite per alkaline cation at 298.15 K (kJ mol1 ). 3R Polytype Mg Zn 1H Polytype Mg Zn

∆G 103.08 102.74

∆H 95.09 93.71

T∆S -7.99 -9.03

109.05 105.45

96.73 94.86

-12.32 -10.59

∆ZPE ∆Eel -1.87 84.60 -1.34 82.93 -0.55 -1.20

85.93 84.22

The results shown in Table 3 indicate that the cation K+ can be thermodynamically replaced by Na+ . Although not favorable thermodynamically, experimental data 12 indicated that using one higher excess of K2 SO4 in the dispersion of Na-phase and stirring the dispersion for longer times, the exchange reaction can proceed until the completeness. The reaction rates of those processes play a major in the ion selectivity, however our study takes into account only thermodynamic potentials, the ion exchange kinetics is not regarded. The cation exchange rate is higher in reactions that are thermodynamically guided. 12 This behavior is observed for both polytypes. Potassium delivery is something interesting, because it is an important major nutrient for plant nutrition. 36 This nutrient is also crucial for structural integrity of cereals, photosynthesis and osmoregulation. 37 The values of ∆ZPE are particularly relevant in the computations of ∆H, the magnitude of this contribution is small for this reactions. The change in the layer divalent cation did not show inuence in those reactions, the thermodynamic potentials are very close.

Electronic Analysis The electronic analyses were done for the 3R polytype structures, by the fact that the dierent stacking, in this case, will not aect the electronic properties, due to the similarity between the structures and nally because the 3R polytype structures are the most stables. The charge density dierence plot was constructed for 3R polytype structures in order to 11



Page 12 of 22

study the intermolecular interactions between the layers and the interlayer region, anion, cation and the hydration water molecules and between the interlayer components. The charge density dierence plots for the Figures 4 and 5 is dened as:

∆ρlayer−interlayer (~r) = ρLDH (~r) − ρlayer (~r) − ρSO4 (~r) − ρX (~r) − ρH2 O (~r) − ρH2 O

r) surrounding X (~ (9)

The ρLDH , ρlayer , ρSO4 , ρX , ρH2 O and ρH2 O

surrounding X

are the charge density of the layer

double hydroxyde, of the layer, of the SO4 , of the X, where X is Na or K, of the interlayer water molecules and of the Na or K hydration water molecules, respectively. In the Figures 4 and 5 the red regions indicate the surplus of charge density, while the blue regions represent the depletion of charge density and the contour spacing is 0.004 electrons/Borh3 . The interactions between the anion and the hydroxyl layer are stronger than those between the water molecules and the hydroxyl layer. The ZnAl hydroxyl layer interacts stronger with the interlayer region than the MgAl hydroxyl layer. It can be seen in the Figure 4, comparing the electronic density in the hydroxyl groups. There is more electronic density transfer in ZnAl hydroxyl groups than MgAl layer, it is also showed in density of states calculation. Figure 5 shows the interaction between water molecules and anions, the stabilization of the cation by its hydration water molecules and the whole hydrogen bound network which stabilize the interlayer region and consequently the layered double hydroxide structure.

12


Page 13 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Figure 4: Charge density dierence of Mg-Al-(NaSO4 ), Mg-Al-(KSO4 ) (upper line) and ZnAl-(NaSO4 ), Zn-Al-(KSO4 ) (bottom line). The contour spacing is of 0.004 electrons/Bohr3 .

Figure 5: Charge density dierence focusing on interlamellar interactions of Mg-Al-(NaSO4 ), Mg-Al-(KSO4 ) (upper line) and Zn-Al-(NaSO4 ), Zn-Al-(KSO4 ) (bottom line). The contour spacing is of 0.004 electrons/Bohr3 The density of states and the projected density of states analysis, have the main objective to identify the basicity or the acidity of the structure sites. The most basic sites are those 13



in the valence band closest to the Fermi level and the most acid sites are those closest to the Fermi level in the conduction band. Figure 6 shows the PDOS of the intercalated SO2− 4 and hexahydrated cation (Na+ and K+ ). The intercalated anions and cations have the same behavior in each structure. For Mg/Al layer composition the anion electronic density in the valence band is closer to the Fermi level than for Zn/Al layer. However for Zn/Al layer composition the most basic sites are the hydroxyl groups in the layer, thus explaining the observed dierence in the electronic density transfer between the two layers composition.

Figure 6: Density of States (DOS) and projected density of states (PDOS) of the calculated 3R layered double hydroxydes with dierent lamellar composition and intercalated cation. The Bader charges analysis was performed to verify the atomic charge distribution in the Motukoreaite structure. The charge calculated using the Bader method is a good approximation to the total electronic charge of an atom. 38 The Bader charges are shown in Table 4.

14


Page 14 of 22

Page 15 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


Table 4: Bader charges of the calculated 3R polytype layer double hydroxide(e). Mg Zn Al O (layer) H (layer) SO2− 4 Na K H2 O (cation)

Mg-Al-(NaSO4 ) 1.75 3.00 -1.95 1.00 -1.59 0.89 -0.02

Zn-Al-(NaSO4 ) 1.35 3.00 -1.82 1.00 -1.59 0.89 -0.01

Mg-Al-(KSO4 ) 1.75 3.00 -1.95 1.00 -1.59 0.88 -0.02

Zn-Al-(KSO4 ) 1.35 3.00 -1.82 1.00 -1.59 0.88 -0.01

The anion charge is the same for all structures and the cations charges are similar. The main charge dierence is in the hydroxyl layer between Mg and Zn charges. This dierence agrees with the analysis made before. The hydroxyl groups are the most basic sites in the Zn/Al layer double hydroxides showed in the PDOS analysis and have greater electronic density transfer. Therefore the lower charge of the Zn atoms comparing to Mg atoms explain why this hydroxyl groups are more basic. The water molecules surrounding the cation have a slightly negative behavior, caused by the cation eletrostatic interaction.

Conclusions DFT calculations were performed to investigate the structure of the Mg-Al-(NaSO4 ), Mg-Al(KSO4 ) Zn-Al-(NaSO4 ) and Zn-Al-(KSO4 ) hydrotalcite-like compounds that have exchangeable interlayer hexahydrated monovalent alkaline cations. The simulation of the natural Motukoreaite (Mg-Al-(NaSO4 )) and Natroglaucocerinite (Zn-Al-(KSO4 )) showed a good agreement with the experimental parameters. The thermodynamic analyses indicate the polytype 3R for all four tested compositions. The main dierence in the electronic density was observed when the layer composition was changed. The hydroxyl groups in the Zn/Al layer showed a basic behavior. This is evident when we compare the Bader charges of the divalent cations.

15



As far as the authors know, this is the rst time that DFT calculations were performed using LDH as cationic exchangers. The thermodynamic properties for the exchange reaction indicated that the cations can be exchanged without removing the intercalated sulfate, as also indicated experimentally. The proposed reactions is just performed by replacing the alkaline metal in the triple charges anions between the layers [Na(SO4 )2 (H2 O)12 ]3− by other cations like K+ . The thermodynamic potentials for the exchange of Na+ by K+ show that the exchange does not occur, but the reverse reaction does, although experimentally both reactions has been performed. 12 The change in the hydroxyl layer composition did not aect these exchange reactions. The same behavior is observed with a dierent polytype. Calculations and experimental data are under way to replace both the cations and sulfate by the carbonate/bicarbonate system and will be objective of forthcoming publications.

Acknowledgement This study was nanced in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior Brasil (CAPES) - Finance Code 001 (fellowship for PIRM), CNPq (research grant for AAL and process 309729/2017-3, for FW and process 303846/2014-3 and 400117/20169), FAPEMIG (CEX APQ 02191/2017), Brazilian agencies and the enterprise PETROBRAS S/A - CENPES (AAL). We also acknowledge the CENAPAD-SP computational center for the use of its facilities by performed calculations.

References (1) Cavani, F.; Trirò, F.; Vaccari, A. Hydrotalcite-type Anionic Clays: Preparation, Properties and Applications.

Catal. Today

1991, 11, 173301.

(2) Rodgers, K. A.; Chisholm, J. E.; Davis, R. J.; Nelson, C. S. Motukoreaite, a New Hy-

16


Page 16 of 22

Page 17 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


drated Carbonate, Sulfate, and Hydroxide of Magnesium and Aluminum from Auckland, New Zealand. (3) Witzke,

T.;

Mineral. Mag.

Pöllmann,

1977, 41, 389390.

H.;

Vogel,

[Zn8x Alx (OH)16 ][(SO4 )x/2+y/2 .Nay (H2 O)6 ]. (4) Cooper,

M.

A.;

Hawthorne,

A.

Zeit. Kristal.

F.

C.

Shigaite,[AlMn+2 2 (OH)6 ]3 (SO4 )2 .Na(H2 O)6 (H2 O)6 , Can. Miner.

Struktur

synthese

von

1995, 9, 252.

The a

und

Crystal

Structure

Hydrotalcite-group

of

Mineral.

1996, 34, 9197.

(5) Huminicki, D. M. C.; Hawthornes, F. C. The Crystal Structure of Nikischerite, NaFe2+ 6 Al3 (SO4 )2 .(OH)18 .(H2 O)12 , a Mineral of the Shigaite Group. Can. Miner. 2003, , 7982.

41

(6) Moore, P. B. Wermlandite, a New Mineral from Langban, Sweden.

Lithos

1971, 4,

213217. (7) Rius,

J.;

Allmann,

[Ca(H2 O)6 .2SO4 .6H2 O]2 .

R.

Structure

Fort. Miner.

of

Wermlandites,[Mg7 (Al,Fe)2 (OH)18 ]2+

1978, 56, 113114.

(8) Mills, S. J.; Christy, A. G.; Génin, J.-M. R.; Kameda, T.; Colombo, F. Nomenclature of the Hydrotalcite Supergroup: Natural Layered Double Hydroxides.

Mineral. Mag.

2012, 76, 12891336. (9) Britvin, S. N.; Chukanov, N. V.; Bekenova, G. K.; Yagovkina, M. A.; Antonov, A. V.; Bogdanova, A. N.;

Krasnova, N. I. Karchevskyite, [Mg18 Al9 (OH)54 ][Sr2 (CO3 ,

PO4 )9 (H2 O, H3 O)11 ], a New Mineral Species of the Layered Double Hydroxide Family. Geol. Ore Deposit.

2008, 50, 556564.

(10) Christiansen, B. C.; Balic-Zunic, T.; Petit, P.-O.; Frandsen, C.; Mørup, S.; Geckeis, H.; Katerinopoulou, A.; Stipp, S. S. Composition and Structure of an Iron-bearing, Layered

17



Page 18 of 22

Double Hydroxide (LDH) Green Rust Sodium Sulphate.

Geochim. Cosmochim. Acta

2009, 73, 35793592. (11) Hansen, H. C. B. In

Layered Double Hydroxides: Present and Future

; Rives, V., Ed.;

Nova Science Publishers: New York, 2001; pp 413434. (12) Sotiles, A. R.; Baika, L. M.; Grassi, M. T.; Wypych, F. Cation Exchange Reactions in Layered Double Hydroxides Intercalated with Sulfate and Alkaline Cations 2+ (A(H2 O)6 )[M2+ = Mn, Mg, Zn; A+ = Li, Na, K). J. 6 Al3 (OH)18 (SO4 )2 ]·6H2 O (M Chem. Soc.

Am.

2019, 141, 531540.

(13) Costa, D. G.; Rocha, A. B.; Diniz, R.; Souza, W. F.; Chiaro, S. S. X.; Leitão, A. A. Structural Model Proposition and Thermodynamic and Vibrational Analysis of Hydrotalcite-Like Compounds by DFT Calculations.

J. Phys. Chem. C

2010,

,

114

1413314140. (14) Xu, S.; Yan, H.; Wei, M. Band Structure Engineering of TransitionMetalBased Layered Double Hydroxides Toward Photocatalytic Oxygen Evolution from Water: A TheoreticalExperimental Combination Study.

J. Phys. Chem. C

2017, 121, 26832695.

(15) Fonseca, C. G.; de Carvalho, G. S.; Wypych, F.; Diniz, R.; Leitão, A. A. Na+ as a Probe to Structural Investigation of Dehydrated Smectites Using NMR Spectra Calculated by DFT.

Appl. Clay Sci.

2016, 126, 132140.

(16) Nangoi, I. M.; Faro, A. C.; Souza, W. F.; Chiaro, S. S.; Leitão, A. A. Direct Comparison Among the Formation of Terephthalate and Carbonateintercalated Mg-Al-LDH: The Inuence of the High Aluminum Content.

Appl. Clay Sci.

2018, 151, 194200.

(17) Tavares, S. R.; Vaiss, V. S.; Wypych, F.; ; Leitão, A. A. Similarities Between Zinc Hydroxide Chloride Monohydrate and Its Dehydrated Form: A Theoretical Study of Their Structures and Anionic Exchange Properties.

J. Phys. Chem. C

19113. 18


2014,

, 19106

118

Page 19 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


(18) Tavares, S. R.; Vaiss, V. S.; Antunes, F. P. N.; Fonseca, C. G.; Nangoi, I. M.; Moraes, P. I. R.; Soares, C. V.; Haddad, J. F. S.; Lima, L. L.; Silva, B. N. N.; Leitão, A. A. DFT Calculations for Structural Prediction and Applications of Intercalated Lamellar Compounds.

Dalton Trans.

2018, 47, 28522866.

(19) Giannozzi, P.; Baroni, S.; Bonini, N.; Calandra, M.; Car, R.; Cavazzoni, C.; Ceresoli, D.; Chiarotti, G. L.; Cococcioni, M.; Dabo, I.; et al., QUANTUM ESPRESSO: A Modular and Open-source Software Project for Quantum Simulations of Materials. Condens. Matter

J. Phys.:

2009, 21, 395502395521.

(20) Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas.

Phys. Rev.

1964,

, B864

136

B871. (21) Kohn, W.; Sham, L. J. Self-consistent Equations Including Exchange and Correlation Eects.

Phys. Rev.

1965, 140, A1133A1138.

(22) Perdew, J. P.; Wang, Y. Accurate and Simple Analytic Representation of the Electrongas Correlation Energy.

Phys. Rev. B

1992, 45, 1324413249.

(23) Vanderbilt, D. Soft Self-Consistent Pseudopotentials in a Generalized Eingenvalue Formalism.

Phys. Rev. B

1990, 41, 78927895.

(24) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-zone Integrations. Phys. B

Rev.

1976, 13, 51885192.

(25) Kokalj, A. Computer Graphics and Graphical User Interfaces as Tools in Simulations of Matter at the Atomic Scale.

Comp. Mater. Sci.

2003, 28, 155168.

(26) Baroni, S.; Giannozzi, P.; Testa, A. Green's-Function Approach to Linear Response in Solids.

Phys. Rev. Lett.

1987, 58, 18611864.

(27) Baroni, S.; Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and Related Crystal

19



Properties from Density-Functional Perturbation Theory.

Page 20 of 22

Rev. Mod. Phys.

2001,

,

73

515562. (28) Costa, D. G.; Rocha, A. B.; Souza, W. F.; Chiaro, S. S. X.; Leitão, A. A. Comparative Structural, Thermodynamic and Electronic Analyses of Zn-Al-An− Hydrotalcite-Like − Compounds (An− = Cl− , F− , Br− , OH− , CO2− 3 or NO3 ): An Ab Initio Study. Clay Sci.

Appl.

2012, 56, 1622.

(29) Tavares, S. R.; Vaiss, V. S.; Wypych, F.; Leitão, A. A. Theoretical Study of the Anion Exchange Properties and the Thermal Decomposition of Zn5 (OH)8 (NO3 )2 ·2H2O and Zn5 (OH)8 (NO3 )2 ·2NH3 .

Appl. Clay Sci.

2015, 114, 103111.

(30) Moraes, P. I. R.; Tavares, S. R.; Vaiss, V. S.; Leitão, A. A. Ab Initio Study of Layered Double Hydroxides Containing Iron and Its Potential Use as Fertilizer. J. C

Phys. Chem.

2016, 120, 99659974.

(31) Moraes, P. I. R.; Tavares, S. R.; Vaiss, V. S.; Leitão, A. A. Investigation on Sustainable Phosphate Release in Agriculture: Structural and Thermodynamic Study of Stability, Dehydration and Anionic Exchange of Mg-Al-HPO4 Layered Double Hydroxide by DFT Calculations.

Appl. Clay Sci.

2018, 162, 428434.

(32) Fonseca, C. G.; Tavares, S. R.; Soares, C. V.; daFonseca, B. G.; Henrique, F. J.; Vaiss, V. S.; Souza, W. F.; Chiaro, S. S.; Diniz, R.; Leitão, A. A. The Role of Zn2+ Dopants in the Acid-basic Catalysis on MgO(001) Surface: Ab initio Simulations of the Dissociative Chemisorption of R-O-R` and R-S-R` (R, R`=H, CH3 , C2 H5 ). Surf.

Sci.

2017, 661, 6068. (33) Makov, G.; Payne, M. C. Periodic Boundary Conditions in Ab initio Calculations. Phys. Rev. B

1995, 51, 40144022.

(34) Lide, D. R.

CRC Handbook of Chemistry and Physics

20


; CRC Press, 2003; pp 12472.

Page 21 of 22 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60


(35) Rius, J.; Plana, F. Contribution to the Superstructure Resolution of the Double Layer Mineral Motukoreaite.

Neues Jb. Miner. Monat.

1986, 6, 263273.

(36) Timilsena, Y. P.; Adhikari, R.; Casey, P.; Muster, T.; Gill, H.; Adhikari, B. Enhanced Eciency Fertilisers: A Review of Formulation and Nutrient Release Patterns. Food Agric.

J. Sci.

2014, 95, 11311142.

(37) Pettigrew, W. T. Potassium Inuences on Yield and Quality Production for Maize, Wheat, Soybean and Cotton. (38) Bader, R. F. W.

Physiol. Plant.

2008, 133, 670681.

Atoms in Molecules: A Quantum Theory

USA, 1990.

21


; Oxford University Press,


Graphical TOC Entry

Graphic.png

22


Page 22 of 22

6H2O (M2+= Mg, Zn and A+= Na - ACS Publications

Recommend Documents