Abstract:Quantum and anharmonic effects are investigated in (H2)2-Li(+)-benzene, a model for hydrogen adsorption in metal-organic frameworks and carbon-based materials, using rigid-body diffusion Monte Carlo (RBDMC) simulations. The potential-energy surface (PES) is calculated as a modified Shepard interpolation of M05-2X/6-311+G(2df,p) electronic structure data. The RBDMC simulations yield zero-point energies (ZPE) and probability density histograms that describe the ground-state nuclear wavefunction. Binding a second H2 molecule to the H2-Li(+)-benzene complex increases the ZPE of the system by 5.6 kJ mol(-1) to 17.6 kJ mol(-1). This ZPE is 42% of the total electronic binding energy of (H2)2-Li(+)-benzene and cannot be neglected. Our best estimate of the 0 K binding enthalpy of the second H2 to H2-Li(+)-benzene is 7.7 kJ mol(-1), compared to 12.4 kJ mol(-1) for the first H2 molecule. Anharmonicity is found to be even more important when a second (and subsequent) H2 molecule is adsorbed; use of harmonic ZPEs results in significant error in the 0 K binding enthalpy. Probability density histograms reveal that the two H2 molecules are found at larger distance from the Li(+) ion and are more confined in the θ coordinate than in H2-Li(+)-benzene. They also show that both H2 molecules are delocalized in the azimuthal coordinate, ϕ. That is, adding a second H2 molecule is insufficient to localize the wavefunction in ϕ. Two fragment-based (H2)2-Li(+)-benzene PESs are developed. These use a modified Shepard interpolation for the Li(+)-benzene and H2-Li(+)-benzene fragments, and either modified Shepard interpolation or a cubic spline to model the H2-H2 interaction. Because of the neglect of three-body H2, H2, Li(+) terms, both fragment PESs lead to overbinding of the second H2 molecule by 1.5 kJ mol(-1). Probability density histograms, however, indicate that the wavefunctions for the two H2 molecules are effectively identical on the "full" and fragment PESs. This suggests that the 1.5 kJ mol(-1) error is systematic over the regions of configuration space explored by our simulations. Notwithstanding this, modified Shepard interpolation of the weak H2-H2 interaction is problematic and we obtain more accurate results, at considerably lower computational cost, using a cubic spline interpolation. Indeed, the ZPE of the fragment-with-spline PES is identical, within error, to the ZPE of the full PES. This fragmentation scheme therefore provides an accurate and inexpensive method to study higher hydrogen loading in this and similar systems.

From Clusters to Liquid: What Are the Preferred Ways for Benzene and Pyrrole to Interact?

A Theoretical Study of N–H ··· Π H-bond Interaction of Pyrrole: from Clusters to the Liquid

Natures of Benzene-Water and Pyrrole-Water Interactions in the Forms of Σ and Π Types: Theoretical Studies from Clusters to Liquid Mixture

Simultaneous interactions of pyridine with substituted benzene ring and H–F in X-ben⊥pyr···H–F complexes: a cooperative study

Ab Initio Study of the Π-Π Interactions Between CO2 and Benzene, Pyridine, and Pyrrole

The study on interactions between 1‐ethyl‐3‐methylimidazolium chloride and benzene/pyridine /pyrrole/thiophene

Ab Initio Study of the Π–π Interaction in the Pyridine Dimer: Effects of Bond Functions

Molecular Interactions Between Pyrazine and N-Propanol, Chloroform, or Tetrahydrofuran.

Intermolecular Interactions Between Σ- and Π-Holes of Bromopentafluorobenzene and Pyridine: Computational and Experimental Investigations

Intermolecular Π/π and H/π Interactions in Dimers Researched by Different Computational Methods

Differentiation of Cation−π Bonding from Cation−π Intermolecular Interactions: A Quantum Chemistry Study Using Density-Functional Theory and Morokuma Decomposition Methods

Vapor-liquid and vapor-solid phase equilibria for united-atom benzene models near their triple points: the importance of quadrupolar interactions.

Exploring the nature of interactions among thiophene, thiophene sulfone, dibenzothiophene, dibenzothiophene sulfone and a pyridinium-based ionic liquid.

Theoretical Insight into the Interactions of TMA-Benzene and TMA-Pyrrole with B3LYP Density-Functional Theory (DFT) and ab Initio Second Order Møller−Plesset Perturbation Theory (MP2) Calculations

Effect of methyl group on the cooperativity between cation–π interaction and NH···O hydrogen bonding

A Study of Pi-Pi Interaction in High Performance Liquid Chromatography by Molecular Modeling

Quantitative study of interactions between oxygen lone pair and aromatic rings: substituent effect and the importance of closeness of contact

Assessment of Some Density Functional Theory Methods and Force Field Models in Describing Various Interaction Modes of Benzene Dimer

Coarse-grained interaction potentials for polyaromatic hydrocarbons

"Plug-and-Play" potentials: Investigating quantum effects in (H2)2-Li(+)-benzene

Theoretical Study on Interactions Between Pyridinium-Based Ionic Liquids and Thiophenic Compounds