No abstract

The phase diagram and dynamical properties of systems of particles interacting through a repulsive screened Coulomb (Yukawa) potential have been calculated using molecular and lattice dynamics techniques. The phase diagram contains both a melting transition and a transition from fcc to bcc crystalline phases. These phase transitions have been studied as a function of potential shape (screening length) and compared to phenomenological criteria for transition temperatures such as those of Lindemann and of Hansen and Verlet. The transition from fcc to bcc with increasing temperature is shown to result from a higher entropy in the bcc phase because of its softer shear modes. Even when the stable solid phase below the melting temperature is fcc, bcc-like local order is found in the liquid phase. This may substantially slow crystallization. The calculated phase diagram and shear modulus are in good agreement with experiments on colloidal suspensions of polystyrene spheres. The single particle dynamics of Yukawa systems show several unusual features. There is a pronounced subdiffusive regime in liquids near and below the melting temperature. This regime reflects the existence of two time scales: a typical phonon period, and the time for a particle to feel a new environment. The second time scale becomes longer as the temperature is lowered or the range of interaction (screening length) increases.

Forces acting within the area of atomic contact between surfaces play a central role in friction and adhesion. Such forces are traditionally calculated using continuum contact mechanics, which is known to break down as the contact radius approaches atomic dimensions. Yet contact mechanics is being applied at ever smaller lengths, driven by interest in shrinking devices to nanometre scales, creating nanostructured materials with optimized mechanical properties, and understanding the molecular origins of macroscopic friction and adhesion. Here we use molecular simulations to test the limits of contact mechanics under ideal conditions. Our findings indicate that atomic discreteness within the bulk of the solids does not have a significant effect, but that the atomic-scale surface roughness that is always produced by discrete atoms leads to dramatic deviations from continuum theory. Contact areas and stresses may be changed by a factor of two, whereas friction and lateral contact stiffness change by an order of magnitude. These variations are likely to affect continuum predictions for many macroscopic rough surfaces, where studies show that the total contact area is broken up into many separate regions with very small mean radius.

Finite element methods are used to study non-adhesive, frictionless contact between elastic solids with self-affine surfaces. We find that the total contact area rises linearly with load at small loads. The mean pressure in the contact regions is independent of load and proportional to the rms slope of the surface. The constant of proportionality is nearly independent of Poisson ratio and roughness exponent and lies between previous analytic predictions. The contact morphology is also analyzed. Connected contact regions have a fractal area and perimeter. The probability of finding a cluster of area ac drops as a −τ c where τ increases with decreasing roughness exponent. The distribution of pressures shows an exponential tail that is also found in many jammed systems. These results are contrasted to simpler models and experiment.