A DFT Study of Mn(Ⅱ) Adsorption on Illite Surfaces
-
摘要: 为探究Mn(Ⅱ)在伊利石表面的吸附机理，采用密度泛函理论(DFT)模拟Mn(Ⅱ)在伊利石(001)面和(010)面的吸附，研究了活性位、吸附构型、电荷和态密度。结果表明：在(001)面，Mn(Ⅱ)优先吸附于硅氧环空穴处，且与活性氧OS1形成1个共价键，吸附能为-262.55 kJ/mol；在(010)面，Mn(Ⅱ)与羟基基团的氧原子形成1~3个共价键，随着共价键数量的增加，吸附能增大，吸附的最稳定构型为Mn(Ⅱ)吸附于3个≡Al—OH基团之间的空穴处，吸附能为-533.62 kJ/mol；Mn(Ⅱ)与(001)面和(010)面均存在共价键作用和静电作用，在(001)面的吸附能小于(010)面，且与(001)面以静电作用为主，与(010)面以共价键作用为主；Mn(Ⅱ)与伊利石表面共价键的形成主要是Mn(Ⅱ)的4s轨道与表面OS的2p轨道间的相互作用。研究结果可为黏土吸附材料的开发和污染土壤的净化提供理论基础。Abstract: In order to explore the adsorption mechanism of Mn(Ⅱ) on illite surfaces, the density functional theory (DFT) was used to simulate the adsorption of Mn(Ⅱ) on the (001) and (010) surfaces. The active site, adsorption configuration, charge and state density were studied, and the following conclusions were obtained. On the (001) surface, Mn(Ⅱ) was adsorbed preferentially on the hole of silicon oxygen ring and formed a covalent bond with the active oxygen OS1, with an adsorption energy of -262.55 kJ/mol. On the (010) surface, Mn(Ⅱ) formed 1 to 3 covalent bonds with OS from hydroxyl groups, and the adsorption energy increased with the increase of covalent bonds. The most stable adsorption configuration was Mn(Ⅱ) adsorbed on the holes between three ≡Al—OH groups, and the adsorption energy was -533.62 kJ/mol. There were covalent bond interactions and electrostatic interactions between Mn(Ⅱ) and the (001) and (010) surfaces. The adsorption energy of Mn(Ⅱ) on (001) surface was less than that on the (010) surface. The adsorption of Mn(Ⅱ) on the (001) surface was mainly electrostatic interaction and that on the (010) surface was mainly covalent interaction. The formation of the surface covalent bond between Mn(Ⅱ) and illite was mainly due to the interaction between the 4s orbital of Mn and the 2p orbital of OS. The results can provide a theoretical basis for the development of clay adsorption materials or the purification of contaminated soil.
-
Keywords:
- Mn(Ⅱ) /
- illite /
- adsorption /
- density functional theory
-
重金属对水体和土壤的污染一直是环境治理的难题，黏土矿物对环境中重金属的吸附主要集中在：(1)为开发黏土矿物材料，对矿物表面进行改性以提高吸附性能[1]；(2)通过一定的物理化学方法促使重金属离子从黏土矿物表面解吸，实现土壤的净化或黏土吸附材料的回收利用[2-3]。研究内容主要探究宏观规律，然而对于更深层次的吸附机理的探讨较为匮乏。随着计算化学的发展，密度泛函理论被广泛应用于吸附微观机理的研究。郑刘春等[4]通过密度泛函理论多角度的阐述了螯合纤维对水体中Cd(Ⅱ)的吸附机理。CHEN等[5]通过DFT研究Hg(Ⅱ)在高岭石(001)面的吸附，结果表明Hg(Ⅱ)与表面O形成共价键，且主要是Hg 5d与O 2p轨道间的相互作用而成键。SANCHEZ等[6]借助量子力学理论研究金属Fe、Cu和Zn在伊利石的分布，结果表明离子的高化合价态更容易与伊利石晶体结合，形成稳定结构。可见，密度泛函理论能深入原子层面，为吸附微观机理的研究提供了可行的手段。伊利石作为一种黏土矿物，在地壳中广泛分布。为探究Mn(Ⅱ)在伊利石表面的吸附微观机理。本文以伊利石的(010)面和(001)面为研究对象，采用密度泛函理论研究Mn(Ⅱ)在伊利石表面的吸附行为，通过对吸附活性位、吸附能、吸附构型、态密度和电荷进行分析，以揭示吸附机理，研究结果可为吸附新材料的开发及土壤污染的治理提供理论基础。
1. 研究方法
1.1 表面模型建模
伊利石晶胞采用DRITS模型[7]，为满足伊利石的结构要求，构建2×1×1的超晶胞，并设置晶格取代，最终晶胞为K0.5Al2(Si4Al0.5)O10(OH)2。为研究Mn(Ⅱ)在伊利石不同晶面的微观吸附机理，构建(001)面和(010)面的表面模型，并在表面上方添加1.5 nm厚的真空层，以避免周期性的影响。对表面进行结构优化，优化后的伊利石表面模型见图 1，在图 1A伊利石(001)面中，表面结构由上下2层SiO4四面体以及中间AlO6八面体组成，即TOT结构。(001)面的Si—O环内存在Al对Si的晶格取代而使表面荷负电，因而以K+作为平衡离子、平衡电荷；在图 1B的(010)面中，裸露的Si和Al在环境中形成羟基基团，如≡Si—OH、≡Al—OH[8]。通过模型的构建，形成82个原子的(001)面和拥有94个原子的(010)面。
1.2 结构优化
采用基于第一性原理的Dmol3模块[9]对伊利石超晶胞进行结构优化，并选用广义梯度(GGA)的PBE交换相关泛函[10]，设置Grimme色散校正。采用全电子体系计算，选择双数值正极化(DNP)原子中心基组，运用BFGS优化算法，将自洽场收敛精度设定为1×10-6 eV/atom。几何优化的收敛标准为：原子的最大位移为5×10-4 nm，原子间的最大作用力为2×10-2 Ha/nm，体系总能量的变化为1×10 -5 Ha，所有计算均在倒易空间中进行，布里渊区[11]k点采样选择1×2×1。表 1是伊利石的晶格常数，模拟值与DRITS模型的实验值基本一致，表明密度泛函理论计算的参数设置具有合理性。
表 1 伊利石晶格常数Table 1. The lattice constant of illite晶格常数 模拟值 DRITS[7] a/nm 1.045 1.040 b/nm 0.898 0.895 c/nm 1.017 1.018 α/(°) 90.01 90.00 β/(°) 101.87 101.68 γ/(°) 90.00 90.00 表面模型及吸附体系优化采用与体相一致的交换关联函数和收敛标准。为避免吸附体系在真空环境下运行，选择COSMO模式，吸附环境在水介质中，介电常数为78.54，且通过晶格取代或增减平衡离子实现吸附体系的荷电中性。k点限制在г点。Mn(Ⅱ)在1.5 nm×1.5 nm ×1.5 nm的空间，k点采样限制在г点，优化标准与表面模型一致。
1.3 吸附能计算
通过吸附能的计算分析Mn(Ⅱ)在表面不同吸附位的强弱，其在伊利石表面的吸附能计算公式：
Eads =E[Mn−illite]−E[illite ]−E[Mn], (1) 其中，E[Mn-illite]表示Mn(Ⅱ)在表面吸附后体系的总能量；E[illite]为伊利石表面模型的能量；E[Mn]为Mn(Ⅱ)的能量；Eads为Mn(Ⅱ)在表面的吸附能，其中正、负分别表示吸热反应和放热反应，在放热反应中，绝对值越大表示吸附越稳定；在吸热反应中，吸附不能自发进行，吸附不稳定。
2. 结果与讨论
2.1 吸附点分布
伊利石在(001)面的硅氧环有Al取代，能吸附平衡离子，静电作用较强，分析(001)面静电势的分布规律，能查清吸附活性位。图 2为Mn(Ⅱ)与(001)面的静电势和初始吸附位分布，颜色的差异代表静电势的变化，由蓝色→白色→红色，表示静电势由负电性到正电性的转变。图 2A为(001)面的静电势分布，在硅氧环的氧原子被蓝色包裹，存在负电性，且在环中央有蓝色区域的重叠，存在较强的负电性。这与蒙脱石在(001)面的静电势基本相似，在硅氧环内存在较强的负电势[12]。图 2B为Mn(Ⅱ)的静电势，Mn(Ⅱ)被红色包裹，存在较强的正电势。根据电势的分布特征，在(001)面设置了不同Mn(Ⅱ)的吸附位，图 2C为Mn(Ⅱ)初始吸附位的分布，包括：空穴位H1~H3、表面氧原子顶位T1~T3、氧原子之间或氧与硅氧环空穴位间的桥位B1~B4。将Mn(Ⅱ)置于以上氧原子空穴位、顶位、桥位处进行模拟计算，以获得吸附稳定构型。
![]()
图 2 Mn(Ⅱ)与(001)面的静电势和初始吸附位分布Figure 2. The electrostatic potential and initial adsorption site distribution of Mn(Ⅱ) and the (001) surface伊利石(010)面有羟基基团暴露，易与重金属离子发生鳌合反应[13]。前线轨道理论表明最高占据轨道(HOMO)和最低未占据轨道(LUMO)是决定体系化学反应的关键[14]。因而，通过前线轨道分析能获得吸附位。图 3为Mn(Ⅱ)与(010)面的前线轨道和初始吸附位分布，绿色为HOMO轨道，紫色为LUMO轨道。图 3A为(010)面的前线轨道分布，羟基的氧原子被绿色包裹，主要为HOMO轨道，图 3B为Mn(Ⅱ)的前线轨道分布，Mn(Ⅱ)被紫色包裹，主要为LUMO轨道。根据前线轨道理论，HOMO轨道与LUMO轨道的结合为最佳吸附组合，因此，对(010)面的吸附位进行设置。图 3C为(010)面的Mn(Ⅱ)初始吸附位，其中羟基间的空穴位有H1~H3，羟基氧顶位T1~T5，羟基氧桥位B1~B7。模拟过程中将Mn(Ⅱ)布置在吸附位并进行结构优化，计算吸附能，以确定不同的稳定构型。
![]()
图 3 Mn(Ⅱ)与(010)面的前线轨道和初始吸附位分布Figure 3. The frontier orbital and initial adsorption site distribution of Mn(Ⅱ) and the (010) surface2.2 吸附构型分析
通过在初始吸附位进行吸附体系的模拟计算，获得吸附能及吸附构型。图 4为Mn(Ⅱ)在伊利石表面的吸附能，吸附能均为负，表明Mn(Ⅱ)在表面吸附为放热反应，吸附可以自发进行，吸附能绝对值越大，吸附越稳定。
在伊利石(001)面的吸附能均为-262.55 kJ/mol，其绝对值低于在(010)面的吸附能(-436.47~-533.62 kJ/mol)的绝对值，结合图 5中Mn(Ⅱ)在伊利石表面的吸附构型进行进一步分析。图 5A中(001)-A为在(001)面的吸附构型，虽然有不同初始吸附位，但是吸附的稳定构型只有一种，位于伊利石(001)面的硅氧环内，且与活性氧原子存在共价键。在伊利石(010)面的吸附能随着初始吸附位的不同而不同。具体表现为：(1)在图 5的D、E、F、K和M构型中，即(010)-B3、(010)-B4、(010)-B5、(010)-H3、(010)-T2吸附体系，Mn(Ⅱ)的吸附位于≡Al—OH或≡Si—OH构成的3羟基基团的空穴处，与羟基中的氧原子形成3条共价键，吸附能在-507.36~-533.62 kJ/mol范围，其中3个羟基为≡Al—OH时，基团空穴处的Mn(Ⅱ)吸附能最大，即(010)-B5构型；(2)在图 5的I、J、O和P构型中，即(010)-H1、(010)-H2、(010)-T4和(010)-T5吸附体系，Mn(Ⅱ)吸附于≡Al—OH或≡Si—OH构成的2个羟基基团的氧桥位处，且形成2条共价键，吸附能在-475.91~-496.86 kJ/mol范围，其中吸附点位于2个≡Si—OH氧桥位的吸附能最低，即(010)-T4构型；(3)在图 5的B、C、G、H、L和N构型中，即(010)-B1、(010)-B2、(010)-B6、(010)-B7、(010)-T1和(010)-T3吸附体系，Mn(Ⅱ)吸附于≡Al—OH或≡Si—OH单个羟基基团上方，且形成1个共价键，吸附能为-436.47~-475.86 kJ/mol，且在吸附位于≡Si—OH上方的吸附能最低，即(010)-B1。综上所述，在伊利石(001)面的吸附能小于在伊利石(010)面的吸附能，(010)面的吸附能随着共价键数目的增多而增大，≡Si—OH对Mn(Ⅱ)的吸附能小于≡Al—OH。
![]()
图 5 Mn(Ⅱ)在伊利石表面的吸附构型Figure 5. The adsorption configuration of Mn on the illite surface为进一步探究吸附构型，选取不同吸附面的最稳定构型进行深入分析。图 6为Mn(Ⅱ)在伊利石表面吸附的最稳定构型，表 2为吸附稳定构型参数，构型结合参数进行分析。图 6A为(001)-A构型，Mn(Ⅱ)吸附于空穴中央，与OS1形成键长为0.218 9 nm的共价键，这是由于荷电的(001)面产生极化诱导，使得活性OS1与Mn(Ⅱ)形成共价键。LI等[15]的研究说明这种共价键的强度受溶液介质的pH、电解质类型及浓度的影响。Mn(Ⅱ)与其余的氧原子OS2~OS6的距离为0.267 0~0.354 8 nm。LATA等[16]研究Mn(Ⅱ)在水溶液中形成的水化膜，Mn(Ⅱ)与O的间距为0.219~0.375 nm，均存在配位作用。因而，表面氧原子对Mn(Ⅱ)存在一定的相互作用。吸附后平衡离子K+与(001)面的距离由0.155 0 nm增加到0.181 6 nm，表明K+有远离(001)面的趋势，这主要是由于在黏土矿物的基面，存在高价离子对低价离子的离子交换吸附[17]。图 6B为(010)-B5构型，Mn(Ⅱ)吸附于3个≡Al—OH基团的空穴处，Mn(Ⅱ)与OS1、OS2、OS3分别形成键长为0.216 1、0.209 6、0.210 3 nm的共价键。根据LIU等[18]对2∶1型的黏土矿物端面羟基基团酸度的计算，(010)面的≡Al—OH和≡Si—OH的pHa分别为15.1、7.0。根据酸碱理论，≡Al—OH获得质子的能力强于≡Si—OH，SiO4四面体层的Si被Al取代，形成的≡Al—OH强化了对Mn(Ⅱ)的吸附。
![]()
图 6 Mn(Ⅱ)在伊利石表面的吸附稳定构型Figure 6. The stable adsorption configurations of Mn(Ⅱ) on the illite surface表 2 稳定吸附构型参数Table 2. The parameters of stable adsorption configurations参数 (001)-A (010)-B5 Mn—OS1的键长/nm 0.218 9 0.216 1 Mn—OS2的键长/nm 0.329 3 0.209 6 Mn—OS3的键长/nm 0.267 0 0.210 3 Mn—OS4的键长/nm 0.354 8 — Mn—OS5的键长/nm 0.292 3 — Mn—OS6的键长/nm 0.337 9 — 吸附前K+与(001)面的距离/nm 0.155 0 — 吸附后K+与(001)面的距离/nm 0.181 6 — 吸附能/(kJ·mol-1) -262.55 -533.62 Mn(Ⅱ)在(001)面和(010)面的吸附能存在较大差异，在(010)面的吸附能明显大于在(001)面的吸附能，为进一步探究吸附机理，需对稳定构型(001)-A、(010)-B5进行电荷分析和态密度分析。
2.3 电荷分析
对吸附最优构型进行电荷分析，进一步探究Mn(Ⅱ)与表面活性位间的相互作用机理。图 7为吸附稳定构型的电子密度图，颜色变化表示电子密度的大小，由红色→白色→蓝色表示电子密度由小到大。图 7A为(001)-A构型，电子密度图沿Mn(Ⅱ)与最近距离的2个表面OS原子的剖面，即OS1—Mn—OS3剖面，Mn(Ⅱ)的电子密度与OS1间有部分重叠，与OS3无重叠，表明Mn(Ⅱ)与OS1存在共价键作用；同样，在图 7B的(010)-B5构型中，Mn—OS1—OS2—OS3的电子密度剖面显示Mn(Ⅱ)与OS1、OS2和OS3均存在电子密度的重叠，表明Mn(Ⅱ)与OS1、OS2和OS3均存在共价键作用。总之，Mn(Ⅱ)与伊利石的(001)、(010)面氧原子的电子密度均存在一定程度的重叠，从而形成共价键，这与之前的成键分析结果相吻合。
表 3为吸附构型的Mulliken布居及电荷数，在(001)面Mn—OS1的布居为0.03，Mn—OS3和Mn—OS5的布居均为-0.01，在(010)面Mn—OS1、Mn—OS2和Mn—OS3的布居分别为0.18、0.27和0.24，布居越大表面共价键作用越强，因此在(001)面，Mn(Ⅱ)与表面氧间存在弱的共价键作用，在(010)面存在较强的共价键作用；Mn(Ⅱ)在(001)面的电荷数为2.32 e，在(010)面的电荷数为0.48 e，由于吸附体系为电中性体系，因而Mn(Ⅱ)与(001)面的静电作用明显强于(010)面。
表 3 稳定吸附构型的Mulliken布居及电荷Table 3. The Mulliken layout and charge of stable adsorption configurations参数 (001)-A (010)-B5 Mn—OS1的布居 0.03 0.18 Mn—OS2的布居 — 0.27 Mn—OS3的布居 -0.01 0.24 Mn—OS4的布居 — — Mn—OS5的布居 -0.01 — Mn—OS6的布居 — — Mn电荷数/e 2.32 0.48 2.4 态密度分析
态密度在相同的能量下轨道间的相互作用可表征共价键的形成[19]。图 8为在表面吸附的Mn—OS的态密度，图 8A为在(001)面Mn—OS1的态密度，OS1的2p轨道与Mn的4s轨道在-0.1~-0.48 eV成键，以OS1的2p轨道为主，在0~0.17 eV形成反键，以Mn的4s轨道为主。图 8B为在(010)面Mn—OS2的态密度，OS2的2p轨道与Mn的4s轨道在-3.3~-11 eV成键，以2p轨道为主，在0~4.2 eV形成反键，以Mn的4s轨道为主。对比吸附在(001)面和(010)面的态密度，在(001)面Mn—OS1间的态密度峰跨度小，峰形尖锐，局域性较强，在(010)面Mn—OS2间的态密度峰跨度大，峰形平缓，离域性较强，表明在(010)面Mn与OS2间的作用强于在(001)面Mn与OS1间的作用。这与键长和布居的规律一致。Mn的4s、3d轨道均穿过费米能级，表明反应活性较强，这为解吸药剂对重金属离子的捕获提供捕捉位。
3. 结论
借助密度泛函理论对Mn(Ⅱ)在伊利石表面的吸附进行模拟计算与分析，从原子层面揭示吸附机理，归纳如下：
(1) 在(001)面，Mn(Ⅱ)优先吸附于硅氧环空穴处，且与活性氧OS1形成1个共价键，吸附能为-262.55 kJ/mol；在(010)面，Mn(Ⅱ)与羟基基团形成1~3个共价键，随着共价键的增多，吸附能越大，吸附的最稳定构型为Mn(Ⅱ)吸附于3个≡Al—OH之间的空穴处，最优吸附能为-533.62 kJ/mol。
(2) Mn(Ⅱ)与(001)面和(010)面均存在共价键作用和静电作用，在(001)面的吸附能小于(010)面，且与(001)面以静电作用为主，与(010)面以共价键作用为主。
(3) Mn(Ⅱ)与伊利石表面共价键的形成主要是Mn的4s轨道与表面O的2p轨道间的相互作用。
(4) 从吸附微观机理出发，可以通过对伊利石剥片以暴露更多的吸附位，增加对Mn(Ⅱ)的吸附；可以添加电解质以降低静电作用或通过强酸以削弱共价键作用，来完成对Mn(Ⅱ)的解吸附。
-
![]()
图 2 Mn(Ⅱ)与(001)面的静电势和初始吸附位分布
Figure 2. The electrostatic potential and initial adsorption site distribution of Mn(Ⅱ) and the (001) surface
![]()
图 3 Mn(Ⅱ)与(010)面的前线轨道和初始吸附位分布
Figure 3. The frontier orbital and initial adsorption site distribution of Mn(Ⅱ) and the (010) surface
![]()
图 5 Mn(Ⅱ)在伊利石表面的吸附构型
Figure 5. The adsorption configuration of Mn on the illite surface
![]()
图 6 Mn(Ⅱ)在伊利石表面的吸附稳定构型
Figure 6. The stable adsorption configurations of Mn(Ⅱ) on the illite surface
表 1 伊利石晶格常数
Table 1 The lattice constant of illite
晶格常数 模拟值 DRITS[7] a/nm 1.045 1.040 b/nm 0.898 0.895 c/nm 1.017 1.018 α/(°) 90.01 90.00 β/(°) 101.87 101.68 γ/(°) 90.00 90.00 
下载: 导出CSV
表 2 稳定吸附构型参数
Table 2 The parameters of stable adsorption configurations
参数 (001)-A (010)-B5 Mn—OS1的键长/nm 0.218 9 0.216 1 Mn—OS2的键长/nm 0.329 3 0.209 6 Mn—OS3的键长/nm 0.267 0 0.210 3 Mn—OS4的键长/nm 0.354 8 — Mn—OS5的键长/nm 0.292 3 — Mn—OS6的键长/nm 0.337 9 — 吸附前K+与(001)面的距离/nm 0.155 0 — 吸附后K+与(001)面的距离/nm 0.181 6 — 吸附能/(kJ·mol-1) -262.55 -533.62
下载: 导出CSV
表 3 稳定吸附构型的Mulliken布居及电荷
Table 3 The Mulliken layout and charge of stable adsorption configurations
参数 (001)-A (010)-B5 Mn—OS1的布居 0.03 0.18 Mn—OS2的布居 — 0.27 Mn—OS3的布居 -0.01 0.24 Mn—OS4的布居 — — Mn—OS5的布居 -0.01 — Mn—OS6的布居 — — Mn电荷数/e 2.32 0.48
下载: 导出CSV
-
[1] UDDIN M K. A review on the adsorption of heavy metals by clay minerals with special focus on the past decade[J]. Chemical Engineering Journal, 2017, 308: 438-462. doi: 10.1016/j.cej.2016.09.029
[2] YOU X, LIU S, DAI C, et al. Effects of EDTA on adsorption of Cd(Ⅱ) and Pb(Ⅱ) by soil minerals in low-permeability layers: batch experiments and microscopic characterization[J]. Environmental Science and Pollution Research, 2020, 27: 41623-41638. doi: 10.1007/s11356-020-10149-9
[3] POCIECHA M, LESTAN D. Using electrocoagulation for metal and chelant separationfrom washing solution after EDTA leaching of Pb, Zn, and Cd contaminated soil[J]. Journal of Hazardous Materials, 2010, 174: 670-678. doi: 10.1016/j.jhazmat.2009.09.103
[4] 郑刘春, 张涛, 孟佩佩, 等. 硫化螯合秸秆纤维吸附水体中Cd(Ⅱ)行为的DFT计算[J]. 华南师范大学学报(自然科学版), 2021, 53(2): 14-20. doi: 10.6054/j.jscnun.2021021 ZHENG L C, ZHANG T, MENG P P. The DFT calculation of Cd(Ⅱ) adsorption on Sulfur-chelating-based straw cellulose[J]. Journal of South China Normal University (Natural Science Edition), 2021, 53(2): 14-20. doi: 10.6054/j.jscnun.2021021
[5] CHEN G, LI X, ZHOU L, et al. Mechanism insights into Hg(Ⅱ) adsorption on kaolinite(001) surface: a density functional study[J]. Applied Surface Science, 2019, 488: 494-502. doi: 10.1016/j.apsusc.2019.05.227
[6] SANCHEZ-CORONILLA A, MARTIN E I, FERNANDEZ F J, et al. A theoretical study on the inclusion of Fe, Cu, and Zn in illite clays[J]. Journal of Nanomaterials, 2019, 2019: 1-14.
[7] DRITS V A, ZVIAGINA B B, MCCARTY D K, et al. Factors responsible for crystal-chemical variations in the solid solutions from illite to aluminoceladonite and from glauconite to celadonite[J]. American Mineralogist, 2010, 95: 348-361. doi: 10.2138/am.2010.3300
[8] LIU X, LU X, WANG R, et al. Atomic scale structures of interfaces between kaolinite edges and water[J]. Geochimica et Cosmochimica Acta, 2012, 92: 233-242. doi: 10.1016/j.gca.2012.06.008
[9] DELLEY B. From molecules to solids with the DMol3 approach[J]. Journal of Chemical Physics, 2000, 113(18): 7756-7764. doi: 10.1063/1.1316015
[10] PERDEW J P, BURKE K, ERNZERHOF M. Generalized gradient approximation made simple[J]. Physical Review Letters, 1996, 77(18): 3865-3868. doi: 10.1103/PhysRevLett.77.3865
[11] MONKHORST H J, PACK J D. Special points for Brillon-zone integrations[J]. Physical Review B, 1976, 13(12): 5188-5192. doi: 10.1103/PhysRevB.13.5188
[12] ULIAN G N, MORO D, VALDRE G. DFT simulation of the water molecule interaction with the (001) surface of montmorillonite[J]. Minerals, 2021, 11(5): 501-511. doi: 10.3390/min11050501
[13] BERGAYA F, LAGALY G. Hand book of clay sicence[M]. London: Elsevier Science Ltd, 2013: 587-598.
[14] FUKUI K. Role of frontier orbitals in chemical reactions[J]. Science, 1982, 218: 747-754. doi: 10.1126/science.218.4574.747
[15] LI Q, SHI W, YANG Q. Polarization induced covalent bonding: a new force of heavy metal adsorption on charged particle surface[J]. Journal of Hazardous Materials, 2021, 412: 1-10.
[16] LATA K, RAMANⅡ A, SAHOO N, et al. Nature of second hydration shell in Mn2+ aquoion system Influence on proton relaxivity in nuclear magnetic resonance[J]. The Journal of Chemical Physics, 1991, 94(5): 3715-3721. doi: 10.1063/1.459742
[17] BERGAYA F, LAGALY G. Hand book of clay sicence[M]. London: Elsevier Science Ltd, 2013: 333-351.
[18] LIU X, CHENG J, SPRIK M, et al. Surface acidity of 2∶1-type dioctahedral clay minerals from first principles molecular dynamics simulations[J]. Geochimica et Cosmochimica Acta, 2014, 140: 410-417. doi: 10.1016/j.gca.2014.05.044
[19] PENG C, ZHONG Y, WANG G, et al. Atomic-level insights into the adsorption of rare earth Y(OH)n3- (n=1-3) ion[J]. Applied Surface Science, 2019, 469(1): 357-367.
-
期刊类型引用(0)
其他类型引用(1)
计量
- 文章访问数: 408
- HTML全文浏览量: 99
- PDF下载量: 75
- 被引次数: 1
