X-ray reflectivity from curved liquid interfaces

X-ray reflectivity studies of the structure of liquid–vapour and liquid–liquid interfaces at modern sources, such as free-electron lasers, are currently impeded by the lack of dedicated liquid surface diffractometers. It is shown that this obstacle can be overcome by an alternative experimental approach that uses the natural curvature of a liquid drop for variation of the angle of incidence. Two modes of operation are shown: (i) sequential reflectivity measurements by a nanometre beam and (ii) parallel acquisition of large ranges of a reflectivity curve by micrometre beams. The feasibility of the two methods is demonstrated by studies of the Hg/vapour, H2O/vapour and Hg/0.1 M NaF interface. The obtained reflectivity curves match the data obtained by conventional techniques up to 5 c in micro-beam mode and up to 35 c in nano-beam mode, allowing observation of the Hg layering peak.


Introduction
X-ray reflectivity (XRR) techniques enable structural studies with molecular resolution on the electron-density profiles along interface normals (Als-Nielsen & McMorrow, 2011;Seeck & Murphy, 2015). Their application on liquid/liquid systems is of particular interest as such systems play a major role in biological systems and provide opportunities for technical applications. Conventional XRR methods require both an extended flat fluidic interface and the means to change the angle of the incident X-ray beam. At synchrotron X-ray sources, beam tilting is currently achieved by highly specialized liquid surface diffractometers. Available instrument designs employ either Bragg reflection at one (Als-Nielsen & Pershan, 1983;Schlossman et al., 1997) or at two subsequent crystals (Honkimä ki et al., 2006;Murphy et al., 2014;Seeck & Murphy, 2015) for tilting the beam down on the liquid surface. However, these setups and others outlined by Pershan & Schlossman (2012) are available at only a few synchrotron beamlines worldwide and do not exist at all at current or upcoming free-electron laser (FEL) sources, limiting the application of XRR methods to liquid interfaces. Alternative techniques such as angular-dispersive (Naudon et al., 1989) and energy-dispersive (Metzger et al., 1994;Pietsch et al., 2001) XRR allow a fixed setup, but necessitate specialized optics and energy-sensitive detectors.
Here, we present a method that circumvents the experimental challenges by utilizing the curved surface of a liquid drop to reflect an X-ray beam impinging horizontally. If a sessile drop is used in combination with a narrow beam (of the order of 100 nm high), the surface curvature enables reflection angles from 0 to 180 , selected by a simple vertical translation of the drop relative to the beam. Use of a broader (of the order of 10 mm) beam gives access to a large range of the reflectivity curve in one measurement, thus enabling fast acquisition for low reflection angles similar to the RAPID XRR (Yano et al., 2009) method. The sole experimental prerequisites for both modes are a vertical translation stage for the sample, a two-dimensional detector and the ability to vertically focus the X-ray beam down to the required height at the sample. These requirements are nowadays met at most modern synchrotron and FEL beamlines.
Extended discussions of the influence of the interface curvature on the reflected beam and the difficulties arising exist for solid/liquid (Briscoe et al., 2012) and liquid/vapour interfaces (Kawamoto et al., 1993;Regan et al., 1997). The application of the broad-beam geometry as a useful method was proposed by Bosio et al. in 1989(Bosio et al., 1989, but abandoned due to experimental difficulties at that time.
We demonstrate the feasibility of this reflection geometry on a polar liquid (water), a liquid metal (mercury), and the interface between a liquid metal and an aqueous electrolyte solution. Specifically, the nano-beam mode was used to record the XRR curve of the mercury/vapour interface, while the micro-beam mode was used to record the XRR curves of the water/vapour interface, the mercury/vapour interface and the mercury/electrolyte interface. These samples were chosen because of their high surface tension and high X-ray reflectance, which results in intense and sharp reflections. In addition, liquid metals such as mercury exhibit a strong stratification of the atomic layers close to the interface, resulting in a pseudo-Bragg peak (layering peak) at q z ¼ 2=ð3 1=2 rÞ, where r is the atomic radius (Magnussen et al., 1995;. For mercury, the peak is located at q z ' 2:2 Å À1 . Here, we pay special attention to the layering peak's visibility as an indicator for the reflectivity method's performance. A high-brilliance synchrotron beam is required as the layering peak intensity is not only six to seven orders of magnitude lower than the incident-beam's intensity, but is also superimposed on the liquid scattering peak.

Principle of measurement
We consider a sessile drop placed on the x/y plane centred around the z axis. The angle between the x/y plane and a local tangent to the surface then varies from 0 at the drop's apex to the contact angle Â C at its base (see Fig. 1 inset). A horizontal ray with negligible width and height travelling parallel to the x axis and impinging on the drop's surface is reflected at the local tangent plane. For drops that are symmetrical with respect to the x/z plane, the reflected ray is located in that plane. The angle of reflection measured relative to the x/y plane depends on the vertical position of the impingement point. The following discussion is restricted to the case of a drop with negligible curvature in the x/y plane and a ray impinging perpendicular to the drop's edge. This can be experimentally achieved by a very large puddle or confinement of the drop inside a rectangular container. If a highly collimated beam with finite height and width is used instead of a point-like ray, the reflection from the drop is spread over an angular range (see Fig. 1) which depends on the height of the beam, the drop's shape and the vertical position of the impingement point on the drop. We now derive the relationship between the vertical position z s of a ray with respect to the drop's apex and the reflection angle. Since a universal expression for the shape of a drop does not exist, we assume a semi-infinite puddle with its edge at x ¼ 0. It can then be shown that for non-wetting systems (Â C > 90 ) the angle between the local tangent and the horizontal i is given by (Adamson & Gast, 1997;de Gennes et al., 2004) where Ã c ¼ ð=gÞ 1=2 is the capillary length of the liquid given by its surface tension , the liquid's density and the gravitational acceleration g. For small incident angles, equation (1) can be approximated by a linear expression, This approximation deviates by less than 0.55% from the exact solution for i < 10 . A ray impinging on the drop is reflected by i þ f ¼ 2 i . Consequently, a perfectly collimated beam with finite height h is spread over an angular range of Sketch of the proposed reflection geometry. A horizontal X-ray beam (along the x direction) is reflected from the surface of a liquid drop and detected by a two-dimensional detector. Incident angle i and reflection angle f depend solely on the local curvature.
corresponding to a range of wavevector transfer values of approximately In the following paragraph, we will discuss the effect of vertical divergence (in the x/z plane) and horizontal divergence (in the x/y plane) on the reflection. Divergence Â D in the vertical direction leads to a distribution of incident angles around i with width Â C . If a collimated beam is spread into an angular range Áð i þ f Þ, the divergent beam is spread into the range Áð i þ f Þ þ Â D . The vertical divergence thus effectively increases the angular range covered by the beam, but lowers the q z resolution to which is most relevant at small angles. Let us now consider the horizontal divergence Â D only. While the angle between the local tangent and a non-divergent beam would be i , the effective angle eff in the divergent case ( Fig. 2a) can be calculated given the vector of an incoming ray r ¼ ðcos Â D ; sin Â D ; 0Þ and the tangent plane's normal vector n ¼ ðsin i ; 0; cos i Þ, which effectively lowers the incident angle and consequently the effective momentum transfer q eff z . Fig. 2(b) shows the relative deviation of q eff z from q z for divergence values of 10 mrad (Seeck et al., 2012) (blue) and 1 mrad (Riekel et al., 2010) (green). Even for high angles ! 2 Å À1 and a divergence of 1 mrad, the relative deviation is below 10 À6 . Fig. 2(c) shows the relative deviation of the effective momentum transfer from q z for a momentum transfer of 2 Å À1 over a wide range of divergence values. The deviation is negligible for divergence values 10 À2 rad. These considerations show that divergence effects are negligible at undulator beamlines of modern synchrotron sources.
To utilize this geometry for XRR measurements with a horizontal X-ray beam (wavelength ), the droplet has to be placed on a vertical translation stage to change the reflection angle by variation of z s . A two-dimensional detector, placed at a distance d from the impingement point ( Fig. 1), is required for measuring the reflected beam as well as the background. The reflection angle can be calculated from the vertical position s of the reflected beam on the detector relative to the primary beam position via i þ f ¼ arctanðs=dÞ ( Fig. 1). From this, the surface normal momentum transfer is easily calculated as q z ¼ ð4=Þ sin ½ð i þ f Þ=2. It should be noted that the position of the impingement point in the horizontal direction (along the primary beam) changes with z s as the point moves along the surface. With a fixed detector, this leads to an error in d of Ád L=2, where L is the size of the droplet in the horizontal direction. The observed beam position on the detector s therefore is systematically shifted by Ás=s ¼ Ád=d L=ð2dÞ. Ás can thus be minimized by increasing the detector-sample distance and reducing the size of the drop. To correct for this shift, it is necessary to know the inverse of the function describing the shape of the drop, xðz s Þ. Even for the simple case of a semi-infinite sample, this is not possible analytically (de Gennes et al., 2004) and will not be discussed here. However, near the apex the differential equation ðdz s Þ=ðdxÞ ¼ tan i ' 2 1=2 ðz s À aÞ=Ã c can be solved within the linear approximations [equations (3) and (4)]. The resulting offset in q z can then be estimated for a drop with contact angle Â C as where a ¼ 2Ã c sin ðÂ C =2Þ is the height of the sessile drop. At typical X-ray wavelengths, e.g. 1 Å , Áq off z is 10 À6 Å À1 for the mercury/vapour interface and the water/vapour interface at q z positions close to the critical angle and therefore negligible. Even though Áq off z increases with q z , it does not exceed 3 Â 10 À3 Å À1 at q z > 2:0 Å À1 for both systems.
A beam with a height h of the order of 10 to 100 mm can be utilized to probe a large range of the XRR curve with a single detector image. We will refer to this mode of operation as micro-beam mode throughout this paper. Because the primary beam in this case is scattered into a range of q z values, the intensity at each q z is lower compared with conventional XRR methods, leading to worse counting statistics. This is particularly problematic if diffuse scattering from the surface or Horizontal divergence. (a) Because of the angle between the horizontal (dash-dotted line) and a ray of a beam (orange) with non-zero horizontal divergence Â D , the incident angle is effectively lowered from i to eff . The effective angle can be calculated from the ray vector r and the tangent plane normal n. (b) Deviation of q eff z from q z for two different horizontal divergence values. The effect is most pronounced for large q z . (c) Deviation of q eff z from q z depending on divergence for a relatively high value of q z ¼ 2 Å À1 . scattering from the liquid bulk become comparable with the specular intensity (Braslau et al., 1988). For liquid surfaces, the diffuse scattering and the reflectivity are inherently correlated. Separation of these components requires complex mathematical models (Shpyrko et al., 2003;Runge et al., 2016) that take the beam profile and drop shape into account and are beyond the scope of this paper. For the systems studied here, no increased intensity due to diffuse scattering was experimentally observed (see below). The micro-beam mode is thus limited to low q z values where the reflected intensity is at least an order of magnitude higher than those of other scattering contributions.
In this mode, the q z resolution is determined by the detector pixel size and the detector's distance from the sample. Table 1 lists the calculated resolution for three combinations of pixel size and detector distance used in the experiments described in this paper. The table shows that these parameters have to be chosen carefully to obtain sufficient resolution at low angles. At a wavelength of = 1 Å , a detector with 50 mm pixel size placed at a distance of 1.2 m from the sample allows a resolution of 0.0005 Å À1 , sufficiently small to resolve the critical angle of, for example, water (0.0217 Å À1 ).
Micro-beam mode requires a highly collimated beam, otherwise the recorded intensity is convoluted with the angular distribution. Furthermore, the intensity depends on the vertical beam profile. To obtain a constant scaling factor, a flat-top beam profile is preferred for this mode.
For a small vertical beam size h of the order of nanometres, the spreading of the reflected beam Áð i þ f Þ and Áq z becomes small enough to resolve the critical angle of the studied liquids at typical X-ray wavelengths. Assuming a beam with h 200 nm and ¼ 1 Å , equation (4) yields Áq z 0.0013 Å À1 for the water/vapour and Áq z 0.0019 Å À1 for the mercury/vapour interface. Both values are considerably smaller than the respective critical angles. By translation of the sample in the vertical direction, the XRR curve can be recorded step by step, similar to classical techniques. In practice, a beam with h ' 100 nm is sufficiently small to apply this technique. We will therefore refer to this mode of operation as nano-beam mode. Assuming a beam of height h with a flat-top profile, the Fresnel reflectivity R F within the measured q z range Áq z given by equation (4) is largest near the critical angle, with variations of up to 15% for the mercury/ vapour interface. This variation biases the measured q z to lower values at q z positions where the reflectivity exhibits strong variations. For a beam with h ¼ 170 nm and ¼ 1 Å , the deviation of q z observed on the mercury/vapour interface is below 3% at the critical angle. The bias can be reduced further by using a beam with a Gaussian profile. Further uncertainties in the q z position arise from the vertical beam divergence Â D as described by equation (5).

Experimental details
Measurements in nano-beam mode were carried out at the nanofocus endstation of ESRF beamline ID13, which provides a vertical beam size of 170 nm at an energy of 15.2 keV and a vertical divergence in the order of 1 mrad (Riekel et al., 2010). These experiments were performed on clean mercury surfaces in N 2 (purity N5.7) atmosphere [Ã c = 1.91 mm (Adamson & Gast, 1997;Dean & Lange, 1999), drop height a = 3.71 mm]. Equation (5) yields a Áq z =q z of 12.5% at the critical angle due to the high beam divergence, dropping to 7% at q z ¼ 0:1 Å À1 and becoming less than 1% above q z ¼ 0:7 Å À1 . Instead of a free sessile drop, we used Hg samples contained in a purposebuilt KEL-F cell of rectangular shape (1.0 Â 0.8 cm, beam impinging on the long edge) to ensure reproducible drop size and orientation. A two-dimensional detector (Dectris Ltd, Eiger 4M, 75 Â 75 mm pixel size, distance d = 0.3 m) was used for data acquisition. XRR curves were recorded up to the liquid metal layering peak at approximately 2.2 Å À1 , corresponding to i þ f ¼ 16.42 at the given photon energy. Starting with the incoming beam half-shaded by the sample, the vertical sample position was changed stepwise and a single detector image was recorded for each point of the XRR curve. The detector images were monitor-normalized and corrected by subtraction of the background intensity measured 15 pixels left and right of the specular reflection. The specular intensity was extracted by integration of a 12 Â 150 pixel area. The reflection's vertical position s on the detector was obtained by a Gaussian fit to the vertical profile of the same area, then mapped to an angle i and converted to a momentum transfer value q z . Measurements were carried out on two mercury samples prepared under identical conditions. XRR data were recorded for q z > 1:3 Å À1 on the first sample and for q z < 1:2 Å À1 on the second sample with exposure times from 1 to 240 s per step.
The micro-beam mode was used to record XRR curves of the water/vapour interface [in air, Ã c = 2.71 mm (Adamson & Gast, 1997;Dean & Lange, 1999)] at a photon energy of 15 keV using a beam of 100 mm vertical size and a Dectris Ltd Eiger 1M detector (75 Â 75 mm pixel size, distance d = 1.196 m) with total exposure times of 1-100 s. Identical energy and beam size were used to record XRR curves of the mercury/vapour interface (in N 2 atmosphere) with a Dectris Ltd Pilatus 100k detector (172 Â 172 mm pixel size, distance d = 1.394 m) and a total exposure time of 1 s. The XRR curves of the mercury/electrolyte interface (in N 2 atmosphere) were recorded at an energy of 25 keV using a Dectris Ltd Eiger 1M detector (distance d = 1.22 m) with 1 s total exposure time. The water sample was contained in a sample cell with a rectangular shape (4.0 Â 19.5 cm, PTFE), mercury samples were contained in circular sample cells (diameter 4 cm, KEL-F). The recorded two-dimensional data sets were integrated in  Table 1 In micro-beam mode, angular resolution Áð i þ f Þ and momentum transfer resolution Áq z depend on detector pixel size and distance between sample and detector. the horizontal direction (ten pixels wide for water/vapour and mercury/electrolyte and five pixels wide for mercury/vapour to obtain similar integration regions with different detector pixel size) and the vertical pixel coordinates were mapped to q z values. All micro-beam experiments were carried out at room temperature using the sample and detector stage of the LISA liquid diffractometer (Murphy et al., 2014) at the PETRA III high-resolution diffraction beamline P08, where the vertical divergence does not exceed 10 mrad in collimation mode (Seeck et al., 2012). The vertical spread of the beam due to the divergence is below the size of two detector pixels and is therefore neglected. The sufficiently large incoming beam was cut by a vertical slit system to obtain a quasi-flat-top beam profile. At both sources, Gaussian beam profiles were confirmed by knife-edge scans and slit settings were optimized to avoid slit diffraction. All sample cell parts that came in contact with the sample liquids were cleaned in a mixture of H 2 SO 4 and H 2 O 2 (water/vapour and mercury/electrolyte experiments) or HNO 3 (mercury/vapour experiments) to avoid contamination of the interface and ensure cleanliness of the sample cell. Milli-Q water was used for the water/vapour sample and to prepare the 0.1 M NaF (99.995%, Alfa Aesar) aqueous solution for the mercury/electrolyte sample. High-purity mercury (99.99+%, Chempur) was used in all studies of the mercury interface.

Results and discussion
Examples of reflections from the mercury/vapour interface recorded in nano-beam mode are shown in Fig. 3. The presented images were recorded at q z = 0.073, 1.18 and 2.18 Å À1 and are background-corrected. The data shown in Fig. 3(a) (exposure time 10 s) include a part of the primary beam that is transmitted over the sample due to the very small vertical translation (visible at the bottom of the plot). As the sample is further translated upwards, the primary beam intensity is completely cut off by the sample. Data recorded at q z = 1.18 Å À1 (Fig. 3b) exhibit the expected Gaussian-shaped reflection. The reflection shown in Fig. 3(c) was recorded at the mercury layering peak position. The intensity at this position is expected to be a factor of approximately 10 À7 reduced compared with the primary beam intensity. Still, a distinct peak was observed with an exposure time of 480 s. Fig. 4 shows the XRR curves of the first (green triangles) and the second mercury sample (blue circles). The observed integrated intensity below the critical angle was 10 5 counts s À1 . Compared with the Fresnel reflectivity of a perfectly sharp interface (dotted black line), the data measured on the first sample exhibit the expected maximum at 2.2 Å À1 , arising from the stratification of the near-surface Hg into atomic layers. The data are in good agreement with previous studies of the free liquid Hg surface at room temperature, measured in conventional geometry on liquid surface diffractometers (Magnussen et al., 1995).    XRR curves of the mercury/vapour interface measured on two samples using a 170 nm beam. Data on the first sample (green triangles) exhibit the characteristic layering peak and are in good agreement with literature values (Magnussen et al., 1995). Oscillations in the XRR curve of the second sample (blue circles) suggest a mercury oxide layer at the interface. The dotted black line indicates the Fresnel reflectivity. Inset: incident angle versus sample translation measured on the second sample (blue circles) with fit of equation (1)  is commonly used to describe the liquid metal layering (solid green line), with a fixed layer spacing of 2.76 Å yielded an intrinsic roughness of (0.982AE0.013) Å and layer broadening of (0.459AE0.027) Å , which are in good agreement with the values reported previously. Data measured at low angles on the second sample revealed a distinct intensity minimum at 1 Å À1 , indicating an additional layer at the interface. A fit of a DCM modified by an adlayer (dashed blue line) yields a layer thickness of (5.82AE0.42) Å , which suggests a mercury oxide adlayer beginning to form at the interface Downs & Hall-Wallace, 2003;Gražulis et al., 2009Gražulis et al., , 2012, most likely a result of oxygen contamination due to a faulty sample cell seal. The inset in Fig. 4 shows the measured incident angle i in dependence of the vertical sample translation z s (blue circles, measured on the second sample). A fit of equation (1) to the data for i þ f > 1 is shown in black. The fit yields Ã c = (3.12AE0.05) mm, 62.5% larger than the value for a semi-infinite mercury drop. The larger Ã c as well as the notable deviations in i from equation (1) are most likely a result of the deviation of the surface shape from that of the semi-infinite case due to confinement by the sample cell. However, this does not affect the extraction of the reflectivity curve.
Micro-beam XRR curves were collected at all interfaces around the critical angle. In this range in particular, rapid changes in Rðq z Þ occur, providing a stringent test of this approach. The XRR curves obtained from the mercury samples are shown in Fig. 5(a) as red squares and green triangles (mercury/vapour) and blue circles (mercury/electrolyte). The solid lines indicate the theoretical Fresnel reflectivity including absorption effects. For both samples, the observed critical q z value q c matches the theoretical values very well [q c ¼ 0:065 Å À1 (Henke et al., 1993), see Fig. 5(a) inset]. To avoid damaging the detector, the mercury/electrolyte curve and the first part of the mercury/vapour curve (red squares) were recorded with the incident-beam intensity of 3 Â 10 11 photons s À1 (Seeck et al., 2012) attenuated by a factor of 10 5 and 10 4 , respectively. The maximum q z value of the mercury/electrolyte curve is therefore limited to 0.12 Å À1 or 1:8q c by the minimum count rate of one photon per exposure over the 1 s counting time (indicated by the dashed blue line). This q z range is comparable with that reported by Bosio et al. (1989). For the second part of the mercury/vapour curve (green triangles), the attenuation was reduced to a factor of ten, thus extending the accessible q z range to 0.24 Å À1 or 3:7q c . For q z > 0:2 Å À1 , the curve shows deviations from the Fresnel curve, most likely due to oxidation of the sample beginning to occur.
Another way to overcome the limitation of the q z range is to increase the exposure time. Fig. 5(b) shows the XRR curves of the water/vapour interface with total exposure times of 1 s (red circles), 10 s (green squares) and 100 s (blue triangles). As before, the horizontal dashed lines indicate the minimum absolute count rate of one photon per pixel over the exposure period. While this threshold decreases with increasing exposure time, the maximum q z value increases from 0.06 Å À1 for 1 s exposure time to 0.18 Å À1 for an exposure time of 100 s.
The accessible q z range can be extended further by binning of the raw data. Considering the high resolution of the raw data (approximately 0.0005 Å À1 ), the reduction of resolution from averaging over ten points does not impose serious limitations on the usability of the data. The effect of binning is shown in Fig. 5(c), in which the water/vapour XRR curve recorded with 100 s exposure time [as shown in Fig. 5(b) as blue triangles] was binned over q z intervals of 0.005 Å À1 width. In comparison with the original data, the lowest accessible intensity is more than a factor of ten lower (dashed blue line), thus extending the accessible q z range from 0.18 to 0.35 Å À1 or approximately 5q c .

Conclusion and outlook
We successfully implemented a method that enables the measurement of liquid-liquid and liquid-vapour XRR curves using a fixed, horizontal source beam, a vertical translation stage and a two-dimensional detector. Using a nanometre beam, it is possible to probe one single angle at a time (a) XRR curves recorded with the micro-beam method: mercury/vapour interface measured in two parts with 1 s exposure time each (red squares and green triangles, 100 mm beam) and mercury/electrolyte interface (blue circles, 20 mm beam, 1 s total exposure time). The solid black lines denote the Fresnel reflectivity for these interfaces. The dashed horizontal line marks the value corresponding to a total count of one photon. Data for mercury/vapour are offset by a factor of 100. Inset: region close to the critical angle of the Hg reflectivity. (b) XRR curves of the water/vapour interface recorded with 1 s (red circles), 10 s (green squares) and 100 s (blue triangles) integration time. Dashed lines indicate the level equivalent to a total flux of one photon. The Fresnel reflectivity is shown as a solid black line. (c) Water/vapour reflectivity measured in 100 s. The raw data were binned to extend the accessible q z range from 0.18 to 0.35 Å À1 (blue triangles). analogous to classical XRR techniques while maintaining the accessible q z range. Alternatively, a larger beam (of the order of micrometres) can be used that is reflected into a wider angular range, enabling parallel measurements of a larger range of an XRR curve up to 5 c . The feasibility of both modes was demonstrated by recording the XRR curves of the mercury/vapour, mercury/electrolyte and water/vapour interfaces. The obtained XRR curves match the results known from previous studies by conventional XRR techniques.
The presented method is experimentally simple and does not require a specific liquid diffractometer, but requires tradeoffs regarding accessible q z range, signal-to-noise ratio and most notably q z resolution. Which of these limitations apply depends on the employed measurement mode. As low divergence is crucial to the feasibility of the method, its application is limited to high-brilliance X-ray sources such as thirdgeneration synchrotrons and FELs.
XRR measurement on drop samples therefore will not render specialized diffractometers obsolete. It will, however, enable studies of liquid interfaces at sources lacking liquid diffractometers, especially ultra-fast XRR measurements via pump-probe techniques at FELs. A future combination of the method with drop shape analysis methods (Hoorfar & Neumann, 2011) should allow independent determination of z s ðxÞ, which can be used for improving the quality of the extracted data. In addition, the method may be extended to more complex experiments, for example the study of free falling or levitated drops.