# Adsorption-induced slip inhibition for polymer melts on ideal substrates

Mark Ilton, Thomas Salez, Paul D Fowler, Marco Rivetti, Mohammed Aly, Michael Benzaquen, Joshua D Mcgraw, Elie Raphaël, Kari Dalnoki-Veress, Oliver Bäumchen

► To cite this version:

Mark Ilton, Thomas Salez, Paul D Fowler, Marco Rivetti, Mohammed Aly, et al., Adsorption-induced slip inhibition for polymer melts on ideal substrates. Nature Communications, 2018, 9, pp.1172. 10.1038/s41467-018-03610-4. hal-01803716

HAL Id: hal-01803716

https://hal.science/hal-01803716

Submitted on 30 May 2018

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Adsorption-induced slip inhibition for polymer melts on ideal substrates

Mark Ilton<fnref n="1" />,², Thomas Salez ID ³,⁴,⁵, Paul D. Fowler<fnref n="1" />,⁶, Marco Rivetti⁶, Mohammed Aly⁷, Michael Benzaquen⁵,⁸, Joshua D. McGraw<fnref n="1" />,⁷, Elie Raphaël ID ⁵, Kari Dalnoki-Veress ID <fnref n="1" />,⁵ & Oliver Bäumchen ID ⁶

Hydrodynamic slip, the motion of a liquid along a solid surface, represents a fundamental phenomenon in fluid dynamics that governs liquid transport at small scales. For polymeric liquids, de Gennes predicted that the Navier boundary condition together with polymer reptation implies extraordinarily large interfacial slip for entangled polymer melts on ideal surfaces; this Navier-de Gennes model was confirmed using dewetting experiments on ultra-smooth, low-energy substrates. Here, we use capillary leveling—surface tension driven flow of films with initially non-uniform thickness—of polymeric films on these same substrates. Measurement of the slip length from a robust one parameter fit to a lubrication model is achieved. We show that at the low shear rates involved in leveling experiments as compared to dewetting ones, the employed substrates can no longer be considered ideal. The data is instead consistent with a model that includes physical adsorption of polymer chains at the solid/liquid interface.

<footnote n="1">¹Department of Physics & Astronomy, McMaster University, Hamilton, ON L8S 4M1, Canada. ²Polymer Science & Engineering Department, University of Massachusetts Amherst, Amherst, MA 01003, USA. ³Univ. Bordeaux, CNRS, LOMA, UMR 5798, 33405 Talence, France. ⁴Global Station for Soft Matter, Global Institution for Collaborative Research and Education, Hokkaido University, Sapporo, Hokkaido 060-0808, Japan. ⁵Laboratoire de Physico-Chimie Théorique, UMR CNRS Gulliver 7083, ESPCI Paris, PSL Research University, 75005 Paris, France. ⁶Max Planck Institute for Dynamics and Self-Organization (MPIDS), Am Faßberg 17, 37077 Göttingen, Germany. ⁷Département de Physique, Ecole Normale Supérieure/PSL Research University, CNRS, 24 Rue Lhomond, 75005 Paris, France. ⁸Ladhyx, UMR CNRS 7646, Ecole Polytechnique, 91128 Palaiseau Cedex, France. Correspondence and requests for materials should be addressed to O.Bäu. (email: oliver.baeumchen@ds.mpg.de)</footnote>
When a liquid flows along a solid surface, molecular friction at the solid/liquid interface can have a large effect on the overall dynamics. For a sufficiently high solid/liquid interfacial friction, the fluid velocity parallel to the interface goes to zero at the boundary. This “no-slip” boundary condition is a standard approximation for describing fluid flow at macroscopic length scales. In the past few decades, there have been many experiments measuring deviations from the no-slip boundary condition at microscopic length scales¹⁻⁵. These measurements have stimulated the interest in hydrodynamic slip for both the fundamental understanding of the molecular mechanisms involved, as well as the impact of slip on technological applications⁶⁻²⁰.

Hydrodynamic slip was first modeled by Navier²¹ using a solid/liquid stress balance at the substrate. This balance can be used to define the slip length as $b = [u_x/\partial_z u_x]|_{z=0}$. As shown schematically in Fig. 1a, the linear extrapolation length of the horizontal fluid velocity profile to zero is the slip length, *b*. Experimental techniques used to quantitatively measure the slip length can be broadly classified into four main categories²²: hydrodynamic drainage experiments, where the pressure is measured as fluid is squeezed out of (or drawn into) a small gap between two solid surfaces²³⁻²⁵; a direct measurement of the velocity profile near the interface, by either using tracer particles²⁶⁻²⁸ or fluorescence recovery²⁹; pressure-drop experiments for flows driven through microchannels³⁰⁻³², and dewetting experiments, which measure the retraction of a thin layer of fluid from a substrate³³⁻³⁶. Reported slip lengths from these methods vary from nm to mm³⁷,³⁸, while several studies also indicate the validity of the no-slip boundary condition ($b=0$)³⁹⁻⁴². From these experiments, some of the parameters affecting the magnitude of hydrodynamic slip have been elucidated which include the interfacial properties⁴³⁻⁴⁵, surface roughness²⁹,⁴⁶, shear rate¹,²³,⁴⁷, and molecular weight in polymer fluids¹³,³⁸,⁴⁸. Thus, *b* is highly sensitive to the solid/liquid combination.

With polymeric fluids flowing across ultra-smooth, low-energy surfaces, dewetting experiments have been used to measure large (>1 µm) slip lengths⁴⁸⁻⁵⁰. These results confirmed the scaling of the slip length with molecular weight for polymeric fluids, as originally predicted by de Gennes based on connecting the Navier boundary condition to polymer theory⁵¹; the experimental results furthermore suggested that the surfaces used were “ideal” as defined in⁵¹ to exhibit no chain adsorption. Specifically, Navier assumed that the substrate/liquid stress balance can be expressed as $\kappa u_x = \eta \partial_z u_x$, where $\kappa$, $u_x(x, z)$ and $\eta$ are the linear friction coefficient, the flow velocity parallel to the substrate and the viscosity, and $\partial_z$ denotes differentiation with respect to the surface-normal. De Gennes assumed a constant $\kappa$ and used the classic result of a strongly varying viscosity with molecular weight⁵², $M_w$. These essential ingredients allowed the prediction of the Navier-de Gennes model and subsequent experimental verification⁴⁸ of the scaling law $b = \eta/\kappa \sim M_w^2$.

A different experimental technique with shear rates orders of magnitude smaller than those characteristic of dewetting has recently been developed⁵³,⁵⁴. This capillary leveling technique is an experimental approach that invokes a film with an initially non-uniform thickness profile, resulting in an excess surface area compared to a flat film, the latter representing a metastable equilibrium. To suppress the energy cost of the excess surface area, surface tension causes the film to flow, driving it towards a uniform thickness. For an initially stepped film, and after a transient regime, the surface profile evolves in a self-similar fashion—that is, flow causes the profile to broaden, but the characteristic shape remains fixed and the broadening is determined by a power law in time. By monitoring the self-similar profile and fitting it to a lubrication model, quantitative nano-

<figure><img src="image_3.png" /><figcaption>**Fig. 1** Capillary leveling experiments on different substrates. **a** Schematic of the experimental procedure in which identical PS samples on two different substrates are annealed simultaneously. **b, c** Temporal series of experimentally measured atomic force microscopy profiles, normalized to demonstrate self-similarity of the film thickness profiles. Interfacial slip causes a faster broadening of the film thickness profile</figcaption></figure>

rheological information about the film can be obtained. Previous works on capillary leveling measured the viscosity of thin polymer films supported by a substrate with a no-slip ($b=0$) boundary condition⁵⁴, and addressed the infinite-slip ($b \to \infty$) limit using freestanding polymer films where the two free interfaces provide no resistance to flow⁵⁵.

In between the no-slip and infinite-slip extremes, the current work demonstrates the utility of capillary leveling as a quantitative probe of finite interfacial slip. We find that the measured slip length is independent of both temperature and sample geometry over the ranges studied. However, the slip length increases with the molecular weight of the polymer, and eventually saturates to a plateau at large molecular weight despite the apparently ideal character of the substrate. This latter fact contrasts drastically with the previous high-shear-rate dewetting studies using identical materials and conditions³⁶,⁴⁸⁻⁵⁰. Inspired by the case of permanent, chemically-grafted chains⁵⁶, we propose a low-shear-rate description based on the presence of a dilute number of physically adsorbed polymer chains, thus reconciling the two sets of experiments.

## Results

**Experimental approach.** Capillary leveling of stepped films is used to measure the slip length of polystyrene (PS) on Teflon™ fluoropolymer (AF) substrates (see Methods), a combination of materials which has been previously shown to exhibit interfacial slip³⁶,⁴⁸⁻⁵⁰. As a calibration, identical films of PS are prepared on silicon (Si) substrates since the Si/PS interface has no interfacial
slip⁵⁴. Both types of samples are annealed simultaneously and side-by-side on the same heating stage as outlined schematically in Fig. 1a. Surface profiles are obtained with atomic force microscopy (AFM). The self-similar profile of the Si/PS sample provides a calibration measurement of the PS capillary velocity, $v_c = \gamma/\eta$, where $\gamma$ and $\eta$ are the surface tension and viscosity, respectively. The value of $v_c$ depends on temperature and molecular weight, which are identical for the simultaneously studied Si/PS and AF/PS samples. The protocol thus allows the unambiguous and quantitative determination of the slip length of the solid/liquid (AF/PS) interface, the only differing quantity between the two simultaneously annealed samples. The measured film thickness profiles are self-similar in the reduced variable $x/t^{1/4}$, where $x$ is the horizontal coordinate and $t$ is the annealing time, for PS stepped films leveling on both substrates shown in Fig. 1b. We find there that PS films broaden more rapidly on the AF substrates than on the Si calibration substrates for all investigated molecular weights, and this faster broadening is more significant at higher PS molecular weight. In order to demonstrate that we can resolve slip lengths as small as a few tens of nanometers, we show a zoom on the dip region of the lower molecular weight film as an inset of Fig. 1b, where it can be clearly seen that the film on AF has also leveled further than that on Si. Capillary leveling thus provides one advantage over dewetting, for which small slip lengths have comparatively larger measurement error.

**Theoretical approach.** To extract quantitatively the slip length at the solid/liquid interface, we employ a continuum hydrodynamic model for the thin liquid film. Using the incompressible Stokes' equations in the lubrication approximation⁵⁷, and allowing for weak slip (slip length much smaller than the characteristic film thickness) at the solid/liquid interface, leads to a partial differential equation describing the evolution of the film thickness profile $h(x, t)^{58}$:

$$ \frac{\partial h}{\partial t} = -\frac{v_c}{3} \frac{\partial}{\partial x} \left[ (h^3 + 3bh^2) \frac{\partial^3 h}{\partial x^3} \right]. \quad (1) $$

We have also analyzed our experiments using the intermediateslip thin-film equation outlined in ref.⁵⁸. In the worst case, this discrepancy between the models results in a small (on the order of 30%) increase in the measured slip length, which does not affect any of the conclusions of this work. Besides, we stress that strongslip⁵⁸ or infinite-slip⁵⁵ descriptions would be incompatible with the observed $t^{1/4}$ self-similarity (see Fig. 1b).

By introducing an arbitrary reference length scale $h_0 = h_1 + h_2/2$, and the associated time scale $3h_0/v_c$, one can normalize Eq. (1) into a nondimensional form. Furthermore, for a given stepped initial profile (Fig. 1a), the rescaled solution $(h - h_1)/h_2$ of Eq. (1) is self-similar in the variable⁵⁸,⁵⁹:

$$ U_0 = \left( \frac{3x^4}{h_0^3 v_c t} \right)^{1/4}, \quad (2) $$

but depends intrinsically on two parametric ratios, $h_2/h_1$ and $b/h_1$. As a particular case, for a stepped initial profile with $h_2/h_1 \ll 1$ one can linearize Eq. (1). Nondimensionalizing the obtained equation by introducing the previous length scale $h_0$, but a different time scale $3h_0/[v_c(1 + 3b/h_0)]$, one obtains the result that the rescaled solution $(h - h_1)/h_2$ is now a single universal function of only the following

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 2** Capillary leveling as a probe to measure slip length. **a–c** For three different sample geometries, the rescaled self-similar theoretical profiles (dashed lines) fit the experimentally measured ones (solid lines) with one free parameter, the slip length $b$. Parameters are indicated in legends, and theoretical details are provided in main text. **c** The position $x$ is replaced by the radial coordinate $r$</figcaption></figure>

generalized variable:

$$ U_b = \left[ \frac{3x^4}{h_0^3 v_c (1 + 3b/h_0) t} \right]^{1/4} \quad (3) $$

In general, numerical solutions of Eq. (1)⁵⁹ can be used to fit the data (Fig. 2a). For the particular case of $h_2 \ll h_1$, analytical solutions of the linearized version of Eq. (1)⁶⁰ can also be used to fit the data (Fig. 2b). Since $v_c$ is fixed by the simultaneous no-slip calibration experiment, and the sample geometry is directly measured using AFM, the slip length $b$ is the only free parameter in fitting the theory to experimentally measured AF/PS profiles (Fig. 2). The slip length is found to be independent of temperature (Fig. 2a) in the considered range between 120 and
<figure><img src="image_7.png" /><figcaption>**Fig. 3** Slip-inhibition on ideal substrates. **a** Results from PS leveling experiments (blue circles) on AF substrates. Each data point consists of 2-12 individual measurements and the error bars represent the standard deviation for each measurement. For comparison, results from PS dewetting experiments (orange squares, data from$^{48}$) on AF substrates are also shown. Two equations with one free parameter (solid lines) describe both sets of experiments: it assumes adsorption of chains in the low-shear-rate leveling experiments and no chain adsorption in the high-shear-rate dewetting experiments. **b** The difference in the measured slip length between leveling (blue) and dewetting (orange) experiments is confirmed using a different substrate (SAM) for PS (9 kg mol$^{-1}$ at 110 °C and 65 kg mol$^{-1}$ at 135 °C)</figcaption></figure>

160 °C, and is not sensitive to changes in the sample geometry through $h_1$ and $h_2$ (Fig. 2b).

Complementary experiments were performed in a different geometry, in which the PS film was created with a cylindrical hole at the top$^{61}$ (see Methods) instead of a step. The result is shown in Fig. 2c. The slip length is determined by fitting the radially averaged normalized profile to the analytical asymptotic solution of the linearized axisymmetric thin-film equation$^{61}$—including weak slip through the variable $U_b$ above, with $x$ replaced by the radial coordinate.

**Effect of molecular weight on slip.** The effect of chain length on interfacial slip was studied using a series of 13 different PS molecular weights $8 \le M_w \le 373 \text{ kg mol}^{-1}$ spanning the range between unentangled and well-entangled PS$^{52,62}$. Results are shown in Fig. 3a (blue circles). At low molecular weight, the slip length increases with increasing PS molecular weight, but becomes approximately constant for molecular weights greater than ~100 kg mol$^{-1}$.

The molecular weight dependence of the AF/PS slip length has previously been found in dewetting studies (Fig. 3a, orange squares) to increase sharply at large molecular weight$^{48–50}$. The dewetting dependence contrasts markedly with the leveling results (Fig. 3a, blue circles) of the current work. Although the results from the two techniques agree at low PS molecular weight, the leveling results exhibit slip lengths that are reduced by two orders of magnitude at the highest molecular weights.

To determine if the difference in slip length at high molecular weight is specific to AF/PS, we performed a set of experiments with PS on self-assembled monolayer (SAM) substrates which are known to provide a slip boundary condition for PS$^{34,50,63}$ in dewetting experiments. In the SAM/PS experiments, both leveling

and dewetting measurements were performed for two different molecular weights (9 and 65 kg mol$^{-1}$). Results are shown in Fig. 3b. As for AF/PS, both molecular weights show a discrepancy between the slip length accessed with leveling and dewetting. Furthermore, the difference grows with molecular weight. Therefore, the observed difference in slip length between dewetting and leveling experiments exists also in the SAM/PS system, and is thus not specific to AF/PS.

## Discussion

Using lubrication theory$^{57}$ applied to the leveling experiments, typical shear rates at the substrate can be estimated through $\partial_z u|_{z=0} = v_c h \partial_x^3 h$. Using the experimental data (Fig. 1), we find strain rates of order $10^{-5}$ to $10^{-6} \text{ s}^{-1}$ for the molecular weights used. This range is three orders of magnitude lower than the average shear rates calculated for dewetting with the same molecular weights$^{49}$, and even lower if the maximum shear rate at the dewetting contact line is used.

A quantitative analysis of the residence time for polymer molecules under flow illustrates the effects of the different shear rates between leveling and dewetting experiments (see Supplementary Note 1 for details). The energy associated with the external force acting on an adsorbed polymer chain is estimated as $R_g^4 \partial_x P$, where $P$ denotes the (Laplace) pressure and $R_g$ the radius of gyration of the chain. Typical values from the leveling experiments ($\partial_x P \approx 0.5 \text{ kPa µm}^{-1}$), with $R_g = 29 \text{ nm}^{52}$, provide $R_g^4 \partial_x P \approx 3.5 \times 10^{-22} \text{ J}$, which is substantially smaller than thermal energy $k_B T \approx 5 \times 10^{-21} \text{ J}$ at $T = 150 \text{ °C}$, or the van der Waals interaction energies. Conversely, the energy associated with dewetting experiments is larger than $k_B T$ or van der Waals interaction energies by about one order of magnitude, thus crossing the threshold energy scale for chain desorption. The low shear rates for leveling thus lead to residence times for polymer chains at the solid/liquid interface that are long compared to the typical polymer relaxation time, and these adsorbed chains are expected to significantly contribute to the solid/liquid friction. This analysis is supported by studies, that find a shear dependence of polymer adsorption$^{64,65}$ and a recent work demonstrating that dewetting processes are faster when chain adsorption becomes weaker$^{66}$. Assuming that such an adsorption scenario is operative in leveling but not in dewetting, the large difference in measured slip lengths at high molecular weights between the low-shear-rate leveling experiments and high-shear-rate dewetting experiments (Fig. 3) can be rationalized, as detailed below.

The Navier-de Gennes model$^{51}$ predicts that under ideal conditions of no adsorption, where the polymer melt slips along a smooth passive surface, the slip length follows the form:

$$b_{\text{ideal}} = a \frac{\eta}{\eta_0} = a \frac{M_w}{M_0} \left[ 1 + \left( \frac{M_w}{M_e} \right)^2 \right], \quad (4)$$

where $a$ is the monomer size, $\eta$ is the polymer-melt viscosity, $\eta_0$ is the viscosity of a melt of monomers, $M_0$ is the monomeric molecular weight, and $M_e$ is the entanglement molecular weight. The right-hand side of Eq. (4) corresponds to a smooth interpolation between the Rouse and mean-field reptation regimes for the polymer-melt viscosity$^{52}$. To be consistent with the reference dewetting measurements$^{48,67}$, we have chosen the scaling of Eq. (4) to match the large-molecular-weight limit ($b_{\text{ideal}} \sim M_w^3$) found previously. We stress that the choice of a different power law (namely, the empirical reptation scaling $b_{\text{ideal}} \sim M_w^{3.4}$) does not alter the conclusions of the present work. Using the parameters and data from ref.$^{48}$ (i.e., $a = 0.3 \text{ nm}$, $M_0 = 104 \text{ g mol}^{-1}$, $M_e/M_0 = 517$), we recall on Fig. 3a (orange line) that Eq. (4) agrees with the dewetting data over the entire molecular-weight range used.
We now turn to the case of low-shear-rate experiments, and we describe the influence of transient physically adsorbed chains in an analogous fashion to the case of permanent chemically grafted chains⁵⁶. In the dilute-adsorption regime, adding the adsorption-induced frictional stress of Eq. (10b) from ref.⁵⁶ to the previous ideal frictional stress $\eta u|_{\perp=0}/b_{\text{ideal}}$, within the Navier-de Gennes construction⁵¹, leads to the dilute-adsorption prediction for the slip length:

$$b_{\text{ads}} = \frac{b_{\text{ideal}}}{1 + b_{\text{ideal}}/b^*} \quad (5)$$

Here, $b^* = aM_e/(nM_0)$ and $n$ is the number of adsorbed chains in one cross-sectional chain area. Invoking the parameters from ref.⁴⁸ as above, the dimensionless number $n$ is thus the only unknown quantity, and we now make the assumption (justified from the fit below) that $n$ does not vary (or varies weakly) with $M_w$. Stated differently, the density of physically adsorbed chains per unit surface is assumed to scale inversely with the cross-sectional chain area.

By fitting Eq. (5) to the leveling experimental data in Fig. 3a, we find an excellent agreement (blue line) and extract $n = 0.45 \pm 0.02$. Therefore, with a single free parameter we are able to reconcile the two very different experimental measurements of the slip length on the same AF/PS system. The saturation value of the slip length at high molecular weights and for low shear rates appears to be set by $b^* \sim N_e a$, where $N_e$ is the number of monomers between entanglements (omitting the numerical prefactor $1/n$ in front). The prefactor $1/n$ is expected to increase with shear rate, and to eventually diverge, thus allowing for a continuum of curves in between the two shown in Fig. 3a. In addition, the leveling data appears to have a sufficient resolution to observe for the first time the low-$M_w$ Rouse limit of the Navier-de Gennes prediction. We add two remarks: (i) we self-consistently find $n < 1$ which validates the dilute-adsorption hypothesis; (ii) $n$ is indeed nearly constant, as having a variation of $n$ with $M_w$ would correspond to not having a plateau for $b$ at large $M_w$.

Although the substrates we use are very smooth (see Methods), it is reasonable to expect that chain adsorption may occur at least temporarily at low shear stresses. First, even ultra-smooth surfaces show contact-angle hysteresis: if a contact line can be pinned on atomic-scale roughness, then so too can polymer chains. Secondly, unfavorable wetting does not imply repulsive interaction between the solid and the liquid, as wetting is rather controlled by a balance between this interaction and the solid-air interaction. Finally, molecular dynamics simulations have shown that adsorbed groups of connected monomers can occur at unfavorable interfaces, and the length of these adsorbed groups increases with molecular weight⁶⁸. Chains which are adsorbed for long enough durations of time to affect the interfacial fluid dynamics are likely to have multiple attached monomers. Therefore, larger adsorbed chains should exclude other chains from adsorbing to the substrate. The fact that $n$ is a constant smaller than 1 could be a signature of this exclusion mechanism.

In conclusion, we have demonstrated that capillary leveling can quantitatively probe interfacial dynamics at low shear rates. The use of simultaneously-annealed measurement samples on AF substrates and calibration samples on no-slip Si substrates, combined with weak-slip lubrication theory, allow for a robust one-parameter-fit of the slip length to the experimental data. For the case of PS films on an AF substrate, we find the slip length to increase with PS molecular weight before reaching a plateau value at large molecular weights. This contrasts with previous dewetting measurements on the same AF/PS system, which showed a strong increase in slip length at large PS molecular weights, consistent with the Navier-de Gennes model. Inspired by previous results

for grafted chains, we propose an extension of the Navier-de Gennes model which takes into account a dilute physical adsorption of polymer chains in the low-shear-rate leveling experiments, and no adsorption in the high-shear-rate dewetting experiments. With one free parameter, the proposed extension of the Navier-de Gennes model is able to capture the molecular-weight dependence of the slip length for both sets of experiments. Beyond providing new fundamental insights on the actively-studied problem of hydrodynamic slip, these results demonstrate that even ultra-smooth low-energy surfaces such as Teflon cannot always be considered as ideal substrates.

## Methods

**Substrate preparation and characterization.** Silicon (Si) wafers (obtained from University Wafer and Si-Mat) were cleaved into 1 cm × 1 cm squares. To create the calibration samples, the silicon wafers were rinsed with ultra-pure water (18.2 MΩ cm, Pall), methanol, and toluene (Fisher Scientific, Optima grade). To create a slip substrate, the wafers were coated with a thin film of the amorphous fluoropolymer AF (AF1600/AF2400, obtained from Sigma Aldrich) by dip coating from a dilute solution (solvent FC-72, obtained from Acros Organics, 0.5% w/w concentration solution, 0.5 mm/s retraction speed). Following the manufacturer's recommended procedure, the AF substrates were annealed for 20 min at 5 °C above the glass-transition temperature of AF (160 °C for AF1600 or 240 °C for AF2400) to remove residual solvent. The AF film thickness was 10–15 nm, measured using ellipsometry (EP3, Accurion). AFM (Caliber, Veeco; Dimension and Multimode, Bruker) measurements showed that the AF substrates have a 0.3 nm root mean square (RMS) surface roughness, and that PS droplets have a Young's contact angle of 88° on these substrates. As a second set of ultra-smooth, low-energy substrates, we decorated Si wafers with a dense SAM of octadecyltrichlorosilane (OTS) and dodecyltrichlorosilane (DTS, both purchased from Sigma-Aldrich). The SAM was composed of a mixture of equal parts OTS and DTS, providing the largest slip length for low-molecular-weight PS, see ref.⁶³ for details. Silane molecules cova-lently bind to the native oxide layer of the Si wafer during the established procedure⁶³,⁶⁹ for fabrication. These substrates have an RMS roughness of 0.2 nm and PS droplets have a long-time, receding contact angle of 63° on these mixed OTS/ DTS SAMs⁶³.

**Polymer film preparation.** Polystyrene (PS) with molecular weight ($M_w$) ranging between 8 and 373 kg mol⁻¹ and polydispersity less than 1.1 was obtained from polymer source and polymer standards service (PSS). Films with an initially stepped thickness profile (as in Fig. 1a) were made using a previously-described technique⁵⁴, with only minor modification. A bottom PS film (thickness range $h_1 = 100–800$ nm) and a top PS film (thickness range $h_2 = 40–400$ nm) were spin cast from a dilute or semi-dilute toluene solutions (liquid chromatography grade) onto freshly cleaved mica substrates (Ted Pella). The PS films were pre-annealed on mica in a home-built vacuum oven for at least ten times longer than the calculated longest relaxation time of the PS⁵² (pre-annealing temperature 140–180 °C, pre-annealing time 4–72 h; depending on the PS molecular weight). After annealing, the bottom PS film was floated onto an ultra-pure water bath (18.2 MΩ cm, Pall), and picked up onto either a silicon substrate (calibration) or AF substrate (measurement). The bottom film was then allowed to dry for at least 2 h before undergoing another annealing (annealing for at least two times the calculated longest relaxation time) to relax residual stress. The top PS film was then floated off its mica substrate onto the water bath. Sharp edges in the top PS film were created by the floating process for low $M_w$ PS⁷⁰, or for high $M_w$ PS by a procedure which involved floating onto Si, cleaving, and refloating onto the water bath⁵⁴. The sharp-edged top film was then picked up off the water bath with the previously-prepared bottom PS film on a substrate. A final drying of the film at room temperature concluded the sample preparation procedure. Identical procedures were applied for the experiments on the SAM substrates. Additional samples where the second film was prepared with a hole (as in Fig. 2c) were made in the same manner as the stepped films described above, except for the creation of sharp edges. For the hole geometry, a top film was floated onto the water bath and picked up using a metal washer with a millimetric circular hole to create a freestanding film. The top film was then heated above the PS glass-transition temperature in the freestanding state until small holes were nucleated with a diameter between ~3 and 10 µm. After quenching to room temperature, the top film was transferred onto the bottom film supported by a solid substrate. Full details on the hole-geometry sample preparation are presented in ref.⁶¹.

**Experimental setup.** Pairs of otherwise identically-prepared samples were used with only the substrate being different (Si or AF). The film thickness profile of each sample was determined by measuring the surface topography of the film using AFM, and averaging the 3D topography along the direction of translational quasi-invariance of the sample to obtain a 2D thickness profile. The pairs of samples were then placed side-by-side for simultaneous annealing in either the home-built vacuum oven or on a hot stage (Linkam, UK). For a given pair of samples, the
annealing temperature was held constant (between 120 and 160 °C), and chosen such that the PS was in its liquid melt state inducing the capillary-driven leveling of the thickness profiles. After a chosen duration of annealing *t*, the samples were rapidly quenched to room temperature, deep into the glassy state of PS, where the leveling process was temporarily halted. The broadening of the thickness profiles were measured using AFM. The samples were then further annealed, quenched, and measured again using AFM. The process of alternate annealing and AFM measurements was repeated until the measured thickness profiles became self-similar (Fig. 1b).

**Data availability.** All relevant data are available from the authors upon request.

Received: 11 August 2017 Accepted: 28 February 2018

Published online: 21 March 2018

## References

1. Léger, L., Raphaël, E. & Hervet, H. Surface-anchored polymer chains: Their role in adhesion and friction. *Adv. Polym. Sci.* **138**, 185–225 (1999).

2. Neto, C., Evans, D. R., Bonaccurso, E., Butt, H.-J. & Craig, V. S. J. Boundary slip in Newtonian liquids: a review of experimental studies. *Rep. Prog. Phys.* **68**, 2859–2897 (2005).

3. Lauga, E., Brenner, M. & Stone, H. Microfluidics: The no-slip boundary condition In: C., Tropea, A. L., Yarin. & J. F., Foss. (eds.) Springer handbook of experimental fluid mechanics. 1219–1240 (Springer, Berlin, 2007).

4. Bäumchen, O. & Jacobs, K. Can liquids slide? Linking stability and dynamics of thin liquid films to microscopic material properties. *Soft Matter* **6**, 6028 (2010).

5. Rothstein, J. P. Slip on superhydrophobic surfaces. *Annu. Rev. Fluid Mech.* **42**, 89–109 (2010).

6. Ajdari, A. & Bocquet, L. Giant amplification of interfacial driven transport by hydrodynamic slip: Diffusio-osmosis and beyond. *Phys. Rev. Lett.* **96**, 186102 (2006).

7. Heryudono, A. et al. Single-equation models for the tear film in a blink cycle: Realistic lid motion. *Math. Med. Biol.* **24**, 347–377 (2007).

8. Ren, Y. & Stein, D. Slip-enhanced electrokinetic energy conversion in nanofluidic channels. *Nanotechnology* **19**, 195707 (2008).

9. Vermesh, U. et al. Fast nonlinear ion transport via field-induced hydrodynamic slip in sub-20-nm hydrophilic nanofluidic transistors. *Nano Lett.* **9**, 1315–1319 (2009).

10. Ilton, M., Bäumchen, O. & Dalnoki-Veress, K. Onset of area-dependent dissipation in droplet spreading. *Phys. Rev. Lett.* **115**, 046103 (2015).

11. Chen, F., Peng, D., Lam, C.-H. & Tsui, O. K. C. Viscosity and surface-promoted slippage of thin polymer films supported by a solid substrate. *Macromolecules* **48**, 5034–5039 (2015).

12. Secchi, E. et al. Massive radius-dependent flow slippage in single carbon nanotubes. *Nature* **537**, 210–213 (2016).

13. Inn, Y. & Wang, S.-Q. Hydrodynamic slip: Polymer adsorption and desorption at melt/solid interfaces. *Phys. Rev. Lett.* **76**, 467–470 (1996).

14. Pryamitsyn, V. & Ganesan, V. Mechanisms of steady-shear rheology in polymer-nanoparticle composites. *J. Rheol.* **50**, 655 (2006).

15. Begam, N., Chandran, S., Sprung, M. & Basu, J. K. Anomalous viscosity reduction and hydrodynamic interactions of polymeric nanocolloids in polymers. *Macromolecules* **48**, 6646–6651 (2015).

16. Blake, T. D. The physics of moving wetting lines. *J. Coll. Interface Sci.* **299**, 1–13 (2006).

17. Klein, J., Kumacheva, E., Perahia, D., Mahalu, D. & Warburg, S. Interfacial sliding of polymer-bearing surfaces. *Faraday Discuss.* **98**, 173 (1994).

18. Klein, J., Perahia, D. & Warburg, S. Forces between polymer-bearing surfaces undergoing shear. *Nature* **352**, 143–145 (1991).

19. Yang, X., Wang, S. Q., Halasa, A. & Ishida, H. Fast flow behavior of highly entangled monodisperse polymers: 1. Interfacial stick-slip transition of polybutadiene melts. *Rheol. Acta* **37**, 415–423 (1998).

20. Boukany, P. E., Tapadia, P. & Wang, S.-Q. Interfacial stick-slip transition in simple shear of entangled melts. *J. Rheol.* **50**, 641–654 (2006).

21. Navier, C.-L. Mémoire sur les lois du mouvement des fluides. *Mémoires De l'Académie R. Des. Sci. De. l'Institut De. Fr.* **6**, 389–440 (1823).

22. Bouziques, C. et al. Using surface force apparatus, diffusion and velocimetry to measure slip lengths. *Philos. Trans. R. Soc. A* **366**, 1455–1468 (2008).

23. Craig, V. S. J., Neto, C. & Williams, D. R. M. Shear-dependent boundary slip in an aqueous Newtonian liquid. *Phys. Rev. Lett.* **87**, 054504 (2001).

24. Restagno, F., Crassous, J., Charlaix, É., Cottin-Bizonne, C. & Monchanin, M. A new surface forces apparatus for nanorheology. *Rev. Sci. Instrum.* **73**, 2292 (2002).

25. Garcia, L. et al. A micro-nano-rheometer for the mechanics of soft matter at interfaces. *Rev. Sci. Instrum.* **87**, 113906 (2016).

26. Ou, J. & Rothstein, J. P. Direct velocity measurements of the flow past drag-reducing ultrahydrophobic surfaces. *Phys. Fluids* **17**, 103606 (2005).

27. Joly, L., Ybert, C. & Bocquet, L. Probing the nanohydrodynamics at liquid-solid interfaces using thermal motion. *Phys. Rev. Lett.* **96**, 046101 (2006).

28. Li, Z. et al. Near-wall nanovelocimetry based on total internal reflection fluorescence with continuous tracking. *J. Fluid Mech.* **766**, 24 (2014).

29. Pit, R., Hervet, H. & Léger, L. Direct experimental evidence of slip in hexadecane: Solid interfaces. *Phys. Rev. Lett.* **85**, 980–983 (2000).

30. Cheng, J. T. & Giordano, N. Fluid flow through nanometer-scale channels. *Phys. Rev. E* **65**, 031206 (2002).

31. Cheikh, C. & Koper, G. Stick-slip transition at the nanometer scale. *Phys. Rev. Lett.* **91**, 156102 (2003).

32. Choi, C. H., Westin, K. J. A. & Breuer, K. S. Apparent slip flows in hydrophilic and hydrophobic microchannels. *Phys. Fluids* **15**, 2897–2902 (2003).

33. Redon, C., Brzoska, J. B. & Brochard-Wyart, F. Dewetting and slippage of microscopic polymer films. *Macromolecules* **27**, 468–471 (1994).

34. Fetzer, R., Jacobs, K., Wagner, B. & Witelski, T. P. New slip regimes and the shape of dewetting thin liquid films. *Phys. Rev. Lett.* **95**, 127801 (2005).

35. Münch, A. & Wagner, B. Impact of slippage on the morphology and stability of a dewetting rim. *J. Phys. Condens. Matter* **23**, 184101 (2011).

36. Haefner, S. et al. Influence of slip on the Plateau-Rayleigh instability on a fibre. *Nat. Commun.* **6**, 7409 (2015).

37. Granick, S., Zhu, Y. & Lee, H. Slippery questions about complex fluids flowing past solids. *Nat. Mater.* **2**, 221–227 (2003).

38. Wang, S.-Q. & Drda, P. A. Superfluid-like stick-slip transition in capillary flow of linear polyethylene melts. 1. General features. *Macromolecules* **29**, 2627–2632 (1996).

39. Zettner, C. M. & Yoda, M. Particle velocity field measurements in a near-wall flow using evanescent wave illumination. *Exp. Fluids* **34**, 115–121 (2003).

40. Honig, C. D. F. & Ducker, W. A. No-slip hydrodynamic boundary condition for hydrophilic particles. *Phys. Rev. Lett.* **98**, 028305 (2007).

41. Schaeffel, D. et al. Hydrodynamic boundary condition of water on hydrophobic surfaces. *Phys. Rev. E* **87**, 051001(R) (2013).

42. Ahmad, K., Zhao, X., Pan, Y., Wang, W. & Huang, Y. Atomic force microscopy measurement of slip on smooth hydrophobic surfaces and possible artifacts. *J. Phys. Chem. C* **119**, 12531–12537 (2015).

43. Barrat, J.-L. & Bocquet, L. Large slip effect at a nonwetting fluid-solid interface. *Phys. Rev. Lett.* **82**, 4671–4674 (1999).

44. Cho, J.-H., Law, B. & Rieutord, F. Dipole-dependent slip of Newtonian liquids at smooth solid hydrophobic surfaces. *Phys. Rev. Lett.* **92**, 166102 (2004).

45. Cottin-Bizonne, C., Cross, B., Steinberger, A. & Charlaix, E. Boundary slip on smooth hydrophobic surfaces: Intrinsic effects and possible artifacts. *Phys. Rev. Lett.* **94**, 056102 (2005).

46. Zhu, Y. & Granick, S. Limits of the hydrodynamic no-slip boundary condition. *Phys. Rev. Lett.* **88**, 106102 (2002).

47. Mhetar, V. & Archer, L. A. Slip in entangled polymer solutions. *Macromolecules* **31**, 6639–6649 (1998).

48. Bäumchen, O., Fetzer, R. & Jacobs, K. Reduced interfacial entanglement density affects the boundary conditions of polymer flow. *Phys. Rev. Lett.* **103**, 247801 (2009).

49. Bäumchen, O. et al. Slippage and nanorheology of thin liquid polymer films. *J. Phys. Condens. Matter* **24**, 325102 (2012).

50. Bäumchen, O. et al. Influence of slip on the Rayleigh-Plateau rim instability in dewetting viscous films. *Phys. Rev. Lett.* **113**, 014501 (2014).

51. de Gennes, P.-G. Ecoulement viscosimétriques de polymeres enchevetres. *Comptes rendus Hebd. Des. séances De l'Académie Des. Série B* **288**, 219–220 (1979).

52. Rubinstein, M. & Colby, R. H. *Polymer physics* (Oxford University Press, Oxford, 2003).

53. Buck, E., Petersen, K., Hund, M., Krausch, G. & Johannsmann, D. Decay kinetics of nanoscale corrugation gratings on polymer surface: Evidence for polymer flow below the glass temperature. *Macromolecules* **37**, 8647–8652 (2004).

54. McGraw, J. D., Salez, T., Bäumchen, O., Raphaël, E. & Dalnoki-Veress, K. Self-similarity and energy dissipation in stepped polymer films. *Phys. Rev. Lett.* **109**, 128303 (2012).

55. Ilton, M. et al. Capillary leveling of freestanding liquid nanofilms. *Phys. Rev. Lett.* **117**, 167801 (2016).

56. Brochard-Wyart, F., Gay, C. & de Gennes, P.-G. Slippage of polymer melts on grafted surfaces. *Macromolecules* **29**, 377–382 (1996).

57. Oron, A., Davis, S. & Bankoff, S. Long-scale evolution of thin liquid films. *Rev. Mod. Phys.* **69**, 931–980 (1997).

58. Münch, A., Wagner, B. & Witelski, T. P. Lubrication models with small to large slip lengths. *J. Eng. Math.* **53**, 359–383 (2005).

59. Salez, T. et al. Numerical solutions of thin-film equations for polymer flows. *Eur. Phys. J. E* **35**, 114 (2012).
60. Salez, T. et al. Capillary-driven flow induced by a stepped perturbation atop a viscous film. *Phys. Fluids* **24**, 102111 (2012).

61. Backholm, M., Benzaquen, M., Salez, T., Raphaël, E. & Dalnoki-Veress, K. Capillary levelling of a cylindrical hole in a viscous film. *Soft Matter* **10**, 2550–2558 (2014).

62. Fetters, L. J., Lohse, D. J., Milner, S. T. & Graessley, W. W. Packing length influence in linear polymer melts on the entanglement, critical, and reptation molecular weights. *Macromolecules* **32**, 6847–6851 (1999).

63. McGraw, J. D. et al. Influence of bidisperse self-assembled monolayer structure on the slip boundary condition of thin polymer films. *J. Chem. Phys.* **146**, 203326 (2017).

64. Cohen, Y. & Metzner, A. B. Adsorption effects in the flow of polymer solutions through capillaries. *Macromolecules* **15**, 1425–1429 (1982).

65. McGlinn, T., Kuzmenka, D. & Granick, S. Influence of shear on polymer adsorption kinetics. *Phys. Rev. Lett.* **60**, 805–808 (1988).

66. Wang, L. et al. The influence of polymer architectures on the dewetting behavior of thin polymer films: from linear chains to ring chains. *Soft Matter* **13**, 3091–3098 (2017).

67. Hénot, M., Drockenmuller, E., Léger, L. & Restagno, F. A proof of the Navier interfacial equation through quantitative comparison between solid-polymer melt and solid-elastomer friction. *ACS Macro Lett.* **7**, 112–115 (2018).

68. Smith, K. A., Vladkov, M. & Barrat, J. L. Polymer melt near a solid surface: A molecular dynamics study of chain conformations and desorption dynamics. *Macromolecules* **38**, 571–580 (2005).

69. Lessel, M. et al. Self-assembled silane monolayers: An efficient step-by-step recipe for high-quality, low energy surfaces. *Surf. Interface Anal.* **47**, 557–564 (2015).

70. Bäumchen, O. et al. Relaxation and intermediate asymptotics of a rectangular trench in a viscous film. *Phys. Rev. E* **88**, 035001 (2013).

## Acknowledgements

The authors thank Vincent Bertin, Pascal Damman, Frédéric Restagno, Barbara Wagner, Andreas Münch and Dirk Peschka for interesting discussions. They acknowledge financial support from the German Research Foundation (DFG) under grant BA3406/2 and NSERC (Canada). T.S. acknowledges financial support from the Global Station for Soft Matter, a project of Global Institution for Collaborative Research and Education at Hokkaido University. O.B. acknowledges financial support from the Joliot ESPCI Paris Chair and the Total-ESPCI Paris Chair. J.D.M. and M.A. were supported by LabEX ENS-ICFP: No. ANR-10- LABX-0010/ANR-10-IDEX-0001-02 PSL.

## Author contributions

M.I., T.S., J.D.M., E.R., K.D.-V. and O.B. conceived the project and designed research. M. I., P.D.F., M.R., M.A., J.D.M., and O.B. performed experiments. T.S., M.B., J.D.M., and E. R. developed the theory. All authors contributed to the analysis and interpretation of the results. M.I., T.S., J.D.M., E.R., K.D.-V., and O.B. wrote the manuscript.

## Additional information

Supplementary Information accompanies this paper at https://doi.org/10.1038/s41467-018-03610-4.

**Competing interests:** The authors declare no competing interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions/

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018