Stability, adsorption, and diffusion of CH4, CO2, and H2 in clathrate hydrates.

We present a study of the adsorption and diffusion of CH4, CO2, and H2 molecules in clathrate hydrates using ab initio van der Waals density functional formalism [M. Dion, Phys. Rev. Lett. 92, 246401 (2004)10.1103/PhysRevLett.92.246401]. We find that the adsorption energy is dominated by van der Waals interactions and that, without them, gas hydrates would not be stable. We calculate the maximum adsorption capacity as well as the maximum hydrocarbon size that can be adsorbed. The relaxation of the host lattice is essential for a good description of the diffusion activation energies, which are estimated to be of the order of 0.2, 0.4, and 1.0 eV for H2, CO2, and CH4, respectively.

Efficient implementation of a van der Waals density functional: application to double-wall carbon nanotubes.

We present an efficient implementation of the van der Waals density functional of Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)], which expresses the nonlocal correlation energy as a double spatial integral. We factorize the integration kernel and use fast Fourier transforms to evaluate the self-consistent potential, total energy, and atomic forces, in O(NlogN) operations. The resulting overhead, for medium and large systems, is a small fraction of the total computational cost, representing a dramatic speedup over the O(N(2)) evaluation of the double integral. This opens the realm of first-principles simulations to the large systems of interest in soft matter and biomolecular problems. We apply the method to calculate the binding energies and the barriers for relative translation and rotation in double-wall carbon nanotubes.

Energetics and dynamics of H(2) adsorbed in a nanoporous material at low temperature.

Molecular hydrogen adsorption in a nanoporous metal-organic framework structure (MOF-74) is studied via van der Waals density-functional calculations. The primary and secondary binding sites for H(2) are confirmed. The low-lying rotational and translational energy levels are calculated, based on the orientation and position dependent potential energy surface at the two binding sites. A consistent picture is obtained between the calculated rotational-translational transitions for different H(2) loadings and those measured by inelastic neutron scattering exciting the singlet to triplet (para to ortho) transition in H(2). The H(2) binding energy after zero-point energy correction due to the rotational and translational motions is predicted to be approximately 100 meV in good agreement with the experimental value of approximately 90 meV.