The effect of self-affine roughness on solid contact is examined with molecular dynamics and continuum calculations. The contact area and normal and lateral stiffnesses rise linearly with the applied load, and the load rises exponentially with decreasing separation between surfaces. Results for a wide range of roughnesses, system sizes and Poisson ratios can be collapsed using Persson's contact theory for continuous elastic media. The atomic scale response at the interface between solids has little affect on the area or normal stiffness, but can greatly reduce the lateral stiffness. The scaling of this effect with system size and roughness is discussed. PACS numbers: 46.55.+d, 62.20.Qp, 81.40.Pq The presence of roughness on a wide range of length scales has profound effects on contact and friction between experimental surfaces. Under a broad range of conditions [1][2][3][4][5][6], the area of intimate contact between rough surfaces A c is orders of magnitude smaller than the apparent surface area A 0 . As discussed below, this provides the most common explanation for Amontons' laws that friction is proportional to load and independent of A 0 . Because A c is small, the interfacial region is very compliant. In a range of applications the interfacial compliance can significantly reduce the stiffness of macroscopic joints formed by holding two components together under pressure [1,7].In this paper, we examine the effect of surface roughness on the normal and lateral stiffness of contacts between elastic solids using molecular dynamics (MD) and continuum calculations. The results provide a numerical test of recent continuum theories [8,9] and their applicability to real solids. The contact area and normal stiffness approach continuum predictions rapidly as system size increases. Continuum theory also captures the internal deformations in the solid under tangential load, but the total lateral stiffness may be greatly reduced by atomic scale displacements between contacting atoms on the opposing surfaces. This makes it a sensitive probe of the forces underlying friction and may help to explain unexpectedly small experimental results [10].The topography of many surfaces can be described as a self-affine fractal [2,11]. Over a wide range of lengths, the root mean squared (rms) change in height dh over a lateral distance ℓ scales as a power law: dh ∼ ℓ H , where the roughness or Hurst exponent H is typically between 0.5 and 0.9. Greenwood and Williamson (GW) considered the peaks of rough landscapes as independent asperities and found that A c rose linearly with normal load F N for nonadhesive surfaces [2]. This explains Amontons's laws if there is a constant shear stress at the interface. A linear scaling of area with load is also obtained from Persson's scaling theory, which includes elastic coupling between contacts approximately [12,13].Dimensional analysis implies that the linear relation between load and area must have the form( 1) where a modulus like the contact modulus E ′ is the only dimensional quantity char...