The adsorption energy of MgmHn is
defined by taking H2 and Mgm clusters as
references,
Ea(MgmHn) =
E(MgmHn) - E(Mgm) - n/2
E(H2) (1)
where E(MgmHn), E(Mgm)
represent the energy of MgmHn,
Mgm clusters, respectively. E(H2) is the
energy of H2 molecule.
For hydrogen desorption reaction, the average desorption energy per
H2 mole of MgmHn cluster
is given by
<Ed(MgmHn)>
= [E(Mgm) + n/2 E(H2) -
E(MgmHn)]/n (2)
The stepwise desorption energy (∆Ed) for
MgmHn cluster is described in relation
to the energy of MgmHn-2 of the previous
step (E(MgmHn-2)),
∆Ed(MgmHn) =
E(MgmHn-2) + E(H2) -
E(MgmHn) (3)
Note that the zero-point energy correction (ZPE) has been included in
the calculation of adsorption energy, average desorption and stepwise
desorption energies.
At each step of the hydrogen desorption reactions of
MgmHn (n = 2m) clusters, the pair of
hydrogen atoms with lowest dissociation barrier is detached from the
cluster. To compute transition state (TS), we first performed relaxed
potential energy scan using M062X/6-31G(d) method for the hydrogen pairs
with the pair distance less than 4.0 Å. Only the hydrogen pair distance
was kept constant while all other degrees of freedom were relaxed. The
hydrogen pair distance was reduced with a step of 0.1 Å until the
distance is less than 1 Å. The structure with highest energy was used as
the initial guess for the TS optimization at a high level of
M062X/def2TZVP. The frequency analysis was performed to determine the
saddle point of potential energy surface.
Results and discussion
3.1 MgmHn Cluster Structures
Figure 1 displays the structures of
MgmHn (m = 1-6, n ≥ 2m) clusters
obtained from the global configuration search. It can be seen that the
skeleton of Mgm cluster is dramatically reconstructed as
more hydrogen atoms are added. Most hydrogen atoms locate in bridging
positions between two or more Mg atoms. The saturated
Mg3H6,
Mg4H8,
Mg5H10,
Mg6H12 clusters are identical to the
structures reported in a recent study by Emmanuel [4], which
confirms the accuracy of the global search method. In addition to
saturate structures, the oversaturated
Mg3H7,
Mg4H9,
Mg5H11,
Mg6H13 (stoichiometric composition of
Mg:H > 1:2) clusters are found from the global search.
These hydrogen-enriched structures may be the intermediates in hydrogen
adsorption/desorption reaction. Note that the structures of
MgmHn (n > 2m) obtained by
geometry optimization using other functionals (B3PW91, HSE06 and PBE0)
are identical with M062X structures. Adding more hydrogen atom leads to
the decomposition of the nanoclusters
(Mg3H8,
Mg4H10,
Mg5H12,
Mg6H14) into a hydrogen molecule and
MgmH2m clusters.