Methods
To determine the reliable and accurate computational method, we examined
the performance of a variety of computational methods, including
B3LYP-D3, ωB97XD, M062X and MP2, in the calculations of the bond
lengths, frequencies, and binding energies of MgH, H2,
and Mg2. The comparison with experimental
results [31, 32] is shown in
Table 1. The bond lengths of MgH, H2 and
Mg2 calculated by B3LYP-D3, ωB97XD and M062X functionals
are in good agreements with experimental values, while the MP2 method
strongly overestimates the bond lengths of Mg2. The
B3LYP-D3 overestimates the stretching frequency of Mg2.
For the binding energies of MgH, M062X reproduces the experimental value
quite well, while ωB97XD gives a higher binding energy. Since M062X
functional provides a good description of the Mg-Mg and Mg-H bonding
interaction, it will be used in the study of hydrogen desorption
reactions of MgmHn clusters.
Initially, 400 configurations of MgmHn(m = 1-6) nanoclusters with the hydrogen-hydrogen distance less than 1.0
Å were generated from 3000 configurations sampled using Molclus program
[33].
These
structures were used as initial structures for the geometry optimization
using M062X/6-311G (d, p) method implemented in Gaussian 16 [34].
For each type of MgmHn cluster, the
structure with minimum electronic energy was employed for the geometry
optimization using a higher level of M062X/def2TZVP. Frequency
calculations were performed to determine the local minimum in potential
energy surfaces. The electron localized function (ELF) were calculated
using Multiwfn 3.9 program [35]. The ab initio molecular
dynamics (AIMD) simulations for the hydrogen dissociation reactions of
MgmHn (n ≥ 2m) nanoclusters were carried
out using M062X/def2TZVP method implemented in ORCA program [36].
The simulation time was 5000 fs with a time step of 1.0 fs. The
Berendsen thermostat [37] was used to obtain a correct canonical
ensemble at 300K.