We discuss the origin of the breakloose (or static) friction force when an ice block is slid on a hard randomly rough substrate surface. If the substrate has roughness with small enough amplitude (of order a 1 nm or less), the breakloose force may be due to interfacial slip and is determined by the elastic energy per unit area, Uel/A0, stored at the interface after the block has been displaced a short distance from its original position. The theory assumes complete contact between the solids at the interface and that there is no elastic deformation energy at the interface in the original state before the application of the tangential force. The breakloose force depends on the surface roughness power spectrum of the substrate and is found to be in good agreement with experimental observations. We show that as the temperature decreases, there is a transition from interfacial sliding (mode II crack propagation, where the crack propagation energy GII = Uel/A0) to opening crack propagation (mode I crack propagation with GI the energy per unit area to break the ice–substrate bonds in the normal direction).

Friction on ice is a fascinating topic with a long history.1–22 Pioneering work was done by Faraday more than 150 years ago.23 He proposed that the surface of ice is covered by a liquid-like water layer, now understood as a result of premelting, which is a very common phenomenon for crystalline solids.24–33 Thomson suggested instead that the liquid-like layer was due to pressure melting, as observed for materials that at the melting temperature have a bigger volume in the solid state than in the fluid state.34 Pressure melting was often used as an explanation for the low friction on ice for the following 80 years or so, until Bowden and Hughes suggested that frictional heating might result in melting of the ice surface and that this was the main reason for the low friction of ice at high enough sliding speeds.35 Yet, frictional heating cannot explain why the surface of ice remains slippery when standing still or at low sliding speeds. These observations were resolved by reconsidering Faraday’s original premolten layer proposal and confirming its existence using a battery of modern experimental techniques.1,30,31 Why does ice exhibit a nonzero friction force as the sliding speed v → 0 if a liquid-like premelted water film separates the surfaces at the sliding interface? For an elastically hard solid with large surface roughness, a contribution to the friction force on ice comes from plowing of the ice by the hard asperities of the counter (substrate) surface. However, for a perfectly smooth substrate surface, one would expect a much smaller friction coefficient than observed when the solids are separated by a premelted water film with water-like viscosity. Shearing of a water film with the thickness d = 1 nm results in a viscous frictional shear stress τ = ηv/d ≈ 1 kPa when the sliding speed is v = 1 mm s−1. If the ice has large enough surface roughness, the contact pressure in the ice contact regions will be of the same order as the ice penetration hardness, which is typically σP ≈ 10 MPa, giving the friction coefficient μ = ηv/(σPd) ≈ 10−4, which is much smaller than usually observed for the sliding speed v = 1 mm s−1. In reality, the premelted water film may have larger effective viscosity than water and may exhibit shear thinning,50 but one would still expect μ → 0 as v → 0, which is not observed.

Experiments have shown that there is a premelted water film at the interface between (hydrophilic) silica glass and ice.36–38 When a tangential force displaces an ice block frozen on a silica surface, it is usually found that the slip plane is located not at the ice–silica interface but in the ice (cohesive failure). However, for very smooth silica surfaces, the slip occurs at the ice–silica interface (adhesive failure). Thus, for a silica surface with the root-mean-square (rms) roughness hrms ≈ 7.5 nm, only cohesive failure was observed, while on a silica surface with hrms ≈ 0.45 nm (when measured over a 2 × 2 μm2 surface area), interfacial sliding occurred for temperatures between 0 and −20 °C and cohesive failure for T < −25 °C (see Ref. 36). For −25 °C < T < −20 °C, adhesive failure (but no sliding) occurred, i.e., no ice could be detected remaining on the silica surface. When interfacial sliding takes place, the maximum (or breakloose) shear stress τmax increases linearly with the difference between the ice bulk melting temperature and the actual temperature, TmT, with τmax ≈ 0.75 MPa at T = −20 °C. This result was tentatively explained by the assumption of a premelted liquid-like water film at the ice–silica interface for T > −20 °C, as verified in Ref. 36 using several other experimental methods.

Here, we will discuss the origin of the breakloose force for ice in contact with a hard randomly rough substrate surface. We will assume that the substrate surface is so smooth that (because of adhesion) the ice (with the premelted water film) is in complete contact with it at any stage in the sliding act. We also assume that in the initial state, there are no stress and elastic deformations in the ice induced by the roughness of the substrate surface. This is the case if the ice is formed by freezing water on the substrate surface or if the ice is pressed against the substrate and left for sufficiently long time that the stress field initially generated in the vicinity of the hard substrate asperities can relax by creep or regelation to a nearly stress free state. We assume that a thin liquid-like film occurs at the ice–silica interface (due to ice premelting) and that the (viscous) frictional shear stress can be neglected. We show that as the temperature decreases, there is a transition from interfacial sliding (mode II crack propagation, where the crack propagation energy GII = Uel/A0) to opening crack propagation (mode I crack propagation with GI the energy per unit area to break the ice–substrate bonds in the normal direction).

Consider applying a tangential force to an ice block in full contact with a hard randomly rough substrate. We assume first that the interface displaces uniformly laterally in response to the applied tangential force, but later we will show that the breakloose friction involves frictional (type II) crack propagation. We assume no pressure melting or plastic deformation of the ice and that no change occur in the ice surface topography due to flow of the premelted water film from high pressure to low pressure interfacial regions, where it could freeze (regelation resulting from premelting, see Fig. 5).

Under the assumptions above, it is easy to derive an expression for the breakloose force. Let h(x) be the surface roughness profile of the substrate. When the bottom surface of the ice block has moved the lateral distance s=sx̂, the ice surface deforms elastically in response to the change in the surface roughness profile it experiences. The normal displacement field is
We define
so that
(1)
The elastic energy stored in the resulting deformation field is39,
(2)
where ⟨⋯⟩ stands for ensemble averaging and where E* is the effective Young’s modulus, which for a rigid substrate is given by E* = E/(1 − ν2), where E and ν are Young’s modulus and Poisson's ratio of the ice. Using (1) and (2), we get
Since the surface roughness power spectrum is39,40
(3)
where A0 is the surface area, we get
Using the integral representation of the Bessel function of order 0,
and introducing polar coordinates in q-space with d2q = dqqdϕ, we get
(4)
where we have assumed isotropic surface roughness so the power spectrum depends on only |q| = q and not the angle ϕ. The force F = dUel/ds is given by
(5)
where J1(x)=J0(x) is the Bessel function of order 1. Using that for small x, we have J1(x) ≈ x/2, and we get Fks as s → 0, where
Note that the mean square slope and the mean square curvature of the surface are determined by41,42
which both are usually determined mainly by the shortest wavelength surface roughness. This is, therefore, also the case for k and for the maximum Fmax of the F(s) curve.
It is clear from (5) that scaling the power spectrum with a factor g2 will increase the breakloose force with the same factor g2. Scaling C(q) by g2 is equivalent with scaling h(x) with a factor of g. Thus, the breakloose force will scale with the rms surface roughness hrms as hrms2. For a self-affine fractal surface, it is also clear that increasing the fractal dimension Df (which is equivalent to reducing the Hurst exponent H) but keeping the rms roughness amplitude and the roll-off wavenumber unchanged will increase the breakloose force as the short wavelength roughness is increased by this transformation. Similarly keeping hrms and H constant but increasing the roll-off wavenumber of the power spectrum will increase the amplitude of the short wavelength roughness and, hence, the breakloose force. For self-affine fractal surfaces with the power spectrum
for q0 < q < q1 and zero otherwise, the asymptotic value for Uel/A0 for large displacement is given by (see  Appendix A)
To study the dependency of F and Uel on the displacement s quantitatively, we have performed numerical calculations using a realistic surface roughness power spectrum.

We have measured the surface topography of flat fused silica widows obtained from Thorlabs (super polished mirror blank 2 Å), which is similar to the samples studied in Ref. 36 (which were obtained from Laser Components GmbH, Olching, Germany). The topography measurements were obtained using a MultiMode8 atomic force microscope (AFM, Bruker) operated in ScanAsyst imaging mode. The images were recorded with ScanAsyst-air cantilevers with a tip radius of 2 nm and a spring constant of 0.4 N/m. From the measured height profile z = h(x, y), we have calculated the surface roughness power spectra using (3).

Figure 1 shows the surface roughness power spectrum of the silica surface as obtained from AFM topography measurements (1024 × 1024 data points) over square areas with the side 10 µm (red curve) and 1 µm (green and blue curves). The rms-roughness values are 0.59, 0.48, and 0.46 nm and the rms-slopes are 0.05, 0.20 and 0.19 for the red, green, and blue curves, respectively. The slope of the dashed line correspond to the Hurst exponent H ≈ 0.21. Such small Hurst exponents are observed only on very smooth surfaces (small hrms) for reasons related to surface fragility (see Ref. 43). In  Appendix B, we show the height probability distribution Ph and a topography image. The height distribution Ph is nearly Gaussian as expected for a randomly rough surface, but the theory presented above does not require the surface to be randomly rough.

FIG. 1.

The surface roughness power spectrum of the silica substrate as obtained from AFM topography measurements (1024 × 1024 data points) over square areas with the side 10 µm (red curve) and 1 µm (green and blue curves). The rms-roughness values are 0.59, 0.48, and 0.46 nm and the rms-slopes are 0.05, 0.20, and 0.19, for the red, green, and blue curves, respectively. The slope of the dashed line correspond to the Hurst exponent H ≈ 0.21.

FIG. 1.

The surface roughness power spectrum of the silica substrate as obtained from AFM topography measurements (1024 × 1024 data points) over square areas with the side 10 µm (red curve) and 1 µm (green and blue curves). The rms-roughness values are 0.59, 0.48, and 0.46 nm and the rms-slopes are 0.05, 0.20, and 0.19, for the red, green, and blue curves, respectively. The slope of the dashed line correspond to the Hurst exponent H ≈ 0.21.

Close modal
For numerical calculations, it is convenient to write qs = x so that (5) takes the form
We introduce the shear stress τ = F/A0 and get
(6)
Since J1(x) oscillates as a function of x with a period of order 2π, it is necessary to use a uniform distribution of x-integration steps with step length Δx ≈ 0.1 (or less) in order to evaluate the integral (6) numerically.

Figure 2 shows the shear stress τ divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. We have used the surface roughness power spectrum shown in Fig. 1. The red line is using the combined power spectra from the AFM data obtained on 1 and 10 µm scan lengths (green + red line in Fig. 1), while the green line is using only the 1 µm scan length data. Clearly, the shear stress is dominated by the large wavenumber part of C(q), and the two calculations give nearly identical results. For surfaces with smaller fractal dimension (larger H), the long wavelength roughness becomes more important and the maximal stress occurs at larger displacements. Thus, in these cases, Uel will reach its asymptotic value for larger displacements, which implies that the crack tip process zone will be longer (see Sec. V).

FIG. 2.

The shear stress τ divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. The green and red lines are calculated for the surface roughness power spectrum given by the green line and the green + red line in Fig. 1.

FIG. 2.

The shear stress τ divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. The green and red lines are calculated for the surface roughness power spectrum given by the green line and the green + red line in Fig. 1.

Close modal

Another important quantity is the elastic energy stored at the interface44 as a function of the displacement s. Figure 3 shows the elastic energy per unit surface area divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. Note again that using the power spectrum of the 1 µm scan length AFM data gives nearly the same result as the combined power spectrum from the 1 and 10 µm data.

FIG. 3.

The elastic energy per unit surface area divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. The green and red lines are calculated for the surface roughness power spectrum given by the green line and the green + red line in Fig. 1.

FIG. 3.

The elastic energy per unit surface area divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. The green and red lines are calculated for the surface roughness power spectrum given by the green line and the green + red line in Fig. 1.

Close modal

If the shear stress τ given by (6) is very small (see below), then the bottom surface of the ice will displace uniformly due to the applied tangential force. However, this is not the case for the systems studied in Ref. 36 where instead an interfacial mode II (frictional) crack will propagate at the interface during onset of sliding (see Fig. 4).45 Thus, the ice–substrate contact region will consist of two regions: one where slip has occurred and one where no slip has occurred. At the boundary line (crack tip) between these two regions occur a narrow strip (crack tip process zone) of width of order a few nm in the present case (but larger for surfaces with smaller fractal dimension) where the contact change from the commensurate-like no-slip state to the incommensurate-like slip state. The static or breakloose friction force is the determined by the force needed to move the crack tip, which depends on the crack propagation energy G = Uel/A0.

FIG. 4.

A rectangular ice block in frozen contact with a silica substrate. On one side of the ice block, a stress or pressure σ is applied. The contact between the ice block and the substrate is broken by interfacial (type II) crack propagation. We assume that the ice–silica contact in the broken area x < x0 exhibits zero friction (zero frictional shear stress).

FIG. 4.

A rectangular ice block in frozen contact with a silica substrate. On one side of the ice block, a stress or pressure σ is applied. The contact between the ice block and the substrate is broken by interfacial (type II) crack propagation. We assume that the ice–silica contact in the broken area x < x0 exhibits zero friction (zero frictional shear stress).

Close modal
The mode II crack propagation problem has been studied in many publications. If we assume that the kinetic frictional shear stress τk vanishes in the slip region (because of an ice premelted water film), then (see  Appendix C and Ref. 46),
where α is a dimensionless geometrical factor, which is of order unity if all length scales (height and length of the block, and height of the region over which the pressure σ is applied) are of similar magnitude. Thus, the effective frictional shear stress τ = F/A0, where F = σhw and A0 = wL, is given by
(7)
In the present case, we have a cylinder shaped ice block with the height h = 9 mm and the diameter L ≈ 10 mm so that h/L ≈ 1; we expect α to be of order unity.
If the shear stress τmax [given by (6)] is smaller than τs, then the slip at the interface at the onset of sliding will be (nearly) uniform, but if τs < τmax, the onset of sliding will involve the interfacial crack propagation. If we write G = τmaxD, where D is a distance of order the length of the crack-tip process zone, then in the present case, from Figs. 2 and 3, D ≈ 2.3 nm. The condition τs = τmax gives
or (with α = 1)
(8)
This equation is equivalent to σ = 2ED/L. The stress σ acting on a block of length L gives a displacement u given by σ = Eu/L or u = /E = 2D. If the displacement u is smaller than the crack tip process zone (which is of order D), then the breakloose friction cannot be described by the crack propagation theory but is determined by (6). If hL ≈ 1 cm and D ≈ 5 nm, then (8) gives for ice τmax ≈ 104 Pa. In the present case, τmax ≈ 108 Pa and interfacial crack propagation will occur.

Using ice Young’s modulus E = 9 GPa and Poisson ratio ν = 0.3, we get the effective modulus of ice E* = E/(1 − ν2) ≈ 1010 Pa and (from Fig. 3) Uel/(E*A0) = 2.1 × 10−11 m−1, and we get G = 0.21 J/m2. Using (7) with L = 1 cm and α = 1, we get τs = 0.55 MPa, which is close to the measured breakloose shear stress 0.75 MPa.

In Ref. 36, the maximum (breakloose) shear stress τs = Fmax/A0 was measured for ice blocks in frozen contact with a very smooth silica surface with hrms ≈ 0.45 nm. Interfacial slip was observed for −20 °C < T < 0 °C and the shear stress decreased linearly with TmT as the temperature approached the melting temperature Tm = 0 °C. For T < −25 °C, cohesive fracture was observed. The breakloose shear stress for T ≈ −20 °C was found to be τs ≈ 0.75 MPa, which we interpret as the shear stress when the premelted film is first formed and which may be of monolayer thickness. This shear stress can be compared to the theory prediction τs ≈ 0.55 MPa, assuming the surface roughness power spectrum shown in Fig. 1.

We note that the mode II (sliding) crack propagation energy found above, GII = G ≈ 0.2 J/m2, is similar to the expected crack propagation energy for mode I (opening) crack propagation caused by the adhesive ice–countersurface interaction. For ice–ice adhesion, the (adiabatic) opening crack propagation energy GI = 2γice is twice the surface energy per unit surface area of ice. The latter is47  0.1J/m2, so for the ice–ice interface, GI ≈ 0.2 J/m2. We expect a similar interfacial binding energy for the ice–silica interface because the silica surface may be covered by a ∼monolayer thick water film. The energy to propagate an opening crack may be larger than given above because of nonadiabatic effects. In addition, the normal stresses acting on the ice at the silica–ice interface are so high that some plastic deformations may occur in the ice at the sliding interface (see  Appendix D).

Why does the breakloose shear stress decrease as the temperature T approaches the bulk melting temperature Tm? Here we propose the following explanation. Figure 5(a) shows ice in full adhesive contact with a hard substrate with sinus-like corrugation. If the ice is displaced laterally with a distance s equal to half the wavelength of the roughness profile, a positive pressure will result at the top of asperities and negative pressure in the valleys [see (b)]. As a result, water will flow in the premelted film from top of asperities to the bottom of the valleys. Since the premelted water film thickness must be approximately constant in time, ice must melt at the top of asperities and refreeze at the bottom of the valleys. Thus, ice is removed from the top of asperities and is added at the bottom of the wells. After some relaxation time t*, the stress free state (c) is obtained. The relaxation time t* depends on the temperature T and is shorter the closer the T is to the ice bulk melting temperature Tm because of the increased thickness of the premelted film and because of the increased mobility of the water molecules in the film as T increases.

FIG. 5.

(a) Ice in full adhesive contact with a hard substrate with a sinus-like corrugation. (b) The ice is displaced laterally with a distance s equal to half the wavelength of the roughness profile; it will result in positive pressure at the top of asperities and negative pressure in the valleys. As a result, water will flow in the premelted film from top of asperities to the bottom of the valleys. Since the premelted water film thickness must be approximately constant in time, ice is removed from the top of asperities and is added at the bottom of the wells. After some relaxation time t*, the stress free state (c) is obtained. The relaxation time t* depends on the temperature T and is shorter the closer the T is to the ice bulk melting temperature Tm because of the increased thickness of the premelted film and because of the increased mobility of the water molecules in the film as T increases.

FIG. 5.

(a) Ice in full adhesive contact with a hard substrate with a sinus-like corrugation. (b) The ice is displaced laterally with a distance s equal to half the wavelength of the roughness profile; it will result in positive pressure at the top of asperities and negative pressure in the valleys. As a result, water will flow in the premelted film from top of asperities to the bottom of the valleys. Since the premelted water film thickness must be approximately constant in time, ice is removed from the top of asperities and is added at the bottom of the wells. After some relaxation time t*, the stress free state (c) is obtained. The relaxation time t* depends on the temperature T and is shorter the closer the T is to the ice bulk melting temperature Tm because of the increased thickness of the premelted film and because of the increased mobility of the water molecules in the film as T increases.

Close modal

If the time for the displacement s is much longer than t*, the ice can nearly adiabatically adjust to the displacement and the breakloose friction force will nearly vanish. We believe this is the case close to the melting temperature Tm. However, as TmT increases, the thickness of the premelted film decreases and the mobility of the water molecules decreases, and at some point, there is no time for the water to move from the top of asperities to the bottom of the valley resulting in the maximum in the breakloose friction force.

Experimentally, for very smooth surfaces (here a fused silica plate with the rms roughness 0.5nm) for −20 °C < T < 0 °C, interfacial sliding is observed with the ice in complete contact with the substrate, with a friction force that decreases toward zero as the temperature approaches the ice melting temperature. We have proposed that this is due to changes in the thickness and molecule mobility of the water molecules in the premelted water film and to regelation. This reduces the elastic energy Uel/A0 (the energy for type II crack propagation) with increasing temperature and finally results in a negligible breakloose force at the ice melting temperature.

A premelted film of water and regelation may occur also for −25 °C < T < −20 °C, but we have found that at T ≈ −20 °C, the elastic energy GII = Uel/A0 associated with the mode II (frictional) crack propagation becomes similar to the energy per unit area, GI, to break the ice–substrate bonds in the normal direction as involved in type I crack propagation. This explains why the breakloose friction switches from interfacial sliding (type II crack propagation) to (type I) opening crack propagation at this temperature.

Here, it is interesting to note that an ice block on a silica surface with hrms ≈ 0.8 nm (when measured over a 2 × 2 μm2 surface area) did (at T = −10 °C) slide in adhesive contact, while this was not the case for any temperatures when hrms ≈ 7.5 nm where the contact broke by adhesive or cohesive failure.36 When comparing these results to the theory, we need to take into account that because both the substrate and the ice surfaces have the same roughness, the effective hrms is a factor of 2 bigger than hrms of the substrate. This gives 1.1 and 10.6 nm for the two cases considered above.

The study presented above may be relevant for other types of systems. For example, viscoelastic solids have frequency dependent modulus E(ω), where the low-frequency modulus E0 = E(0) may be much smaller than the high-frequency modulus E1 = E(∞). In this case, when the surface roughness is so small that complete contact occur at the interface (because of adhesion) in both the static and the sliding state, a contribution to the breakloose friction force will come from the effect described above. That is, after long enough time of stationary contact, the substrate roughness induced elastic energy stored at the interface is determined by E0, while at the onset of sliding, the viscoelastic solid behaves as a much stiffer material because of the high frequency perturbing deformations induced by the (short wavelength) substrate roughness (typical frequencies ωqv, where q is the wavenumber of a roughness component). However, in this case, the frictional shear stress is not negligible as assumed for ice (where a premelted water film is assumed to occur at the interface) and must be included in this study. In particular, a non-uniform shear stress would contribute to the energy stored at the interface and would change the definition of GII. In addition, if the frictional stress in the sliding state is too big, elastic instabilities may occur (e.g., Schallamach waves), and the breakloose process will not involve the type of slip of the viscoelastic solid at the interface as assumed in the theory presented above.

Consider now a rubber block with a perfectly smooth surface in complete adhesive contact with a rigid substrate with at most nanoscale surface roughness. We assume that the temperature is so high and contact time is so long that the rubber effectively is in the rubbery state with the elastic modulus E0. Next, the temperature is lowered so that the rubber reach its glassy state where the elastic modulus is E1. The deformation field induced in the rubber is “frozen-in” when the temperature is lowered to the glassy state. Thus, the force to displace the rubber block will again be given by the theory above with the elastic modulus given by the rubber modulus in the glassy state, which typically is similar to that of ice. Of course, if the roughness amplitude is larger than 1nm, the contact may no longer be complete in the displaced state, in which case the breakloose force would instead be determined mainly by the change in the interfacial binding energy.

People living in countries with cold winters know that it is not easy to remove ice formed on automobile windshields.49 Window glass, e.g., windshields, are produced from float glass. Float glass is a sheet of glass made by floating molten glass on a bed of molten metal, typically tin. This results in very smooth surfaces, which in ideal cases would have surface roughness resulting only from frozen capillary waves. Frozen capillary waves have the surface roughness power spectrum (see, e.g., Refs. 51 and 52)
(9)
Here, we have assumed that when reaching the temperature T*, the diffusive movement of the atoms in the glass becomes so slow that the capillary waves that exist at this temperature are frozen. Note that the gravity term ρg becomes important only for very small wavenumber roughness, namely, for q < (ρg/γ)1/2 ≈ 300 m−1, and the gravity term can be neglected in most cases giving C(q) ∼ q−2, which corresponds to a self-affine fractal surface with the Hurst exponent H = 0 (and fractal dimension Df = 3). For such surfaces, the rms slope is determined mainly by the shortest (nanometer) wavelength roughness. Real glass surfaces have larger roughness than predicted by (9) because of effects such as mechanical vibrations during cooling from the melted state or scratches from mechanical interactions or weathering.53 

Using the power spectrum given by (9) with T* = 600 K, γ = 0.3 J/m2, and ρ = 2200 kg/m2 and including the roughness wavenumber range q0 < q < q1 with q0 = 300 m−1 and q1 = 2 × 1010 m−1 give the rms roughness amplitude hrms = 0.27 nm and the rms slope ξ = 0.95. In Fig. 6, we show the elastic energy per unit surface area (4) divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s obtained for the surface roughness power spectrum given by (9). The asymptotic value for G = Uel/A0 for large displacement (here only s ≈ 10 nm) is given by (A2), GE*ξ2/q1. Because of the large rms slope for frozen capillary waves, the elastic energy per unit surface area produced when the ice displaced away laterally from its original position is G ≈ 0.45 J/m2. This is even larger than found above for the silica surface (see Fig. 1) because the surface roughness of the glass surface with capillary waves has a larger fractal dimension. Thus, we expect cohesive failure rather than adhesive failure at the ice–glass interface. This explains why it is hard to remove ice from the windshields of cars even for temperatures (say T = −5 °C) where one expects a very thin water film of premelted ice.

FIG. 6.

The elastic energy per unit surface area divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. For the surface roughness, the power spectrum is given by (9).

FIG. 6.

The elastic energy per unit surface area divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. For the surface roughness, the power spectrum is given by (9).

Close modal

To summarize, we have proposed a new mechanism for the ice breakloose friction force on very smooth substrate surfaces. The ice in the initial state is assumed to adhere fully in a stress free state on the substrate, and the breakloose force is associated with the roughness induced elastic energy stored at the interface as the ice is displaced laterally relative to the substrate. The temperature dependency of the breakloose friction force is attributed to regelation resulting from premelting.

Experimentally (see Ref. 36), for ice frozen on top of a silica surface with the rms roughness 0.5nm with decreasing temperature, we observe the following three different breakloose scenarios:

  • Sliding breakloose: the ice is sliding in full (adhesive) contact with the substrate.

  • Adhesive breakloose: ice–substrate bonds are broken, resulting in a transition to partial contact but with no ice left on the countersurface.

  • Cohesive breakloose: the shear plane occur inside the ice. In this case, ice is left on the countersurface.

It is clear that if the roughness is big enough, the breakloose due to a tangential force will always occur by process (c), and during slip, the hard substrate asperities will plow the ice surface, resulting in ice wear and a contribution to the sliding friction from the plowing.

Experimentally, for very smooth surfaces (here a fused silica window with the rms roughness 0.5nm), all three processes (a)–(c) occur depending on the temperature. For −20 °C < T < 0 °C, process (a) occurs with a breakloose friction force, which decreases toward zero as the temperature approaches the ice melting temperature. We have proposed that this is due to changes in the thickness and the mobility of the water molecules in the premelted water film and to regelation. This reduces the elastic energy Uel/A0 (the energy for type II crack propagation) with increasing temperature and finally results in a negligible breakloose force at the ice melting temperature.

A premelted film of water and regelation may occur also for −25 °C < T < −20 °C, but we have found that at T ≈ −20 °C, the elastic energy GII = Uel/A0 associated with the mode II (frictional) crack propagation becomes similar to the energy per unit area, GI, to break the ice–substrate bonds in the normal direction as involved in type I crack propagation. This explains why the nature of the breakloose friction switches from process (a) (type II crack propagation) to (b) (type I crack propagation) at this temperature.

E.C.T. acknowledges support from the Swedish Research Council (VR). We thank Trinh Nguyen and Matthew Fielden for their help when performing the AFM topography measurements.

The authors have no conflicts to disclose.

B. N. J. Persson: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal). E. C. Tyrode: Conceptualization (equal); Data curation (equal); Formal analysis (equal); Funding acquisition (equal); Investigation (equal); Methodology (equal); Project administration (equal); Resources (equal); Software (equal); Supervision (equal); Validation (equal); Visualization (equal); Writing – original draft (equal); Writing – review & editing (equal).

The data that support the findings of this study are available from the corresponding author upon reasonable request.

For q1s ≫ 1, we can neglect the J0(qs) term in (4) since the amplitude of J0(x) for large x varies as x−1/2. Thus, for large displacement s, we get
Assume that for q0 < q < q1,
and C(q) = 0 for q > q1 and q < q0. We get
Writing q = q0x, we get
or
(A1)
For H = 0, this gives
(A2)
where we have used that q1/q0 ≫ 1 in all cases of interest. As H → 1, we get
(A3)
Using that the rms roughness amplitude hrms is obtained as
and we get

Figure 7 shows an AFM topography map from a 1 µm × 1 µm surface region of the polished silica surface used in calculating the green line in Fig. 1. Note the polishing tracks that can be observed also on other surface area (not shown).

FIG. 7.

The surface height of the silica substrate as obtained from AFM topography measurements (1024 × 1024 data points) over a square area with the side 1 µm. The topography was used when calculating the green C(q) curve in Fig. 1.

FIG. 7.

The surface height of the silica substrate as obtained from AFM topography measurements (1024 × 1024 data points) over a square area with the side 1 µm. The topography was used when calculating the green C(q) curve in Fig. 1.

Close modal

Figure 8 shows the surface height probability distribution Ph over square areas with the side 10 µm (red curve) and 1 µm (green and blue curves). Note that Ph is nearly Gaussian as expected for a randomly rough surface.

For another study of the surface roughness of polished fused quartz, see Ref. 48.

FIG. 8.

The surface height probability distribution Ph of the silica substrate as obtained from AFM topography measurements (1024 × 1024 data points) over square areas with the side 10 µm (red curve) and 1 µm (green and blue curves). Ph is nearly Gaussian as expected for a randomly rough surface.

FIG. 8.

The surface height probability distribution Ph of the silica substrate as obtained from AFM topography measurements (1024 × 1024 data points) over square areas with the side 10 µm (red curve) and 1 µm (green and blue curves). Ph is nearly Gaussian as expected for a randomly rough surface.

Close modal

We assume that the onset of sliding occur by interfacial mode II (frictional) crack propagation at the interface. Thus, the ice–substrate contact region will consist of two regions: one where slip has occurs and one where no slip has occurred. At the boundary line (crack tip) between these two regions occurs a narrow strip (crack tip process zone, of width a few nm in the present applications) where the contact changes from the commensurate-like no-slip state to the incommensurate-like slip state. The static or breakloose friction force is the force needed to move the crack tip, which depend on the crack propagation energy G = Uel/A0.

The mode II interfacial crack propagation problem has been studied in many publications, but here we give for the readers convenience a simple summary of the basic theory for the simplest case where the kinetic frictional stress τk in the slip region can be neglected (for a treatment including τk, see Ref. 46), and we neglect edge effects. Consider a rectangular ice block in frozen contact with a silica substrate (see Fig. 4). On the one side of the ice block, a pressure or stress σ is applied. The contact between the ice block and the substrate is broken by interfacial crack propagation. We assume that the ice–silica contact in the broken area x < x0 exhibits zero friction (zero frictional shear stress). To increase the length of the frictional crack (mode II crack) by Δx requires an energy GwΔx, where w is the width of the block in the y-direction and where G = Uel/A0 is given in Fig. 3. This energy is supplied by the external stress σ. Thus, if the crack length increases with Δx, the surface where the pressure acts will move with the distance u determined by σ = Eux, where ux is the strain in a Δx wide slab of material on the cracked side of the ice block. Thus, the external work is
This external work is partly used at the interface to move the ice over the roughness barrier, GwΔx, and partly stored as bulk elastic energy in the added segment Δ x of compressed solid,
Using Uσ = Uel + GwΔx gives
or
Thus, the effective frictional shear stress τ = F/A0, where F = σhw and A0 = wL, is given by
In the present case, we have a cylinder shaped ice block. For this case, we denote the maximum effective shear stress by
(C1)
where α is a dimensionless geometrical factor that depends on the ratio h/L (where h ≈ 9 mm is the height of the ice cylinder and L ≈ 10 mm the diameter of the ice cylinder), and since h/L ≈ 1, we expect α to be of order unity.
Here, we derive the normal rms stress at the interface. We have39 
Using that40 
we get
Using (1) and (3), this gives
(D1)
We will denote σ2=σrms2. Note that for large displacement s,
(D2)
or σrms/E* ≈ ξ/2.

Figure 9 shows the average shear stress τ (red line) and the average normal stress σrms (blue) divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s for the surface roughness power spectrum shown by the green + red line in Fig. 1. Note that when the shear stress is maximal, the rms-average normal stress is about 15 times bigger than the shear stress, and that σrmsξE*/2 as the displacement s → ∞, while the tangential stress τ → 0 as s → ∞. The stress needed to plastically deform ice depends on the ice temperature and the deformation rate and is typically rather low under most practical conditions of order 10 MPa for T = −10 °C and the strain rates 1 s−1. At very short (nanometer) length scale, which is relevant for the applications in this paper where the crack-tip process zone is of order a few nm, the plastic yield stress may be much higher. Short wavelength roughness generates very high deformation rates even at low sliding speeds, and under this condition, thermal activation is strongly reduced, and since dislocations and other imperfections usually involved in plastic deformation of ice may be less important at the nanoscale, the plastic yield stress may be close to its ideal value αE, where α ≈ 0.1 (see Ref. 54). We conclude that nanoscale roughness may result in plastic deformation if the rms slope ξ is larger than 0.1, which is close to the rms slope we found on the polished fused silica plate. However, molecular dynamics studies of ice sliding on a corrugated substrate are needed to better understand the role of plastic deformations.

FIG. 9.

The average shear stress τ (red line) and the average normal stress σrms (blue) divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. The surface roughness power spectrum is shown by the green + red line in Fig. 1.

FIG. 9.

The average shear stress τ (red line) and the average normal stress σrms (blue) divided by the effective Young’s modulus E* as a function of the logarithm of the displacement s. The surface roughness power spectrum is shown by the green + red line in Fig. 1.

Close modal
1.
R.
Rosenberg
, “
Why is ice slippery?
,”
Phys. Today
58
(
12
),
50
(
2005
).
2.
A.-M.
Kietzig
,
S. G.
Hatzikiriakos
, and
P.
Englezos
, “
Physics of ice friction
,”
J. Appl. Phys.
107
,
081101
(
2010
).
3.
J. F.
Nye
, “
Glacier sliding without cavitation in a linear viscous approximation
,”
Proc. R. Soc. A
315
,
381
(
1970
).
4.
B. N. J.
Persson
, “
Ice friction: Glacier sliding on hard randomly rough bed surface
,”
J. Chem. Phys.
149
,
234701
(
2018
).
5.
L.
Bäurle
,
T. U.
Kaempfer
,
D.
Szabó
, and
N. D.
Spencer
, “
Sliding friction of polyethylene on snow and ice: Contact area and modeling
,”
Cold Reg. Sci. Technol.
47
,
276
(
2007
).
6.
R. W.
Liefferink
,
F. C.
Hsia
,
B.
Weber
, and
D.
Bonn
, “
Friction on ice: How temperature, pressure, and speed control the slipperiness of ice
,”
Phys. Rev. X
11
,
011025
(
2021
).
7.
L.
Canale
,
J.
Comtet
,
A.
Nigues
,
C.
Cohen
,
C.
Clanet
,
A.
Siria
, and
L.
Bocquet
, “
Nanorheology of interfacial water during ice gliding
,”
Phys. Rev. X
9
,
041025
(
2019
).
8.
M. E. H.
van Dongen
and
D. M. J.
Smeulders
, “
Ice speed skating: Onset of lubrication by frictional heating
,”
Europhys. Lett.
134
,
34005
(
2021
).
9.
A. D.
Roberts
, “
Rubber-ice adhesion and friction
,”
J. Adhes.
13
,
77
(
1981
).
10.
A. D.
Roberts
and
J. C.
Richardson
, “
Interface study of rubber-ice friction
,”
Wear
67
,
55
(
1981
).
11.
D. D.
Higgins
,
B. A.
Marmo
,
C. E.
Jeffree
,
V.
Koutsos
, and
J. R.
Blackford
, “
Morphology of ice wear from rubber-ice friction tests and its dependence on temperature and sliding velocity
,”
Wear
265
,
634
(
2008
).
12.
O.
Lahayne
,
B.
Pichler
,
R.
Reihsner
,
J.
Eberhardsteiner
,
J.
Suh
,
D.
Kim
,
S.
Nam
,
H.
Paek
,
B.
Lorenz
, and
B. N. J.
Persson
, “
Rubber friction on ice: Experiments and modeling
,”
Tribol. Lett.
62
,
17
(
2016
).
13.
N.
Orndorf
,
S.
Singla
, and
A.
Dhinojwala
, “
Transition in the acid-base component of surface free energy of ice upon the premelting of its second molecular bilayer
,”
J. Phys. Chem. C
124
,
19588
(
2020
).
14.
Ł.
Baran
,
P.
Llombart
,
W.
Rżysko
, and
L. G.
MacDowell
, “
Ice friction at the nanoscale
,”
Proc. Natl. Acad. Sci. U. S. A.
119
(
49
),
e2209545119
(
2022
).
15.
A. J.
Tuononen
,
A.
Kriston
, and
B.
Persson
, “
Multiscale physics of rubber-ice friction
,”
J. Chem. Phys.
145
,
114703
(
2016
).
16.
C.
Klapproth
,
T. M.
Kessel
,
K.
Wiese
, and
B.
Wies
, “
An advanced viscous model for rubber-ice-friction
,”
Tribol. Int.
99
,
169
(
2016
).
17.
F. E.
Kennedy
,
E. M.
Schulson
, and
D. E.
Jones
, “
The friction of ice on ice at low sliding velocities
,”
Philos. Mag. A
80
,
1093
(
2000
).
18.
B. N. J.
Persson
, “
Ice friction: Role of non-uniform frictional heating and ice premelting
,”
Chem. Phys.
143
,
224701
(
2015
).
19.
S.
Hemette
,
J.
Cayer-Barrioz
, and
D.
Mazuyer
, “
Thermal effects versus viscoelasticity in ice-rubber friction mechanisms
,”
Tribol. Int.
162
,
107129
(
2021
).
20.
F. P.
Bowden
, “
Friction on snow and ice
,”
Proc. R. Soc. London, Ser. A
217
,
462
(
1953
).
21.
A.
Almqvist
,
B.
Pellegrini
,
N.
Lintzén
,
N.
Emami
,
H.-C.
Holmberg
, and
R.
Larsson
, “
A scientific perspective on reducing ski-snow friction to improve performance in Olympic cross-country skiing, the biathlon and Nordic combined
,”
Front. Sports Act. Living
4
,
844883
(
2022
).
22.
K.
Kalliorinne
,
B. N. J.
Persson
,
J.
Sandberg
,
G.
Hindér
,
R.
Larsson
,
H.-C.
Holmberg
, and
A.
Almqvist
, “
Snow contact characterisation of cross-country skis: A macro-scale study of the apparent contact
,”
Lubricants
11
,
225
(
2023
).
23.
M.
Faraday
, “
On regelation, and on the conservation of force
,”
Philos. Mag.
17
,
162
(
1859
).
24.
J. W. M.
Frenken
and
J. F.
Van Der Veen
, “
Observation of surface melting
,”
Phys. Rev. Lett.
54
,
134
(
1985
).
25.
B.
Slater
and
A.
Michaelides
, “
Surface premelting of water ice
,”
Nat. Rev. Chem.
3
,
172
(
2019
).
26.
R.
Lipowski
, “
Critical surface phenomena at first-order bulk transitions
,”
Phys. Rev. Lett.
49
,
1575
(
1982
).
27.
U.
Tartaglino
,
T.
Zykova-Timan
,
F.
Ercolessi
, and
E.
Tosatti
, “
Melting and nonmelting of solid surfaces and nanosystems
,”
Phys. Rep.
411
,
291
(
2005
).
28.
D. T.
Limmer
Closer look at the surface of ice
,”
Proc. Natl. Acad. Sci. U. S. A.
113
,
12347
(
2016
).
29.
Y.
Qiu
and
V.
Molinero
, “
Why is it so difficult to identify the onset of ice premelting?
,”
J. Phys. Chem. Lett.
9
,
5179
(
2018
).
30.
Y.
Li
and
G. A.
Somorjai
, “
Surface premelting of ice
,”
J. Phys. Chem. C
111
,
9631
(
2007
).
31.
E. S.
Thomson
,
H.
Hansen-Goos
,
J. S.
Wettlaufer
, and
L. A.
Wilen
, “
Grain boundary melting in ice
,”
J. Chem. Phys.
138
,
124707
(
2013
).
32.
M.
Torabi Rad
,
G.
Boussinot
, and
M.
Apel
, “
Dynamics of grain boundary premelting
,”
Sci. Rep.
10
,
21074
(
2020
).
33.
M.
Demmenie
,
P.
Kolpakov
,
Y.
Nagata
,
S.
Woutersen
, and
D.
Bonn
, “
Scratch-healing behavior of ice by local sublimation and condensation
,”
J. Phys. Chem. C
126
,
2179
(
2022
).
34.
J.
Thomson
, “
On recent theories and experiments regarding ice at or near its melting-point
,”
Proc. R. Soc. London
10
,
151
(
1859
).
35.
F. P.
Bowden
and
T. P.
Hughes
, “
The mechanism of sliding on ice and snow
,”
Proc. R. Soc. London, Ser. A
172
,
280
(
1939
).
36.
J. F. D.
Liljeblad
,
I.
Furó
, and
E. C.
Tyrode
, “
The premolten layer of ice next to a hydrophilic solid surface: Correlating adhesion with molecular properties
,”
Phys. Chem. Chem. Phys.
19
,
305
(
2017
).
37.
H. H. G.
Jellinek
, “
Ice adhesion
,”
Can. J. Phys.
40
,
1294
(
1962
).
38.
S.
Yamaguchi
,
Y.
Suzuki
, and
T.
Otosu
, “
Perspective on sum frequency generation spectroscopy of ice surfaces and interfaces
,”
Chem. Phys.
522
,
199
(
2019
).
39.
B. N. J.
Persson
, “
On the elastic energy and stress correlation in the contact between elastic solids with randomly rough surfaces
,”
J. Phys.: Condens. Matter
20
,
312001
(
2008
).
40.
B. N. J.
Persson
, “
Theory of rubber friction and contact mechanics
,”
J. Chem. Phys.
115
,
3840
(
2001
).
41.
P. R.
Nayak
, “
Random process model of rough surfaces
,”
ASME J. Lubr. Technol.
93
,
398
(
1971
).
42.
B. N. J.
Persson
,
O.
Albohr
,
U.
Tartaglino
,
A. I.
Volokitin
, and
E.
Tosatti
, “
On the nature of surface roughness with application to contact mechanics, sealing, rubber friction and adhesion
,”
J. Phys.: Condens. Matter
17
,
R1
(
2004
).
43.
B. N. J.
Persson
, “
On the fractal dimension of rough surfaces
,”
Tribol. Lett.
54
,
99
(
2014
).
44.
B. N. J.
Persson
, “
Adhesion between an elastic body and a randomly rough hard surface
,”
Eur. Phys. J. E
8
,
385
(
2002
).
45.
T. L.
Anderson
,
Fracture Mechanics, Fundamentals and Applications
,
3rd ed.
(
CRC Press
,
2004
).
46.
B.
Lorenz
and
B. N. J.
Persson
, “
On the origin of why static or breakloose friction is larger than kinetic friction, and how to reduce it: The role of aging, elasticity and sequential interfacial slip
,”
J. Phys.: Condens. Matter
24
,
225008
(
2012
).
47.
L. B.
Boinovich
and
A. M.
Emelyanenko
, “
Experimental determination of the surface energy of polycrystalline ice
,”
Dokl. Phys. Chem.
459
,
198
(
2014
).
48.
N. I.
Chkhalo
,
S. A.
Churin
,
A. E.
Pestov
,
N. N.
Salashchenko
,
Y. A.
Vainer
, and
M. V.
Zorina
, “
Roughness measurement and ion-beam polishing of super-smooth optical surfaces of fused quartz and optical ceramics
,”
Opt. Express
22
,
20094
(
2014
).
49.
Y.
Zhang
,
Y.
Zhang
,
R.
Guo
, and
B.
Cui
, “
Method for testing shear and tensile strengths of freshwater/seawater ice
,”
Water
14
,
1363
(
2022
).
50.
Y.
Zhao
,
Y.
Wu
,
L.
Bao
,
F.
Zhou
, and
W.
Liu
, “
A new mechanism of the interfacial water film dominating low ice friction
,”
J. Chem. Phys.
157
,
234703
(
2022
).
51.
B. N. J.
Persson
Influence of frozen capillary waves on contact mechanics
,”
Wear
264
,
746
(
2008
).
52.
B. N. J.
Persson
, “
Contact mechanics for randomly rough surfaces
,”
Surf. Sci. Rep.
61
,
201
(
2006
).
53.
N.
Miyashita
,
A. E.
Yakini
,
W.
Pyckhout-Hintzen
, and
B. N. J.
Persson
, “
Sliding friction on ice
,”
J. Chem. Phys.
158
,
174702
(
2023
).
54.
B. N. J.
Persson
, “
Surface roughness induced stress concentration
,”
Tribol. Lett.
71
,
66
(
2023
); arXiv:2304.02159.