Patent Application
20110038776
This invention pertains to compounds useful for solid-state storage of hydrogen. More specifically, this invention pertains to a family of new hydride compounds, XLi2MgHn, where X may be Al, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, Ge, Pd, Ag, Cd, In, Sn, Sb, Pt, Au, Hg, Tl, Pb, or Bi. Many of these hydride compounds have known non-hydrogen-containing analog compounds, XLi2Mg, and the combinations of hydrides and non-hydrogen containing compounds may provide the basis for useful hydrogen storage systems.
Considerable development effort is currently being expended on the development of hydrogen and oxygen consuming fuel cells, and there is also interest in hydrogen burning engines. Such power systems require means for storage of hydrogen fuel which hold hydrogen in a safe form at ambient conditions and which are capable of quickly receiving and releasing hydrogen. In the case of automotive vehicles, fuel storage is required to be on-board the vehicle, and storage of hydrogen gas at high pressure is generally not acceptable for such applications.
These requirements have led to the study and development of solid-state compounds for temporary storage of hydrogen, often as hydrides. For example, sodium alanate, NaAlH4, can be heated to release hydrogen gas, and a mixture of lithium amide, LiNH2, and lithium hydride, LiH, can be heated and reacted with the same effect. Despite such progress, however, no known solid-state system currently satisfies targets for on-board vehicular hydrogen storage.
Certain compounds, XLi2Mg, where X is Al, Zn, Ga, Ge, Ag, Cd, In, Sn, Sb, Au, Hg, Tl, Pb, and Bi are known to exist. It is an object of this invention to determine whether corresponding hydrides, XLi2MgHn, may be formed that could serve with such compounds in hydrogen storage combinations of hydrogen releasing and hydrogen accepting compounds.
This invention involves the use of state-of-the-art density functional theory to examine the possible formation of XLi2Mg materials and their hydrides with X being any of the elements Al, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, Ge, Pd, Ag, Cd, In, Sn, Sb, Pt, Au, Hg, Tl, Pb, and Bi (listed in order of increasing atomic number). The results indicate that previously unknown XLi2Mg compounds exist with X═Ni, Cu, Pd, Pt and that previously unknown XLi2MgHn hydrides may form for all elements X considered. This discovery provides the basis for a hydrogen storage system comprising an easy-to-use combination of a non-hydrogen-containing base compound and its hydride, its hydrogen-containing analog. Hydrogen is released from the hydride by application of heat or the like to yield a non-hydride compound in which hydrogen may later be re-absorbed.
Computational methods in chemistry, coupled with advances in affordable computing power, are now able to compute, with reasonable precision, the electronic total energies of elements and compounds. In turn, these electronic total energies may be combined to derive the enthalpy of formation of compounds from their constituent elements. Hence the feasibility of forming previously-unknown compounds, such as those with potential for controlled uptake and release of hydrogen, may be investigated computationally.
The Vienna ab initio simulation package (VASP), a state-of-the-art method implementing density functional theory, was employed with projector-augmented wave potentials constructed using the generalized gradient approximation for the exchange-correlation energy functional. Given a crystal structure, VASP computes the electronic structure, including the total electronic energy Eel.
Two template structures, 1 and
2 (
designating parent) were selected for the XLi2Mg compounds.
1 is described by the fcc F
2 is the BiF3-type (Fm
Eight XLi2MgHn templates i (
designating hydride) were constructed from two known hydride structures. Seven of these were derived from the disordered tetragonal (P4/mmm; No. 123) PdSr2LiH5 structure. Enthalpies of formation ΔH were obtained for each template structure from differences of electronic total energies:
ΔH(XLi2Mg)=Eel(XLi2Mg)−Eel(X)−2Eel(Li)−Eel(Mg) (1)
for the parent compounds, and
ΔH(XLi2MgHn)=(2/n)[Eel(XLi2MgHn)−Eel(X)−2Eel(Li)−Eel(Mg)−(n/2)Eel(H2)] (2)
for the hydrides, where n is the number of H atoms in a given configuration. Each ΔH, specified per XLi2Mg formula unit (f. u.) in Eq. (1) and per H2 molecule in Eq. (2), is the standard enthalpy of formation at zero temperature in the absence of zero point energy contributions. A negative ΔH indicates stability of the material relative to its elemental metal and molecular H2 constituents.
Thus, a group of new hydrides are provided as compounds capable of releasing hydrogen for a hydrogen-consuming device. These new compounds are XLi2MgHn, where X will typically be any one of the elements Al, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, Ge, Pd, Ag, Cd, In, Sn, Sb, Pt, Au, Hg, Tl, Pb, and Bi, and n is an integer having a value of 1 to 8 and, preferably, 4-7. In preferred hydrogen storage systems, these new hydrides, typically stored as a body of particles, release their hydrogen upon heating to yield a solid de-hydrogenated compound to which hydrogen may subsequently be restored. The original hydride may be restored by contacting the base compounds with hydrogen under suitable pressure and temperature conditions.
Other objects and advantages of the invention will be apparent from the following description of preferred embodiments.
1 template.
3 template.
State-of-the-art computational electronic structure methods implementing density functional theory (DFT) have been employed with substantial success to model hydride properties, including the crucial enthalpies of formation. That success encourages the development of strategies for harnessing these calculation tools to guide the discovery of novel hydrides. The approach in this case is to choose a related compound having a known crystal structure and calculate enthalpies of formation for isostructural, hypothetical compounds constructed by elemental replacement. In a further step, a parallel process is followed for a hypothetical hydride derived from the hypothetical compound through the addition of hydrogen to the original lattice.
The goal is to identify compounds which may take up and release hydrogen in a reversible manner. Thus in pursuing computational approaches using isostructural templates to guide the calculation, at least one template, for the compound itself, is required. This represents the simplest outcome and physically corresponds to a situation where the template structure of the compound is sufficiently open to accommodate hydrogen without distorting the structure. If this situation does not arise, then one template structure will be required for the compound and another for the hydride.
In practice of this invention it has been found that even this situation is inadequate to fully capture the complexities of the reaction. First, more than one compound template is required and these will be designated as 1 and
2 and then eight template structures, designated as
1-
8, are required for the hydride. Details are provided in subsequent sections but it should be emphasized that the structural choices are not arbitrary but are guided and informed by the known behavior of either representative examples of the family of compounds or by knowledge of analogous compounds.
For the compounds, XLi2Mg, two template structures are considered, both face-centered cubic. The choice is based on crystallographic information available for known XLi2Mg ternary compounds which include Ag, Al, Au, Bi, Cd, Ga, Ge, Hg, In, Pb, Sb, Sn, Tl and Zn. Of these, all are disordered with the exception of PbLi2Mg and an ordered variant of TlLi2Mg. The disordered structures possess 4 Li and 4 Mg, or 4 X and 4 Mg atoms on an 8-fold site in the conventional cell of either the BiF3-type (Fm
These four ordered structures can all be described by the fcc F 1 whose structure is shown in
1 template form. The computed X-ray powder diffraction patterns are virtually indistinguishable, demonstrating that the
1 template form is an excellent approximation to the observed disordered structure.
The ordered structures represented by PbLi2Mg and ordered TlLi2Mg are both BiF3-type (Fm 2.
There are multiple available sites in the 1 template structure for incorporation of hydrogen atoms. To identify a structural template for the hydride, preliminary computations were performed to assess the stability of hydrides based on the structure of AgLi2Mg. Various numbers of hydrogen atoms were inserted at various locations in the
1 lattice and the stability of the resulting hydrides AgLi2MgHn evaluated. The results indicated that stable hydrides would not form and thus that this structure is not a suitable template for the hydride.
Consideration was then given to structures of lower symmetry. Eight XLi2MgHn templates i (
designating hydride) were constructed from two known hydride structures having lower lattice symmetries. Seven of these were derived from the disordered tetragonal (P4/mmm; No. 123) PdSr2LiH5 structure and satisfied the cases when n, the number of hydrogen atoms, was 4, 5 or 6. For the case where the case of 7 hydrogen atoms the template structure was derived from the ordered hexagonal (P63/mmc; No. 194) RuMg2LiH7 structure.
To progress from the known PdSr2LiH5 to the desired XLi2MgHn structure the following atom substitutions are first made to transform the parent PdSr2Li atom placement into an atom placement for the desired parent XLi2Mg. X is substituted for Pd; Li is substituted for Sr; and Mg is substituted for Li. This locates X, Li, and Mg on the 1a, 2h, and 1b sites in P4/mmm, respectively. To progress from the known RuMg2Li atom placement to the desired XLi2Mg atom placement, X is substituted for Ru; Mg is substituted for Li; and Li is substituted for Mg.
In the hydride PdSr2LiH5 the H atoms occupy the 2e, 2g, and one of the 2f sites. Seven XLi2MgHn templates with n=4, 5, and 6 were generated by various fillings of these, with the eighth template based on RuMg2LiH7. The templates are:
1-XLi2MgH4 with four H atoms on the 2e and 2f sites. The P4/mmm space group symmetry (of PdSr2LiH5) is preserved.
2-XLi2MgH4 with four H atoms on the 2e and 2g sites. The P4/mmm space group symmetry (of PdSr2LiH5) is maintained.
3-XLi2MgH4 with four H atoms on the 2f and 2g sites. The structure remains P4/mmm. This structure is shown in
4-XLi2MgH5 with five H atoms occupying one of the 2e and all 2f, 2g sites. This structure is equivalent to an ordered orthorhombic Pmmm (No. 47) lattice with X (1a), Li (2l), Mg (1b), and H (1c, 1e, 1f, and 2i) sites.
5-XLi2MgH5 with five H atoms occupying one of the 2f and all 2e, 2g sites. This structure can also be described as ordered orthorhombic Pmmm with X (1a), Li (2l), Mg (1b), and H (1c, 1d, 1f, and 2i) sites.
6-XLi2MgH5 with five H atoms occupying one of the 2g and all 2e, 2f sites. This structure is equivalent to an ordered tetragonal P4m (No. 99) lattice with X (1a1), Li (1b1, 1b2), Mg (1a2), and H (1a3, 2c1, and 2c2) sites.
7-XLi2MgH6 with six H atoms on the 2e, 2f, and 2g sites. The space group symmetry is again P4/mmm.
8-XLi2MgH7 ordered hexagonal (P63/mmc; No. 194). The populated sites are X (2a), Li (4f), Mg (2b), and H (2c, 12k).
Since the purpose of the calculations is to identify previously unknown compounds and hydrides, there can be no foreknowledge of which compounds or hydrides are stable, and if stable, what structure they will adopt. Hence calculations must be made for both parent structures and for all eight hydride structures.
Electronic total energies Eel were calculated with the Vienna ab initio simulation package (VASP), which implements DFT using a plane wave basis set. Projector-augmented wave potentials were employed for the elemental constituents, and the generalized gradient approximation (GGA) of Perdew and Wang was used for the exchange-correlation energy functional μxc. Non-magnetic calculations were performed for all materials with the exception of those with X═Fe, Co, Ni, and Pd, for which spin-polarized calculations were done to assess the possibility of magnetic states In all cases a 900 eV plane wave cutoff energy was imposed. The number of points in the irreducible Brillouin zone for the k-space meshes utilized was 60 (1,
2), 125 (
4,
5), 126 (
1,
2,
3,
6,
7), and 133 (
8). At least two simultaneous relaxations of the lattice constants and nuclear coordinates not fixed by the space group were carried out. The electronic total energies and forces were converged to 10−6 eV/cell and 10−3 eV/Å, respectively. Calculations for the H2 molecule and the elemental metals were performed with the same computational machinery to the same levels of precision.
Enthalpies of formation ΔH were obtained from differences of electronic total energies:
ΔH(XLi2Mg)=Eel(XLi2Mg)−Eel(X)−2Eel(Li)−Eel(Mg)
for the parent compounds, and
ΔH(XLi2MgHn)=(2/n)[Eel(XLi2MgHn)−Eel(X)−2Eel(Li)−Eel(Mg)−(n/2)Eel(H2)]
for the hydrides, where Eel(Y) is the electronic total energy of constituent Y and n is the number of H atoms in a given configuration. Each ΔH, specified per XLi2Mg formula unit (f. u.) in Eq. (1) and per H2 molecule in Eq. (2), is the standard enthalpy of formation at zero temperature in the absence of zero point energy contributions.
1
2
AlLi
2
Mg
ZnLi
2
Mg
GaLi
2
Mg
GeLi
2
Mg
AgLi
2
Mg
CdLi
2
Mg
SnLi
2
Mg
SbLi
2
Mg
AuLi
2
Mg
HgLi
2
Mg
TlLi
2
Mg
PbLi
2
Mg
BiLi
2
Mg
A negative ΔH indicates stability of the material relative to its elemental metal and molecular H2 constituents. Table I lists ΔH(XLi2Mg) calculated according to Eq. (1) for all known XLi2Mg compounds as well as for those with X═Sc—Cu, Pd, and Pt for structures based on both the 1 and
2 templates. In qualitative agreement with experiment, ΔH is negative for all the materials reported to form (shown in underlined bold in Table I) using either template. In these cases the lower ΔH value corresponds to the observed structure except for PbLi2Mg, TlLi2Mg and BiLi2Mg. For PbLi2Mg and TlLi2Mg the energy differences are modest but the energy difference for the BiLi2Mg is relatively large and may merit additional review of the experimental diffraction analyses.
For the group of compounds not previously reported, ΔH is large and positive for XLi2Mg with X═Sc, Ti, V, Cr, Mn, Fe, and Co, indicating that the compounds likely do not exist. For NiLi2Mg, CuLi2Mg, PdLi2Mg, and PtLi2Mg, however, negative ΔH values are obtained, raising the possibility that these materials form. ΔH is especially large and negative for the Pd and Pt compounds in either the 1 or
2 template.
Turning now to the hydrides, enthalpies of formation relative to the elemental metals and H2 for the XLi2MgHn hydrides in the eight i template structures are presented in Table II. At least one, and in most cases several, negative ΔH values are obtained for every X, raising the possibility of hydride formation in each case. For each X the lowest ΔH value is shown in bold underline and Table III summarizes the minimum ΔH results from Table II and includes the hydrogen mass percentage for each hydride.
i (n)
1 (4)
−20
2 (4)
3 (4)
−33
−62
−87
4 (5)
5 (5)
6 (5)
−32
−86
−51
−24
−15
7 (6)
−51
−27
−12
8 (7)
−89
i (n)
1 (4)
−98
−112
2 (4)
−14
−26
−17
−64
−28
−30
−32
−22
3 (4)
4 (5)
5 (5)
6 (5)
−35
−28
7 (6)
8 (7)
i)
6)
6)
7)
7)
7)
3)
3)
8)
3)
6)
6)
6)
1)
1)
6)
2)
6)
2)
2)
1)
2)
2)
2)
2)
2)
According to the van't Hoff relation
ln p/p0=ΔH/RT−ΔS/R,
where ΔS is the entropy of formation and R the gas constant, the configuration having the most negative ΔH is that which is stable at the lowest H2 pressure p.
For controlled and repeated hydrogen uptake and release, it is more desirable to have XLi2MgHn release hydrogen while reverting to XLi2Mg than to have XLi2MgHn revert to its constituent elements while releasing hydrogen. The preferred reversible reaction is then:
XLi2MgHnXLi2Mg+n/2H2,
and the formation enthalpy of the hydride with respect to its parent compound, ΔH*(XLi2MgHn), is given by:
ΔH*(XLi2MgHn)≡(2/n)[Eel(XLi2MgHn)−Eel(XLi2Mg)−(n/2)Eel(H2)].
Listings of ΔH*(XLi2MgHn) are also given in Table III. For AlLi2MgH5, NiLi2MgH4, CuLi2MgH5, ZnLi2MgH5, PdLi2MgH4, AgLi2MgH5, and PtLi2MgH4 ΔH*(XLi2MgHn)<0 in Table III and a stable parent is either known (AlLi2Mg, ZnLi2Mg, AgLi2Mg) or predicted (NiLi2Mg, CuLi2Mg, PdLi2Mg, PtLi2Mg) to exist (cf. Table I) suggesting the desirable possibility of cycling between the two according to the reaction given above.
AgLi2MgH5 (6) is the only Ag hydride with ΔH*<0. For X═Al, Ni, Cu, Zn, Pd, and Pt, all of which also feature stable parent compounds in Table I, it is important to observe that there are at least two hydrides having ΔH*(XLi2MgHn)<0. The hydride with the minimum ΔH* (and hence the most stable with respect to its parent) is listed in Table IV. Only the X═Ni, Cu, and Ag entries are common to Table III and Table IV; for the others (X═Al, Zn, Pd, Pt) the hydride most stable with respect to the elemental components differs from that which is most stable relative to the parent, with the latter having a larger hydrogen content n. The materials in Table IV are the preferred members of the XLi2Mg/XLi2MgHn class of compounds for reversible hydrogen storage.
i)
8)
3)
6)
7)
6)
6)
7)
Thus, a family of hydrides is provided as follows: AlLi2MgHn, ScLi2MgHn, TiLi2MgHn, VLi2MgHn, CrLi2MgHn, MnLi2MgHn, FeLi2MgHn, CoLi2MgHn, NiLi2MgHn, CuLi2MgHn, ZnLi2MgHn, GaLi2MgHn, GeLi2MgHn, PdLi2MgHn, AgLi2MgHn, CdLi2MgHn, InLi2MgHn, SnLi2MgHn, SbLi2MgHn, PtLi2MgHn, AuLi2MgHn, HgLi2MgHn, TlLi2MgHn, PbLi2MgHn, and BiLi2MgHn.
In these hydride formulas, n is an integer having a value from 4 to 7.
Particles of one or more of these compositions may be used in a suitable, predetermined mass to release a desired quantity of hydrogen upon heating for delivery to a hydrogen-consuming device. Upon release of hydrogen they may revert to like hydrogen-depleted analogs. Hydrogen may be reacted with the analogs to restore the hydrides for re-use.
