2.1. Fabrication
The pristine soil composition was listed in Table 1, which was measured by XRF (Axios, Netherlands). To remove adsorptive water, the soil (200 mesh) and Nd2O3 powder (Tianjin Kermel Co. Ltd., purity ≥ 99.99%) were preheated at 100 °C for 12 h. Four kinds of Nd2O3 addition contents were selected as 10 wt %, 15 wt %, 20 wt % and 25 wt %, respectively. The mixed soil samples were weighed at 5.0 g. Then, the pristine soil and the Nd2O3 powder were sufficiently mixed in a mortar by grinding with ethyl alcohol (AR grade) until the mixed powder dried absolutely. All mixed soil samples were sintered by microwave sintering furnace (HAMiLab-M1500) at 1000 °C、1100 °C and 1200 °C for 30 min. The samples were heated at the rate of 30 °C/min, and then cooled down to room temperature.
Table 1 Composition of the pristine soil.
|
SiO2
|
Al2O3
|
CaO
|
Fe2O3
|
MgO
|
K2O
|
Na2O
|
TiO2
|
Content (wt %)
|
66.32
|
16.57
|
4.67
|
5.87
|
1.66
|
2.86
|
0.81
|
0.74
|
2.2. Characterization
To confirm the phase structure of samples, X-ray diffraction (XRD, X’Per PRO, Netherlands) was used with Cu Kα radiation (λ = 1.5406 Å). The scanning ranged from 10o to 80o with the scanning rate of 8o/min. In addition, Fourier Transform Infrared Spectroscopy (FT-IR, PE, Spectrum One) was employed to obtain the detailed structure information. The micro-morphology and element distribution of samples were analyzed by using Scanning Electron Microscope (SEM, TM-4000, Hitachi) and energy dispersive spectrometry (EDS), respectively.
2.3. Computational methods
In MD simulations, the General Utility Lattice Program (GULP) code [20] was chosen to analyze the micro-structure of system. The 25 wt % Nd2O3 doped soil was selected. The Buckingham potential was employed to describe the inter-atom interaction, which can be expressed as follow [21]
where rij is the distance between atoms i and j, Zi and Zj are the effective charges of atoms i and j, Aij, ρij and Cij are the parameters of atom i-j interaction. The related parameters [21, 22] of Buckingham potential were shown in Table 2.
Table 2 Parameters of the inter-atomic interaction potential.
|
Aij (eV)
|
ρij (Å)
|
Cij( eV∙Å6)
|
O-O
|
9022.82
|
0.265
|
85.0924
|
Al-O
|
28538.48
|
0.172
|
34.5779
|
Si-O
|
50306.25
|
0.161
|
46.2979
|
Ca-O
|
155668.01
|
0.178
|
42.2598
|
Fe-O
|
8020.29
|
0.19
|
0.0000
|
Mg-O
|
32652.71
|
0.178
|
27.2810
|
K-O
|
2284.78
|
0.29
|
0.0000
|
Na-O
|
120304.17
|
0.17
|
0.0000
|
Ti-O
|
50126.74
|
0.178
|
46.2979
|
Nd-O
|
13084.217
|
0.255
|
0.0000
|
Around 2100 atoms were distributed randomly in simulations box, which were used as the initial modes of MD simulations. The 11 Å cutoff distance and 1 fs time step were applied in the whole simulations process. At first, the canonical (NVT) ensemble was employed to heat the randomly distributed atoms to 5000 K for 100 ps. Then the systems were cooled to 300 K with the 5 K/ps cooling speed, and the NVT ensemble was applied in the cooling process. Lastly, at 300 K, the systems were reduced the inner stress with the microcanonical (NVE) ensemble for 100 ps, and the related trajectory data were collected by the final 800 configurations for structural analysis.