Abstract:We show how we can find the enthalpy of a chemical reaction under non-ideal conditions using the Small System Method to sample molecular dynamics simulation data for fluctuating variables. This method, created with Hill's thermodynamic analysis, is used to find properties in the thermodynamic limit, such as thermodynamic correction factors, partial enthalpies, volumes, heat capacities and compressibility. The values in the thermodynamic limit at (T,V, μj) are then easily transformed into other ensembles, (T,V,Nj) and (T,P,Nj), where the last ensemble gives the partial molar properties which are of interest to chemists. The dissociation of hydrogen from molecules to atoms was used as a convenient model system. Molecular dynamics simulations were performed with three densities; ρ = 0.0052 g cm(-3) (gas), ρ = 0.0191 g cm(-3) (compressed gas) and ρ = 0.0695 g cm(-3) (liquid), and temperatures in the range; T = 3640-20,800 K. The enthalpy of reaction was observed to follow a quadratic trend as a function of temperature for all densities. The enthalpy of reaction was observed to only have a small pressure dependence. With a reference point close to an ideal state (T = 3640 K and ρ = 0.0052 g cm(-3)), we were able to calculate the thermodynamic equilibrium constant, and thus the deviation from ideal conditions for the lowest density. We found the thermodynamic equilibrium constant to increase with increasing temperature, and to have a negligible pressure dependence. Taking the enthalpy variation into account in the calculation of the thermodynamic equilibrium constant, we found the ratio of activity coefficients to be in the order of 0.7-1.0 for the lowest density, indicating repulsive forces between H and H2. This study shows that the compressed gas- and liquid density values at higher temperatures are far from those calculated under ideal conditions. It is important to have a method that can give access to partial molar properties, independent of the ideality of the reacting mixture. Our results show how this can be achieved with the use of the Small System Method.

Molecular Dissociation in Hot, Dense Hydrogen

Dissociation of Fluid Hydrogen at High Pressure and Temperature

Quantum Molecular Dynamics Study on the Proton Exchange, Ionic Structures, and Transport Properties of Warm Dense Hydrogen-Deuterium Mixtures.

Bound state breaking and the importance of thermal exchange-correlation effects in warm dense hydrogen

Distinct Metallization and Atomization Transitions in Dense Liquid Hydrogen

Isotopic and vibrational-level dependence of H$_2$ dissociation by electron impact

Thermodynamic anomalies and three distinct liquid-liquid transitions in warm dense liquid hydrogen

High-pressure Raman Study of Solid Hydrogen Up to 300 GPa

On the Room-Temperature Phase Diagram of High Pressure Hydrogen: an Ab Initio Molecular Dynamics Perspective and A Diffusion Monte Carlo Study

The reaction enthalpy of hydrogen dissociation calculated with the Small System Method from simulation of molecular fluctuations

An Apparent Dissociation Transition in Anharmonically Bound 1D Systems

Simulation of dissociation effect at high temperature and high pressure by REMC method

Two-dimensional molecular para-hydrogen and ortho-deuterium at zero temperature

Dense Molecular Hydrogen: Order, Disorder, and Localization

Black-body radiation induced photodissociation and population redistribution of weakly bound states in H$_2^+$

Quantum Monte Carlo Simulation of the High-Pressure Molecular-Atomic Crossover in Fluid Hydrogen

Ground-state properties of one-dimensional matter and quantum dissociation of a Luttinger liquid

Benchmarking vdW-DF first principle predictions against Coupled Electron-Ion Monte Carlo for high pressure liquid hydrogen

Experimental phase transition mapping for hydrogen above 300 K up to 300 GPa

Molecular dynamics study of CO2 hydrate dissociation: Fluctuation-dissipation and non-equilibrium analysis

Structure and phase boundaries of compressed liquid hydrogen