Abstract:We present an anharmonic approach for molecule-surface vibrations that employs rigid body coordinates, based on the Rodrigues rotation formula, to describe curvilinear displacements (rotations and translations) of the molecule along normal modes. These displacements are used to calculate energy data points from which one-dimensional polynomial potentials are fitted using cubic splines. In these potentials, for each of the six rigid body modes separately, one-dimensional Schrödinger equations are solved with harmonic oscillator or Fourier functions (for pure rotations) as basis sets. The anharmonic vibrational energies obtained are used to calculate partition functions and from them enthalpies, entropies, and Gibbs free energies of adsorption. Our numerical implementation has been successfully tested for Morse and cosine potentials with known analytical solutions. The methods have been applied to adsorption of CH4 on the hydroxyl group of the proton form of the chabazite zeolite (H-CHA), as well as to adsorption of CH4 and CO on the Mg2+ ions of the metal-organic framework (MOF) Mg2(dobdc). To obtain the best estimates for thermodynamic functions, we include the coupling between molecule-surface vibrations and intrasystem vibrations at the harmonic level. The calculated Gibbs free energies show that the coupling is small for CH4/H-CHA and CO/MOF (between -0.7 and +0.1 kJ/mol) but substantial for CH4/MOF (-3.4 kJ/mol). The predicted anharmonic effect on the Gibbs free energy of adsorption for CH4/H-CHA, CH4/MOF, and CO/MOF is -4.7, 0.3 ± 0.7, and -2.4 ± 0.6 kJ/mol, respectively, which results in +4.2, +0.9 ± 0.7, and -0.4 ± 0.6 kJ/mol, respectively, for the deviation from experiment. This is well within chemical accuracy limits (±4.2 kJ/mol) for the adsorption of CH4 and CO in the MOF. The larger deviation for CH4/H-CHA, at the edge of the chemical accuracy range, is most likely due to contributions from soft zeolite modes which are neglected in our approach.

Lyapunov instability of rigid diatomic molecules in three dimensions - a simpler method

Algebraic Chern Simons Theory

Computation of 1-D Shock Structure in A Gas in Rotational Non-Equilibrium Using A New Set of Simplified Burnett Equations

Assessment of random-phase approximation and second order Møller-Plesset perturbation theory for many-body interactions in solid ethane, ethylene, and acetylene

Density ordering instabilities of quasi-two-dimensional fermionic polar molecules in single-layer and multi-layer configurations: exact treatment of exchange interactions

Two-dimensional Electronic Spectra from the Hierarchical Equations of Motion Method: Application to Model Dimers.

Electromagnetism, Optics, Acoustics, Heat Transfer, Classical Mechanics, and Fluid Dynamics Multiple-Relaxation-Time Lattice Boltzmann Approach to Richtmyer---Meshkov Instability

The {\it ab initio} calculation of spectra of open shell diatomic molecules

Molecular-Dynamics Simulation Of Richtmyer-Meshkov Instability On A Li-H-2 Interface At Extreme Compressing Conditions

Multiple-Relaxation-Time Lattice Boltzmann approach to richtmyer-Meshkov instability

Rigid Body Approximation for the Anharmonic Description of Molecule-Surface Vibrations

Rotation–vibration Motion of Pyramidal XY3molecules Described in the Eckart Frame: Theory and Application to NH3

Molecular Structure Theory without the Born-Oppenheimer Approximation: Rotationless Vibrational States of LiH

B 97 , 134513 ( 2018 ) Instability of Fulde-Ferrell-Larkin-Ovchinnikov states in atomic Fermi gases in three and two dimensions

A well-tempered density functional theory of electrons in molecules

Molecular dynamics with quantum forces: Vibrational spectra of localized systems

An energy-modified quantum defect method for the analysis of Rydberg spectra: Application to 2-butyne

Computing linear optical spectra in the presence of nonadiabatic effects on Graphics Processing Units using molecular dynamics and tensor-network approaches

Structure Relaxation Via Long Trajectories Made Stable

ELECTRONIC STRUCTURE OF MOLECULAR RYDBERG STATES OF SOME SMALL MOLECULES AND MOLECULAR ION

Linear Response pCCD-Based Methods: LR-pCCD and LR-pCCD+S Approaches for the Efficient and Reliable Modeling of Excited state Properties