Source author record

William Graham Hoover

William Graham Hoover appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

18 published item(s)

preprint2022arXiv

Computational Implementation of Maxwell's Knudsen-Gas Demon and Its Extension to a Two-Dimensional Soft-Disk Fluid

An interesting preprint by Puru Gujrati, "Maxwell's Demon Must Remain Subservient to Clausius' Statement" [ of the Second Law of Thermodynamics ], traces the development and application of Maxwell's Demon. He argues against the Demon on thermodynamic grounds. Gujrati introduces and uses his own version of a generalized thermodynamics in his criticism of the Demon. The complexity of his paper and the lack of any accompanying numerical work piqued our curiousity. The internet provides well over two million "hits'' on the subject of "Maxwell's Demon". There are also hundreds of images of the Demon, superimposed upon a container of gas or liquid. However, there is not so much along the lines of simulations of the Demonic process. Accordingly, we thought it useful to write and execute relatively simple FORTRAN programs designed to implement Maxwell's low-density model and to develop its replacement with global Nosé-Hoover or local purely-Newtonian thermal controls. These simulations illustrate the entropy decreases associated with all three types of Demons.

preprint2022arXiv

Thermodynamic Entropy from Sadi Carnot's Cycle using Gauss' and Doll's-Tensor Molecular Dynamics

Carnot's four-part ideal-gas cycle includes both isothermal and adiabatic expansions and compressions. Analyzing this cycle provides the fundamental basis for statistical thermodynamics. We explore the cycle here from a pedagogical view in order to promote understanding of the macroscopic thermodynamic entropy, the state function associated with thermal energy changes. From the alternative microscopic viewpoint the Hamiltonian ${\cal H}(q,p)$ is the energy and entropy is the (logarithm of the) phase-space volume $Ω$ associated with a macroscopic state. We apply two novel forms of Hamiltonian mechanics to Carnot's Cycle: [1] Gauss' isokinetic mechanics for the isothermal segments and [2] Doll's Tensor for the isentropic adiabatic segments. We explore the equivalence of the microscopic and macroscopic views of Carnot's cycle for simple fluids here, beginning with the ideal Knudsen gas and extending the analysis to a prototypical simple fluid.

preprint2021arXiv

Nonequilibrium Time Reversibility with Maps and Walks

Time-reversible dynamical simulations of nonequilibrium systems exemplify both Loschmidt's and Zermélo's paradoxes. That is, computational time-reversible simulations invariably produce solutions consistent with the {\it irreversible} Second Law of Thermodynamics (Loschmidt's) as well as {\it periodic} in the time (Zermélo's, illustrating Poincaré recurrence). Understanding these paradoxical aspects of time-reversible systems is enhanced here by studying the simplest pair of such model systems. The first is time-reversible, but nevertheless dissipative and periodic, the piecewise-linear compressible Baker Map. The fractal properties of that two-dimensional map are mirrored by an even simpler example, the one-dimensional random walk, confined to the unit interval. As a further puzzle the two models yield ambiguities in determining the fractals' information dimensions. These puzzles, including the classical paradoxes, are reviewed and explored here. We review our investigations presented in Budapest in 1997 and end with presentday questions posed as the Snook Prize Problems in 2020 and 2021.

preprint2020arXiv

Compressible Baker Maps and Their Inverses. A Memoir for Francis Hayin Ree [ 1936-2020 ]

This memoir is dedicated to the late Francis Hayin Ree, a formative influence shaping my work in statistical mechanics. Between 1963 and 1968 we collaborated on nine papers published in the Journal of Chemical Physics. Those dealt with the virial series, cell models, and computer simulation. All of them were directed toward understanding the statistical thermodynamics of simple model systems. Our last joint work is also the most cited, with over 1000 citations, "Melting Transition and Communal Entropy for Hard Spheres", submitted 3 May 1968 and published that October. Here I summarize my own most recent work on compressible time-reversible two-dimensional maps. These simplest of model systems are amenable to computer simulation and are providing stimulating and surprising results.

preprint2016arXiv

A Tutorial: Adaptive Runge-Kutta Integration for Stiff Systems : Comparing the Nosé and Nosé-Hoover Oscillator Dynamics

"Stiff" differential equations are commonplace in engineering and dynamical systems. To solve them we need flexible integrators that can deal with rapidly-changing righthand sides. This tutorial describes the application of "adaptive" [ variable timestep ] integrators to "stiff" mechanical problems encountered in modern applications of Gibbs' 1902 statistical mechanics. Linear harmonic oscillators subject to nonlinear thermal constraints can exhibit either stiff or smooth dynamics. Two closely-related examples, Nosé's 1984 dynamics and Nosé-Hoover 1985 dynamics, are both based on Hamiltonian mechanics, as was ultimately clarified by Dettmann and Morriss in 1996. Both these dynamics are consistent with Gibbs' canonical ensemble. Nosé's dynamics is "stiff" and can present severe numerical difficulties. Nosé-Hoover dynamics, though it follows exactly the same trajectory, is "smooth" and relatively trouble-free. Our tutorial emphasises the power of adaptive integrators to resolve stiff problems like the Nosé oscillator. The solutions obtained illustrate the power of computer graphics to enrich numerical solutions. Adaptive integration with computer graphics are basic to an understanding of dynamical systems and statistical mechanics. These tools lead naturally into the visualization of intricate fractal structures formed by chaos as well as elaborate knots tied by regular nonchaotic dynamics. This work was invited by the American Journal of Physics.

preprint2016arXiv

From Ann Arbor to Sheffield: Around the World in 80 Years. I

Childhood and graduate school at Ann Arbor Michigan prepared Bill for an interesting and rewarding career in physics. Along the way came Carol and many joint discoveries with our many colleagues to whom we both owe this good life. This summary of Bill's early work prior to their marriage and sabbatical in Japan is Part I, prepared for Bill's 80th Birthday celebration at the University of Sheffield in July 2016.

preprint2016arXiv

Nonequilibrium Systems : Hard Disks and Harmonic Oscillators Near and Far From Equilibrium

We relate progress in statistical mechanics, both at and far from equilibrium, to advances in the theory of dynamical systems. We consider computer simulations of time-reversible deterministic chaos in small systems with three- and four-dimensional phase spaces. These models provide us with a basis for understanding equilibration and thermodynamic irreversibility in terms of Lyapunov instability, fractal distributions, and thermal constraints

preprint2016arXiv

Singly-Thermostated Ergodicity in Gibbs' Canonical Ensemble and the 2016 Ian Snook Prize

For a harmonic oscillator, Nosé's single-thermostat approach to simulating Gibbs' canonical ensemble with dynamics samples only a small fraction of the phase space. Nosé's approach has been improved in a series of three steps: [ 1 ] several two-thermostat sets of motion equations have been found which cover the complete phase space in an ergodic fashion, [ 2 ] sets of single-thermostat motion equations, exerting "weak control" over both forces and momenta, have been shown to be ergodic, and [ 3 ] sets of single-thermostat motion equations exerting weak control over two velocity moments provide ergodic phase-space sampling for the oscillator and for the rigid pendulum, but not for the quartic oscillator or for the Mexican Hat potential. The missing fourth step, motion equations providing ergodic sampling for anharmonic potentials requires a further advance. The 2016 Ian Snook Prize will be awarded to the author(s) of the most interesting original submission addressing the problem of finding ergodic algorithms for Gibbs' canonical ensemble using a single thermostat.

preprint2016arXiv

Yokohama to Ruby Valley : Around the World in 80 Years. II

We two had year-long research leaves in Japan, working together fulltime with several Japanese plus Tony De Groot back in Livermore and Harald Posch in Vienna. We summarize a few of the high spots from that very productive year ( 1989-1990 ), followed by an additional fifteen years' work in Livermore, with extensive travel. Next came our retirement in Nevada in 2005, which has turned out to be a long-term working vacation. Carol narrates this part of our history together.

preprint2015arXiv

Comparison of Very Smooth Cell-Model Trajectories Using Five Symplectic and Two Runge-Kutta Integrators

Time-reversible symplectic methods, which are precisely compatible with Liouville's phase-volume-conservation theorem, are often recommended for computational simulations of Hamiltonian mechanics. Lack of energy drift is an apparent advantage of such methods. But all numerical methods are susceptible to Lyapunov instability, which severely limits the maximum time for which chaotic solutions can be "accurate". The "advantages" of higher-order methods are lost rapidly for typical chaotic Hamiltonians. We illustrate these difficulties for a useful reproducible test case, the two-dimensional one-particle cell model with specially smooth forces. This Hamiltonian problem is chaotic and occurs on a three-dimensional constant-energy shell, the minimum dimension for chaos. We benchmark the problem with quadruple-precision trajectories using the fourth-order Candy-Rozmus, fifth-order Runge-Kutta, and eighth-order Schlier-Seiter-Teloy integrators. We compare the last, most-accurate particle trajectories to those from six double-precision algorithms, four symplectic and two Runge-Kutta.

preprint2015arXiv

Deterministic Time-Reversible Thermostats : Chaos, Ergodicity, and the Zeroth Law of Thermodynamics

The relative stability and ergodicity of deterministic time-reversible thermostats, both singly and in coupled pairs, are assessed through their Lyapunov spectra. Five types of thermostat are coupled to one another through a single Hooke's-Law harmonic spring. The resulting dynamics shows that three specific thermostat types, Hoover-Holian, Ju-Bulgac, and Martyna-Klein-Tuckerman, have very similar Lyapunov spectra in their equilibrium four-dimensional phase spaces and when coupled in equilibrium or nonequilibrium pairs. All three of these oscillator-based thermostats are shown to be ergodic, with smooth analytic Gaussian distributions in their extended phase spaces ( coordinate, momentum, and two control variables ). Evidently these three ergodic and time-reversible thermostat types are particularly useful as statistical-mechanical thermometers and thermostats. Each of them generates Gibbs' universal canonical distribution internally as well as for systems to which they are coupled. Thus they obey the Zeroth Law of Thermodynamics, as a good heat bath should. They also provide dissipative heat flow with relatively small nonlinearity when two or more such bath temperatures interact and provide useful deterministic replacements for the stochastic Langevin equation.

preprint2015arXiv

Ergodic Time-Reversible Chaos for Gibbs' Canonical Oscillator

Nosé's pioneering 1984 work inspired a variety of time-reversible deterministic thermostats. Though several groups have developed successful doubly-thermostated models, single-thermostat models have failed to generate Gibbs' canonical distribution for the one-dimensional harmonic oscillator. Sergi and Ferrario's 2001 doubly-thermostated model, claimed to be ergodic, has a singly-thermostated version. Though neither of these models is ergodic this work has suggested a successful route toward singly-thermostated ergodicity. We illustrate both ergodicity and its lack for these models using phase-space cross sections and Lyapunov instability as diagnostic tools.

preprint2015arXiv

Ergodicity of the Martyna-Klein-Tuckerman Thermostat and the 2014 Snook Prize

Nosé and Hoover's 1984 work showed that although Nosé and Nosé-Hoover dynamics were both consistent with Gibbs' canonical distribution neither dynamics, when applied to the harmonic oscillator, provided Gibbs' Gaussian distribution. Further investigations indicated that two independent thermostat variables are necessary, and often sufficient, to generate Gibbs' canonical distribution for an oscillator. Three successful time-reversible and deterministic sets of two-thermostat motion equations were developed in the 1990s. We analyze one of them here. It was developed by Martyna, Klein, and Tuckerman in 1992. Its ergodicity was called into question by Patra and Bhattacharya in 2014. This question became the subject of the 2014 Snook Prize. Here we summarize the previous work on this problem and elucidate new details of the chaotic dynamics in the neighborhood of the two fixed points. We apply six separate tests for ergodicity and conclude that the MKT equations are fully compatible with all of them, in consonance with our recent work with Clint Sprott and Puneet Patra.

preprint2015arXiv

The Equivalence of Dissipation from Gibbs' Entropy Production with Phase-Volume Loss in Ergodic Heat-Conducting Oscillators

Gibbs' thermodynamic entropy is given by the logarithm of the phase volume, which itself responds to heat transfer to and from thermal reservoirs. We compare the thermodynamic dissipation described by phase-volume loss with heat-transfer entropy production. Their equivalence is documented for computer simulations of the response of an ergodic harmonic oscillator to thermostated temperature gradients. In the simulations one or two thermostat variables control the kinetic energy or the kinetic energy and its fluctuation. All of the motion equations are time-reversible. We consider both strong and weak control variables. In every case the time-averaged dissipative loss of phase-space volume coincides with the entropy produced by heat transfer. Linear response theory nicely reproduces the small-gradient results obtained by computer simulation. The thermostats considered here are ergodic and provide simple dynamical models, some of them with as few as three ordinary differential equations, while remaining capable of reproducing Gibbs' canonical phase-space distribution and precisely consistent with irreversible thermodynamics.

preprint2015arXiv

Time-Reversible Ergodic Maps and the 2015 Ian Snook Prize

The time reversibility characteristic of Hamiltonian mechanics has long been extended to nonHamiltonian dynamical systems modeling nonequilibrium steady states with feedback-based thermostats and ergostats. Typical solutions are multifractal attractor-repellor phase-space pairs with reversed momenta and unchanged coordinates, $(q,p)\longleftrightarrow (q,-p)$ . Weak control of the temperature, $\propto p^2$ and its fluctuation, resulting in ergodicity, has recently been achieved in a three-dimensional time-reversible model of a heat-conducting harmonic oscillator. Two-dimensional cross sections of such nonequilibrium flows can be generated with time-reversible dissipative maps yielding æsthetically interesting attractor-repellor pairs. We challenge the reader to find and explore such time-reversible dissipative maps. This challenge is the 2015 Snook-Prize Problem.

preprint2014arXiv

Ergodicity of a Time-Reversibly Thermostated Harmonic Oscillator and the 2014 Ian Snook Prize

Shuichi Nosé opened up a new world of atomistic simulation in 1984. He formulated a Hamiltonian tailored to generate Gibbs' canonical distribution dynamically. This clever idea bridged the gap between microcanonical molecular dynamics and canonical statistical mechanics. Until then the canonical distribution was explored with Monte Carlo sampling. Nosé's dynamical Hamiltonian bridge requires the "ergodic" support of a space-filling structure in order to reproduce the entire distribution. For sufficiently small systems, such as the harmonic oscillator, Nosé's dynamical approach failed to agree with Gibbs' sampling and instead showed a complex structure, partitioned into a chaotic sea, islands, and chains of islands, that is familiar textbook fare from investigations of Hamiltonian chaos. In trying to enhance small-system ergodicity several more complicated "thermostated" equations of motion were developed. All were consistent with the canonical Gaussian distribution for the oscillator coordinate and momentum. The ergodicity of the various approaches has undergone several investigations, with somewhat inconclusive ( contradictory ) results. Here we illustrate several ways to test ergodicity and challenge the reader to find even more convincing algorithms or an entirely new approach to this problem.

preprint2014arXiv

Heat Conduction, and the Lack Thereof, in Time-Reversible Dynamical Systems: Generalized Nosé-Hoover Oscillators with a Temperature Gradient

We use nonequilibrium molecular dynamics to analyze and illustrate the qualitative differences between the one-thermostat and two-thermostat versions of equilibrium and nonequilibrium (heat-conducting) harmonic oscillators. Conservative nonconducting regions can coexist with dissipative heat conducting regions in phase space with exactly the same imposed temperature field.