Thermodynamic properties of hydrogen dissociation reaction from the small system method and reactive force field ReaxFF

Significance Statement

Thermodynamic attributes are important in chemical engineering and chemistry in general. These properties have been used to describe macroscopic systems, in laboratory chemical reactor scales. The descriptions are insufficient at the nanoscale owing to the fact that thermodynamic variables are not extensive when the particle number drops to a few particles.But the Small System Method can be used to connect properties of the system in the nanoscale with macroscopic limit. It has now been applied to compute thermodynamic correction factors, partial molar enthalpies, derivatives of activity coefficients, and reaction enthalpies of macroscopic systems. The results obtained in all these cases indicate agreement with data obtained from other approaches.

Thermodynamic factors enable us to compute Maxwell-Stefan diffusion coefficients from Fick diffusion coefficients for both binary and ternary systems. Most of the computations have been based on non-reactive systems. Therefore, obtaining thermodynamic data in reacting systems would be groundbreaking for this method.

Researchers led by Professor Signe Kjelstrup at the Norwegian University of Science and Technology demonstrated a combination of the small system method and the Reactive Force Field in a bid to obtain thermodynamic data of the hydrogen dissociation reaction. In their work, they provided a number of properties of the reaction such as heat of reaction, heat capacity, thermodynamic correction factors with the reactive force field. They compared their results with the already obtained values from the three-body interaction potential. Their work is published in peer-reviewed journal, Chemical Physics Letters.

The development of accurate thermodynamic methods and precise potentials for the dissociation of hydrogen gas is the starting point for a number of complex computations. The dissociation normally occur at high temperatures, therefore the need for computer simulations. The authors simulated a cubic box system with periodic boundary conditions in three dimensions.

The researchers varied the box size from 28.47 Å-100 Å and the number of hydrogen molecules from 37-1500. They then performed NVT simulations with 0.1 fs step size. They equilibrated the system with 1 ns run at 300 K.

In general, the authors applied the Small System Method to the hydrogen dissociation reaction, implementing Reactive Force Field to obtain thermodynamic data. The main aim of the study was to compare the results from the Reactive Force Field with the results obtained earlier from the accurate interaction potential. They reported heat capacity, thermodynamic correction factors, activation energy, and reaction enthalpy.

The results of their study indicated that the two approaches of computing thermodynamic data were in agreement. However, simplifying ideal mixture approximation was considered for enthalpies of reaction. Three body potential outcomes appeared more precise, but the Reactive Force Field gave the trends and magnitudes that were expected.

These results open the possibility of combining small system method and reactive force field for more complex reactions. The approach would be helpful in bringing important information about reactions away from equilibrium, and studying transport characteristics in mixtures where chemical reactions occur.

Thermodynamic properties of hydrogen dissociation reaction from the small system method and reactive force field ReaxFF- Advances in Engineering

About The Author

Dr. Thuat T. Trinh

Dr. Trinh is working at Department of Civil and Environmental Engineering, Norwegian University of Science and Technology, NTNU, Norway.  He earned his PhD degree of Chemistry and Chemical Engineering at Eindhoven Technical University, The Netherlands. He also had some postdoctoral experience in Ecole Normale Superieuse de Lyon and Dept. of Chemistry, NTNU. His research in atomistic multi-scale modeling approach (quantum mechanics, molecular dynamics, Monte Carlo and continuum description) to study chemical, physical, thermodynamic and mechanical properties of complex systems.

About The Author

Signe Kjelstrup is a professor of physical chemistry since 1985 at the Norwegian University of Science and Technology, Trondheim. Since 2005, she holds a part time chair on irreversible thermodynamics and sustainable processes at the Technical University of Delft, The Netherlands. Her works in non-equilibrium thermodynamics concern electrochemical cells, membrane systems and entropy production minimization in process equipment. She is honorary doctor of the University of North East China, and has been a guest professor of Kyoto University. Her book on irreversible thermodynamics coauthored with K.S. Forland and T. Forland (Wiley, 1988 and 1994, Tapir 2001) has been translated into Japanese and Chinese. Together with Dick Bedeaux she wrote the World-Scientific bestseller “Non-equilibrium Thermodynamics of Heterogeneous Systems” in 2008.


Thuat T. Trinh, Nora Meling, Dick Bedeaux, Signe Kjelstrup. Thermodynamic properties of hydrogen dissociation reaction from the small system method and reactive force field ReaxFF. Chemical Physics Letters, volume 672 (2017), pages 128–132.

Go To Chemical Physics Letters