Publications

Abstract not provided.

Wetting in a system where the kinetics of drop spreading are controlled by the rate of formation of a precursor film is modeled for the first time at the atomistic scale. Molecular dynamics simulations of Pb(l) wetting Cu(111) and Cu(100) show that a precursor film of atomic thickness evolves and spreads diffusively. This precursor film spreads significantly faster on Cu(111) than on Cu(100). For Cu(100), the kinetics of drop spreading are dramatically decreased by slow advancement of the precursor film. Slow precursor film kinetics on Cu(100) are partly due to the formation of a surface alloy at the solid-liquid interface which does not occur on Cu(111). © 2003 The American Physical Society.

Abstract not provided.

In this study, twelve embedded atom method (EAM) function sets were tested for their ability to predict liquid/vapor surface tension. Testing was carried out in the isochoric-isothermal (NVT) ensemble with a Berendsen thermostat. It was shown that the use of charge gradient corrections in conjunction with appropriate EAM functions provide surface property predictions in excellent agreement with experiment for solid and liquid metals.

The effect of polymer architecture on macroscopic properties were investigated using the self-consistent integral equation theory. Using several types of polyolefin polymers, the results obtained using the self consistent polymer reference interaction site model (PRISM) and molecular dynamics (MD) simulations were compared. The results from the two methods were then compared with experimental X ray scattering data.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In this paper the authors investigate the relationship between glassy and ferromagnetic phases in disordered Ising ferromagnets in the presence of transverse magnetic fields, {Lambda}. Iterative mean field simulations probe the free energy landscape and suggest the existence of a glass transition line in the {Lambda}, temperature T plane well within the ferromagnetic phase. New experimental field-cooled and zero-field-cooled data on LiHo{sub x} Y{sub 1{minus}x}F{sub 4} provide support for our theoretical picture.

The authors present the results of molecular dynamics simulations of n-butane and isobutane in silicalite. They begin with a comparison of the bulk adsorption and diffusion properties for two different parameterizations of the interaction potential between the hydrocarbon species, both of which have been shown to reproduce experimental gas-liquid coexistence curves. They examine diffusion as a function of the loading of the zeolite, as well as the temperature dependence of the diffusion constant at loading and for infinite dilution. They continue with simulations in which interfaces are formed between single component gases and the zeolite. After reaching equilibrium, they examine the dynamics of exchange between the bulk gas and the zeolite. Finally, they calculate the permeability of the zeolite for n-butane and isobutane as a function of pressure. Their simulations are performed for a number of different gas temperatures and pressures, covering a wide range of state points.

The authors report and analyze the results of numerical studies of dense granular flows in two and three dimensions, using both linear damped springs and Hertzian force laws between particles. Chute flow generically produces a constant density profile that satisfies scaling relations suggestive of a Bagnold grain inertia regime. The type for force law has little impact on the behavior of the system. Failure is not initiated at the surface, consistent with the absence of surface flows and different principal stress directions at vs. below the surface.

The non-linear stress-strain relation for crosslinked polymer networks is studied using molecular dynamics simulations. Previously we demonstrated the importance of trapped entanglements in determining the elastic and relaxational properties of networks. Here we present new results for the stress versus strain for both dry and swollen networks. Models which limit the fluctuations of the network strands like the tube model are shown to describe the stress for both elongation and compression. For swollen networks, the total modulus is found to decrease like (V{sub o}/V){sup 2/3} and goes to the phantom model result only for short strand networks.

We present the results of molecular dynamics simulations of very long model polymer chains analyzed by various experimentally relevant techniques. The segment motion of the chains is found to be in very good agreement with the reptation model. We also calculated the plateau modulus G⁰_N. The predictions of the entanglement length N_e from G⁰_N and from the mean square displacement of the chain segments disagree by a factor of about 2.2(2), indicating an error in the prefactor in the standard formula for G⁰_N. We show that recent neutron spin echo measurements were carried out for chain lengths which are too small to allow for a correct determination of N_e.

The effect of excluded volume on the coil size of dilute linear polymers was investigated by off-lattice Monte Carlo simulations. The radius of gyration R{sub g} was evaluated for a wide range of chain lengths at several temperatures and at the athermal condition. The theta temperature and the corresponding theta chain dimensions were established for the system, and the dependence of the size expansion factor, a{sub s} = R{sub g} /(R{sub g}){sub {theta}}, on chain length N and temperature T was examined. For long chains and at high temperatures, a{sub s} is a function of N/N{sub s}{sup 2} alone, where the length scale N{sub s}{sup 2} depends only on T. The form of this simulations-based master function compares favorably with {alpha}{sub s}(M/M{sub s}{sup 2}), an experimental master curve for linear polymers in good solvents, where M{sub s}{sup 2} depends only on polymer-solvent system. Comparisons when N{sub s}{sup 2}(T) and M{sub s}{sup 2}(system) are reduced to common units, numbers of Kuhn steps, strongly indicate that coil expansion in even the best of good solvents is small relative to that expected for truly athermal solutions. An explanation for this behavior is proposed, based on what would appear to be an inherent difference in the equation of state properties for polymeric and monomeric liquids.

Evidence for capillary waves at a liquid-vapor interface are presented from extensive molecular dynamics simulations of a system containing up to 1.24 million Lennard-Jones particles. Careful measurements show that the total interfacial width depends logarithmically on L∥, the length of the simulation cell parallel to the interface, as predicted theoretically. The strength of the divergence of the interfacial width on L∥ depends inversely on the surface tension γ. This allows us to measure γ two ways since γ can also be obtained from the difference in the pressure parallel and perpendicular to the interface. These two independent measures of γ agree provided that the interfacial order parameter profile is fit to an error function and not a hyperbolic tangent, as often assumed. We explore why these two common fitting functions give different results for γ.