Mathematical Analysis of Nanostructured Surfaces: The Period-Scale Transform

*is work has been motivated by the urgent need for accurate and complete characterization of patterns consisting of almost periodic arrangements of specific features (trenches, bumps, holes, spikes, and so on) amply used in the industries of nanotechnology, microelectronics, and photonics. *e quantitative characterization of such surface structures demands mathematical methods able to reveal both periodand feature-scale aspects. Given that the conventional approaches (Fourier or wavelet transform) are limited to either periodicity or feature-scale characterization, our work contributes with the proposal of a transformation which combines Fourier and wavelet merits to quantify simultaneously the period and feature scale of a periodic or almost periodic surface pattern. *e output of our study has been (a) a detailed investigation of the mathematical properties of the proposed period-scale transform (PST) along with its relationship with other well-known transforms, (b) a presentation of some examples of PST of model 1D periodic surfaces to identify its benefits, and (c) first applications of PST in real profiles extracted from experimental polymer surfaces after plasma treatment.


Introduction
From the tiny scales of nanotechnology to the much larger feature sizes of geophysics and ocean studies, the morphology of surfaces plays an important role springing from its critical impact on surface properties and functionalities [1]. In order to understand and get control of these properties, it is necessary to find ways and tools to characterize surface morphology in a quantitative and concise way. A big concern in this endeavor is the presence of randomness and its co-existence with order on the arrangement of peaks and valleys on surfaces. Two extreme cases are the fully random surface where no point correlations exist and the fully periodic surface where a fixed surface unit is repeated endlessly along surface with a specific translation vector. More closely to the first case and usually met in nature is the scale-limited fractal (self-affine) surfaces where an anisotropic scale invariance in surface fluctuations prevails at small scales before being overwhelmed by uncorrelated randomness at larger scales. In technology, on the other side, and more specifically in nanotechnology where surface morphology can be controlled at nanometer scale, it is more common to fabricate surfaces either by top-down techniques or self-assembly processes with almost periodic structures consisting of similar humps or holes or trenches arranged in a periodic manner [2]. e same feature of (almost) periodicity is encountered in biological structures with an outstanding example-the surface morphology of moth eyes which has played a principal role in recent biomimetic technology [3]. Also, quite interestingly, recent studies have shown that periodic behaviour can be derived from stochastic mathematical models such as impulsive piecewise fractional functional differential equations or the stochastic Cohen-Grossberg neural networks with delays [4][5][6].
e most widely used method to detect and characterize periodicity of a surface is the Fourier transform (FT) and the diagram of the square of its amplitude versus frequency called power spectral density (PSD). Indeed, FT and mainly PSD are very sensitive to the presence of periodic repetition of a specific surface feature and display it with a well-defined peak at the frequency of repetition which is the inverse of period. Despite its wide and successful application to the detection and measurement of periodicity, PSD fails at the direct and straightforward extraction of the scale (width) of the repeated feature.
is does not mean that PSD peak pattern is not sensitive to feature-scale modifications. ese changes may affect the relevant PSD and therefore the number of well-defined peaks, but it is not always possible to get the real scales. For example, one can look at the 1D surfaces (Gaussian pulse trains) of Figures 1(a) and 1(b) both having the same period but repeated features with different scales. e PSDs shown in Figures 1(c) and 1(d) demonstrate the inability to map PSD peak changes to feature-scale values in a direct way despite the apparent effect of feature scale on the number of peaks.
FT provides a very accurate representation of the frequency content of a surface (or more generally of a signal). However, it does not have good localization properties. For instance, high jumps or holes in the surface cannot be read off easily from the Fourier transform. In many applications, this localization is desirable and can be first achieved by considering the windowed Fourier transform (WFT) (see [7] or [8]). WFT fixes a scale (the width of the window) and analyzes the surfaces from the point of view of this scale. e wavelet transform (WT) provides a similar description, but in contrast with WFT, the so-called "wavelets" do not have a fixed width but the width is adapted to the frequency. As a result, WT is a significant tool for analyzing short-lived phenomena such as singularities or steep and sudden scale changes. WT has found a plethora of profound applications including (but not limited to) geophysics (analysis of seismic data), engineering, and technology.
For instance, wavelet transform methods have been exploited in conditioning monitoring, that is, the process of monitoring some parameter or condition in machinery in order to detect any significant change that would be evidence for a developing default. is early detection and isolation of developing faults and the prediction of fault propagation allows for preventive maintenance and it can also serve as a countermeasure to the possibility of catastrophic incidence as a result of a failure. Furthermore, machinery and technological devices (for example wind turbines) are often subject to unexpected failures due not only to operational but also to environmental reasons. In such cases, prognosis techniques that can forecast the remaining life of machinery assure reliable production and lower maintenance costs. Several techniques have been developed in recent years towards this direction, which are based on WT, artificial intelligence, stochastic methods, and so on (see [9][10][11][12]). We also refer to the review [13] which investigates condition monitoring and fault diagnosis using various methods like wavelets and intelligent methods.
In summary, FT can be used to analyze the frequency content of a surface (or signal), but it has no localization properties and it cannot be used to detect the scale of the repeated feature of an almost periodic surface. On the other hand, the wavelet transform can capture the scale changes but cannot detect periodicity directly. Consequently, the need for a transform being capable of detecting and quantifying in a straightforward way in both aspects of periodicity (period and feature scale) is emerging. e aim of this paper is to propose and elaborate such a transform, which we name period-scale transform (PST). Its derivation is based on the periodization of a fundamental function satisfying some decay conditions. Periodization of functions has been studied extensively in analysis (see [14]), but, to the best of our knowledge, has not been used to define a transform.
is paper is organized as follows. In the next section, we define the transform rigorously, while in Section 3 we investigate its fundamental properties. In Section 4, we discuss its relevance to Poisson summation formula. e limit behaviour of the PST and the connection with other transforms (Fourier, wavelet, and Laplace) are the issues discussed in Sections 5 and 6, respectively. First examples of transform applications in selected model and real 1D surfaces are reported in Section 7 and the paper closes with summarizing the main conclusions in Section 8.

Definition of the PS Transform
We start with a non-zero function g: R ⟶ C satisfying the decay condition for any x ∈ R, where C, ε are positive constants. e function g will be called the fundamental function. Given (scale) s > 0 and (period) λ > 0, we generate a periodic function G(·; s, λ) as follows. Firstly, we change the scale of the function g, that is, we set g s (x) � g(x/s). en, we apply the periodization operation to the function g s in order to obtain a periodic function with period λ. In other words, the function G is defined as e decay condition guarantees that the function G(·; s, λ) is well defined. Furthermore, for any x ∈ R, the following inequality holds: In some cases, the precise values of the constants are not of great importance. Consequently, instead of using (3), we often appeal to the following coarser inequality: or even to the simple fact that, for 0 < s ≤ λ, the (family of ) function(s) G(·; s, λ) is bounded |G(x; s, λ)| < C, where C denotes some absolute constant (not always the same). Furthermore, we consider the mean value of the function G in the interval [0, λ] which will be denoted by G. at is, An easy calculation shows that 2 Mathematical Problems in Engineering where interchanging the sum with the integral is justified by the dominated convergence theorem. e function G(x; s, λ) − G is the kernel of the periodscale transform which is defined as follows.
Definition 1. Let f: R ⟶ R be a function in L 1 (R). e period-scale transform (PS transform or PST in short) of f with respect to the fundamental function g is defined by Since G(·; s, λ) is bounded, for f ∈ L 1 (R), the function f · G also belongs to L 1 (R) and the above transform is well defined.
Additionally, using the dominated convergence theorem, we can analyze further the defining formula (7) of the PST and we obtain Remark 1. In principle, formula (7) (or, equivalently (8)) makes sense for any s > 0 and λ > 0, that is, the PST can be defined on the domain (0, ∞) × (0, ∞). However, given a function f, detecting and characterizing the periodicity and the scale feature of f means that the scale s is not bigger than the period λ. Hence, without losing important information for our purposes, the domain of PS g f can be restricted to the

Linearity and Boundedness.
It is clear that the PST is linear. at is, for any f, h ∈ L 1 (R) and any scalar α, we have Furthermore, the PST defines a bounded linear operator  Figure 1: (a, b) Two profiles (Gaussian pulse trains) with identical periods (�80) but different peak widths (5 and 10, respectively). (c, d) e power spectra (PS) of the two profiles of (a) and (b) showing the difference in two cases: less peaks in profile of (b) with wider features. Despite the difference in the number of peaks in two power spectra, one cannot extract the real widths of profiles in a quantitative straightforward way.
Indeed, since G is bounded, formula (7) gives that

Scaling Properties.
Suppose that f is a periodic function, is also a periodic function; however, both the period and the scale have been multiplied with the factor 1/2. More generally, it is reasonable to expect that the PST of the function f(tx), t > 0, is connected with the PST of the original function f evaluated at the point (ts, tλ). is is true and is given in the next theorem.
, then the following holds: Proof. By the definition of the PST, it follows that Setting tx � x ′ implies that e second assertion follows immediately by the first one and the linearity of the PST.
Having dealt with dilation, we now turn our attention to translation.
at is, given a function f ∈ L 1 (R) and any y ∈ R, we are interested in finding the PST of the function e PST of the translation τ y f of f is given by where h � τ − y/s (g).
Observe that the fundamental function h in the righthand side of the above equation depends on s. In particular, this implies that the PST is not phase invariant.
Proof. It is not hard to see that the function τ − y/s (g) satisfies a decay condition similar to (1), where the constant C depends on y/s. Hence, it can be used as a fundamental function for the PST. Now, for fixed s > 0 and λ > 0, we observe that Clearly, we have and the desired result follows immediately from equation (8).
In the special case λ � y, since G(·, s, λ) has period λ, it follows that translating f by λ brings no change to the PST.
at is, we have the next result whose proof is straightforward.
3.3. Continuity. e purpose of the present section is to show that the PST of any f ∈ L 1 (R) is a continuous function provided that g is continuous. Firstly, we need the following lemma whose proof depends on elementary calculus arguments and it is omitted. Lemma 1. Assume that the fundamental function g is continuous. en, the periodization G of g is continuous with respect to (s, λ). at is, for any (fixed) x ∈ R and any s 0 > 0, Using the above lemma, we can prove the main result of this section.

Theorem 3. Assume that the fundamental function g is continuous. en, for any
For 0 < s ≤ λ, the family of functions G(·; s, λ) is bounded by some absolute constant C independent of s and λ.
erefore, the dominated convergence theorem and the previous lemma imply that and consequently, By equation (8), it follows now that the PST of f is a continuous function. e aim of this section is to present some other forms of the PST which will be useful in the sequel. We start with the following proposition whose proof is straightforward.

Proposition 2.
Assume that the function f satisfies a decay condition similar to (1). en, the PST of f is given by where Recall that for functions f, g ∈ L 1 (R), the convolution of f and g is defined (for almost every x ∈ R) by Using convolution of functions, formula (8), which gives the PST, can be rewritten in the following form.

Proposition 3.
For any function f ∈ L 1 (R) and any 0 < s ≤ λ, we have 3.5. e Calculus of the PST. In this section, we calculate the PST of the derivative z k f of a function f, for any positive integer k.
, and let f ∈ L 1 (R) and k be a positive integer. We assume that f, g are differentiable up to k times and the derivatives satisfy the following: (ii) z i g satisfies a decay condition similar to (1), for any en, the PST of z k f is given by

Mathematical Problems in Engineering
Proof. By the definition of the PST, we obtain Integrating by parts k times gives that Remark 2. If we are not interested in much generality, we may assume that f, g belong to the Schwartz class of functions. Recall that a function f: R ⟶ C is called a Schwartz function if f is infinitely differentiable and for any pair of nonnegative integers k and n, we have Clearly, if f, g are Schwartz functions, then the conditions of the above theorem are fulfilled.

The PST and Discretization: The Poisson Summation Formula
It is well known that periodization and discretization are dual procedures in the following sense: the periodization of a function is completely determined by discrete values of the original function's Fourier transform. And, conversely, the periodization of a function's Fourier transform is determined by discrete values of the original function. is statement becomes accurate with the Poisson summation formula (see [15] [17][18][19]).

Theorem 5 (Poisson summation formula). Assume that
for some C, ε > 0. en, both f and f are continuous and for all x ∈ R, we have , a ≠ 0, and using the simple fact that h(ω) � 1/af(ω/a), we immediately get the following corollary.

Corollary 1.
If f and f satisfy the conditions of the previous theorem, then for any a ≠ 0 and any x ∈ R, we have that Since the construction of PST is based on the periodization of a fundamental function after changing its scale, it is not surprising that PS g f is determined by discrete samples of the Fourier transforms of the original function f and the function g. More precisely, we prove the following result. Theorem 6. Assume that the fundamental function g satisfies the conditions: where C, ε are positive constants. Let also f ∈ L 1 (R) be a function such that en, the PST of f with respect to g can be written in the form: In particular, if Proof. Firstly, we observe that 6 Mathematical Problems in Engineering where C denotes a constant (not always the same). Furthermore, erefore, the assumptions on the functions f and g show that the convolution product f * g satisfies the condition of eorem 5 and Poisson's formula can be applied. Since for any (fixed) s > 0 the function g s satisfies similar assumptions to the function g, we can also apply Poisson's formula to the product f * g s . Now we recall the form of the PST which exploits the convolution of functions and is given in Proposition 3. Hence, we have Applying Poisson's formula to the function f * g s gives that and the result follows.

Limit Behaviour
is section is entirely devoted to the study of the behaviour of the PST for the limit values of the scale s and the period λ.

Theorem 7.
For any function f ∈ L 1 (R) and any s > 0, the following holds:

Lemma 2.
For any x ∈ R and any s > 0, we have Proof. Without loss of generality, we may assume that λ > 2|x|. en, using the decay condition (1), we can easily obtain and the result follows. □ Proof of eorem 7. By inequality (4), for 0 < s ≤ λ, the family of functions G(·; s, λ) is bounded by some absolute constant C. erefore, f(x)G(x; s, λ) is dominated by the integrable function Cf(x). e dominated convergence theorem and the previous lemma show now that Clearly, and the proof of the theorem is complete.

Theorem 8.
For any function f ∈ L 1 (R) and any λ > 0, the following holds: Lemma 3. For any x ∈ R and any λ > 0, we have Mathematical Problems in Engineering where Proof. Since G is periodic with period λ, it suffices to prove the lemma for any x ∈ [0, λ). We first consider the case x � 0. en, an easy calculation exploiting the decay condition (1) gives that Hence, it follows that lim s⟶0 + G(0; s, λ) � g(0). We next consider the case where x ∈ (0, λ). Since s ⟶ 0 + , we may assume, without loss of generality, that (52) erefore, lim s⟶0 + G(x; s, λ) � 0 and this completes the proof of the lemma. □ Proof of eorem 8. As in the Proof of eorem 7, the dominated convergence theorem and the previous lemma imply that Since G 0 is zero for almost every x ∈ R, we have and the desired result follows. It remains to investigate the behaviour of the PST when the period λ tends to zero. is means that the inverse of period, that is, the frequency, tends to infinity. e well-known Riemann-Lebesgue lemma asserts that the FT, f(ω), of a function f ∈ L 1 (R) converges to zero, when the frequency |ω| ⟶ ∞. Since the FT can be derived from the PST (see the next section), it is reasonable to expect that the PST of an L 1 function converges to 0 as λ goes to 0. is is true; however, we have to impose some further restrictions on the functions f and g. □ Theorem 9. Assume that the fundamental function g satisfies the conditions: (56) Let also f ∈ L 1 (R) be a function such that For example, f, g may be Schwartz functions. en, we have that (58) Proof. By eorem 6, the PST of f is given by discrete samples of the FTs of f and g, namely: By our assumptions on the FT of f, it follows that for any λ ∈ (0, 1) and any n ≠ 0: erefore, the series in (59) converges uniformly for 0 < s ≤ λ < 1 and thus However, for n � 0, we have and for n ≠ 0, It follows that e purpose of this section is to show that the FT can be derived as a special case of the PST.
Consider the fundamental function g(x) � exp(− 2π ix) � e − 2πix , 0 ≤ x < 1. en, g(x/s) � exp(− 2πix/s), for 0 ≤ x < s, and it is not hard to see that G(x; s, s) � exp(− 2πix/s) for any x ∈ R. Furthermore, Consequently, for any function f ∈ L 1 (R), we have that is, the diagonal of PST gives the FT of f for the positive frequencies.
In a similar way, we can produce the inverse Fourier transform. Indeed, we take the fundamental function en, G(x; s, s) � exp (2πix/s) for any x ∈ R and,

Connection with the WT. A wavelet is a function
We point out that the term wavelet does not possess a unique definition. However, the above definition is appropriate and for all practical purposes. Given the wavelet ψ and a function f ∈ L ∞ (R), the non-normalized continuous wavelet transform (WT) of f with respect to ψ is defined as Fundamental works on wavelet theory are [7,8,20,21]. If the fundamental function g of the PST is a wavelet, then equation (8) becomes (69) at is, the PST is an infinite sum of wavelet transforms of f.
Furthermore, it follows, by eorem 7, that erefore, if we replace the initial function g with its translation τ t g(x) � g(x − t), then the above observation implies that Remark 3. Summarizing the results of Sections 6.1 and 6.2, we can say that the PST unifies the two well-known and extensively studied transforms: the FT and the WT. Roughly speaking, these two transforms correspond to the two extreme cases of the PST. e case s � λ corresponds to the FT, while the case λ � ∞ gives the WT.

Connection with the Laplace Transform.
In view of eorem 7, we observe that the limit values of PST as λ tends to infinity can produce many other integral transforms, that is, transforms of the form f(x)K(x, s)dx which are defined by a kernel K(x, s). Take for example the (one-sided) Laplace transform defined by Using the fundamental function g(x) � e − x , if x > 0 and g(x) � 0 elsewhere, it is not hard to conclude that for any s > 0,

The PST of Some Selected Functions
e purpose of the present section is to illustrate the PST to some selected functions and real 1D surfaces.

e PST of the Characteristic Function of an Interval.
We start with the characteristic function of the interval [0, T], T > 0, which is defined by en, for 0 < s ≤ λ, formula (8) gives that In order to complete the above calculation, we will use several fundamental functions. Firstly, we consider the function g(x) � χ [0,1] (x). en, for 0 < s < λ, it follows that where for any real number x, [x] denotes the integer part of x.

Mathematical Problems in Engineering
Secondly, we consider the fundamental function g(x) � e − 2πix , for 0 ≤ x < 1. en, ∞ − ∞ g(x)dx � 0 and for s < λ, it follows that As it was expected, the above formula is close to the FT of f which is given by Finally, we take the fundamental function g(x) � sin x, for 0 ≤ x < 1. en, for s < λ, we have

e PSTof the Gaussian.
We continue with the Gaussian function with mean value μ and standard deviation σ 2 . We also take as a fundamental function the Gaussian with μ � 0 and σ 2 � 1. en, the PS transform of the function f is given by

e PST of a Rectangular Pulse Train.
Consider the rectangular pulse train function defined by its period P � 50 and width (duty cycle) S � 20 shown in Figure 2(a) and described by the following formula: e PST of this function is calculated by using the kernel G(x; s, λ) of a rectangular pulse train with varied period λ and feature width s (see equation (7)): e resulting two-dimensional PST PS g f(s, λ) is shown in Figure 2(b) for 1 < s < 50 and 5 < λ < 120. One can directly notice the presence of a strong peak at λ � P � 50, while less pronounced peaks also appear at the half and double periods, i.e., at λ � P/2 � 25 and λ � 2P � 100. Similar peaks are expected at higher multiples or lower submultiples of the fundamental period (λ � 150, 200, . . .) if these periods were included in the diagram of Figure 2(b). e ordinate value of all peaks is at s very close to 20 as also shown in the inset diagram where the 1D section of PST at λ � 50 (PS g f(s, λ) � 50) is illustrated. erefore, this simple example of periodic function reveals that PST is able to quantify both period and feature scale simultaneously.

e PST of a Gaussian Pulse Train.
A little more complicated example of 1D surface is the Gaussian pulse train function comprising Gaussian peaks with fixed width arranged at periodic positions as shown below: An illustration of such a 1D surface is shown in Figure 3(a) where the width of the repeated Gaussian is S � 10 while the period P � 80.
e PST is calculated using definition (7) with the kernel function G consisting of periodically arranged Gaussian peaks with varying period λ and scale s: (86) Figure 3(b) shows PS g f which is dominated by the presence of two peaks at s � 10 and λ � 80 and λ � 160, respectively, while PST is almost zero elsewhere. e first peak at (s � 10, λ � 80) corresponds to the period and scale of the input pulse train and the second is the first harmonic at λ � 2 · 80 � 160. e inset of Figure 3(b) magnifies the 1D behaviour of PST for λ � 80 versus scale s. e above simple deterministic examples show that (a) PST is a sufficient tool to reveal simultaneously both the period and scale of 1D surfaces and (b) PST is characterized by distinct peaks with asymmetrical structure since they are much steeper along period axis (horizontally) than scale axis (vertically).
One can readily observe the asymmetry in the steepness of the peaks: along period axis they are much steeper than scale axis.

e PSTof an Experimental 1D Surface.
In order to get a more realistic example of PST application, we consider a polymer surface treated with oxygen plasma for 1 minute inside a low temperature plasma reactor. e surface topography is measured by an atomic force microscope (AFM) in a grid of 512 × 512 points while the total measured area is 2 × 2 μm 2 . e experimental background and measurement details can be found elsewhere [22]. A randomly selected profile f(x) of this surface is displayed in Figure 4(a). Its morphology is roughly characterized by the presence of more or less well-defined peaks which are distributed along profile with an apparently random manner. However, the PSD analysis shown in Figure 4(b) reveals the existence of several periodicities manifested by the presence of peaks at specific frequencies. e smooth peak A with the highest frequency at 0.02nm − 1 corresponds to the smaller period of 50nm which roughly quantifies the mean distance between nearby peaks when all peaks (high and low) are considered. However, the peaks are diversified with respect to their scales (height and width) and the higher ones define new approximate periodicities with larger periods on the background of the fundamental period of 50nm. ese periods are manifested by the peaks B-E in PSD and reflect periodicities of approximately 80, 140, 155, and 225nm, respectively. erefore, the profile morphology is characterized by multiperiodicity induced by the multiscaling of peak dimensions and especially their heights. However, the multiscale nature of peaks cannot be fully quantified by the PSD analysis since the information about their widths is missing. On the other hand, Figure 4(c) shows the PST of the same profile using as kernel a Gaussian pulse train similar to those used in the Section 7.4. As expected, PST reproduces the periodicities detected by PSD at 50, 75, 140, 160, and 225nm, but at the same time, it proceeds one step further to characterize the scales (i.e., widths) of the peaks involved in each periodicity. us, for example, the fundamental period of 50nm is created by peaks of average width 8 − 10nm, while the largest periodic structure with period 225nm includes bolder protrusions of 35nm width each.
is example certifies the usefulness of PST analysis to disclose at the same diagram both hidden periodicities and feature scales in apparently random morphologies of profiles and surfaces.

Summary and Conclusions
Fourier and wavelet transforms play a prominent role in several applications, e.g., signal analysis, surface measurement, etc. FT is the standard tool for analyzing periodic surfaces as it provides us with accurate information about the harmonic components of the analyzed surface. On the other hand, WT can detect irregularities of the surface such as cracks or scratches (see [23] for a discussion of the advantages and disadvantages of the use of FT and WT in the analysis of a geometrical surface structure). It is noteworthy that several surfaces encountered in nature and technology are almost periodic and they can be viewed as generated by a feature repeated endlessly. Properties of such surfaces depend not only on period or irregularities but also on the width (or scale) of the repeated feature and, in many cases, on the ratio of period and scale. Despite the success of FT and WT in surface characterization, they are not able to quantify straightforwardly the feature scale and ratio.
Taking into account the above considerations, we have proposed and presented the period-scale transform. Its purpose is to detect simultaneously the harmonic components of the analyzed surface as well as the scale of the repeated feature. e PST is an integral transform, i.e., it is of the form where the kernel G(·; s, λ) comes from the periodization with period λ of a fundamental function g, with spatial scale s. e PST defines a bounded linear operator. As far as scaling properties, the transform behaves very well with respect to dilation ( eorem 1) but it is not phase invariant ( eorem 2). Several other properties of the PST have been detailed in this work. Furthermore, we have addressed the limit behaviour of the PST and we have shown that PST bridges the gap between FT and WT in the sense that these two classical transforms can be realized as special cases of PST.
Finally, we have examined the application of PST to model and real 1D surfaces. It was found that PST provides valuable information for both period and feature scale aspects of surface morphology. For example, in the case of the experimental profile, we compared the outputs of PST and FT and explicitly showed the additional information about the scales of repeated surface features provided by PST in excess of FT.
Future work will explore more applications of PST to real surfaces and focus on the relationship with specific surface functionalities (light reflectance, wetting, friction, and so on). Also, more attention will be paid on the phase dependence of PST and its impact on the final result.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.