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.

Quantum Monte Carlo study of Be2 and group 12 dimers M2 (M = Zn, Cd, Hg)

A Systemic Dft Study on Several 5D-Electron Element Dimers: Hf-2, Ta-2, Re-2, W-2, and Hg-2

Application of the CC(P;Q) Hierarchy of Coupled -Cluster Methods to the Beryllium Dimer

A quantum Monte Carlo study of systems with effective core potentials and node nonlinearities

Structure and Analytic Potential Energy Function of the Molecules Beh, H-2 and Beh2

Structural and electronic properties of medium-sized beryllium doped magnesium BeMg clusters and their anions

Density Functional Theory Study of Structure and Electronic Properties of MgBe N ( N = 2-12 ) Clusters

Spectroscopic/Bond Property Relationship in Group 11 Dihydrides via Relativistic Four-Component Methods

Diffusion quantum Monte Carlo study on magnesium clusters as large as nanoparticles

Properties and Stabilities of MX, MX2, and M2X2 Compounds (M = Zn, Cd, Hg; X = F, Cl, Br, I)

Quantum mechanical study of transition metal hydrides: Comparison of determined molecular properties with experimental data

Ab initio calculations of the spectra and lifetimes of the lead dimer

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

Quantum Chemistry Study on Chemical Bonds of Transition-Metal Dimers

Structural and Electronic Properties of Medium-Sized Beryllium Doped Magnesium BeMgn Clusters and Their Anions

Analyzing and Modeling the Interaction Potential of the Ground-State Beryllium Dimer

A DFT Study on Potential High-Energy Compounds:H2n5mn5h2(m=be,mg,ca,zn,and Cd)

A Globally Accurate Potential Energy Surface and Quantum Dynamics Calculations on the Be(1S) + H2(v0 = 0, J0 = 0) → BeH + H Reaction

Unusual chemical bond and spectrum of beryllium dimer in ground $X^1Σ_g^+$ state

Bond−Energy−Electron Relaxation of BeN Nanoclusters and BeX Alloys

Ground and Excited States of [be–xe]+: A Multireference Configuration Interaction Study