Microcanonical thermodynamics of three and four atoms

Resolving the Debate between Boltzmann and Gibbs Entropy: Relative Energy Window Eliminates Thermodynamic Inconsistencies and Allows Negative Absolute Temperatures

Small systems consisting of a few particles are increasingly technologically relevant. In such systems, an intense debate in microcanonical statistical mechanics has been about the correctness of Boltzmann's surface entropy versus Gibbs' volume entropy. Both entropies have shortcomings─while Boltzmann entropy predicts unphysical negative/infinite absolute temperatures for small systems with an unbounded energy spectrum, Gibbs entropy entirely disallows negative absolute temperatures, in disagreement with experiments. We consider a relative energy window, motivated by the Heisenberg energy-time uncertainty principle and eigenstate thermalization in quantum mechanics. The proposed entropy ensures positive, finite temperatures for systems without a maximum limit on their energy and allows negative absolute temperatures in bounded energy spectrum systems, e.g., with population inversion. It also closely matches canonical ensemble predictions for prototypical systems, thus correctly describing the zero-point energy of an isolated quantum harmonic oscillator. Overall, we enable accurate thermodynamic models for isolated systems with few degrees of freedom.

Classical statistical mechanics in the μVL and μpR ensembles

Molecular expressions for thermodynamic properties and derivatives of the entropy up to third order in the adiabatic grand-isochoric μVL and adiabatic grand-isobaric μpR ensembles are systematically derived using the methodology developed by Lustig for the microcanonical and canonical ensembles [J. Chem. Phys. 100, 3048 (1994)10.1063/1.466446; Mol. Phys. 110, 3041 (2012)10.1080/00268976.2012.695032]. They are expressed by phase-space functions, which represent derivatives of the entropy with respect to the chemical potential, the volume, and the Hill energy L in the μVL ensemble and with respect to the chemical potential, the pressure, and the Ray energy R in the μpR ensemble. The derived expressions are validated for both ensembles by Monte Carlo simulations for the simple Lennard-Jones model fluid at three selected state points.

Rigorous expressions for thermodynamic properties in the NpH ensemble

Molecular expressions for thermodynamic properties of fluids and derivatives of the entropy up to third order in the isoenthalpic-isobaric ensemble are derived by using the methodology developed by Lustig for the microcanonical and canonical ensembles [J. Chem. Phys. 100, 3048 (1994)JCPSA60021-960610.1063/1.466446; Mol. Phys. 110, 3041 (2012)MOPHAM0026-897610.1080/00268976.2012.695032]. They are expressed in a systematic way by phase-space functions, which represent derivatives of the phase-space volume with respect to enthalpy and pressure. The expressions for thermodynamic properties contain only ensemble averages of combinations of the kinetic energy and volume of the system. Thus, the calculation of thermodynamic properties in the isoenthalpic-isobaric ensemble does not require volume derivatives of the potential energy. This is particularly advantageous in Monte Carlo simulations when the interactions between molecules are described by very accurate ab initio pair and nonadditive three-body potentials. The derived expressions are validated by Monte Carlo simulations for the simple Lennard-Jones model fluid as a test case.