No Arabic abstract
We have obtained by Monte Carlo NVT simulations the constant-volume excess heat capacity of square-well fluids for several temperatures, densities and potential widths. Heat capacity is a thermodynamic property much more sensitive to the accuracy of a theory than other thermodynamic quantities, such as the compressibility factor. This is illustrated by comparing the reported simulation data for the heat capacity with the theoretical predictions given by the Barker-Henderson perturbation theory as well as with those given by a non-perturbative theoretical model based on Baxters solution of the Percus-Yevick integral equation for sticky hard spheres. Both theories give accurate predictions for the equation of state. By contrast, it is found that the Barker-Henderson theory strongly underestimates the excess heat capacity for low to moderate temperatures, whereas a much better agreement between theory and simulation is achieved with the non-perturbative theoretical model, particularly for small well widths, although the accuracy of the latter worsens for high densities and low temperatures, as the well width increases.
The chemical potentials of multicomponent fluids are derived in terms of the pair correlation functions for arbitrary number of components, interaction potentials, and dimensionality. The formally exact result is particularized to hard-sphere mixtures with zero or positive nonadditivity. As a simple application, the chemical potentials of three-dimensional additive hard-sphere mixtures are derived from the Percus-Yevick theory and the associated equation of state is obtained. This Percus-Yevick chemical-route equation of state is shown to be more accurate than the virial equation of state. An interpolation between the chemical-potential and compressibility routes exhibits a better performance than the well-known Boublik-Mansoori-Carnahan-Starling-Leland equation of state.
The equilibrium properties of a Janus fluid made of two-face particles confined to a one-dimensional channel are revisited. The exact Gibbs free energy for a finite number of particles $N$ is exactly derived for both quenched and annealed realizations. It is proved that the results for both classes of systems tend in the thermodynamic limit ($Ntoinfty$) to a common expression recently derived (Maestre M A G and Santos A 2020 J Stat Mech 063217). The theoretical finite-size results are particularized to the Kern--Frenkel model and confirmed by Monte Carlo simulations for quenched and (both biased and unbiased) annealed systems.
Heat engines used to output useful work have important practical significance, which, in general, operate between heat baths of infinite size and constant temperature. In this paper we study the efficiency of a heat engine operating between two finite-size heat sources with initial temperature differences. The total output work of such heat engine is limited due to the finite heat capacity of the sources. We investigate the effects of different heat capacity characteristics of the sources on the heat engines efficiency at maximum work (EMW) in the quasi-static limit. In addition, we study the efficiency of the engine working in finite-time with maximum power of each cycle is achieved and find the efficiency follows a simple universality as $eta=eta_{mathrm{C}}/4+Oleft(eta_{mathrm{C}}^{2}right)$. Remarkably, when the heat capacity of the heat source is negative, such as the black holes, we show that the heat engine efficiency during the operation can surpass the Carnot efficiency determined by the initial temperature of the heat sources. It is further argued that the heat engine between two black holes with vanishing initial temperature difference can be driven by the energy fluctuation. The corresponding EMW is proved to be $eta_{mathrm{EMW}}=2-sqrt{2}$, which is two time of the maximum energy release rate $mu=left(2-sqrt{2}right)/2approx0.29$ of two black hole emerging process obtained by S. W. Hawking.
The equilibrium properties of a Janus fluid confined to a one-dimensional channel are exactly derived. The fluid is made of particles with two faces (active and passive), so that the pair interaction is that of hard spheres, except if the two active faces are in front of each other, in which case the interaction has a square-well attractive tail. Our exact solution refers to quenched systems (i.e., each particle has a fixed face orientation), but we argue by means of statistical-mechanical tools that the results also apply to annealed systems (i.e., each particle can flip its orientation) in the thermodynamic limit. Comparison between theoretical results and Monte Carlo simulations for quenched and annealed systems, respectively, shows an excellent agreement.
The supercooled liquid silicon, modeled by Stillinger-Weber potential, shows anomalous increase in heat capacity $C_p$, with a maximum $C_p$ value close to 1060 K at zero pressure. We study equilibration and relaxation of the supercooled SW Si, in the temperature range of 1060 K--1070 K at zero pressure. We find that as the relaxation of the metastable supercooled liquid phase initiates, a straight line region (SLR) is formed in cumulative potential energy distributions. The configurational temperature corresponding to the SLR is close to 1060 K, which was earlier identified as the freezing temperature of 4-coordinated amorphous network. The SLR is found to be tangential to the distribution of the metastable liquid phase and thus influences the broadness of the distribution. As the bath temperature is reduced from 1070 K to 1060 K, the effective temperature approaches the bath temperature which results in broadening of the metastable phase distribution. This, in turn, causes an increase in overall fluctuations of potential energy and hence an increase of heat capacity. We also find that during initial stages of relaxation, 4-coordinated atoms form 6-membered rings with a chair--like structure and other structural units that indicate crystallization. Simultaneously a strong correlation is established between the number of chair-shaped 6-membered rings and the number of 4-coordinated atoms in the system. This shows that all properties related to 4-coordinated particles are highly correlated as the SLR is formed in potential energy distributions and this can be interpreted as a consequence of `freezing of amorphous network formed by 4-coordinated particles.