Publications

Abstract not provided.

Abstract not provided.

Atomistic simulations were performed to investigate high temperature wetting phenomena for metals. A sessile drop configuration was modeled for two systems: Ag(l) on Cu and Pb(l) on Cu. The former case is an eutectic binary and the wetting kinetics were greatly enhanced by the presence of aggressive interdiffusion between Ag and Cu. Wetting kinetics were directly dependent upon dissolution kinetics. The dissolution rate was nearly identical for Ag(l) on Cu(100) compared to Cu(111); as such, the spreading rate was very similar on both surfaces. Pb and Cu are bulk immiscible so spreading of Pb(l) on Cu occurred in the absence of significant substrate dissolution. For Pb(l) on Cu(111) a precursor wetting film of atomic thickness emerged from the partially wetting liquid drop and rapidly covered the surface. For Pb(l) on Cu(100), a foot was also observed to emerge from a partially wetting drop; however, spreading kinetics were dramatically slower for Pb(l) on Cu(100) than on Cu(111). For the former, a surface alloying reaction was observed to occur as the liquid wet the surface. The alloying reaction was associated with dramatically decreased wetting kinetics on Cu(100) versus Cu(111), where no alloying was observed. These two cases demonstrate markedly different atomistic mechanisms of wetting where, for Ag(l) on Cu, the dissolution reaction is associated with increased wetting kinetics while, for Pb(l) on Cu, the surface alloying reaction is associated with decreased wetting kinetics.

The effects of side wall movement on granular packings were investigated. The studies showed that the resultant structure of the pack did not depend strongly on the magnitude of the wall movement, as long as the packing was moved for an equivalent distance. The main effect of wall movement was to drive the particle-wall and particle-particle contacts to the Coulomb criterion. This forced the packing in the high wall velocity case to obey the Janssen form, which took the Coulomb criterion as one of its main assumptions.

The spreading of polymer droplets is studied using molecular dynamics simulations. To study the dynamics of both the precursor foot and the bulk droplet, large hemispherical drops of 200 000 monomers are simulated using a bead-spring model for polymers of chain length 10, 20, and 40 monomers per chain. We compare spreading on flat and atomistic surfaces, chain length effects, and different applications of the Langevin and dissipative particle dynamics thermostats. We find diffusive behavior for the precursor foot and good agreement with the molecular kinetic model of droplet spreading using both flat and atomistic surfaces. Despite the large system size and long simulation time relative to previous simulations, we find that even larger systems are required to observe hydrodynamic behavior in the hemispherical spreading droplet.

Several methods for preparing well equilibrated melts of long chains polymers are studied. We show that the standard method in which one starts with an ensemble of chains with the correct end-to-end distance arranged randomly in the simulation cell and introduces the excluded volume rapidly, leads to deformation on short length scales. This deformation is strongest for long chains and relaxes only after the chains have moved their own size. Two methods are shown to overcome this local deformation of the chains. One method is to first pre-pack the Gaussian chains, which reduces the density fluctuations in the system, followed by a gradual introduction of the excluded volume. The second method is a double-bridging algorithm in which new bonds are formed across a pair of chains, creating two new chains each substantially different from the original. We demonstrate the effectiveness of these methods for a linear bead spring polymer model with both zero and nonzero bending stiffness, however the methods are applicable to more complex architectures such as branched and star polymer.

Abstract not provided.

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.

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.

The interface between liquid hexadecane and the (010) surface of silicalite was studied by molecular dynamics. The structure of molecules in the interracial region is influenced by the presence of pore mouths on the silicalite surface. For this surface, whose pores are the entrances to straight channels, the concentration profile for partially absorbed molecules is peaked around 10 monomers inside the zeolite. No preference to enter or exit the zeolite based on absorption length is observed except for very small or very large absorption lengths. We also found no preferential conformation of the unabsorbed tails for partially absorbed molecules.

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 γ.