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.

The average bond-length of diatomic molecules in thermodynamical equilibrium depends on the volume

An Apparent Dissociation Transition in Anharmonically Bound 1D Systems

Systematic Comparison of Experimental, Quantum Mechanical, and Molecular Mechanical Bond Lengths for Organic Molecules

Deduction of Bond Length Changes of Symmetric Molecules from Experimental Vibrational Progressions, Including A Topological Mass Factor

Electronic correlations influence on the anomalous constitutive law between the structure and spectra of hydrogen bonds

Intermolecular distance and density scaling of dynamics in molecular liquids

On the effect of nuclear motion on the interaction between diatomic molecules

Long-range interactions between polar bialkali ground-state molecules in arbitrary vibrational levels

Upper Bound for the Free Energy of Dilute Bose Gases at Low Temperature

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

Breakdown of the Total Dipole Moments of Diatomic Molecules into Individual Orbital Contributions

Probabilistic Approach to the Length-Scale Dependence of the Effect of Water Hydrogen Bonding on Hydrophobic Hydration

Diatomic molecules of alkali-metal and alkaline-earth-metal atoms: interaction potentials, dipole moments, and polarizabilities

Dynamics of molecular collisions in air and its mean free path

Dynamics of dilute gases: a statistical approach

An efficient memetic, permutation-based evolutionary algorithm for real-world train timetabling

An ab initio study on boundaries for characterizing cooperative effect of hydrogen bonds by intermolecular compression

Threshold detachment of negative ions by electron impact

Ground state energy of dilute Bose gases in 1D

Are polar liquids less simple?

Dipolar Droplets at 3D-1D Crossover