Outlook and challenges for hydrogen storage in nanoporous materials

Considerable progress has been made recently in the use of nanoporous materials for hydrogen storage. In this article, the current status of the ﬁeld and future challenges are discussed, ranging from important open fundamental questions, such as the density and volume of the adsorbed phase and its relationship to overall storage capacity, to the development of new functional materials and complete storage system design. With regard to fun-damentals, the use of neutron scattering to study adsorbed H 2 , suitable adsorption isotherm equations, and the accurate computational modelling and simulation of H 2 adsorption are discussed. The new materials covered include ﬂexible metal–organic frameworks, core–shell materials, and porous organic cage compounds. The article concludes with a discussion of the experimental investigation of real adsorptive hydrogen storage tanks, the improvement in the thermal conductivity of storage beds, and new storage system concepts and designs.


Introduction
Solid-state hydrogen storage offers the promise of improving upon the conventional forms of hydrogen storage technology, namely liquid or compressed gas. Liquid H 2 storage requires temperatures below *30 K, while compressed gas storage requires high pressures, up to 70 MPa, to achieve practical storage densities. In contrast, the use of hydrogen storage materials can lead to high H 2 storage densities well above 30 K and at lower pressures than those required for compressed gas. The aim of hydrogen storage materials research is thus to develop an effective H 2 storage method that can operate at near-ambient temperatures and at pressures below 10 MPa.
The materials currently being considered can be generally separated into hydrides and nanoporous materials [1,2], the latter of which physically adsorb hydrogen molecules (H 2 ) on their surface via van der Waals and electrostatic forces. Nanoporous materials have several advantages. For example, the adsorption process is completely reversible and the kinetics of adsorption are rapid [3,4], in contrast to hydrogen absorption by many metal and complex hydrides. The low enthalpy of adsorption also reduces thermal management issues. Together with its reversibility, the physical adsorption of H 2 does not typically induce crystallographic phase changes in the host material; so material stability during repeated cycling of hydrogen is less of an issue.
On the other hand, the use of adsorption has drawbacks. For example, the physical interaction of H 2 with solid surfaces is quite weak because H 2 has no charge and no dipole moment, a relatively weak quadrupole moment, and a low polarisability [5,6]. High storage capacities can thus only be achieved at low temperatures. Furthermore, the physical adsorption of H 2 is a surface process, so the advantage of the presence of a nanoporous material in a tank is provided only by the enhanced density of the adsorbed H 2 in the pores. This limits volumetric storage densities, regardless of the impressive gravimetric capacities reported recently for some novel adsorbents [7,8].
Nevertheless, significant advances have been made in the synthesis of new adsorbents with very high surface areas, and hence adsorption capacities, including metalorganic frameworks (MOFs), covalent organic frameworks (COFs), amorphous organic polymers, and novel types of porous carbon. Figure 1, for example, shows the trend of increasing gravimetric storage capacity with increasing BET surface area for MOFs, with the best performing zeolites and activated carbons included for comparison.
Nanoporous materials thus remain competitive, although overcoming some of the problems identified above will involve significant materials research and engineering challenges. However, it is worth noting that porous adsorbents such as zeolites, silica gels, activated aluminas, activated carbons, and carbon molecular sieves have been used industrially in vast quantities for decades, in many applications, including gas, vapour and liquid phase separations, and catalysis, so their economic and practical viability is proven. The main challenge in adsorptive hydrogen storage research has thus been the development of new adsorbents that share this viability while also showing sufficiently high hydrogen storage capacities for practical use.
In this article, the current status of the field is discussed and existing challenges are identified. Initially, the measurement of H 2 adsorption is addressed, including the question of the adsorbed phase density and volume, and the use of neutron scattering to study adsorbed H 2 . Suitable adsorption isotherm equations and the accurate computational modelling and simulation of H 2 adsorption by nanoporous materials are then addressed. Some of the interesting new materials that are being developed are then identified, and the development of practical storage tanks is examined.

Measurement
The measurement of H 2 adsorption is essential in order to characterise the hydrogen storage potential of nanoporous materials, while gaining a better understanding of the behaviour of the adsorbed H 2 in different materials is critical to the development of improved adsorbents. In this section, the measurement techniques, the properties of the adsorbed H 2 phase, and the study of adsorbed H 2 using neutrons are discussed.

H 2 adsorption measurement techniques
The amount of H 2 adsorbed is typically measured using either the manometric or gravimetric techniques. In the manometric case, known amounts of H 2 are prepared and introduced to the sample cell step-by-step, using measurements of pressure, temperature, and volume, and an equation of state for H 2 , to calculate the amount of adsorption (see Fig. 2a). At each point, the amount of H 2 in the gas phase is calculated and any missing H 2 is attributed to the adsorbed phase. In order to determine the number of Fig. 1 A plot of excess gravimetric hydrogen adsorption capacity at 77 K (at 2 MPa or above) versus BET specific surface area for a range of metal-organic frameworks measured in different laboratories, hence the inclusion of two values for Cu-BTC, MOF-5, and MOF-177 (modified and updated from [7]) moles in the gas phase, however, the accessible volume must be known. This void or dead volume is the difference between the volume of the empty sample cell and the volume occupied by the sample and the adsorbate. Determination of the combined sample and adsorbate volume is intrinsically problematic [9,10] and can be further complicated by sample swelling.
The gravimetric technique, on the other hand, directly determines the amount adsorbed from the weight measured by a microbalance (see Fig. 2b). In this case, the amount of gas introduced does not need to be calculated, and the error accumulation associated with the manometric method is also avoided [11]. However, the microbalance reading must be corrected for the buoyancy effect of the sample in the surrounding fluid (gas) [12], for which knowledge of the volume of the sample and the adsorbate is required, in a directly analogous way to the manometric method.
Both techniques have advantages and disadvantages [2], but both rely on the determination of the volume of the sample and the adsorbed phase in order to calculate the absolute adsorption or capacity. For materials amenable to helium pycnometry, this reduces to the problem of knowing the volume of the adsorbate as a function of temperature, pressure, and uptake.
If the adsorbate volume, V ad , is ignored in the manometric calculations mentioned above, i.e. it is assumed that V ad = 0, and the two techniques use only the sample volume in the calculation of the adsorbed quantity, then the uptake is underestimated by the amount of gas that would occupy V ad . The calculated quantity is then known as the excess adsorption or capacity.
Although both techniques are conceptually simple, the results obtained on identical samples in different laboratories by Zlotea et al. [13] give cause for alarm about their accuracy and reproducibility. Several laboratories were sent the same carbon molecular sieve with sample preparation instructions [13]. It can be seen in Fig. 3 that the reported results showed considerable disparity, particularly at higher pressure, so this remains a significant challenge for the characterisation of H 2 adsorption by nanoporous materials. It should be noted, however, that the results from a more recent exercise demonstrated better agreement [14].

Capacity definitions
The capacity is a measure of how much H 2 is stored by a material. It can be normalised gravimetrically (the amount stored on a mass basis) or volumetrically (the amount stored in a given volume). However, in the literature, hydrogen storage capacities are often reported without proper definition or sufficient experimental details. For example, the terms absolute and total, with regard to adsorbed quantities and capacities, are sometimes used interchangeably [2,15,16]. This variation in definitions can cause confusion and lead to unrealistic capacity claims.
The common types of capacity are the excess, n ex , i.e. the amount of gas present over and beyond the amount of a non-adsorbing non-interacting gas under the same Fig. 2 Schematic diagrams of a a simplified manometric apparatus [9] and b gravimetric apparatus employing a symmetric microbalance with both the sample (S) and the tare weight (T) suspended in the gas [11]. conditions of pressure and temperature, the absolute, n abs , the quantity of gas molecules in the adsorption volume, V ad (defined as the volume where the density of the adsorbate differs from the free gas density through gas-solid interactions), and the total, n tot , the quantity of gas molecules associated with the sample in both the free and adsorbed state. The gas molecules that are associated with the sample are determined by defining an appropriate volume, such as the bulk volume of the material.
An alternative to the standard definitions above is the net adsorption (or net capacity). This is defined as the amount of H 2 in the volume containing the adsorbent minus the amount of H 2 in the same volume, without the adsorbent, at an identical temperature and pressure [17]. This simply uses a different reference point for the calculation. It has the advantage that the volume of the adsorbent does not need to be known. It also provides a direct measure of the benefit of filling the tank with adsorbent, in comparison with pure gas compression, with the peak in the net adsorbed quantity indicating the optimum storage pressure, in this respect [18]. However, the use of the net adsorption is yet to be widely adopted.
In hydrogen storage studies, n ex is often reported as the amount adsorbed, but in other cases, n abs or n tot is reported instead. It is important, however, to clearly state the quantity being reported and explicitly define the assumptions used to calculate the capacities. The volume, V ad , of a gas adsorbed on an adsorbent has been a topic of considerable discussion because of its direct relevance to the determination of n abs . In addition, isotherm models, which we cover in Sect. 2.3, conventionally describe n abs and thus require the determination of V ad to compare theory with experiment. The adsorbed phase can range from a simple homogenous monolayer on a well-defined surface, to a complex mix of adsorbate-adsorbent interactions on a heterogeneous material, with pores of varying size and geometry, in the presence of different surface functionalities. Therefore, in complex systems, the resultant adsorbate-adsorbent interactions can vary not only as a function of pressure and temperature, but also by location. In the context of investigating new adsorbents, this translates to an ill-defined V ad that makes n abs difficult to quantify. The volume of the sample itself may also be difficult to accurately quantify throughout the sorption experiment. Although the calculation of n tot can involve a well-defined volume, the densification of the sample can drastically influence this volume and the hydrogen storage properties. Therefore, these methods must be explicitly reported as the densification could alter the physical structure of pore volume and interparticle void space.
It is possible to calculate n abs from n ex with an additional assumption. The two main options are to assume that either the volume occupied by the adsorbed H 2 , V ad , or its average density, q av , is known [2,8,19]. For a purely microporous adsorbent, for example, V ad can be assumed to be equal to the pore volume, V pore [18]. Therefore, n ex can be converted to n abs using the following expression [8,18], where q fg is the free gas density. Note that at low pressures the difference between n ex and n abs tends to be minimal, particularly when there is a significant amount of adsorption, as the V pore q fg term can be insignificant compared to n ex . Alternatively, the density of the adsorbed H 2 , q ad , can be assumed to be represented by an average value, q av [2,8,19]. In this case, n abs is given by, Unfortunately, neither V ad nor q ad is known and they are difficult to measure. Furthermore, the properties of the adsorbed H 2 are dependent on a number of factors, particularly in real, heterogeneous materials [20], and those with pores of varying size and geometry.
An additional concern that can lead to measurement error is the use of helium to determine the skeletal volume of the sample needed for the determination of n ex . This assumes that at an appropriate temperature and pressure helium does not adsorb and that it can access the same volume as H 2 . The problem with this approach is that neither assumption is likely to be entirely satisfied in nanoporous materials [8,21]. This remains a significant challenge for the accurate measurement of H 2 adsorption. To conclude, only the standard capacity definitions defined above [22][23][24] should be used, while investigation of the use of net adsorption more widely would be an interesting topic for future research. In general, it may be advantageous to report several of the capacity types to make the comparison between materials easier. Furthermore, it can be argued that while absolute adsorption is required for the direct fitting of isotherm equations and for the determination of the enthalpy or heat of adsorption [25], it is best avoided for the sole metric for comparison of the hydrogen storage performance of different materials. Finally, for a thorough comparison of different measurements, it is recommended to report all assumptions involved in calculating the capacity values (excess, absolute, total, and net) [24].

Neutron scattering
Microscopic information regarding the behaviour of H 2 at an atomic or molecular level can be obtained using neutron scattering techniques [26]; neutron diffraction, in particular, has widely been used to study D 2 adsorbed in porous materials [27][28][29][30]. Difficulties arise from the large incoherent scattering cross section of H 2 that makes the static structure factor, S(Q), difficult to determine from the resultant background. Standard diffraction techniques are also most helpful if the H 2 or D 2 is adsorbed in a structurally periodic manner. Other methodologies are better suited to the study of disordered fluid-like phases, including liquid diffraction or extrapolation from the kinematically accessible range of the fully dynamic structure factor, S(Q,x) [28]. Indeed, inelastic neutron scattering (INS) has widely been applied to probe the potential energy surface experienced by H 2 and to gain insight into H 2 weakly adsorbed in various carbons [29][30][31] and MOFs [32,33], and in the more strongly bound electrostatic environments found in zeolites [34,35] and Kubas-type inorganic compounds [36]. Furthermore, quasi-elastic neutron scattering (QENS) can be used to determine the dynamics of H 2 molecules adsorbed in nanoporous materials by analysing the broadening, in energy, of the elastic scattering peak [37][38][39]. Both INS and QENS, in contrast to diffraction, can exploit the large incoherent scattering cross section of H 2 .
Neutron diffraction is still typically performed using D 2 to reduce the incoherent scattering from H 2 . This isotopic substitution can affect the physics due to differences in zero point energy, as evidenced by differences in the thermodynamic properties [40][41][42]. Few studies detail the structural differences between H 2 and D 2 adsorption; however, with the advent of high-intensity neutron sources and diffraction instruments, this may change. The negative consequences of the incoherent background may be overestimated since most published neutron diffraction work on hydrogen-containing MOFs does not rely on deuterated MOF ligands, even though D 2 is subsequently used to determine their binding characteristics [43][44][45][46][47]. Several hydrogen-containing structures are also well known from powder diffraction [48]. Early evidence suggests there is little distinction between the structural parameters derived from the weak coherent scattering of H 2 compared to the strongly coherent scattering of D 2 adsorbed in one particular MOF, in which the centre of the H 2 nuclear scattering density was found to be 2.26(4) Å from the metal atom compared to a D 2 distance of 2.23(5) Å ; however, the broader impact of this is yet to be thoroughly evaluated [49].
Most measurements are also currently made at liquid helium temperatures, where the thermal motion and diffusion of weakly bound molecules is minimal. Barriers remain to performing detailed measurements on adsorptive systems over desired in situ operational temperature and pressure ranges, although the practical aspects of performing experiments under these conditions are not difficult to overcome at a neutron source [39,50,51]. Another consideration is the trade-off in neutron intensity versus data resolution. Higher accuracy structural information requires a higher-resolution diffractometer, which typically provides a lower count rate. This is mitigated somewhat at spallation neutron sources compared to reactor sources, and the situation will likely improve further with global development of higher-intensity neutron sources. Some fast diffractometers can currently perform this type of work, even if the results provide limited structural information, as evidenced by the report of negative thermal expansion and site occupation factors as a function of temperature for D 2 in a specific MOF that indicates D 2 retention even above 100 K [47]. The low temperature requirements of the scattering technique might also be mitigated by materials that adsorb H 2 at higher temperatures, as required for H 2 storage. In this case, the enthalpy of H 2 adsorption would presumably be large enough for molecules to remain bound at higher temperatures, thus permitting structural characterisation over a wider temperature range.
Neutron total scattering or pair distribution measurements are an alternative for studying more fluid-like adsorbate phases. This technique additionally incorporates the short-range order present in the background of diffraction data [52]. Methodologies to deal with incoherent scattering backgrounds are being developed [53] and, while unlikely to allow extensive study of H 2 adsorption on structures that are poorly defined crystallographically, may be useful to study MOFs with H-containing ligands.
Many assumptions used to analyse macroscopic gas adsorption data rely on the knowledge of the adsorbed phase volume or density, as discussed in Sect. 2.2, for which independent measurements are difficult or impossible. For complex pore networks, attempts have been made to characterise and reconcile pore structures obtained from small-angle scattering data with those from adsorption/intrusion studies (e.g. [54]). However, the data are complicated by the multiple length-scales, closed porosity, and surface roughness factors that inhibit a unique description and depend upon the adsorption models used (see [55,56], for example). With the advent of MOFs that are crystallographically well defined, and the observation that the adsorbed D 2 /H 2 phases can be characterised in a precise crystallographic manner, it is likely that a robust link between assumptions in traditional gas adsorption techniques and independent observation can be achieved, particularly in the case of the heterogeneous adsorption potentials exhibited by unsaturated metals in MOFs. Using high-resolution diffraction, and quantification of the porefilling for other adsorbates, such as N 2 (typically used for surface area determination) and D 2 in MOFs, it should be possible to obtain a real space view of the gas adsorption process at 77 K and thus extract meaningful parameters corresponding to estimated molecular size and adsorbed phase density. Judicious choice of the pore structure to be studied in a series of MOFs would allow the effects of confinement and dimensionality of the pore and adsorbateadsorbate interactions on curved surfaces to be interrogated in detail.

Modelling
In order to gain further insight into the H 2 adsorption process, to provide data for the design of full storage systems, and to screen materials for practical applications, various modelling and simulation techniques can be used. In this section, the fitting of macroscopic H 2 adsorption data to analytical isotherm equations and the use of molecular-level computational methods are discussed.

H 2 adsorption isotherm modelling
Experimental H 2 adsorption isotherm data can be described by theoretical models through the adjustment of their parameters using, for example, least squares minimisation [57,58]. This fitting process allows the prediction of thermodynamic properties, such as the excess adsorption and isosteric enthalpies, over a wide range of temperatures and pressures [59] and can provide an estimate of parameters such as the adsorbed phase volume, saturation pressure, and limiting adsorption capacity, which are otherwise difficult to determine [59,60]. Isotherm fitting can also provide insight into underlying adsorption mechanisms and the pore structure characteristics of adsorbents [60,61], while data needed for the modelling and design of hydrogen storage systems can also be obtained [62][63][64][65][66][67].
Various models are available [20,68], from simple expressions, such as the Langmuir and Freundlich equations, to those of greater complexity, which include the Tóth [69], Unilan [61,70], modified Dubinin-Astakhov (DA) [59] models, and the multicomponent potential theory of adsorption (MPTA) [60,71]. Different assumptions are used; for example, the Langmuir equation describes adsorption on energetically homogeneous surfaces. Each site is occupied by only one molecule, with no adsorbateadsorbate interactions. The Unilan model uses the Langmuir equation, but accounts for heterogeneity using a uniform distribution of adsorption site energies, while the modified DA model is based on Dubinin's theory of volume filling of micropores [72]. It thus assumes that adsorption occurs due to a Polanyi-type adsorption potential that is present throughout the pore volume, so adsorption occurs via pore volume filling rather than the monolayer formation assumed in the Langmuir model [68].
The simpler models can fit data for a wide range of traditional adsorbents over limited ranges of pressure and allow the transformation of measured isotherms into linear forms from which model parameters can easily be obtained using linear regression [68]. However, due to deviation of the predictions by the simpler models for data measured over wide pressure ranges [68] and the lack of temperature dependence, as required in realistic adsorptive storage systems, more rigorous adsorption models, such as the modified DA, Unilan, and MPTA approaches, tend to be favoured by hydrogen storage system developers.
Dundar et al. [73] recently found that different models were needed to fit data for different MOFs. The modified DA, Unilan, and MPTA models were used to fit experimental H 2 adsorption isotherms for three prototypical materials (MOF-5, Cu-BTC, and MIL-101) at different temperatures. Isotherms for Cu-BTC and MIL-101 were better described using the modified DA model, while H 2 adsorption on MOF-5 was better described by the Unilan and Tóth models. The fits for MOF-5 data, measured and reported by Zhou et al. [74], in the temperature range 77-300 K, are shown in Fig. 4 [73]. The use of the modified DA model to fit H 2 adsorption isotherms for MOF-5 was also examined by Purewal et al. [61]. It produced a well-behaved fit, but the predicted isotherms deviated from experimental data between 200 and 300 K, resulting in anomalous negative adsorption. The results also predicted unusually large saturation pressures, free energies of adsorption, and heterogeneity parameters. However, the model better described H 2 adsorption data for Maxsorb and Cu-BTC [61].
Although both the Unilan and Tóth models outperform the modified DA model when used to fit data for MOF-5, the MPTA has been found to fit H 2 adsorption isotherms for all prototypical MOFs better than the other models. It also provided insights into the properties of the adsorbed phase [73]. For example, it indicated an asymptotic increase in H 2 density inside the pores, at lower temperatures (&40 K), analogous to a gaseous to solid-like phase transition [73]. However, system-level models using the MPTA are yet to be reported, probably due to its cumbersome iterative calculations and the lack of an analytical expression for isosteric enthalpies. Nevertheless, despite its known limitations, the modified DA model is largely favoured for H 2 storage system modelling [62-64, 66, 67, 75]. Not only are its analytical expressions for isosteric enthalpies and absolute adsorption simple to implement in system-level models, but it is also a non-iterative method and is therefore less cumbersome.
The models discussed above have so far been used to describe adsorption that is primarily due to van der Waals interactions. However, if functionalised nanoporous materials that interact more strongly with H 2 are developed, it is possible that adsorption may involve both physisorption and weak chemisorption. A future challenge in this area is thus the determination of the limits of the existing models and adapting them to fit H 2 adsorption isotherms on nextgeneration adsorbents. Nevertheless, testing different models on a wider range of datasets would be a valuable near-term goal.

Computational modelling and simulation of H 2 adsorption
Studying and selecting nanoporous materials for H 2 storage applications is challenging due to the rapid growth in the number of new structures and the time-consuming nature of experiments. Computational methods, however, can predict macroscopic adsorption properties and thus offer an efficient alternative. They can also help design better adsorbents by providing insight into structure-property relationships and can yield molecular-level information that is otherwise inaccessible. Adsorption-based H 2 storage studies most commonly involve grand canonical Monte Carlo (GCMC) simulation and quantum chemical (QC) calculations.

Grand canonical Monte Carlo (GCMC) simulations
GCMC simulations can be used to determine both adsorption isotherms and isosteric enthalpies of adsorption since they imitate experiment [76][77][78]. The results-the number of particles, N, in a model porous solid versus the external chemical potential, l, at a temperature, T-are directly comparable to the output of adsorption experiments, although simulations calculate the total amount of adsorbate in the pores, whereas experiments measure the excess (see Sect. 2.2). Apart from technicalities, including the size of the simulation box, choice of periodic boundary conditions, and the number of iterations, the quality of a GCMC simulation depends primarily on how the adsorbent, the adsorbate, and the interactions between them are described.
For crystalline materials, such as MOFs and COFs, for which the coordinates of all solid atoms are known from diffraction experiments, the adsorbent can be described in crystallographic detail. For amorphous materials, like porous carbons or low-crystallinity polymers, a generic surface model structure, such as graphite or graphene, and a pore model, such as a slit, cylinder, or sphere, are usually adopted. Nevertheless, various attempts have been made to more accurately model complex amorphous pore networks [79]. Methods include the use of virtual porous carbon [80] and finite wall thickness models [81], quenched molecular dynamics [82], simulated polymerisation algorithms [83], and the packing of 3D carbon nanostructures [84].
Molecular H 2 is typically modelled either as an uncharged sphere [85] or as a single mass centre, but threepoint-charged, dumbbell [86]. The latter accounts for the molecular quadrupole moment by considering two positive charges, q, at the ends of the dumbbell, and a negative charge, -2q, at the centre of mass. In most cases, a 12-6 Lennard-Jones (LJ) potential-plus electrostatic interactions for charged H 2 models-is used to represent the H 2 -H 2 interactions, with different sets of LJ parameters considered for either the spherical or dumbbell H 2 models [87]. The Morse potential has also been used as it can be adjusted to fit a particular force field (FF). More complex models are also available [88,89]. In addition, the quantum nature of H 2 becomes evident at low temperatures, in small pores and at high densities [39,87,90,91]. There are two main ways of accounting for quantum effects, although alternatives are being developed [91]. The first is to use the Feynman-Hibbs (FH) correction [92], which is an expanded version of the LJ potential [39]. In most cases, particularly above 50 K, this is sufficiently accurate. The second is the path integral MC formalism outlined by Wang and Johnson [90], which is more accurate but is also more computationally demanding.
The H 2 -solid interactions are often described using generic classical FFs, such as the universal force field (UFF) [93], DREIDING [94], and the optimised potential for liquid simulations (OPLS) [95], but this approach is limited by the accuracy of a given FF to describe the interactions of H 2 with all the framework atoms [78]. In some cases, classical approximations may fail, particularly when short-range interactions of H 2 with metal centres, for example, are dominant [96]. When Coulombic interactions are important, FFs should also be supplemented with atomic partial charges. These can be calculated using Mulliken population analysis [97], charge partitioning methods [98] or with the aid of cluster-based electrostatic potential fitting (e.g. CHelpG [99] and RESP [100]). New electrostatic potential-based methods for periodic solids include REPEAT [101,102] and DDEC [103,104], while charge equilibration methods can be used for fast, highthroughput calculations, at the expense of accuracy [105]. Determination of atomic charges can be avoided by using full DFT-based electrostatic potential maps of the material [106], but this is computationally expensive.
The most accurate strategy for developing case-specific H 2 -solid interactions is to perform QC calculations, as discussed next. These methods can be very accurate and are thus recommended for the study of novel systems, in which the adsorption mechanism is unknown, and to probe very specific interactions. However, such multiscale approaches are computationally intensive. Detailed information can be found in recent reviews focussing on MOFs [107][108][109].

Quantum chemical (QC) calculations
There are two main types of QC calculation: wave function theory (WFT) and electronic density functional theory (DFT). They have many similarities since they share a number of mathematical approximations. In both cases, the Schrödinger equation is solved in order determine the properties of a molecular system using different approximations, including the Born-Oppenheimer approximation, the representation of the wave function as a Slater determinant and the basis set approximation.
The Hartree-Fock (HF) theory underpins both WFT and DFT, but it neglects the dynamical correlation of the electrons in the system, which is fundamental to the calculation of several properties, including the energy of the system. Inclusion of the effect of the dynamic motion of electrons is essential when dealing with molecular systems in which weak interactions are dominant, as is the case for H 2 physisorption. In order to include dynamic electron motion, several WFT methods have been developed, including single-reference methods, such as Møller-Plesset (MP) perturbation or coupled cluster (CC) methods, and multireference methods; however, the latter are rarely used because they are computationally demanding. The most common is second-order MP perturbation theory (MP2) [110], in which electron correlation is treated as an additive perturbation to the one-electron operators. CC methods have also been used for H 2 adsorption in nanoporous materials, but they are limited to systems with only a few atoms due to their computational cost. Nevertheless, CC methods serve as a benchmark in order to check the accuracy of other single-reference or DFT methods; CCSD(T) at the basis set limit, in particular, is considered the gold standard in computational chemistry, although care should be taken because results obtained using this method can still contain errors [109].
Accurate calculations depend not only on the selection of the theoretical method but also on the choice of basis set. The most accurate result for a given method can be achieved using an infinite expansion for each molecular orbital. A small loss of accuracy is expected in the calculated energy values when using finite basis sets, which can be retrieved if a complete basis set (CBS) scheme is applied. Another error originating from the incompleteness of basis sets is the basis set superposition error (BSSE), which artificially increases the calculated energy of the system. If a method that can treat this error is applied, such as the counterpoise (CP) method [111], a 20-50 % reduction in the calculated interaction energy for H 2 adsorption on carbon nanotubes or MOFs has been found.
DFT uses the idea that the spatial distribution of electron density uniquely determines the ground-state wave function, and vice versa, while if the universal density functional (i.e. the functional that uses the electron density to calculate the energy) is available, the ground-state energy can be obtained by variation. Following the work of Kohn and Sham [112], the approximate functionals of DFT partition the energy of the system into several terms, most of which are functions of the electron density, except the nuclear repulsion term. The different flavours of DFT arise from the choice of the exchange-correlation term of the energy, which accounts for the remaining (non-classical) terms and can be further separated into two functionals of the electron density: the exchange functional (same-spin electron interaction) and the correlation functional (mixedspin electron interaction). The different forms of these functionals are broadly known as approximations of DFT.
The earliest were based only on local electron densities [local density approximation (LDA)], while gradient approximations [generalised gradient approximation (GGA)] depend additionally on the density gradient [109]. They have widely been used to study H 2 adsorption in carbon-based materials and MOFs, with the most common examples being PBE, PW91, BLYP, and BP86. Although these approximations include terms that treat electron correlation, they can fail to correctly describe the interaction of H 2 with solids and to predict accurate interaction energies. Beyond LDA and GGA, hybrid functionals that include a percentage of non-local exchange obtained from HF calculations have been used to study adsorption, with B3LYP being the best known. Others include the B3PW91 and PBE0 functionals. Other categories include meta-GGA and hyper-GGA functionals, but these are rarely used in H 2 storage studies. For a detailed description of exchange correlation functionals, see Odoh et al. [109].
A drawback of DFT is the rather poor description of weak interactions, which are important in adsorptive gas storage [78]. This is due to the failure of the functional to accurately describe medium-to long-range correlations. However, attempts have been made to develop DFT methods that are able to correctly account for dispersion interactions. These include the use of a non-local dispersion term in the correlation functional (vdw-DF), the use of a non-local term in the functional to account for mediumrange correlation, and the addition of a pairwise contribution to the electronic energy due to dispersion interactions (DFT-D). They have all been successfully applied to gas adsorption in MOFs and outperform the corresponding pure GGA or hybrid functionals without increasing computational effort.
DFT methods are faster than WFT and can treat larger systems. Moreover, DFT can be applied to a fragment of a material or to a periodic cell, whereas WFT can only be applied to a small fragment. However, WFT can be very accurate and the way to increase the accuracy is known. In contrast, DFT methods require an appropriate functional to be found. The main challenge in this field is the development of more efficient methods and algorithms that offer improved accuracy for a given computational cost. For the modelling of H 2 adsorption, in particular, the accurate but efficient description of weak interactions is a clear priority.

Discussion on computational modelling
Despite recent progress, GCMC simulation methods could be further improved because there are problems that have not yet been resolved or fully investigated. There are also inherent shortcomings in GCMC methodology that should always be considered. For instance, a technical problem surfaces at high adsorbate densities. When pores are almost full, or the adsorbate molecules are highly confined (e.g. in ultramicropores), the normal GCMC technique fails due to the low acceptance ratio for insertion. This renders the exploration of phase space extremely slow. Several biased-GCMC techniques have thus been developed to overcome this problem [113]. This effect is particularly relevant to large molecules, such as aromatics, adsorbed in small cavities, but it may also play a significant role in H 2 adsorption simulations when high densities (high pressures) and/or very narrow channels are considered [114].
Another issue is the ability of pairwise interactions to adequately describe the adsorbate-adsorbent system [78]. Many-body interactions can, in principle, be included by correctly parameterising the LJ parameters to fit the bulk phase behaviour of H 2 . It is then assumed that the behaviour of the pore confined phase is the same as the bulk, or that a particular set of LJ parameters can still describe the behaviour under confinement. However, the contributions of many-body interactions may be very different under confinement, so different case-specific LJ parameters should be used for each pore model. This is more than just a geometric consideration because in adsorption a foreign body is also included (the solid) and therefore the third body, in three-body interactions, for example, may be a framework ion. Further details on the importance of this can be found in Kostov et al. [115].
With regard to future challenges, the requirements differ depending on the material type. For example, a regular pore model is normally used to describe amorphous materials, which in reality consist of a very complex pore network. Pore connectivity is seldom considered, and simulations are based on individual pore models. These models may work well for relatively ordered mesoporous materials, but not for disordered microporous materials, especially those exhibiting ultramicroporosity, which are of most interest for H 2 storage.
Another issue is that an energetically homogeneous surface is normally used to construct pore models. However, real carbon surfaces, for instance, have a degree of inhomogeneity induced by surface heteroatoms and functional groups, but also crystal defects, local curvature, geometrical imperfections, and a lack of long-range periodicity. Accounting for surface inhomogeneity and its effects on adsorption has been one of the most challenging and important topics in adsorption science for over 30 years [20]. Several mathematical approaches exist to determine the amount of inhomogeneity, but no satisfactory way of efficiently including this in molecular simulations has emerged. Energetic inhomogeneity is coupled with the inherent inhomogeneity created by the presence of pores of differing size, a common characteristic of non-or poorly crystalline materials. Specific inhomogeneities such as functional groups or lattice defects can be added to carbon surfaces, and their effect on H 2 adsorption studied; currently, however, this can only be done using an explicit atomistic description of the surface [116]. Notable approaches towards the development of a mean-field description of the heterogeneity of adsorbent surfaces are QSDFT (quenched solid density functional theory) [117] and 2D-NLDFT (non-local density functional theory) incorporating energetic heterogeneity and geometric corrugation [118,119]. The future application of such meanfield approaches to the GCMC simulation of H 2 adsorption by carbons can be expected to significantly improve upon the currently used methods of representing heterogeneity in these and related materials.
For crystalline adsorbents, the pore size and shape are known, so irregularity is less relevant; however, simulations are typically performed on idealised, fully desolvated systems, which can result in disagreement with experiment [108]. For example, the boundary conditions used in GCMC simulations assume crystals of infinite size but real materials consist of small crystals, with defects, intercrystalline voids, and fused crystallites, while in many MOFs, the complete removal of solvents is practically impossible. Moreover, the framework atoms are frequently assumed to be fixed, even though the structure is known not to be rigid; in fact, several pertinent phenomena, including structural framework transitions, and breathing or gating mechanisms in MOFs, for example, have been widely observed [120]. In such cases, optimisation of lattice constants and atomic positions, and consideration of the effects of framework flexibility on adsorption is required. The latter can be achieved using molecular dynamics (MD) [39], osmotic thermodynamic ensemble MC simulations or, more efficiently, hybrid MC/MD (hybrid Monte Carlo, or HMC) [121][122][123], in which short MD trajectories are considered as MC moves, allowing better sampling of the host framework flexibility by following its collective motions. More details can be found in recent reviews [120,124].
Consideration of the above problems and the study of their effect on adsorption remains a significant challenge. Furthermore, the development of a more accurate and less case-specific modelling approach would be a convenient tool for MOFs, due to the number of synthesised structures. The development of a ''MOF force field'' built on the basis of accurate quantum-level calculations and validated using available experimental data for families of MOFs sharing similar chemical characteristics would be an invaluable future asset.
The use of high-throughput computational screening is also a notable development. A series of databases containing a large number of porous materials have recently been reported [125][126][127][128]. This offers the potential for screening materials for H 2 storage. The H 2 storage capacity can also be correlated to the characteristics of different materials in order to focus on the most promising and to guide the synthesis and tailoring of new MOFs. Colón et al. [129], for example, considered frameworks functionalised with Mg and correlated the uptake to the degree of Mg functionalisation and the physical properties of the frameworks. Goldsmith et al. [130] also used high-throughput methods to determine promising materials for H 2 storage. In contrast to other reports, however, the screening in this study included real rather than hypothetical frameworks, using data mining techniques on a vast catalogue of existing MOFs in the Cambridge Structural Database (CSD). Results for approximately 20,000 structures showed that the relationship between gravimetric and volumetric H 2 storage density is concave downward (see Fig. 5). The use of these high-throughput methods in the future can be expected to contribute significantly to the search for new nanoporous materials for H 2 storage.

Materials
Many different porous materials have been investigated for adsorptive H 2 storage, including zeolites [133], various types of porous carbon [134][135][136], MOFs [16,[137][138][139], and porous organic polymers [140,141]. However, as already noted, the interaction of H 2 with most surfaces is weak and high H 2 uptakes are limited mostly to low temperature conditions, around 77 K, and to materials containing micropores, in which the adsorption potentials from the opposing pore walls overlap. As shown in Fig. 1, an approximately linear relationship between the gravimetric capacity and BET surface area has been observed [7,142], although the saturation capacity of an adsorbent is ultimately limited by its pore volume [137].
The highest gravimetric capacities have thus been reported for materials with very high surface areas and large pore volumes. MOF-177, for example, has an excess H 2 adsorption capacity of *7.5 wt% at 7.0 MPa and 77 K [143], while an even higher value of 8.6 wt% has been reported for MOF-210, under essentially the same conditions [144]. MOF-177 and MOF-210 have reported BET surface areas of approximately 4700 and 6250 m 2 g -1 , respectively. Increasing surface area therefore offers a route to increasing gravimetric capacity, but this approach has limitations. Firstly, it leads primarily to an increase in gravimetric capacity only at low temperatures, and secondly, there is a physical limit to the achievement of even higher surface areas. Furthermore, materials with the highest reported surface areas also tend to have larger pores. This, in turn, leads to a reduction in the adsorption potential overlap responsible for the stronger H 2 -solid interactions in materials with very narrow pores. Larger pores are likely to lead to lower total volumetric capacities (see Sect. 2.2) because in larger pore materials, particularly those approaching the mesoporous regime ([2 nm), H 2 is more likely to form a gas-like phase in the core of the pores.
It can be seen that achieving high storage capacities at near-ambient temperatures thus requires an approach other than simply increasing surface area. Furthermore, both the volumetric and gravimetric capacities of materials need to be considered. High gravimetric capacities result in lighter tanks, but a poor volumetric capacity increases bulk. These aspects are correlated, and it is important not to focus on one to the detriment of the other. To increase either the volumetric or gravimetric capacity of a material at nearambient temperatures, it is necessary to increase the enthalpy of adsorption by increasing the strength of the interaction between H 2 and the material. This could potentially be achieved by altering the chemical nature of the adsorbent, but it is a challenge due to the properties of molecular hydrogen. Other options include the exploitation of framework flexibility, which can lead to hysteretic H 2 adsorption behaviour, and modification of the H 2 adsorption behaviour of materials using novel concepts, such as core-shell architectures [145]. Each of these challenges will be considered below.

Volumetric versus gravimetric capacity
The correlation between volumetric and gravimetric capacity was illustrated in the study by Goldsmith et al. [130]. As noted in the previous section, a concave relationship between the volumetric and gravimetric capacities of MOFs was found. The materials with the highest gravimetric capacities thus exhibit lower volumetric  [130]. The dashed lines indicate the current 2020 US Department of Energy (DOE) volumetric and gravimetric hydrogen storage system targets. Note that the results plotted in blue are for materials only. The values for complete storage systems incorporating these adsorbents will be considerably lower due to the factors discussed in more detail in Sect. 5; for example, the figures for a recent prototype adsorption system, labelled ''MATI/MOF-5'' (see [131] for more detail), and for a type 4 compressed H 2 storage system [132] are also shown (in red).
(Reproduced and modified with permission from [130] (http://pubs. acs.org/doi/pdf/10.1021/cm401978e). The red symbols and labelling have been added to the original figure) capacities than some materials with higher volumetric capacities. The results of this study are shown in Fig. 5. These two factors clearly need to be balanced.
One strategy to help overcome limitations in volumetric capacity is pelletisation, which can increase the volumetric uptake of a given material in a practical engineered form (see Sect. 5.1). Also, for MOFs, there is potential for framework interpenetration to increase the surface area of a given material per unit volume. The presence of additional skeletal material from the interpenetrated frameworks in a given volume might be expected to result in a concomitant decrease in the gravimetric capacity; however, studies have suggested that interpenetration can increase both the volumetric and gravimetric capacities of some materials [146].
It is also worth noting the importance of considering the usable or deliverable capacity of a material [61,[147][148][149][150][151][152], rather than just its total or absolute capacity (as discussed in Sect. 2.2). This can be defined as the amount of H 2 stored reversibly between the maximum storage pressure and the delivery pressure of a practical storage unit. The latter is determined by the required back pressure, for example the pressure practically required by a fuel cell stack. For a given set of operating conditions, the usable capacity depends upon the H 2 adsorption properties of the adsorbent. It also exhibits a peak when plotted against temperature, so that there is an optimum operating temperature for each material, which tends to be higher for materials with higher enthalpies of adsorption. The usable capacities of different nanoporous materials are discussed in more detail by Schlichtenmayer and Hirscher [153].

Increasing the H 2 interaction strength
In 2006, Bhatia and Myers [147] considered the optimal enthalpy of adsorption for achieving ambient temperature H 2 storage in nanoporous materials. Using a simple model, assuming Langmuir-type behaviour (see Sect. 3.1), they concluded that the optimum enthalpy of adsorption for the delivery of H 2 from a hydrogen store between the pressures of 3.0 and 0.15 MPa was &15 kJ mol -1 H 2 . Other studies reached similar conclusions, although a higher value is sometimes quoted [15,78]. For example, Bae and Snurr [148] suggested a value of approximately 20 kJ mol -1 H 2 based on GCMC simulations (see Sect. 3.2.1) of H 2 adsorption by a series of MOFs over a wider pressure range, up to 12 MPa. To put this in context, H 2 adsorption by carbons has an isosteric enthalpy of adsorption of approximately 6 kJ mol -1 [147].
Amongst the newer materials, MOFs have the potential to achieve a higher enthalpy of adsorption due to the presence of open metal sites and increased local charge densities in the pores of some materials. For example, the M-MOF-74 (M: Mg 2? , Co 2? , Ni 2? , and Zn 2? ) series, which is also known as CPO-27, has a high density of open metal sites. The trend, in this case, for the H 2 -metal interaction strength was found by Pham et al. [154] [155]. This behaviour is surprising because one would expect higher H 2 binding in the case of Mg-MOF-74 since Mg 2? is small and is therefore the hardest cation (high partial charge) in the above series. The observed trend can be explained on the basis of the different polarisabilities of the metal cations. Theoretical calculations showed that the higher the contribution from polarisation, the stronger the H 2 -metal interaction.
With regard to local charge densities, the study of MOFs with either positively or negatively charged frameworks would be an interesting topic for future work. In the rare case of a cationic rht-MOF with NO 3 anions acting as counterions, it has been demonstrated that H 2 adsorption occurs first near the NO 3 ions because of their preferential interactions with H 2 [156]. This finding is interesting considering that this rht-MOF also contains open metal sites (dimeric copper paddlewheel units) that usually show the strongest interaction with H 2 . Anionic MOFs are more abundant; they are usually obtained when metal cations with ?3 oxidation states such as In 3? , Ga 3? , and lanthanides are combined with carboxylate-based organic linkers. In this case, the extra framework charge-balancing cations provide strong H 2 adsorption sites. In addition, these counterions are, in principle, exchangeable, which provides a potential route to H 2 adsorption sites with tunable energetics. An example is the family of anionic zeolite-like MOFs (ZMOFs), which can be constructed from a variety of organic linkers, including 4,5-imidazoledicarboxylic acid (H 3 ImDC) or 4,6-pyrimidinedicarboxylic acid [157]. In these materials, the enhanced binding of H 2 (9 kJ mol -1 H 2 at zero surface coverage) can be attributed to the electrostatic field created by the counterions present in the pore cavities. Notably, there is an observed increase of almost 50 % in the isosteric enthalpy of adsorption compared to corresponding neutral MOFs.
Another approach is the modification or design of the organic linkers in MOFs in order to increase the strength of their interaction with H 2 . There are two main mechanisms: the introduction of additional adsorption sites on the functional groups and the secondary effect of the functionalities on the polarity of the framework [158]. The latter can increase the H 2 affinity of the secondary building unit. An example is provided by a study of H 2 adsorption in Zn-based MOFs containing internally polarised organic units [159]. The use of 2,6-azulenedicarboxylate in MOF-650, instead of the nonpolar 2,6-naphthalenedicarboxylate linker used in IRMOF-8, resulted in a high initial isosteric enthalpy of adsorption, of 6.8 kJ mol -1 H 2 , compared to its nonpolar counterpart. Another example is the introduction of amide functional groups, which has been shown to increase the interaction strength of H 2 [160]. Mixed linker strategies can also be used [161]. From a synthetic point of view, there are various ways of modifying MOFs, including linker exchange, chemical functionalisation, and postsynthetic cation exchange. Cohen [162] covered the postsynthetic chemical modification of MOFs in some detail, while Deria et al. [163] reviewed the use of linker, nonbridging ligand and metal ion exchange to modify these materials. MOF surfaces can also be modified in various ways [164]. Further investigation of the use of these approaches to improve the H 2 storage properties of MOFs would seem both likely and worthwhile.
A relatively new class of materials that are yet to receive a significant amount of attention for H 2 storage are nanoporous molecular crystals. These consist of discrete molecules interconnected via non-covalent interactions, rather than covalent or coordination bonds [165][166][167][168]. They are formed by removing guest molecules from an inclusion compound. This normally results in a non-porous solid, but permanently porous materials with high BET surface areas, up to *3750 m 2 g -1 , have been reported [169]. There are several types of nanoporous molecular crystal, including those with intrinsic and extrinsic porosity [170], but intrinsically porous organic cage compounds [166] seem particularly interesting for H 2 storage applications. It is possible, for example, that such materials can be functionalised with multiple open metal sites inside one cage and that these opposing strong sites may increase the H 2 interaction strength further.

Flexible frameworks and core-shell materials
Hysteretic H 2 adsorption has been reported for a number of MOFs [171][172][173][174][175], although it is relatively rare. For H 2 storage, the interest in this behaviour lies in the potential for nanoporous materials to release H 2 at a different pressure to that used for H 2 charging or to induce H 2 release via temperature changes, from materials in which the H 2 would otherwise remain trapped, either thermodynamically or kinetically. This is in contrast to the behaviour exhibited by most-more rigid-materials that show completely reversible H 2 adsorption that can be accurately described by the isotherm models covered in Sect. 3.1. In addition, if flexibility is present, the expansion and contraction of the framework can be accompanied by a change in heat that may act positively on the thermodynamics of gas uptake and release. This was demonstrated recently for two flexible MOFs, Co(bdp) and Fe(bdp), that exhibit hysteretic methane adsorption [176]. Further development of MOFs showing hysteretic H 2 adsorption behaviour could potentially allow the tuning of materials with properties specifically suited to the requirements of H 2 storage; however, only preliminary measurements currently exist. Further experiments on well-defined model systems will be necessary to better understand this phenomenon.
Core-shell materials, meanwhile, consist of a core of one material encapsulated in a shell of another [145,161,164,177,178]. This offers the possibility of combining the behaviour of the core and the shell, to provide functionality arising from the different properties of the two components. For example, the shell could act as a gate for uptake and release of the gas stored in the core; the idea being that the gate opening of the shell would occur at a well-defined temperature, thus releasing the H 2 stored in the core. This may result in an isotherm of a form closer to that of a classical hydride than the reversible type I behaviour exhibited by most nanoporous materials. H 2 adsorption and desorption measurements have not yet been performed on such core-shell materials, so this would be an interesting avenue for future research. Furthermore, although a number of core-shell MOFs have been reported [145,161,164,177], other nanoporous materials have also been used [178]. It would thus seem possible to form core-shell architectures from many different combinations of materials, which could in turn lead to composites with properties of interest for H 2 storage.
It should be emphasised that, at the time of writing, both flexible frameworks and core-shell materials are a long way from any application, since very few measurements using H 2 are known. Measurements on model systems must be established in order to better understand the microscopic mechanisms, prior to tailoring materials with properties designed specifically for H 2 storage.

Storage systems
The engineering objectives and operating requirements for storage systems are defined by the intended application. This discussion considers automotive use, which has been the focus of much recent attention [179]. Technical targets were established by the U.S. Department of Energy (DOE) together with the automotive industry [180]. The most challenging of these for adsorption-based systems are gravimetric capacity, volumetric capacity, and charging time. The first is limited by adsorbent characteristics. Volumetric capacity is also material specific, but approaches are available to reduce the adsorbent volume, and hence that of the system. Thirdly, since adsorption occurs very rapidly, charging time is not limited by the adsorption kinetics but by heat and mass transfer.
Gas adsorption is exothermic so heat is typically released during charging. The total amount released depends upon the enthalpy of H 2 adsorption. For MOF-5, for example, this is 3.5-4.5 kJ mol -1 H 2 , depending on the loading [181]. In a system storing 5.6 kg of usable H 2 , 10 MJ of heat must therefore be dissipated due to adsorption alone. In addition, the system is refuelled from a heated (empty) state and then cooled to a charged (full) state. For economic reasons, liquid nitrogen (LN 2 ) can be used for cooling, so the tank will be at around 77 K when full. H 2 adsorption isotherms show that over 95 % of the H 2 is released by 160 K. Accounting for the mass of the adsorbent and an Al type 1 vessel, with a working pressure of 6.0 MPa, the stored thermal energy would be another 17.5 MJ. This must also be dissipated to reach the full state. Operating pressure is critical to performance and cost. Vessels operating below 6.0 MPa can be all-metal (type 1) tanks, which are relatively inexpensive. Above this, however, tank weight increases significantly with increasing pressure. Lighter composite wrapped over metal liner (type 3) or composite wrapped over polymer liner (type 4) tanks can be used but at significantly increased cost.
Dissipation of 27.5 MJ of heat is a serious challenge, particularly when tanks must be charged rapidly. For a refuelling time of 3.3 min, for example, the required heat dissipation rate is 90 kW. This problem is compounded by the low thermal conductivities of adsorbents. Charge cycle requirements dominate heat exchanger design because charging is performed in minutes while discharging occurs over hours. Hence, this will be the focus of the discussion in this section.

Material configuration
Adsorbents are usually synthesised in a finely granulated or powder form. Powder is ideal for increasing gas accessibility, but the interparticle void space is inefficient for gas storage, with its capacity limited to compressed gas densities. Some adsorbents pack densely to approximately 60 % of their crystal density, the highest achievable value for randomly packed spheres. Other adsorbents are less dense due to their electrostatic nature and pack to only 20 % of their crystal density [182]. This void space, while not volumetrically efficient, acts as a buffer tank that can supply H 2 through a pressure drop without thermal intervention. Fine powder requires closely spaced heat flow augmentation, as in typical tube-fin heat exchangers, in which Al fins are arrayed in the powder adsorbent and kept in contact with a cooling/heating source that dictates the adsorbent temperature.
Alternatively, adsorbents can be densified by mechanical consolidation. MOF-5, for example, can bind to itself under relatively low compaction pressures [182]. However, caution is required when using mechanical methods because pore collapse can occur, leading to surface area loss at high compaction pressures. This is a particular problem for MOFs because the ligand structures can collapse irreversibly [183,184]. Binders can be used in some cases, but this decreases gravimetric capacity.

Thermophysical properties
Models indicate that the thermal conductivity of an adsorbent bed for automotive storage systems with a refuelling time of 3.3 min must be in the range 1-3 W m -1 K -1 in order to enable removal of the heat of adsorption. Ahluwalia et al. [185] suggested a minimum thermal conductivity target of 1 W m -1 K -1 , which could be achieved by adding up to 20 wt% expanded natural graphite (ENG) in order to minimise the combined weight of the heat exchanger and the thermal conductivity enhancing media. However, Chakraborty and Kumar [186] designed a system that could use an adsorbent with a thermal conductivity of only 0.3-0.5 W m -1 K -1 by integrating a heating coil into the bed and relying on flowthrough cooling, as discussed below.
One of the first reports on the low thermal conductivity of MOFs was by Huang et al. [187]. For single-crystal MOF-5, it was found to be 0.32 W m -1 K -1 at ambient temperature, and 0.22 W m -1 K -1 near 100 K. MOF-5 powder was shown by Ming et al. [182] to exhibit a lower thermal conductivity, depending on the density of the compacted powder. At 300 K, values from 0.091 W m -1 K -1 at 0.35 g cm -3 to 0.16 W m -1 K -1 at 0.69 g cm -3 were found. Schlemminger et al. [188] stressed the importance of knowing the temperature-dependent thermophysical material properties. The effective thermal conductivity of Fe-BTC in H 2 was shown to decrease from 0.3 W m -1 K -1 at 300 K to 0.17 W m -1 K -1 at 80 K. The heat capacity of the adsorbent also decreased from 950 J kg -1 K -1 at 300 K to 250 J kg -1 K -1 at 100 K, indicating that an equivalent amount of energy at 100 K will cause a larger increase in temperature than at 300 K. This is clearly an important consideration in system design.
As an alternative, Han et al. [189] suggested increasing the intrinsic thermal conductivity of MOFs by reducing the mass of the nodes and shortening the linkers, but this met with modest experimental success. Ren et al. [190] suggested increasing the thermal conductivity through the deposition of MOF material in porous Ni foam. This also facilitated the handling of the adsorbent. Thermal conductivity values exceeding 0.3 W m -1 K -1 have been achieved through compaction with ENG [191]. At a compact density of 0.35 g cm -3 , neat MOF-5 powder had a thermal conductivity of 0.07 W m -1 K -1 . This increased to 0.08, 0.15, and 0.39 W m -1 K -1 after the addition of 1, 5, and 10 wt% ENG, respectively.
Directionality of the thermal conductivity due to compaction was reported by Fedchenia et al. [192]. They performed Hot Disk thermal conductivity measurements in three orthogonal directions on compacted MOF-5, incorporating 10 wt% ENG, with a density of 0.604 g cm -3 . A thermal conductivity of 0.286 W m -1 K -1 was found in the pressing direction, with higher values of 1.49-3.45 W m -1 K -1 in the perpendicular direction. The results were similar to those reported for a MgH 2 /ENG composite [193]. The heat transfer rate between the pelletised MOF-5 ? 10 wt% ENG and the Hot Disk sensor was in the range 644-782 W m -2 K -1 . This is an important consideration when integrating adsorbent pellets with a heat exchanger surface. Directionality of the thermal conductivity was also reported by Ming et al. [184]. A MOF-5 composite with 5 wt% ENG at a density of 0.35 g cm -3 was shown to have a thermal conductivity of 0.15 W m -1 K -1 in the pressing direction and 0.6 W m -1 K -1 in the perpendicular direction, due to preferential alignment of the ENG flakes. Figure 6 shows several heat exchanger concepts. In each case, the adsorbent is packed around the heat exchange material allowing H 2 to flow freely into the bed. The conventional tube-fin concept, shown in Fig. 6a, uses a series of tubes connected with Al fins to provide high thermal conductivity paths throughout the bed, which is packed between the plates and tubes. The spacing of the plates and tubes is critical, in order to minimise the mass, volume, and cost of the heat exchange system, while still providing adequate response to transient operating conditions. Optimisation can be performed using materials performance models, combined with finite element models to track mass and thermal flows within competing designs, as outlined by Hardy et al. [62] and Corgnale et al. [64].

Heat exchanger concepts
A variation of the tube-fin design uses an Al honeycomb for fins, as shown in Fig. 6b. The adsorbent is packed in the honeycomb, together with heating/cooling channels. Lightweight Al foil is used for construction, with the hexagonal cells having dimensions approximately 3-9 mm, flat to flat. The honeycomb heat exchanger has a simple design, low cost, and a low volume. Another benefit is the ability to flow gas through long channels of adsorbent, which is advantageous for cooling purposes. For 300 L of adsorbent, this heat exchanger only costs around US$100 and has a low volume of approximately 3.3 L.
Foams, as shown in Fig. 6c, d, in conjunction with heating/cooling tubes, are also being considered. Metal foams have been used successfully in metal hydride-based hydrogen reservoirs [194], but they are unsuitable for automotive applications due to their high cost (*US$280 L -1 ). Carbon foams are cheaper, with a comparable thermal conductivity, a lower density, and a low coefficient of thermal expansion. However, foams, in general, are difficult to impregnate with adsorbent, which is a major drawback.
Additives such as ENG can be used, either in a homogeneous distribution or inhomogeneously, as in layering, in order to increase its thermal conductivity. An example of a layered structure is shown in Fig. 6e. Such an approach has improved the thermal conductivity of MOF-5 by a factor of five at room temperature, with 10 % ENG layered perpendicularly to the pressing direction [184]. Layering can limit design options, since the layers must be aligned perpendicularly to the pressing direction, but this method of incorporating high thermal conductivity paths in highly insulating adsorbent powders may be advantageous.
Another heat exchanger concept uses microchannels etched into metal plates. This technology was originally pursued by Wegeng et al. [195]. The microchannel design uses a series of fins containing microscale channels (&200 lm) etched into them, as shown in Fig. 6f. The heat transfer fluid flows down one header tube, through the channelled fins, and up the return header, resulting in convective thermal energy transfer that keeps the surface temperature of the fins constant. The microchannel design is not limited by metal heat conduction, unlike conventional fins. The main benefits of the microchannel heat exchanger are its ability to work with compacted adsorbents and its relatively low mass and volume.

Flow-through cooling
In a flow-through system, LN 2 -chilled H 2 is used to cool the adsorbent bed while the H 2 is simultaneously being adsorbed [64]. This requires three to four times the volume of H 2 that is adsorbed, with the remainder returned to the fuelling station for re-cooling. This simplifies the on-board adsorption system since only heating lines need to be added to the system; however, there is a trade-off with greater complexity at the fuelling station because the chilled H 2 and its return must be accommodated.
A subscale prototype was demonstrated recently in a fully instrumented 0.5-L vessel containing a honeycomb heat exchanger, the details of which will be published separately. Thermocouples were arranged both on the Al cell walls and in the adsorbent bed. The system included ten stacked honeycomb units with a resistance heating rod inserted down the centre for discharge. To cool during charging, chilled gas was passed through the honeycomb channels and out the end of the vessel. The system can be cooled rapidly from ambient temperature, using a cold gas flow of 80 K H 2 at a rate of 95 L min -1 , chilled to 80 K at 6.0 MPa. The thermocouples mounted close to the central axis recorded initial warming due to adsorption, followed by rapid cooling. The embedded resistance heater can then be used to heat the bed during the discharge phase.

Discussion on storage systems
Using MOF-5, the majority of the US DOE technical targets can be achieved using advanced engineering methods and techniques [131,196]. H 2 can be charged and discharged rapidly, and storage systems are robust enough to operate under all reasonable terrestrial conditions with minimal impact on performance. The remaining targets needing effort include (1) gravimetric density, (2) volumetric density, (3) fuel cost, and (4) loss of useable H 2 .
Gravimetric density targets can only be achieved using new materials with higher gravimetric capacities or materials that do not require cooling, for example those having an adsorption enthalpy of 20-30 kJ mol -1 H 2 . However, this in turn could lead to fuelling station issues associated with dissipation of this heat during charging. In contrast, volumetric capacity can be partially addressed using more effective adsorbent consolidation methods. However, even at spherical powder packing densities, 40 % of the volume is free space and able to only accommodate H 2 at its gaseous density. Higher volumetric density adsorbents must therefore be developed. In automobile design, volume is more important than weight for several reasons, including the required changes to the vehicle frame design and its effects on range.
Fuel cost is directly associated with the use of LN 2 . H 2 cost could be reduced by using lower H 2 pressures of 6 MPa, rather than 70 MPa, which is becoming the standard for fuel cell vehicles. However, the addition of LN 2 infrastructure at the fuelling station is likely prohibitive. The need for cryogenic temperatures also affects H 2 loss. Even using multilayer vacuum insulation with a heat leak of less than 5 W, venting of a fail-safe full tank begins after 3 days. Calculations indicate significant H 2 will remain even under the harshest terrestrial conditions after 9 months; however, in many instances, the loss of H 2 in a static situation may be unacceptable. The only viable solution to the problem of high fuel cost and the loss of useable H 2 is the identification of adsorbents with higher enthalpies of adsorption.

Discussion and conclusion
This article has covered a number of the important aspects of research into hydrogen storage in nanoporous materials. The discovery and development of new materials undoubtedly plays a central role, with hydrogen storage being only one of many applications currently driving advances in materials chemistry [197]. However, the different experimental and computational approaches to their study for H 2 storage also require attention. Firstly, the measurement of hydrogen uptake by materials, as discussed in Sect. 2.1, can be subject to errors that have led to concerns regarding the accuracy and reproducibility of hydrogen uptake data. No formal guidelines are currently available to help ensure accurate hydrogen sorption measurements are taken by following accepted protocols. This can lead to inconsistencies in the data reported in the literature. Furthermore, once data are measured and the excess adsorption calculated, there are ambiguities in the definition of the absolute or total capacity of a material that can be extracted using different assumptions, as discussed in Sect. 2.2; although efforts are underway to develop a more consistent approach. The use of net adsorption, for example, could provide an interesting alternative. Nevertheless, absolute adsorption isotherms measured at different temperatures can be analytically represented by the models mentioned in Sect. 3.1. Studies in this area, reporting detailed isotherm fitting of H 2 adsorption isotherms measured at different temperatures, are relatively scarce compared to those reporting the H 2 storage capacities of a material or groups of materials at one or two temperatures, 77 and 298 K, for example. However, further studies validating the models for different adsorbents, particularly those exhibiting novel adsorption behaviour, would be valuable.
On a microscopic level, information can be obtained experimentally using, for example, neutron scattering techniques, as discussed in Sect. 2.3. There is clear scope for expanding the range of experimental conditions probed using neutrons-to higher temperatures, in particular-but also the wider use of H 2 instead of D 2 in diffraction studies.
The global development of higher-intensity neutron sources, with the concomitant improvement in neutron scattering instrumentation that can be expected, will certainly aid this task. Computational techniques are another option. Improvements in GCMC methodology to allow a better representation of amorphous or disordered and defective materials, for example, would be particularly valuable, as would the development of a transferable ''MOF force field'' that could accurately describe the H 2 -solid interactions for different groups of MOFs. With regard to quantum chemical calculations, the development of new, efficient methods that can accurately describe dispersion interactions, but with linear scaling with system size, for example, in terms of computational expense, would be particularly welcome. High-throughput computational screening of materials for H 2 storage is another area in which important developments can be expected in the near future.
With regard to materials, there is a need to carefully consider both their volumetric and gravimetric capacities, while the usable capacity is practically important and must not be overlooked. Methods of increasing the interaction strength of H 2 are also a key consideration for increasing the operational temperatures of nanoporous materials for H 2 storage. Possible approaches include the use of open metal sites in MOFs and functionalisation of their organic linkers using, for example, post-synthetic modification. Nanoporous molecular crystals, including intrinsically porous organic cage compounds, are an example of a potentially interesting new material type that should be investigated further for H 2 storage applications. Furthermore, MOFs exhibiting hysteretic H 2 adsorption behaviour and core-shell materials are also interesting targets for future research because they have the potential to significantly increase the usable capacity.
Beyond the discovery of new materials, the further investigation and development of H 2 storage tanks is also important. In Sect. 5, a number of relevant aspects were discussed, including methods of incorporating powders into large adsorbent beds, the thermal conductivity of materials, different heat exchange concepts, and the use of flow-through cooling. Both the gravimetric and volumetric capacities of practical tanks need to be increased. The volumetric capacity, which is a critical consideration for mobile applications, can be increased to a certain extent using alternative adsorbent consolidation methods and more efficient heat exchangers, but, ultimately, new highcapacity adsorbents with higher adsorption enthalpies are required. The use of such materials would also help ameliorate another problem with current storage units, which is the loss of useable H 2 due to the required use of cryogenic temperatures. Nevertheless, any further improvements in storage tank design that can be introduced would be valuable in order to maximise the performance of the current state-of-the-art materials in real scenarios. All of the topics covered in this article would clearly be valuable targets for future research in the field.
Acknowledgments Open access funding provided by Max-Planck-Institut für Intelligente Systeme. The authors acknowledge the contribution of the International Energy Agency (IEA) Hydrogen Implementing Agreement (HIA) from which this paper results, specifically the activities of Task 32: Hydrogen-based energy storage. Part of this research has been co-financed by the European Union (European Social Fund-ESF) and Greek national funds through the Operational Program ''Education and Lifelong Learning'' of the National Strategic Reference Framework (NSRF)-Research Funding Program: THALES. Part of the paper is also based upon work supported by the U. S. Department of Energy (National Nuclear Security Administration) under Award Number DE-FC36-09GO19006. Neither the United States Government nor any agency thereof, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favouring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://crea tivecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.