Date of Award

12-2007

Degree Type

Thesis

Degree Name

Master of Science

Major

Chemical Engineering

Major Professor

Brian J. Edwards

Committee Members

David J. Keffer, William V. Steele

Abstract

Using a recently developed, methodical, Hamiltonian-based procedure, we derive rigorous algorithms for nonequilibrium molecular dynamic (NEMD) simulation in the NVE and NVT ensembles. We demonstrate that the equivalence of the kinetic temperature and configurational temperatures that exists at equilibrium is maintained in nonequilibrium states, given the proper nonequilibrium expression (in the microcanonical ensemble). Specifically, we apply the procedure to the p-SLLOD algorithm, which allows for rigorous NEMD simulation in the presence of an arbitrary, externally imposed flow field. The resulting algorithms are general in that they apply in the presence or absence of external fields in addition to the imposed flow field. Of particular note is the resulting algorithm for the canonical ensemble. The use of the Nosé-Hoover thermostat in the SLLOD and p-SLLOD algorithms has not been rigorously correct to date. We follow a methodical procedure to obtain a rigorous Hamiltonian-based NEMD algorithm using a reformulated Nosé-Hoover and Nosé-Poincaré thermostat. Although the resulting algorithms were unstable, it provided proof that it is the boundary conditions that drive the flow, contrary to the conventional belief that an external flow field has to be introduced into the equations of motion to simulate nonequilbrium flows.

Files over 3MB may be slow to open. For best results, right-click and select "save as..."

Share

COinS