Dynamics of prolate spheroids in the vicinity of an air–water interface

In this article, we present the mobilities of prolate ellipsoidal micrometric particles close to an air–water interface measured by dual wave reflection interference microscopy. Particle’s position and orientation with respect to the interface are simultaneously measured as a function of time. From the measured mean square displacement, five particle mobilities (3 translational and 2 rotational) and two translational–rotational cross-correlations are extracted. The fluid dynamics governing equations are solved by the finite element method to numerically evaluate the same mobilities, imposing either slip and no-slip boundary conditions to the flow at the air–water interface. The comparison between experiments and simulations reveals an agreement with no-slip boundary conditions prediction for the translation normal to the interface and the out-of-plane rotation, and with slip ones for parallel translations and in-plane rotation. We rationalize these evidences in the framework of surface incompressibility at the interface.


Introduction
The dynamics of spherical particles dispersed in a Newtonian fluid near an interface has been widely studied in the past both experimentally and theoretically. Although for a fluid-fluid interface the particle dynamics is expected to be governed by slip flow boundary conditions (BC), 1-3 experiments often deviate from expectations. Indeed, while all experimental results relative to the movement parallel to the interface approach predictions for full-slip BC, [4][5][6] for dynamics normal to the interface the measured mobilities range from slip to no-slip BC. [6][7][8] ‡ Deviation from slip boundary condition has been recently rationalized considering the presence of small concentrations of interfacial active agents determining surface incompressibility. 8 In brief, particle motion induces an interfacial flow perturbing the otherwise homogeneous concentration of surface active agents, which in turn builds up a Marangoni stress opposing the flow. This results in the surface incompressibility condition, i.e. zero divergence of the interfacial flow.
In this perspective, our group has shown that surface incompressibility at air-water interface peculiarly affects the Brownian dynamics of spherical particles, with movement normal to the interface compatible with no-slip BC predictions, while the one parallel approaching full-slip BC ones. 6 This can be understood considering that the movement of a sphere normal to the interface generates a purely radial interfacial flow. As a consequence, incompressibility condition at the interface prevents interfacial flow and the particle dynamics is indistinguishable from the one normal to a solid wall 8 (noslip BC). Conversely, for a movement parallel to the interface, for which the radial component of the flow at the interface is weak, the additional condition of zero flow divergence at the interface only weakly affects the full-slip BC. 9 With high-speed camera measurements we have shown that even a small concentration of surface active agents (9 Â 10 À4 molecules per nm 2 ) ensures a surface incompressibility BC at the air-water interface for the typical timescales probed during a Brownian dynamics experiment. Such small concentrations are below the detection precision of available surface tension measurements and are most likely due to molecular species resulting from the chemical equilibria with the atmosphere, i.e. hydroxide or bicarbonate ions. They are therefore impossible to avoid even if all precautions are taken to keep the interface clean. 8 Differently from spherical particles, the dynamics of colloids with complex morphology is still almost unexplored. [10][11][12][13][14] Earlier studies have mostly focused on spheroidal particles close to solid boundaries, [15][16][17] as this morphology is the simplest after the spherical one. Because of anisotropy, particle hydrodynamic interactions with the confining interface result in a complex tensorial mobility, which demands the knowledge not only of the position but also of the orientation and the aspect ratio (i.e., the ratio between the axes lengths) of the spheroid. Studies of prolate spheroids in the presence of a fluid-fluid interface concerned mainly the absorption 18,19 and the 2D confined dynamics at the interface. 20 Very little is known on the dynamics of an ellipsoidal particle near an air-liquid interface. This lack of knowledge contrasts with the essential role played by the mobility of micrometric non-spherical entities close to fluid interfaces. Many systems of practical interest, like the motion of bacteria in biofilm formation, 21,22 the draining dynamics in flotation processes, 23 and the fabrication of Pickering emulsions, 24 demand indeed the knowledge of the near interface particle mobilities.
In this paper, we study the dynamic behavior of micrometric prolate ellipsoids near an air-water interface. Using Dual Wave Reflection Interference Microscopy (DW-RIM), 25 we are able to measure the dynamics of five (3 translational and 2 rotational) degrees of freedom of a particle, along with its aspect ratio and its distance from the interface. In particular, we here address how the surface incompressibility affects the ellipsoid rotational degrees of freedom and the translational-rotational cross-correlations. The latter result from the presence of the interface, which breaks the homogeneity of space. 26 The experimental results are compared with the outcomes of direct numerical simulations.

Sample preparation and experimental set-up
Prolate ellipsoids are produced at the desired aspect ratio A, i.e. the ratio between the major (a) and the minor (b) semi-axes, in the range 6 to 12 by stretching with a custom developed stretching device sulfate latex beads of diameter 8.7 AE 0.9 mm (MolecularProbest) embedded in a poly(vinyl alcohol) matrix and heated at 100 1C. 20 Such a large size has been chosen to reduce particles surface curvature. A low curvature for a sufficiently large surface area is indeed required to reduce the distortion of the interference fringes due to the finite numerical aperture of the microscope objective (see ESI †).
Samples are prepared dispersing ellipsoids in a cell made of a hollow glass cylinder enclosed on the top by a microscope slide and exposed to air on the bottom (Fig. 1). Pinning of the three phase (air-water-glass) contact line at the lower cylinder edges prevents water falling due to gravity. This experimental geometry introduces a finite curvature of the air-water interface which nonetheless has been evaluated to negligibly affect the experimental results (see ESI †).
Because of gravity, ellipsoids sediment towards the air-water interface until they reach an equilibrium gap distance at few hundreds nanometers where gravity and DLVO interactions balance. 25 If a large fluctuation normal to the interface occurs overcoming the gap distance, ellipsoids occasionally breach the interface and start the dewetting process. 25,27 Due to the DLVO repulsion, however, in most cases the particle remains at the equilibrium gap distance for the entire duration of the experiment. § In the present work we only consider the subset of ellipsoids that are not adsorbed at the interface throughout the duration of the experiment.
Ellipsoids number concentration is chosen low enough to allow neglecting inter-particles hydrodynamic couplings (mutual distance at the interface vicinity larger than 100 mm). Sample cell is located on the sample stage of the DW-RIM setup which, in brief, consists of a custom inverted microscope working in reflection with a 50Â long working distance objective. ¶ The sample cell and the objectives are enclosed in a flexible chamber in order to reduce sample evaporation and to protect the air-water interface from air displacements that would introduce additional particle drifts. The set-up is designed to simultaneously illuminate with a red (wavelength l r = 625 nm) and a blue (wavelength l b = 505 nm) low coherence monochromatic light sources8 and collect the interference signal from the reflections at the particle-water and at the air-water interfaces with an RGB camera.** The choice of an RGB camera is motivated by the requirement of simultaneously recording in different channels the light incoming from the two light sources. Fig. 1 Scheme of the experimental cell. The sample is contained in a hollow glass cylinder (2.), of height and inner radius of 4 mm, glued on a microscope slide (1.) positioned on the microscope sample stage. Micrometric ellipsoids (4.) are dispersed in the water solution where they sediment (black arrow) because of gravity towards the air-water interface (5.). The cell and the 50Â objective (6.) are enclosed in a flexible chamber (3.) designed to reduce water evaporation and particle drift at the air-liquid interface. § Since the electrostatic contribution to the DLVO potential is repulsive (both the particles and the air-water interface have a negative surface potential), the breaching probability significantly increases when the electrostatic interaction is screened. This can be done by adding salt to the solution to decrease the Debye screening length. 25 In principle, salt addition can be therefore exploited for finely tuning the distance of the ellipsoids from the interface in order to characterize the space dependence of the particle dynamics. ¶ Mitutoyo Apochromatic objective MY50X-805, 50Â, WD 3, NA 0.55. 8 Thorlabs collimated LEDs M625L3-C2 and M505L3-C2. ** Basler dart camera daA1280-54uc.
For each particle, we measure the time evolution of the particle center of mass together with the zenithal and azimuthal angles (see Fig. 2). An example of recorded interference pattern can be seen in Fig. 3a. Information on the particle movement parallel to the interface is obtained by tracking the position and the orientation on the image plane of the interference pattern, while information on the movement orthogonal to the interface is accessed from the interference pattern analysis. The superposition of the two wavelengths allows to unequivocally determine the particle-interface distance. 6,25 As a matter of fact, a single channel interference pattern contains a phase ambiguity by an additive integer multiple of 2p. Consequently, from the analysis of each monochromatic interference pattern the distance z between a point of the ellipsoid surface and the air-water interface is given by the expression with m l being an integer number, l = r, b for wavelengths l r and l b respectively, and z l,0 being the minimum possible value of z l obtained from the relative phase A [0,2p) measured from the interference pattern. Among all possible values given by eqn (1), the real distance is found by imposing z r = z b = z. Light sources of low coherence length (in water equal to 13 mm and 5 mm for red and blue LEDs respectively) have been chosen in order to prevent residual ambiguity on the distance due to matching periodicity of z r and z b (within the experimental incertitude) for increasing values of m r and m b . The interference pattern for particles at distances larger than 1 mm is therefore sensibly less contrasted than the one reported in Fig. 3a, as it has been verified by following colloids from the bulk to the air-water interface during sedimentation. Once sedimented, at equilibrium an ellipsoid lies with its major axis on average parallel to the interface, but because of thermal agitation y fluctuates in a range of a few degrees (less than 51) around zero. Once the fluctuating equilibrium position is reached, data are acquired for a few minutes at 30 fps in order to recover trajectories of 5000-10 000 frames length.
A more detailed description of the optical system and of experimental technique can be found in ref. 25.

Data analysis
Dedicated algorithms have been designed specifically for the tracking of the DW-RIM interference pattern of ellipsoidal particles. With reference to Fig. 2, we denote with (x 0 , y 0 , z 0 ) a fixed reference frame, with x 0 y 0 the air-water interface plane and z 0 the axis normal to the interface. Hence, y is the angle between the ellipsoid major axis and the interface plane, and f is the angle formed by the projected major axis on the interface plane. The distance z 0 from the interface and the zenithal angle y are collected analyzing, for each frame, the interference pattern. A threshold is first imposed to the acquired image in order to locate the interference maxima on the binarized image (Fig. 3b). The symmetry axis of the interference pattern, corresponding to the projection of the ellipsoid major axis on the interface plane, is then identified. The movement in the plane x 0 y 0 and the azimuthal angle f are accessed by tracking the interference pattern center of mass and its orientation. A local Cartesian coordinates system xyz is also defined, where x is the direction parallel to the ellipsoids major axis projection on the interface plane, y is the orthogonal direction in the interface plane and z z 0 is the direction orthogonal to the interface. From the in-plane coordinates, the displacements along the major axis projection (x) and the ones in the in-plane orthogonal direction (y) are obtained. A detailed description of the image analysis pipeline is reported in the ESI. † With this procedure the obtained resolutions are 100 nm for the ellipsoid x and y positions, a few nanometers for z and 0.31 for f. Zenithal angle y, finally, is measured with a relative error of about 10%.
The data reported in the present work refer to 18 different ellipsoids with aspect ratios in the range A A [6,12] and mean distance between the air-water interface and the ellipsoid Fig. 2 Geometry of the system where measured and simulated parameters are reported: reference systems fixed with respect to the laboratory (x 0 , y 0 , z 0 ) and to the ellipsoid projection on the x 0 y 0 plane (x, y, z), ellipsoid's center of mass coordinates (x 0 , y 0 , z 0 ), zenithal angle y between ellipsoid major axis and x 0 y 0 plane, azimuthal angle f in the x 0 y 0 plane. The simulation is performed within a cube of lateral size L and with boundary surfaces G p (particle surface), G w (all the external boundaries of the computational domain except the air-water interface) and G s (the airwater interface).

Fig. 3
Typical interference pattern of an ellipsoidal particle of aspect ratio 6 as it is recorded from the camera (a) and after binarization (b). Data relative to the red channel. Please note that the high intensity lateral wings are not part of the interference pattern but consequences of multiple reflection and refraction inside the ellipsoid. In (a), the red dashed line outlines the projection of the ellipsoid on the plane, as it is recovered from the interference pattern analysis.
center of mass in the range z 0 A [2.55,3.92] mm, corresponding to a gap distance of few hundreds nanometers.
In Fig. 4, typical time evolution of the measured positions (Fig. 4a) and orientations (Fig. 4b) is shown for one representative ellipsoid.
From the tracking, the mean square displacement (MSD) (Fig. 4c) and the mean square angular displacement (MSAD) (Fig. 4d) are evaluated: MSD k (t) = h[k(t + t) À k(t)] 2 i t , where k is x, y, z, y and f (for notation simplicity in equations we name MSD k both the MSD and the MSAD). From the MSD and MSAD the diffusion coefficient D k along the freely diffusive degrees of freedom x, y and f are obtained from a linear fit 2D k t. As expected, the in-plane translational diffusivity is larger along the major axis direction (x) than along the minor axis one (y). As it can be seen in Fig. 4a and b, the movement along z and the rotation along y are not free but confined in the DLVO potential due to the colloid and the air-water interface interactions. 25 As a consequence, MSD z and MSD y are both characterized by sub-diffusive behavior with a plateau for large lag times. They are therefore fitted with the analytic solution of the Langevin equation assuming a harmonic potential, 28,29 whose free fitting parameters are the diffusion coefficient and the potential stiffness. Translational and rotational mobilities are obtained as D k /k B T, where k B is the Boltzmann constant and T the temperature.
The fitting procedure decouples the confinement and the viscous drag contributions present in the MSD. Consequently, inter-particle differences in the DLVO and in the gravitational potentials due to the equivalent radius † † and aspect ratio polydispersities only affect the measured mobilities by changing the particle-interface average gap distances. The mobilities measured for different ellipsoids can be therefore duly compared once they are normalized through ad hoc simulations considering same ellipsoids geometries and gap distances as the ones measured. A list of the measured geometrical parameters relative to all the considered ellipsoids together with the corresponding measured mobilities is reported in a table in the ESI. †

Simulations
Because of the lack of exact models for ellipsoids at close distance to the interface, we compare experimental results with numerical simulations. The latter are designed considering a free-buoyant rigid ellipsoid suspended in an incompressible Newtonian liquid in proximity of an interface, where both fullslip (v > = 0, q > v 8 = 0) and no-slip (v > = 0, v 8 = 0) BC are simulated. ‡ ‡ Such BC have been chosen based on previous results on the Brownian diffusion of spherical particles in the same experimental conditions, 6 where the surface incompressibility condition results in an effective no-slip BC for particle diffusion along the z-axis and quasi-slip BC for the parallel one. The same no-slip BC is expected to hold for an ellipsoid diffusing perpendicularly to an incompressible interface. On the ellipsoid surface, no-slip boundary conditions are imposed. The system of partial differential equations along with the BC and the numerical details are reported in the ESI. † We solve the governing equations by applying a single nonzero component of the combined force/torque vector (F,T) and compute the six components of the translational and rotational velocity vector (u,x). By varying the imposed non-zero component of the force/torque vector, in six iterations we can compute the full mobility matrix M: 30 (2) that is symmetric and positive definite and can be partitioned as: with both the submatrices M a and M c being symmetric. The submatrices M a and M c are translational and rotational mobility matrices respectively, whereas M b and the transpose matrix M T b account for the translation-rotation coupling. where M c,kk is the rotational self mobility around the k-axis and M c,kk 0 accounts for the rotational cross-mobility around the axes k and k 0 . Due to the symmetry, for any particle-interface distance z, all the configurations obtained for different azimuthal angles f (rotations around the z-axis) are equivalent. Consequently, the particle-interface hydrodynamic interactions and, in turn, the mobility matrix are functions of z and y only, namely, M = M(z,y). Because of the invariance for rotations around z-axis, it is sufficient to compute the mobility matrix only for configurations with the major axis of the spheroid lying in the x 0 zplane; the mobility or friction matrices for orientations outside the x 0 z-plane are simply obtained by rotations around z-axis.

Results and discussions
Mobilities of particles measured in our experiments depend on the particle properties (size and aspect ratio), on the average distance from the interface and on the boundary conditions at the interface. Indeed, the presence of a boundary hinders the normal dynamics and either hinders or promotes the parallel one, depending on the particular BC. Since all the relevant parameters are experimentally accessible, a benchmarking with predictions can be made at the single particle level.
We therefore compare the measured mobilities M meas i;kk 0 (where i = a, c and k, k 0 = x, y, z according to the definitions of Section 2.3) with simulated mobilities M sim i;kk 0 obtained using slip and no-slip BC. To this aim, for each measured ellipsoid we run a simulation having as input the measured aspect ratio and average z for the particle at hand, with y kept fixed to zero. § § For the five mobilities we define the mobility relative difference as: Fig. 5a-c, we report DM i;k 0 k . M sim i;kk 0 for the translational movement, in the case of slip (full circles) and no-slip (open circles) BC. The reported values with their error bars represent the averages and standard deviations on the 18 ellipsoidal particles. The experimental translational mobilities parallel to the interface ( Fig. 5a and b) are close to the numerical predictions obtained with the slip BC, while the translational mobility orthogonal to the interface (Fig. 5c) approaches the result with no-slip BC.
These results are in agreement with the findings Villa et al. 6 for translational movement of spherical particles. There, the authors argue that, at low frequencies ( f { f c B 10 2 -10 3 Hz), the presence of even tiny quantity of surface active agents at the interface imposes a null divergence of the interfacial velocity field (div 2D v = 0). Such condition strongly affects the particle movement along the z-axis, that is thus equal to the one expected close to a solid wall, while the movement parallel to the interface is almost unaffected if compared to the full-slip BC case. 9 As shown in ref. 6 and 8, surface incompressibility condition is the limiting case for timescales much larger than 1/f c , when a Marangoni stress induced by surface active concentration gradients hinders the interfacial radial flow. In the opposite limit, for fast dynamics at low timescales the concentration gradient at the interface has no time to build up and the interface behaves as a full-slip one. Importantly, even a tiny surface coverage of surface active species (concentration range of c 0 B 10 À2 -10 À3 molecules per nm 2 ), far below the concentrations detectable by standard surface tension instruments, 31 rises the value of f c = 8c 0 k B T/(pZr), where r is the beads radius In the corresponding insets, sketches depict the movement direction or the rotation with respect to the interface (laying on the xy plane). Error bars represent the standard deviation evaluated upon different ellipsoids. § § We checked that the experimental fluctuations in z around its mean value and of y around 01 only slightly affect the simulated mobilities. and Z water viscosity, up to hundreds of Hz. This is a value much higher than typical acquisition frequencies of Brownian dynamics experiments. For the same experimental system, we recorded the dynamics of un-stretched 2r = 9 mm diameter sulphate latex beads using a fast camera (4000 fps). The high frame rate allowed us to measure the transition frequency from full-slip to surface incompressibility BC, resulting in f c B 2000 Hz and corresponding to 9 Â 10 À4 molecules per nm 2 . 6 The range of sampled Fourier components of the particle dynamics probed in the present work, determined by the acquisition frame rate of 30 fps, lies far below f c . Consequently, the Brownian dynamics can be rightfully described considering surface incompressibility at the interface. This explains the apparent discrepancy between the results retrieved for the translational mobilities. Indeed, BC are not different for dynamics normal and parallel to the interface: the unifying condition rationalizing both the directions is surface incompressibility.
Moving from spherical particles to the ellipsoids, translation normal to the interface still induces a purely radial flow at a free interface (see ESI †). As a consequence, surface incompressibility condition for the movement normal to the interface coincides with a no-slip BC, as for the spheres.
In Fig. 5d and e, we also report the in-plane and out-of-plane rotational mobility relative difference, again for slip (full circles) and no-slip (open circles) BC. As for the translations, rotation in the plane (Fig. 5d) is closer to simulations with full-slip BC at the interface. This can also be rationalized in the framework of surface incompressibility as the rotation of a prolate ellipsoid in the plane can be regarded as the translation of its tips, with therefore a weak radial component of the flow at the interface. On the other hand, for the rotation out of the plane (Fig. 5e) the two boundary conditions give comparable results within the experimental error. For small y, such a rotation indeed mainly consists of the vertical upward and downward movements of the two ellipsoid tips. In such an idealized picture, the interfacial velocity field would therefore have the symmetry of a dipolar field with both radial and non-radial relevant components. Please note that for both BC the significant measured deviations from the simulated values can result from the large experimental noise on y. ¶ ¶ We also evaluated the cross-correlations between different degrees of freedom, comparing their trend with simulated cross-mobilities. To this aim we calculate the rotational-translational cross-correlations C v z ,o and C v x ,o (at time delay 0) between the angular velocity of the out-of-plane rotations o and two components of the particle velocity, i.e. v z , along the interface normal and v x , along the ellipsoid major axis: where v i is either v z or v x and hÁ Á Ái t indicates the time average on t. Translational and angular velocities are defined as the variation of the corresponding degree of freedom between two consecutive frames multiplied by the frame rate. Since the dynamics is Brownian, these are average velocities lower than the instantaneous ones. As a consequence, strictly speaking, the cross-correlation in eqn (7) must be intended as the mean mutual displacement between the translational and the rotational degrees of freedom normalized by the product of the corresponding MSDs, evaluated at time delay 1. Angular dependency is obtained splitting the trajectories in 5 different subsets depending on the instantaneous value of y and evaluating eqn (7) for each interval of y. In Fig. 6a the cross-correlations averaged over all the measured ellipsoids are reported, where the abscissa represents the average angle y in each subset of data.
In the case of C v z ,o for negative y we observe that v z and o are anti-correlated, while for positive y they are correlated. To rationalize the observed cross-correlation, note that, as soon as y deviates from zero, one tip of the ellipsoid is farther from the interface than the other. In this condition, if the particle approaches the interface (negative v z ), the ellipsoid tip closer to the interface will experience a greater drag, independently from the assumed BC. Consequently the ellipsoid is subjected to a torque restoring the condition y It is worth noticing that the observed cross-correlation is not a consequence of the ellipsoids-interface potential well, which also acts to restore the y = 0 condition. Indeed, in the limit of a harmonic potential, a rotation of the ellipsoid along y would generate only a restoring torque, with no force in the direction of the z axis.
For C v x ,o , the negative sign observed indicates that the behaviour is anti-correlated for both y 4 0 and y o 0, as in systems where parallel drag decreases approaching the boundary. 26 In order to verify whether the cross-mobilities have the same trend as the cross-correlations, we simulated and reported in  Fig. 6b the cross-mobilities M b,zy (blue curves) and M b,xy (red curves) versus y, where M b,zy corresponds to C v z ,o , while M b,xy corresponds to C v x ,o . For each mobility both no-slip (dashed line) and slip (continuous line) BC are considered. The crossmobilities are simulated for an ellipsoid of aspect ratio A = 8.6 at a particle-interface distance z 0 /2b = 1.167, corresponding to the average aspect ratio and distance of the experiments. Shadowed regions in Fig. 6b account for the inter-sample aspect ratio variance88 DA = AE0.6. The trend of the simulated mobilities is the same of the experimental cross-correlations, for both full-slip and no-slip BC, thus implying that for the considered ellipsoids the sign of the correlation is only dictated by the presence of an interface, regardless of the specific BC applied at the interface. Since cross-mobilities are qualitatively similar for both full-slip and no-slip BC, subtle effects due to surface incompressibility condition are difficult to explore from the analysis of the cross-correlations.
This result is surprising as it is significantly different from the limit case (A = 1) of a sphere: as a consequence of a translation parallel to the interface, a bead is expected to rotate along the axis in the plane orthogonal to the translating direction with a positive correlation (M b,xy 4 0) for no-slip BC and with a negative correlation (M b,xy o 0) for full-slip BC, 26 whereas for the considered ellipsoid both mobilities are negative. This may hint at a crossover in the no-slip case from positive to negative correlation as a function of the aspect ratio which deserves further study.

Conclusions
In this work we first report a study on the mobilities of a micrometric ellipsoidal particle in the close vicinity of an airwater interface. Five different ellipsoidal mobilities (3 translational and 2 rotational) are measured and compared with dedicated simulations for slip and no-slip boundary conditions at the air-water interface. Moreover, the out-of-plane and the in-plane roto-translations cross-correlations are calculated and compared with simulated cross-mobilities. Analogously to the spherical particles, 6 surface incompressibility rules the translational mobilities: parallel ones approach predictions for slip BC, while the orthogonal one is close to no-slip BC. The rotational degrees of freedom strengthen the importance of surface incompressibility on the mobility of the ellipsoid, as the in-plane and out-of-plane rotations relate to the corresponding translational degrees of freedom. The roto-translational crosscorrelations, finally, are found in a qualitative agreement with the simulated cross-mobilities.
In the current work we measured the mobilities in a limited range of particle sizes, aspect ratios and distances from the interface, but the described setup enables a wider characterization as a function of those parameters, by changing the dimensions of the produced ellipsoids and by finely tuning the particle-interface distance. Gap distance can indeed be tuned by changing the DLVO interactions with the addition of salt 25 or through the use of optical tweezers.
We believe that our findings will be relevant to predict the movement of micro and nano entities close to liquid interfaces and membranes in a variety of fields ranging from biophysics to material processing. In particular, the relevance of the present results can be perceived considering the strong geometrical analogy of the investigated system with the one of a living bacteria moving close to an air water interface. Indeed, recent papers report how anomalies on the motion of Escherichia coli 13 and Pseudomonas aeruginosa 32 result from surface incompressibility, thus making the typical assumption of free air-water interface far from reality in such systems. By pointing out the additional effect of surface incompressibility on the rotational dynamics, we would like to take a first step towards a systematic study of surface incompressibility effects on the dynamics of micrometric organisms in water, stimulating further investigations in the field.

Conflicts of interest
There are no conflicts to declare.