Path integral Monte Carlo simulations of H{sub 2} adsorbed to lithium-doped benzene: A model for hydrogen storage materials
- School of Chemistry, The University of Sydney, Sydney (Australia)
- Department of Chemistry, University of Canterbury, Christchurch (New Zealand)
Finite temperature quantum and anharmonic effects are studied in H{sub 2}–Li{sup +}-benzene, a model hydrogen storage material, using path integral Monte Carlo (PIMC) simulations on an interpolated potential energy surface refined over the eight intermolecular degrees of freedom based upon M05-2X/6-311+G(2df,p) density functional theory calculations. Rigid-body PIMC simulations are performed at temperatures ranging from 77 K to 150 K, producing both quantum and classical probability density histograms describing the adsorbed H{sub 2}. Quantum effects broaden the histograms with respect to their classical analogues and increase the expectation values of the radial and angular polar coordinates describing the location of the center-of-mass of the H{sub 2} molecule. The rigid-body PIMC simulations also provide estimates of the change in internal energy, ΔU{sub ads}, and enthalpy, ΔH{sub ads}, for H{sub 2} adsorption onto Li{sup +}-benzene, as a function of temperature. These estimates indicate that quantum effects are important even at room temperature and classical results should be interpreted with caution. Our results also show that anharmonicity is more important in the calculation of U and H than coupling—coupling between the intermolecular degrees of freedom becomes less important as temperature increases whereas anharmonicity becomes more important. The most anharmonic motions in H{sub 2}–Li{sup +}-benzene are the “helicopter” and “ferris wheel” H{sub 2} rotations. Treating these motions as one-dimensional free and hindered rotors, respectively, provides simple corrections to standard harmonic oscillator, rigid rotor thermochemical expressions for internal energy and enthalpy that encapsulate the majority of the anharmonicity. At 150 K, our best rigid-body PIMC estimates for ΔU{sub ads} and ΔH{sub ads} are −13.3 ± 0.1 and −14.5 ± 0.1 kJ mol{sup −1}, respectively.
- OSTI ID:
- 22493251
- Journal Information:
- Journal of Chemical Physics, Vol. 143, Issue 19; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
- Country of Publication:
- United States
- Language:
- English
Similar Records
“Plug-and-Play” potentials: Investigating quantum effects in (H{sub 2}){sub 2}–Li{sup +}–benzene
Investigation of the effect of local atomic polarization on field-enhanced ion transport in insulating binary oxides: The case of
Related Subjects
ORGANIC
PHYSICAL AND ANALYTICAL CHEMISTRY
ADSORPTION
BENZENE
COMPUTERIZED SIMULATION
CORRECTIONS
DEGREES OF FREEDOM
DENSITY FUNCTIONAL METHOD
DOPED MATERIALS
ENTHALPY
EXPECTATION VALUE
HARMONIC OSCILLATORS
HYDROGEN
HYDROGEN STORAGE
LITHIUM
LITHIUM IONS
MOLECULES
MONTE CARLO METHOD
PATH INTEGRALS
POTENTIAL ENERGY
SURFACES
TEMPERATURE DEPENDENCE