We present a universal relation for crack surface cohesion including surface relaxation. Specifically, we analyze how N atomic planes respond to an opening displacement at its boundary, producing structurally relaxed surfaces. Via density-functional theory, we verify universality for metals ͑Al͒, ceramics (␣-Al 2 O 3 ), and semiconductors ͑Si͒. When the energy and opening displacement are scaled appropriately with respect to N, the uniaxial elastic constant, the relaxed surface energy, and the equilibrium interlayer spacing, all energydisplacement curves collapse onto a single universal curve. DOI: 10.1103/PhysRevB.69.172104 PACS number͑s͒: 61.50.Lt, 68.35.Ϫp, 71.15.Mb Macroscopic cohesive theories of fracture often involve empirical postulates on the shape and form of the cohesive law.1-3 While first principles simulations might be preferable, the typical size of engineering finite element models prohibits their direct application. In order to obtain converged finite element results, the cohesive zone size must be resolved by the mesh; for brittle materials, the cohesive zone size is atomistic, making the calculation prohibitively expensive.4 Nanometer scale quantum mechanical calculations have provided insight into cracking at the atomic level, 5-8 but their extrapolation to the macroscopic scale is fraught with difficulty. Indeed, orders-of-magnitude mismatch exist between atomistic predictions of cohesive strengths and critical opening displacements 9-11 and measurements of tensile strength in brittle materials obtained from spallation tests, 12 the latter of which are often employed in engineering simulations. The widely used universal binding energy relation ͑UBER͒ of Rose et al.13 describes cohesion between rigid surfaces based on atomic scale calculations, but application of the UBER to crack propagation simulations is hampered by its inability to capture the shape and absolute energies of cohesive laws for structurally relaxed surfaces.6,14 -16 Here, we address these difficulties by deriving a coarse-grained cohesive energy relation that accounts for structural relaxation of surfaces and exhibits a material-independent universal form.Nguyen and Ortiz recently suggested rescaling interlayer potentials to yield macroscopic cohesive laws. 17 Here we extend their work to account for surface relaxation and reconstruction. Specifically, we consider a perfect crystal acted upon by tensile stresses normal to a cleavage plane. The length scales under consideration range from mesoscopic ͑the dislocation free zone of a metal 18 ͒ to possibly macroscopic ͑brittle materials͒; we conservatively denote both scales as mesoscopic. We assume that atomic layers remain planar after deformation, so that the relative displacement of crystallographic layer i can be described by ␦ i ͑the interlayer spacing minus the equilibrium interlayer spacing, d͒ and that the crystal is periodic with a unit cell containing N atomic layers. We express the total energy per unit area of cleavage plane asis the local energy per layer of a ...