Guillermo Roman-Perez

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.

Density, structure, and dynamics of water: The effect of van der Waals interactions

It is known that ab initio molecular dynamics (AIMD) simulations of liquid water at ambient conditions, based on the generalized gradient approximation (GGA) to density functional theory (DFT), with commonly used functionals fail to produce structural and diffusive properties in reasonable agreement with experiment. This is true for canonical, constant temperature simulations where the density of the liquid is fixed to the experimental density. The equilibrium density, at ambient conditions, of DFT water has recently been shown by Schmidt et al. [J. Phys. Chem. B, 113, 11959 (2009)] to be underestimated by different GGA functionals for exchange and correlation, and corrected by the addition of interatomic pair potentials to describe van der Waals (vdW) interactions. In this contribution we present a DFT-AIMD study of liquid water using several GGA functionals as well as the van der Waals density functional (vdW-DF) of Dion et al. [Phys. Rev. Lett. 92, 246401 (2004)]. As expected, we find that the density of water is grossly underestimated by GGA functionals. When a vdW-DF is used, the density improves drastically and the experimental diffusivity is reproduced without the need of thermal corrections. We analyze the origin of the density differences between all the functionals. We show that the vdW-DF increases the population of non-H-bonded interstitial sites, at distances between the first and second coordination shells. However, it excessively weakens the H-bond network, collapsing the second coordination shell. This structural problem is partially associated to the choice of GGA exchange in the vdW-DF. We show that a different choice for the exchange functional is enough to achieve an overall improvement both in structure and diffusivity.

Flexibility in a Metal–Organic Framework Material Controlled by Weak Dispersion Forces: The Bistability of MIL‐53(Al)

Breathtaking MOFs: DFT calculations reveal that the exceptional, thermally induced density change of the metal-organic framework MIL53(Al) is controlled by a competition between shortand long-range interactions and entropic factors. As shown in the picture (C green, Al cyan, O red, H white), dispersive interactions between the phenyl rings are responsible for stabilizing a narrow-pore form at low temperature. At 325-375 K, vibrational entropy causes the structure to expand markedly, permitting large volumes of light gases to be adsorbed.

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.