Production of Useful Work in an Isothermal Cycle: A Monte Carlo Simulation of a Locally Nonchaotic Energy Barrier
Y. Qiao,1,2,* Z. Shang1 1 Program of Materials Science and Engineering, University of California – San Diego, La Jolla, CA 92093, U.S.A. 2 Department of Structural Engineering, University of California – San Diego, La Jolla, CA 92093-0085, U.S.A. * Corresponding author (Email: [email protected]; phone: 858-534-3388)
Abstract: A Monte Carlo simulation is performed on a billiard-type model system, which contains a locally nonchaotic energy barrier. Without extensive particle collision across the energy barrier, the system steady state is nonequilibrium; that is, the particles follow a non-Boltzmann distribution. Remarkably, as the energy barrier varies in an isothermal cycle, the total produced work is greater than the total consumed work, because of the asymmetry in the cross-influence of the thermally correlated thermodynamic driving forces. Such a phenomenon cannot be explained by the second law of thermodynamics. Similar anomalous effects may be achieved if the barrier is switchable or asymmetric. In essence, the energy barrier is a spontaneously nonequilibrium dimension. It is fundamentally different from Maxwell’s demon, unrelated to the physical nature of information.
KEYWORDS: Nonequilibrium; Nonchaotic; The second law of thermodynamics; Energy barrier; Monte Carlo simulation
1. Introduction
Randomness is a fundamental concept in statistical mechanics [1,2]. An ergodic and chaotic system always reaches thermodynamic equilibrium, while a nonergodic or nonchaotic system may not [e.g., 3-5]. The latter should be analyzed in the framework of nonequilibrium stochastic thermodynamics [6].
1
Figure 1. (A) In a vertical 푦 − 푧 plane in a gravitational field (푔), if 푧̂ ≪ , at the steady state the particle density does not follow the Maxwell-Boltzmann distribution. (B) The billiard-type model system, wherein elastic particles randomly move in the horizontal 푥 − 푦 dimension in between a upper “plateau” and a lower “plain”, across a narrow transition step. If 푧̂ ≪ , the steady-state ⁄ particle density ratio ( = ) is inherently nonequilibrium.
One example is shown in Figure 1(A), in which billiard-like elastic particles randomly move in a vertical 푦 − 푧 plane in a gravitational field (푔). The 푦 dimension is infinitely large. The upper boundary (푧 = 푧̂) is a specular wall, with 푧̂ being the plane height. The lower boundary (푧 = 0) is a diffusive wall at constant temperature (푇), where the reflected particle velocity (푣) is governed by the two-dimensional (2D) Maxwell-Boltzmann distribution [7]: 푝(푣) = ∙ exp − , with 푝 being the probability density, 푚 the particle mass, 퐾 = 푘 푇 the average particle kinetic energy (퐾), and 푘 the Boltzmann constant. If the mean free path of the particles
( ) is small compared to the plane height (푧̂), the system is ergodic and chaotic. At 푧 = 푧̂, the local particle density, , is equal to [8], where is the particle density ratio, is the local ̂ particle density at the bottom of the plane, =푒 is the Boltzmann factor, and 훽 = .
If ≫ 푧̂, the characteristic interaction time 푡 ̅ = ⁄푣̅ is much longer than the characteristic event duration 푡 ̅ =2푧̂⁄푣 ̅ , where 푣̅ is the average particle velocity and 푣 ̅ is the average 푧-dimension particle velocity component (푣 ). On the time scale of 푡 ̅ , few particle collisions take place and thus, the particle trajectories tend to be nonchaotic. The probability for a particle to travel from 푧 =0 to 푧 = 푧̂ is mainly determined by 푣 , relatively unrelated to the 푦- dimension momentum. As a first-order approximation, the particle density ratio at 푧̂ can be
2