## Line tension from dual-geometry sessile droplet measurements: Combining contact-angle size-dependence data for axisymmetric and cylindrical droplets to determine the line tension

Dmitry V. Tatyanenko ∗ and Konstantin D. Apitsin

Department of Statistical Physics, Faculty of Physics, Saint Petersburg State University 7-9 Universitetskaya nab., St. Petersburg, 199034, Russia

(Dated: March 31, 2025)

We perform a thermodynamic analysis of various contributions to the size dependence of the contact-angle cosine for both axisymmetric and cylindrical sessile droplets. This shows that a commonly used method to determine the line tension from the slope of the contact-angle cosine dependence on the three-phase contactline curvature (for axisymmetric droplets) provides a certain combination of the line tension, the adsorptions at the three interfaces, and the macroscopic contact angle. To extract the contribution related to the contact line in the leading order and determine the line tension, we propose a simple technique using droplet-size dependences of the contact angle for axisymmetric and cylindrical droplets under the same conditions.

Phys. Rev. E 111 , 035503 (2025). DOI: 10.1103/PhysRevE.111.035503

©2025 American Physical Society

## I. INTRODUCTION

The contact angle of a large sessile droplet can be found from the well-known classical Young equation [1]

$$\sigma ^ { \alpha \beta } \cos \theta _ { \infty } = \sigma ^ { \beta \gamma } - \sigma ^ { \alpha \gamma } \quad \quad ( 1 ) \quad \text{has}$$

with σ the surface tension of the corresponding interface and θ ∞ the macroscopic contact angle; the superscripts α β , , and γ mark the liquid, gaseous, and solid phases, respectively, and double Greek superscripts mark the corresponding interfaces. The surface tensions here are thermodynamic, defined as the surface excesses of the grand potential per unit area of the interface. For small droplets, the thermodynamic line tension κ , defined as the line excess of the grand potential per unit length of the three-phase contact line [2, 3], affects the contact angle [3-9]:

$$\sigma ^ { \alpha \beta } \cos \theta _ { a } = \sigma ^ { \beta \gamma } - \sigma ^ { \alpha \gamma } - \frac { \kappa } { r _ { a } }. \quad \quad ( 2 ) \quad \text{res} \\ \text{Ger}$$

Here, θ is the contact angle, r is the radius of the three-phase contact line, and the subscript 'a' refers to an axisymmetric droplet (see Fig. 1(a)); the line tension κ is assumed to be constant. Eq. (2) is often called the modified Young equation.

However, some authors have observed that Eq. (3) is generally not correct [21-24], since it assumes the corresponding surface tensions in Eqs. (1) and (2) to be equal, whereas they correspond to different equilibrium sessile droplet sizes, which implies different values of the chemical potential(s). It has been argued that the adsorption at the solid-liquid interface and its effect on the surface tension σ αγ affect cos θ a in the same manner and with a comparable magnitude to that observed in experiments [25, 26]; thus, the line tension itself is not necessary to explain the dependence of cos θ a on 1 / r a. As was demonstrated within an interface displacement model, the contributions of the line tension and the adsorption at the solid-gas interface (and its impact on the surface tension σ βγ ) to the dependence of cos θ a on 1 / r a can be comparable in the presence of a precursor film [21]. Some authors [18, 27] assume that only σ αβ depends on the droplet size and Eq. (3) is correct, although the slope of the graph cos θ a vs 1 / r a yields a quantity different from -κ / σ αβ . Some authors refer to the resulting 'line tension' as the apparent one [20, 22, 27-31]. General thermodynamic analysis demonstrates [21, 24] that the value of this slope is influenced by contributions from both the line tension and the adsorptions at all three interfaces (see also [32] for similar arguments).

By subtracting the modified Young equation (2) from the classical one (1) and dividing both parts of the resulting equation by σ αβ , one usually arrives at [3, 5, 7, 8]

$$\cos \theta _ { \infty } - \cos \theta _ { a } = \frac { \kappa } { \sigma ^ { \alpha \beta } r _ { a } }.$$

This equation is commonly employed for contact-angle-based measurements of the line tension in laboratory experiments (see, e. g., [10]) and computer modeling (e. g., molecular dynamics simulations [11-20]). According to Eq. (3), the slope of the graph cos θ a vs 1 / r a must be equal to -κ / σ αβ . The surface tension σ αβ is typically measured/calculated separately.

∗ d.tatyanenko@spbu.ru

How does one 'extract' the line-tension contribution and determine the value of the line tension from the measure-

Figure 1. Axisymmetric (a) and cylindrical (b) sessile droplets. R is the curvature radius of the liquid-gas interface, θ is the contact angle, r is the radius (half-width) of the droplet on the substrate. Subscript 'a' marks quantities for axisymmetric, 'c' - for cylindrical droplets. Phases: α -liquid, β -gas (vapor), γ -solid substrate.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00547v2-figure-64369425-8806-4cde-b41a-6fe73367de7c.jpg)

ments or calculations of cos θ vs 1 / r ? The thermodynamic analysis [24] demonstrates that the equation corresponding to Eq. (2) for cylindrical droplets lacks a κ / r term (such droplets are long liquid 'channels', 'filaments', or 'ridges' attached to the substrate surface, as illustrated in Fig. 1(b); corresponding quantities will be marked with the subscript 'c'). This means that dependence of the contact angle θ c on the size of a cylindrical droplet (e. g., its half-width r c on the substrate) is governed mainly by dependences of the surface tensions on the droplet size via the chemical potential(s). This suggests an idea to compare or combine data on the size dependence of cos θ for axisymmetric and cylindrical droplets under the same conditions to exclude the surface-tension-related contributions and determine the line tension.

In Sec. II we present a thermodynamic analysis of various contributions to droplet-size dependence of the contact angle for both axisymmetric and cylindrical droplets in a system with a single-component fluid and an ideal, smooth, homogeneous solid substrate in the absence of external fields. Based on this analysis, in Sec. III we propose a simple technique that allows for the determination of the line tension from a combination of contact-angle size dependences for axisymmetric and cylindrical droplets, similarly to the use of Eq. (3).

Cylindrical droplets are often used in modeling, e. g., in molecular dynamics (MD) simulations [14, 15, 18, 20, 22, 3339] due to the convenient use of periodic boundary conditions along their longitudinal axes and anticipated weak dependence of the contact angle on the droplet width. However, they are hardly available for experimental measurements due to Plateau-Rayleigh-type instability [40]. Therefore, the technique we propose is primarily oriented towards modeling, especially with methods that do not permit direct calculation of the line excess of the grand potential (e. g., MD simulations). The dimensions of such simulated axisymmetric and cylindrical droplets, i. e., R r , , and θ , can be 'measured', and these data can be used to determine the line tension.

## II. THERMODYNAMICS AND CONTACT ANGLES

## A. Equilibrium conditions for sessile droplets

In order to describe the thermodynamic properties of a small sessile droplet on a partially wettable substrate, we consider an open system (that can exchange both energy and matter with its surroundings) of a fixed volume including the droplet itself, the substrate, and the surrounding gas (vapor). The liquid and gaseous phases are assumed to constitute a single-component fluid at given values of chemical potential µ and temperature T . These conditions correspond to the grand canonical statistical ensemble. The grand thermodynamic potential of the system can be decomposed into bulk, surface, and line parts:

$$\Omega & = - \, p ^ { \alpha } V ^ { \alpha } - p ^ { \beta } V ^ { \beta } + \omega ^ { \gamma } V ^ { \gamma } \\ & \quad + \sigma ^ { \alpha \beta } A ^ { \alpha \beta } + \sigma ^ { \alpha \gamma } A ^ { \alpha \gamma } + \sigma ^ { \beta \gamma } A ^ { \beta \gamma } + \kappa L \quad \quad \quad ( 4 ) \quad \quad \quad$$

with V the volume, A the area of the interface, L the length of the three-phase contact line, p the pressure in the bulk phase, and ω the bulk density of the grand thermodynamic potential (equals -p in fluid phases). Other quantities as well as notations of phases and interfaces are explained above.

Let us consider two types of droplet geometry: (1) an axisymmetric droplet, shaped as a spherical segment in our 'macroscopic' thermodynamic description, see Fig. 1(a), and (2) a fixed-length piece of a cylindrical droplet, shaped as a segment of a circular cylinder in our 'macroscopic' thermodynamic description, see Fig. 1(b). In both geometries, the cross-section of the droplet in our 'macroscopic' thermodynamic description is a circular segment, thus,

$$r = R \sin \theta.$$

Therefore, the shape and size of the droplet can be described by two independent variables, e. g., ( R r , ) , ( R , θ ) , or ( r , θ ) .

In both geometries, the equilibrium conditions for sessile droplets can be obtained by expressing all the volumes, areas, and the length of the three-phase contact line in (4), and solving one of the following sets of two equations:

$$\inf \quad \partial \Omega ( R, r, T, \mu ) / \partial R = 0 \ \text{ and } \partial \Omega ( R, r, T, \mu ) / \partial r = 0, \ ( 6 a )$$

$$| 3 - \quad \partial \Omega ( R, \theta, T, \mu ) / \partial R = 0 \text{ and } \partial \Omega ( R, \theta, T, \mu ) / \partial \theta = 0, \, ( 6 b )$$

$$\text{di-} _ {. } \quad \partial \Omega ( r, \theta, T, \mu ) / \partial r = 0 \ \text{ and } \partial \Omega ( r, \theta, T, \mu ) / \partial \theta = 0, \ ( 6 c )$$

$$( 6 a ) \\ ( 6 b ) \\ ( 6 c )$$

and then taking into account the geometric relation (5). Since T and µ are constant, the pressures and the surface tensions in (4) will be constant at such variations (for σ αβ , this also requires ∂σ αβ / ∂ R = 0, which is true for the surface of tension chosen as the αβ dividing surface).

The droplet-size-dependent surface tensions σ αβ a and σ αβ c for spherical and cylindrical interfaces may differ, although their difference is expected to be very small at the same value of µ (see Sec. II B).

Let us chose R and θ as independent variables to describe the sessile droplet. Introducing nonvariable quantities for the total volume occupied by the fluid V t ≡ V α + V β and the total area of the substrate A t ≡ A αγ + A βγ and using relation (5), one can express the volumes, areas, and the contact-line length for an axisymmetric droplet as

$$\inf _ { a } \frac { \pi R _ { a } ^ { 3 } } { 3 } ( 2 + \cos \theta _ { a } ) ( 1 - \cos \theta _ { a } ) ^ { 2 }, \quad L _ { a } = 2 \pi R _ { a } \sin \theta _ { a }, \\ \inf _ { \real } \frac { A _ { a } ^ { \alpha \beta } = 2 \pi R _ { a } ^ { 2 } ( 1 - \cos \theta _ { a } ), \quad A _ { a } ^ { \alpha \gamma } = \pi R _ { a } ^ { 2 } \sin ^ { 2 } \theta _ { a }, \\ \intertext { v }. \frac { V _ { a } ^ { \beta } = V _ { t } - V _ { a } ^ { \alpha }, \quad A _ { a } ^ { \beta \gamma } = A _ { t } - A _ { a } ^ { \alpha \gamma }. } { \circ a }$$

Introducing also ∆ σ γ ≡ σ αγ -σ βγ , we can now rewrite decomposition (4) of the grand potential as

$$\mathfrak { c } \, \text{p} \text{-} \\ \mathfrak { rface}, \quad & \Omega _ { \mathfrak { a } } ( R _ { \mathfrak { a } }, \theta _ { \mathfrak { a } } ) = - \left ( p ^ { \alpha } - p ^ { \beta } \right ) \frac { \pi R _ { \mathfrak { a } } ^ { 3 } } { 3 } ( 2 + \cos \theta _ { \mathfrak { a } } ) ( 1 - \cos \theta _ { \mathfrak { a } } ) ^ { 2 } \\ & + 2 \pi \sigma _ { \mathfrak { a } } ^ { \alpha \beta } R _ { \mathfrak { a } } ^ { 2 } ( 1 - \cos \theta _ { \mathfrak { a } } ) + \Delta \sigma ^ { \gamma } \pi R _ { \mathfrak { a } } ^ { 2 } \sin ^ { 2 } \theta _ { \mathfrak { a } } + 2 \pi \kappa _ { \mathfrak { a } } R _ { \mathfrak { a } } \sin \theta _ { \mathfrak { a } } \\ ( 4 ) & - p ^ { \beta } V _ { \mathfrak { t } } - \sigma ^ { \beta \gamma } A _ { \mathfrak { t } } + \omega ^ { \gamma } V ^ { \gamma }$$

We assume the line tension κ = κ ( T , µ , r ) in both geometries. This gives, with use of Eq. (5),

$$\left ( \frac { \partial \kappa } { \partial R } \right ) _ { \theta, T, \mu } = \left ( \frac { \partial \kappa } { \partial r } \right ) _ { T, \mu } \left ( \frac { \partial r } { \partial R } \right ) _ { \theta } = \left ( \frac { \partial \kappa } { \partial r } \right ) _ { T, \mu } \sin \theta,$$

$$\left ( \frac { \partial \kappa } { \partial \theta } \right ) _ { R, T, \mu } = \left ( \frac { \partial \kappa } { \partial r } \right ) _ { T, \mu } \left ( \frac { \partial r } { \partial \theta } \right ) _ { R } = \left ( \frac { \partial \kappa } { \partial r } \right ) _ { T, \mu } R \cos \theta.$$

Application of the equilibrium conditions (6b) then yields

$$- \pi R _ { a } ^ { 2 } ( p ^ { \alpha } - p ^ { \beta } ) ( 2 + \cos \theta _ { a } ) ( 1 - \cos \theta _ { a } ) ^ { 2 }$$

$$+ 4 \pi R _ { a } \sigma _ { a } ^ { \alpha \beta } ( 1 - \cos \theta _ { a } ) + 2 \pi R _ { a } \Delta \sigma ^ { \gamma } \sin ^ { 2 } \theta _ { a }$$

$$+ 2 \pi \kappa _ { a } \sin \theta _ { a } + 2 \pi R _ { a } \frac { \partial \kappa _ { a } } { \partial r _ { a } } \sin ^ { 2 } \theta _ { a } = 0,$$

(

-

π

R

3

a

p

-

p

)

sin

θ

a

+

2

πσ αβ

a

R

α

β

3

2

a

sin

θ

a

+

2

πκ

a

R

a

cos

θ

a

$$+ 2 \pi R _ { a } ^ { 2 } \Delta \sigma ^ { \gamma } \sin \theta _ { a } \cos \theta _ { a } + 2 \pi R _ { a } ^ { 2 } \frac { \partial \kappa _ { a } } { \partial r _ { a } } \sin \theta _ { a } \cos \theta _ { a } = 0. \quad \begin{pmatrix} \gamma & 1 \\ \text{its } c \\ \text{the } \colon \end{pmatrix}$$

In order to exclude the pressure difference p α -p β , we multiply the first equation by ( 2 π R a ) -1 sin 2 θ a ( 1 -cos θ a ) -2 and then subtract the second equation, multiplied by ( 2 + cos θ a )( 2 π R 2 a sin θ a ) -1 [41]. After some trigonometric simplifications and the substitution of the result into one of the aforementioned conditions, we arrive at the Laplace equation and the generalized Young equation as equilibrium conditions for an axisymmetric sessile droplet [3, 7, 8, 21, 24, 42]:

$$p ^ { \alpha } - p ^ { \beta } = \frac { 2 \sigma _ { \mathbf a } ^ { \alpha \beta } } { R _ { \mathbf a } }, \quad \quad \ ( 7 ) \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
 \quad p ^ { \alpha } - p ^ { \beta } = \frac { 2 \theta } { R _ { \mathbf a } }$$

$$\sigma _ { a } ^ { \alpha \beta } \cos \theta _ { a } = - \Delta \sigma ^ { \gamma } - \frac { \kappa _ { a } } { r _ { a } } - \frac { \partial \kappa _ { a } } { \partial r _ { a } } \quad \quad ( 8 ) \quad \ e q \colon \\ \quad \ \text{the}$$

For a cylindrical droplet,

$$V _ { c } ^ { \alpha } = \frac { L _ { c } R _ { c } ^ { 2 } } { 4 } ( 2 \theta _ { c } - \sin 2 \theta _ { c } ),$$

$$V _ { c } ^ { \alpha } = \frac { L _ { c } R _ { c } ^ { 2 } } { 4 } ( 2 \theta _ { c } - \sin 2 \theta _ { c } ), \\ A _ { c } ^ { \alpha \beta } = L _ { c } R _ { c } \theta _ { c }, \quad A _ { c } ^ { \alpha \gamma } = L _ { c } R _ { c } \sin \theta _ { c }, \\ V _ { c } ^ { \beta } = V _ { t } - V _ { c } ^ { \alpha }, \quad A _ { c } ^ { \beta \gamma } = A _ { t } - A _ { c } ^ { \alpha \gamma }, \\ \text{the name to the } T \quad \text{of the three above one at line in} }$$

where the length L c of the three-phase contact line is independent of R c and θ c (see Fig. 1(b)) and is assumed to be constant. The contact angle θ c here is expressed in radians. Decomposition (4) of the grand potential then takes the form

$$\Omega _ { c } & ( R _ { c }, \theta _ { c } ) = - \left ( p ^ { \alpha } - p ^ { \beta } \right ) \frac { L _ { c } R _ { c } ^ { 2 } } { 4 } ( 2 \theta _ { c } - \sin 2 \theta _ { c } ) + \sigma _ { c } ^ { \alpha \beta } L _ { c } R _ { c } \theta _ { c } \quad \text{ asions} \\ & + \Delta \sigma ^ { \gamma } L _ { c } R _ { c } \sin \theta _ { c } + \kappa _ { c } L _ { c } - p ^ { \beta } V _ { t } - \sigma ^ { \beta \gamma } A _ { t } + \omega ^ { \gamma } V ^ { \gamma }. \\ & \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \. \quad \.$$

Application of the equilibrium conditions (6b) yields

$$- \left ( p ^ { \alpha } - p ^ { \beta } \right ) L _ { \text{c} } R _ { \text{c} } ( \theta _ { \text{c} } - \sin \theta _ { \text{c} } \cos \theta _ { \text{c} } ) + \sigma _ { \text{c} } ^ { \alpha \beta } L _ { \text{c} } \theta _ { \text{c} }$$

$$+ \Delta \sigma ^ { \gamma } L _ { \text{c} } \sin \theta _ { \text{c} } + L _ { \text{c} } \frac { \partial \kappa _ { \text{c} } } { \partial r _ { \text{c} } } \sin \theta _ { \text{c} } = 0,$$

$$- \left ( p ^ { \alpha } - p ^ { \beta } \right ) L _ { c } R _ { c } ^ { 2 } \sin ^ { 2 } \theta _ { c } + \sigma _ { c } ^ { \alpha \beta } L _ { c } R _ { c } \quad \begin{pmatrix} & \dots \\ & \text{spo} \\ & + \Delta \sigma ^ { \gamma } L _ { c } R _ { c } \cos \theta _ { c } + L _ { c } R _ { c } \frac { \partial \kappa _ { c } } { \partial r _ { c } } \cos \theta _ { c } = 0. \end{pmatrix} \quad \begin{pmatrix} & \dots \\ & \text{$\text{$\text{$\text{$\text{$\text{$delta}$} } } } \\ & & \text{$\text{$\text{$delta}$} } } \end{pmatrix}$$

Here we have employed trigonometric identities sin2 θ c = 2sin θ c cos θ c and 1 -cos2 θ c = 2sin 2 θ c. To derive a more concise equation for p α -p β , we multiply the first equation by R c cos θ c and subtract the second equation, multiplied by sin θ c . After some simplifications and the substitution of the result into one of the aforementioned conditions, we arrive at the Laplace equation and the generalized Young equation as equilibrium conditions for a cylindrical sessile droplet, slightly different from Eqs. (7) and (8):

$$p ^ { \alpha } - p ^ { \beta } = \frac { \sigma _ { c } ^ { \alpha \beta } } { R _ { c } },$$

$$\sigma _ { c } ^ { \alpha \beta } \cos \theta _ { c } = - \Delta \sigma ^ { \gamma } - \frac { \partial \kappa _ { c } } { \partial \eta _ { c } }. \quad \quad ( 1 0 )$$

The absence of the coefficient 2 in the Laplace equation (9) for the cylindrical droplet is explained by the fact that its αβ interface is curved along a single direction, whereas the spherical droplet has equal curvatures along two directions. The line-tension terms in the generalized Young equations (8) and (10) originate from variations of the line term κ L in (4). For an axisymmetric droplet, L a = 2 π r a and ∂ κ ( a L a ) / ∂ r a = 2 πκ a + 2 π r a ∂κ a / ∂ r a yield the two terms κ a / r a + ∂κ a / ∂ r a in Eq. (8). For a cylindrical droplet, L c is fixed, and ∂ κ ( c L c ) / ∂ r c = L c ∂κ c / ∂ r c yields the only term ∂κ c / ∂ r c in Eq. (10).

For both geometries, we suppose the line tension to depend not only on T and µ , but also on the droplet size. This dependence is typically overlooked, and in this case, Eq. (8) is transformed into the modified Young equation (2), whereas Eq. (10) is transformed into the formally classical Young equation (1) (with surface tensions, however, depending on the chemical potential, i. e., on the equilibrium droplet size, which makes the contact angle θ c different from θ ∞ , as discussed in Sec. II B). We suppose the contact-line radius/halfwidth r to be the only appropriate size variable the line tension can be considered to depend on, as detailed in appendix A. The elimination of the size-dependent term ∂κ / ∂ r from the generalized Young equations (8) and (10) renders them inexact, if κ is assumed to be dependent on the droplet size. It can be seen in a model system if the line tension is calculated directly [21]. Furthermore, the generalized Young equation is the most straightforward approach to determining ∂κ / ∂ r , provided that θ and r are found and the surface and line tensions (the latter for axisymmetric droplet only) are given or calculated directly [21].

## B. Size dependence of contact angle and relations between pressures, surface tensions, droplet size, and chemical potential

Acomparison of the generalized Young equation (8) or (10) with the classical Young equation (1) must take into account the fact that these two equations are written for droplets corresponding to different values of the chemical potential µ . The classical Young equation is applicable to a macroscopically large (strictly speaking, infinite) drop, i. e., the liquid-gas bulk phase coexistence at the binodal. Let us mark the values of the

quantities corresponding to the binodal with the subscript ' ∞ '. The surface tensions at the binodal then should be denoted as σ αβ ∞ , σ αγ ∞ , σ βγ ∞ , and thus, ∆ σ γ ∞ ≡ σ αγ ∞ -σ βγ ∞ . Consequently, the classical Young equation (1) can be rewritten as

$$\sigma _ { \infty } ^ { \alpha \beta } \cos \theta _ { \infty } = - \Delta \sigma _ { \infty } ^ { \gamma }. \quad \quad ( 1 1 ) \quad \text{arr}$$

By subtracting the generalized Young equations (8) and (10) from Eq. (11), we obtain two similar equations: one for an axisymmetric droplet

$$\sigma _ { \infty } ^ { \alpha \beta } \cos \theta _ { \infty } - \sigma _ { a } ^ { \alpha \beta } \cos \theta _ { a } = \delta \Delta \sigma ^ { \gamma } + \frac { \kappa _ { a } } { r _ { a } } + \frac { \partial \kappa _ { a } } { \partial r } \quad ( 1 2 ) \quad \text{$this q$} \\ \dots \quad. \quad. \quad. \quad. \quad.$$

and the other for a cylindrical droplet

$$\sigma _ { \infty } ^ { \alpha \beta } \cos \theta _ { \infty } - \sigma _ { c } ^ { \alpha \beta } \cos \theta _ { c } = \delta \Delta \sigma ^ { \gamma } + \frac { \partial \kappa _ { c } } { \partial r _ { c } }. \quad ( 1 3 ) \quad \text{phas}$$

Here δ ∆ σ γ ≡ ∆ σ γ -∆ σ γ ∞ . To make these equations more similar to Eq. (3), let us introduce δσ αβ a,c ≡ σ αβ a,c -σ αβ ∞ . Then, Eqs. (12) and (13) can be rewritten as

$$\sigma _ { \infty } ^ { \alpha \beta } ( \cos \theta _ { \infty } - \cos \theta _ { a } ) = \delta \Delta \sigma ^ { \gamma } + \delta \sigma _ { a } ^ { \alpha \beta } \cos \theta _ { a } + \frac { \kappa _ { a } } { r _ { a } } + \frac { \partial \kappa _ { a } } { \partial r _ { a } }, \, ( 1 4 ) \quad \stackrel { \dots } { \left ( n ^ { \beta } \ll } } { \text{liquid} \right ) }$$

$$\sigma _ { \infty } ^ { \alpha \beta } ( \cos \theta _ { \infty } - \cos \theta _ { c } ) = \delta \Delta \sigma ^ { \gamma } + \delta \sigma _ { c } ^ { \alpha \beta } \cos \theta _ { c } + \frac { \partial \kappa _ { c } } { \partial r _ { c } }. \quad ( 1 5 ) \quad \text{$\text{$one$}$}$$

It can be observed that Eq. (14), the exact counterpart of Eq. (3) for the difference cos θ ∞ -cos θ , contains four terms for an axisymmetric droplet (instead of one) and three terms for a cylindrical droplet (instead of none). The question thus arises as to how these terms depend on the droplet size (the radius/half-width of the droplet on the substrate).

For the αβ liquid-gas interface, the isothermal Gibbs adsorption equations can be written as

$$d \sigma _ { a, c } ^ { \alpha \beta } = - \Gamma _ { a, c } ^ { \alpha \beta } d \mu, \quad \quad \quad ( 1 6 ) \quad \quad \cdots$$

where Γ αβ is the adsorption (the surface excess of matter per unit surface area) at the αβ interface for corresponding geometry. As the surface of tension is chosen as the dividing surface in both geometries, these equations do not contain ( ∂σ αβ / ∂ R ) d R terms.

For each of the solid-fluid interfaces, similar equations can be written. In general, these equations are more complicated (see, e. g., Sec. 6.3 in [8] and Sec. 3.4.3 in [3]). However, in certain simple cases (e. g., for nondeformable substrate and such a choice of the solid-fluid dividing surfaces that the amounts of matter of the immobile component of solid do not change in the solid phase at any change of µ ), they take [21, 24] the same simple form as (16)

$$d \Delta \sigma ^ { \gamma } \equiv d \left ( \sigma ^ { \alpha \gamma } - \sigma ^ { \beta \gamma } \right ) = \left ( \Gamma ^ { \beta \gamma } - \Gamma ^ { \alpha \gamma } \right ) d \mu. \quad ( 1 7 )$$

By integrating Eqs. (16) and (17) from µ ∞ to µ , we obtain, introducing δµ ≡ µ -µ ∞ ,

$$\delta \sigma _ { a, c } ^ { \alpha \beta } = - \int _ { \mu _ { \infty } } ^ { \mu } \Gamma _ { a, c } ^ { \alpha \beta } ( \mu ^ { \prime } ) \, d \mu ^ { \prime } \simeq - \Gamma _ { \infty } ^ { \alpha \beta } \delta \mu, \quad ( 1 8 ) \quad \text{ this,}$$

$$\delta \Delta \sigma ^ { \gamma } = \int _ { \mu _ { \infty } } ^ { \mu } ( \Gamma ^ { \beta \gamma } - \Gamma ^ { \alpha \gamma } ) \, \mathrm d \mu ^ { \prime } \simeq ( \Gamma _ { \infty } ^ { \beta \gamma } - \Gamma _ { \infty } ^ { \alpha \gamma } ) \, \delta \mu. \quad ( 1 9 )$$

̸

Since Γ αβ a ( µ ∞ ) = Γ αβ c ( µ ∞ ) = Γ∞ αβ , the asymptotic expression in (18) is identical for spherical and cylindrical surfaces at the same value of µ , but the exact integral is generally not the same if Γ αβ a ( µ ) = Γ αβ c ( µ ) . Thus, σ αβ a and σ αβ c generally differ, but at least in the second order in δµ ,

$$\sigma _ { \mathbf a } ^ { \alpha \beta } - \sigma _ { \mathbf c } ^ { \alpha \beta } = \delta \sigma _ { \mathbf a } ^ { \alpha \beta } - \delta \sigma _ { \mathbf c } ^ { \alpha \beta } = O ( ( \delta \mu ) ^ { 2 } ). \quad ( 2 0 )$$

The last asymptotic expressions in Eqs. (18) and (19) are valid for small values of δµ . What is the relationship between this quantity and the droplet size? This can be estimated using the Laplace equation (7) or (9) for the axisymmetric or cylindrical droplet, respectively. By applying the isothermal Gibbs-Duhem equations d p = n d µ to the liquid and gaseous phases, we obtain

$$p ^ { \alpha } - p ^ { \beta } \simeq n _ { \infty } ^ { \alpha } \left ( 1 + \chi _ { \infty } ^ { \alpha } \delta p ^ { \alpha } \right ) \delta \mu \simeq n _ { \infty } ^ { \alpha } \delta \mu + O \left ( ( \delta \mu ) ^ { 2 } \right ) \ \ ( 2 1 )$$

with n α the number density of molecules and χ α the isothermal compressibility in the liquid phase, δ p α ≡ p α -p α ∞ . Here we have neglected the number density of molecules in the gas ( n β ≪ n α ). We will further neglect the compressibility of the liquid since the resulting correction is only of the second order in δµ . Using the Laplace equations (7) and (9), along with the geometric relation (5), we obtain

$$\frac { 2 \sigma _ { a } ^ { \alpha \beta } } { R _ { a } } = \frac { 2 \sigma _ { a } ^ { \alpha \beta } \sin \theta _ { a } } { r _ { a } } \simeq n _ { \infty } ^ { \alpha } \delta \mu, \text{ \quad \ \ } ( 2 2 )$$

$$\frac { \sigma _ { c } ^ { \alpha \beta } } { R _ { c } } = \frac { \sigma _ { c } ^ { \alpha \beta } \sin \theta _ { c } } { r _ { c } } \simeq n _ { \infty } ^ { \alpha } \delta \mu. \text{ \quad \ \ } ( 2 3 )$$

It can now be seen that δµ ∝ 1 / R = O 1 ( / r ) , i. e., small values of δµ correspond to large droplets, and δµ → 0 expectedly corresponds to R → ∞ and r → ∞ .

Accordingly, the correction terms in Eqs. (14) and (15) are given by

$$\Gamma \text{per} \\ \text{1}, \text{ge} _ { \text{using} } \quad \delta \Delta \sigma ^ { \gamma } \simeq \frac { 2 \sigma _ { a } ^ { \alpha \beta } \sin \theta _ { a } } { n _ { \infty } ^ { \alpha } r _ { a } } \left ( \Gamma _ { \infty } ^ { \beta \gamma } - \Gamma _ { \infty } ^ { \alpha \gamma } \right ) = O ( \delta \mu ) = O ( 1 / r ), \ ( 2 4 )$$

$$\delta \sigma _ { \mathfrak { a } } ^ { \alpha \beta } \simeq - \frac { 2 \sigma _ { \mathfrak { a } } ^ { \alpha \beta } \Gamma _ { \infty } ^ { \alpha \beta } \sin \theta _ { \mathfrak { a } } } { n _ { \infty } ^ { \alpha } r _ { \mathfrak { a } } } = O ( \delta \mu ) = O ( 1 / r ), \quad ( 2 5 )$$

$$\kappa _ { \mathfrak { a } } / r _ { \mathfrak { a } } \simeq \kappa _ { \infty } / r _ { \mathfrak { a } } = O ( \delta \mu ) = O ( 1 / r ), \quad \ \ ( 2 6 )$$

where κ ∞ = κ a at µ = µ ∞ ( r a → ∞ ) is the line tension of the macroscopic (infinite) droplet, and

$$\text{id} \text{they} \quad \delta \Delta \sigma ^ { \gamma } \simeq \frac { \sigma _ { c } ^ { \alpha \beta } \sin \theta _ { c } } { n _ { \infty } ^ { \alpha } r _ { c } } \left ( \Gamma _ { \infty } ^ { \beta \gamma } - \Gamma _ { \infty } ^ { \alpha \gamma } \right ) = O ( \delta \mu ) = O ( 1 / r ), \ \ ( 2 7 )$$

$$\mathcal { \delta } _ { \mathbf c } \sigma _ { \mathbf c } ^ { \alpha \beta } \simeq - \frac { \sigma _ { \mathbf c } ^ { \alpha \beta } \Gamma _ { \infty } ^ { \alpha \beta } \sin \theta _ { \mathbf c } } { n _ { \infty } ^ { \alpha } r _ { \mathbf c } } = O ( \delta \mu ) = O ( 1 / r ). \quad ( 2 8 )$$

With regard to the terms ∂κ / ∂ r , it can be shown [3, 21, 24, 42] that they are O ( ( δµ ) 2 ) = O 1 ( / r 2 ) . To demonstrate this, the line adsorption equation [43] (a line counterpart of the Gibbs adsorption equation and the Gibbs-Duhem equation) can be considered [3, 24, 42] at a constant temperature:

$$d \kappa = - \Lambda d \mu + ( \partial \kappa / \partial r ) _ { \mu, T } d r$$

with Λ the line excess of the fluid matter (the number of molecules of the fluid) per unit length of the three-phase contact line (a line counterpart of the adsorption Γ ) in the corresponding geometry. Assuming then

$$\kappa = \kappa _ { \infty } + O ( \delta \mu ), \quad \Lambda = \Lambda _ { \infty } + O ( \delta \mu ) \quad ( 2 9 ) \quad \text{mad}$$

along the equilibrium curve r ( µ ) with the surface of tension chosen as the liquid-gas dividing surface at each value of µ , we arrive at the estimation

$$\left ( \frac { \partial \kappa } { \partial r } \right ) _ { \mu, T } = \frac { d \kappa } { d r } + \Lambda \frac { d \mu } { d r } = \left [ \frac { d \kappa } { d \mu } + \Lambda \right ] \left ( \frac { d r } { d \mu } \right ) ^ { - 1 } = O ( ( \delta \mu ) ^ { 2 } ). \ ( 3 0 ) \quad \begin{array} { c } \text{as $nan$} \\ \text{a film} \\ \text{if the $n$}. \end{array}.$$

The factor in square brackets is O 1 ( ) since (29), while d r / d µ = d ( R sin θ ) / d µ ≃ ( d R / d µ ) sin θ ∞ = O ( ( δµ ) -2 ) . This estimation is supported, in the case of ∂κ a / ∂ r a, by direct calculations within an interface displacement model [21].

This consideration is somewhat simplified; it does not take into account the influence of the stress tensor distribution in the solid substrate on the line tension. We assume that the estimation ∂κ / ∂ r = O ( ( δµ ) 2 ) = O 1 ( / r 2 ) remains valid if the corresponding contribution is taken into account.

For sufficiently large droplets, sin θ can be approximated with sin θ ∞ , with a relative error O ( δµ ) = O 1 ( / r ) . This makes Eqs. (24), (25), (27), and (28) valid with sin θ a,c replaced by sin θ ∞ . The same approximation can be applied to cos θ in the correction terms δσ αβ cos θ in Eqs. (14) and (15).

In light of these considerations, we can now proceed to derive the linear approximations in 1 / r for Eqs. (14) and (15):

$$\cos \theta _ { \infty } - \cos \theta _ { a } \simeq \frac { \varkappa _ { a } } { \sigma _ { \infty } ^ { \alpha \beta } r _ { a } }, \quad \quad ( 3 1 ) \quad \inf _ { \text{for} } ^ { 1 }$$

$$\cos \theta _ { \infty } - \cos \theta _ { c } \simeq \frac { \varkappa _ { c } } { \sigma _ { \infty } ^ { \alpha \beta } r _ { c } }, \quad \quad ( 3 2 ) \quad \text{neg} \\ \det _ { \dots }$$

where

$$\varkappa _ { a } \equiv \kappa _ { \infty } + \frac { 2 \sigma _ { \infty } ^ { \alpha \beta } \sin \theta _ { \infty } } { n _ { \infty } ^ { \alpha } } \left ( \Gamma _ { \infty } ^ { \beta \gamma } - \Gamma _ { \infty } ^ { \alpha \gamma } - \Gamma _ { \infty } ^ { \alpha \beta } \cos \theta _ { \infty } \right ) \quad ( 3 3 ) \quad \text{angle,} \\ \text{subseq} \quad \text{It is}$$

and

$$\varkappa _ { c } \equiv \frac { \sigma _ { \infty } ^ { \alpha \beta } \sin \theta _ { \infty } } { n _ { \infty } ^ { \alpha } } \left ( \Gamma _ { \infty } ^ { \beta \gamma } - \Gamma _ { \infty } ^ { \alpha \gamma } - \Gamma _ { \infty } ^ { \alpha \beta } \cos \theta _ { \infty } \right ) \quad \quad ( 3 4 ) \quad \begin{array} { c } \text{into $\alpha$} \\ \sigma _ { \infty } ^ { \alpha \beta } \end{array} \right )$$

are the quantities sometimes referred to as the apparent line tensions [20, 22, 27-31] of axisymmetric and cylindrical droplets, respectively (see below).

Equations (31) and (33) demonstrate that the dependence of cos θ a on 1 / r a, even in the first order of approximation, is determined not only by the line tension, but also by the adsorptions at the three interfaces. Indeed, Eq. (31) should be treated as Eq. (3) typically is. This implies that the slope of the graph cos θ a vs 1 / r a at 1 / r a → 0 does not directly yield the line tension; rather, it represents a combination of the line tension, the surface tension of the droplet free surface, adsorptions at the three interfaces, and the macroscopic contact angle. Some authors use the term 'apparent line tension' to refer to the calculated value of the line tension from the cos θ ( 1 / r ) dependence at 1 / r → 0, as if Eq. (3) were the correct linear approximation. Then Eq. (33) yields the value of the apparent line tension of axisymmetric droplets.

In previous studies, the impact of individual adsorptions on the contact angle has been considered and estimated.

It is likely that the most substantial contribution can be made by the solid-gas adsorption Γ βγ for a droplet formed on the top of a precursor film. The formation of such a film can, e. g., precede dropwise condensation on a partially wettable substrate [44]. Such precursor films can be as thick as nanometers (e. g., 3-11 nm for alcohols and water on a glass [45] and 11-12 nm for water on fused quartz [46] for a film in equilibrium with saturated vapors). In such a case, if the contact angles of sessile droplets are measured at the substrate surface (not at the top of the precursor film), this film should be considered as an equilibrium state of the solid substrate and treated as an adsorption film resulting in a high adsorption Γ∞ βγ ≈ f ∞ n α ∞ (with f ∞ the precursor film thickness at equilibrium with an infinitely large droplet). Therefore, its contribution to the apparent line tension κ a can be estimated as 2 σ αβ ∞ f ∞ sin θ a ≲ 10 -10 -10 -9 N.

The line tension and the solid-gas adsorption contribution have been directly calculated within an interface displacement model with a model short-range interface potential [21]. The results indicated that they can be comparable in magnitude, with the solid-gas adsorption contribution being even several times more in magnitude than the line tension κ ∞ . It should be noted that other adsorptions are zero within this model.

The impact of the adsorption Γ αγ at the solid-liquid interface on the contact angle was examined by Ward and Wu [25] in the context of a discussion of cos θ a ( 1 / r a ) measurements for dodecane sessile droplets on silanized silicon wafers [47]. The adsorption Γ βγ at the solid-gas interface was deemed negligible compared to Γ αγ due to the low volatility of dodecane [25]. The calculations demonstrated that the effect of the adsorption at the solid-liquid interface was sufficient to explain the observed droplet-size-dependence of the contact angle, even without the line tension. These arguments were subsequently developed in [26].

It is quite common to consider solely the size dependence of the liquid-gas surface tension, which corresponds to taking into account the adsorption Γ∞ αβ . The Tolman equation σ αβ ≃ σ αβ ∞ ( 1 -2 δ T / R a ) with the Tolman length [48] δ T ≈ Γ∞ αβ / n α ∞ is the most commonly employed and provides the contribution to the apparent line tension κ a equal to -2 σ αβ ∞ δ T sin θ a cos θ a. With typical estimates of δ T ∼ 10 -10 mand σ αβ ∞ ∼ 10 -2 N/m, this contribution can be estimated as 10 -12 N, which is smaller or comparable with the values of the apparent line tension κ a observed in simulations [22] and some experiments [47].

Of greater interest is the fact that Eqs. (32) and (34) demonstrate that, for cylindrical droplets, cos θ c also depends on 1 / r c in the first order in 1 / r c. This has been observed in simulations [20, 22] and attributed to the size dependence of the surface tension σ αβ expressed using the Tolman length [48] δ T ≈ Γ∞ αβ / n α ∞ . As it is seen from Eq. (34), the dependence of other surface tensions on the droplet size via the chemical potential can result in their contributions to the apparent line tension of cylindrical droplets.

## III. LINE TENSION FROM DUAL-GEOMETRY CONTACT-ANGLE DEPENDENCES ON DROPLET SIZE

One can easily notice that the apparent line tension (34) of cylindrical droplets, which is proportional to the slope of the cos θ c ( 1 / r c ) graph at 1 / r c → 0, is exactly half of the second, adsorption-related term in the apparent line tension (33) of axisymmetric droplets, which is proportional to the slope of the cos θ a ( 1 / r a ) graph at 1 / r a → 0. Thus, by combining data for the dependences cos θ a ( 1 / r a ) and cos θ c ( 1 / r c ) , one can determine the line tension of a macroscopic droplet as

$$\kappa _ { \infty } = \varkappa _ { a } - 2 \varkappa _ { c }. \text{ \quad \quad \ \ } ( 3 5 )$$

By employing this relation, we have calculated the line tension κ ∞ of a macroscopic droplet from the apparent line tensions κ a and κ c of sessile droplets simulated by Kanduˇ c et al. [20, 22] using data from Fig. 2(c) in [22]. As illustrated in Fig. 2, the calculated values of κ ∞ using Eq. (35) are in good agreement with those obtained in [22] by taking into account the Tolman-length and 'line-tension-stiffness' ∂κ / ∂θ corrections. In [20, 22], the solid-liquid and solidgas surface tensions are assumed to be constant. However, the equilibrium equations for the contact angles θ a,c derived there are different from those presented here. The discrepancies between these works and the present study in the initial assumptions and the resulting equations will be discussed in Sec. IV.

Nevertheless, the Tolman-length correction can be accounted for through the use of Eq. (33). If we consider solely the correction to σ αβ ∞ = 55 mN/m with the Tolman length δ T ≈ Γ∞ αβ / n α ∞ = -0 05 nm as determined from simulations [20], this . leads to σ αβ ∞ δ T = -2 75 pN, and calculate . κ ∞ using Eq. (33) with Γ∞ αγ = Γ∞ βγ = 0, the resulting value of κ ∞ will be quite close to κ a (shown with hollow points and a dashed line in Fig. 2). Therefore, it can be concluded that the Tolman-length correction alone underestimates the adsorption-related contributions to the apparent line tension of axisymmetric droplets.

It is clearly seen that the values of κ a and κ c are comparable in these simulations, which makes κ a visibly different from κ ∞ as calculated in [22] and with the use of Eq. (35).

Equations (31)-(34) suggest another idea. In the context of dual-geometry sessile droplet measurements, there is no need to employ the classical Young equation at all. By subtracting the generalized Young equation for axisymmetric droplets (8) from the one for cylindrical droplets (10), we arrive at the following equation:

$$\sigma _ { c } ^ { \alpha \beta } \cos \theta _ { c } - \sigma _ { a } ^ { \alpha \beta } \cos \theta _ { a } = \frac { \kappa _ { a } } { r _ { a } } + \frac { \partial \kappa _ { a } } { \partial r _ { a } } - \frac { \partial \kappa _ { c } } { \partial r _ { c } }. \quad ( 3 6 ) \quad \stackrel { \dots } { \text{of $1/}$} }$$

In consideration of Eqs. (25) and (28) with estimation (30), the first-order approximation to the combined equation (36) can be expressed as follows:

$$\cos \theta _ { c } - \cos \theta _ { a } \simeq \frac { \kappa _ { \infty } } { \sigma _ { \infty } ^ { \alpha \beta } r _ { a } }.$$

The dependence of κ a on µ and r a, as well as the dependence of σ αβ on µ , are not taken into account since they result in corrections of order O ( 1 / r 2 a ) .

Figure 2. Apparent line tensions of axisymmetric ( κ a ) and cylindrical ( κ c ) MD-simulated droplets [20, 22] (data taken from Fig. 2(c) in [22]) and corresponding line tension κ ∞ of macroscopic droplets calculated with Eq. (35) and by the authors of [22] using additional corrections (data taken from Fig. 3 therein) vs the macroscopic contactangle cosine (tuned by surface polarity). The hollow points and the dashed line show estimations of κ ∞ from the values of κ a and the Tolman-length correction with Eq. (33), corresponding to the αβ (liquid-gas) adsorption correction. An error for κ ∞ from Eq. (35) is estimated as the sum of an error for κ a and a doubled error for κ c . The lines are weighted cubic splines, plotted using gnuplot's [49] smooth acsplines method with the weights w / ∆ , where ∆ is the absolute error of the corresponding value, and w = 500 pN.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00547v2-figure-fbaa669c-219a-4b65-9d08-132ad356473a.jpg)

In order to exclude the terms affected by changes in surface tension, it is necessary to take all the quantities in Eqs. (36) and (37) at the same values of the chemical potential and temperature. The chemical potential is often unknown in simulations or measurements. However, as can be seen from Eqs. (22) and (23), and taking into account Eq. (20), the same value of δµ corresponds to the relation R c ≃ R a / 2 or, in the same order of δµ , to r c ≃ r a / 2.

Thus, to determine the line tension from the measured (in the experiment or modeling) dependences of the cos θ on 1 / r for both axisymmetric and cylindrical droplets under the same conditions, it is sufficient to calculate the difference

$$\underset { \text{ts} } { \text{ts} } ( 8 ) \quad \cos \theta _ { c } \Big | _ { R _ { c } = R _ { a } / 2 } - \cos \theta _ { a } \Big | _ { R _ { a } } \text{ or } \ \cos \theta _ { c } \Big | _ { r _ { c } = r _ { a } / 2 } - \cos \theta _ { a } \Big | _ { r _ { a } } \ ( 3 8 )$$

and then take a linear fit of this quantity vs 1 / r a at small values of 1 / r a. The slope of the resulting graph will yield the value of κ ∞ / σ αβ ∞ :

$$\cos \theta _ { c } \Big | _ { R _ { c } = R _ { a } / 2 } - \cos \theta _ { a } \Big | _ { R _ { a } } \simeq \frac { \kappa _ { \infty } } { \sigma _ { \infty } ^ { \alpha \beta } r _ { a } },$$

$$\cos \theta _ { c } \Big | _ { r _ { c } = r _ { a } / 2 } - \cos \theta _ { a } \Big | _ { r _ { a } } \simeq \frac { \kappa _ { \infty } } { \sigma _ { \infty } ^ { \alpha \beta } r _ { a } }. \quad \quad ( 4 0 )$$

An additional advantage of this technique is that the graph of the difference cos θ c -cos θ a vs 1 / r a must start from the

known point (0 at 1 / r a = 0). Therefore, the slope of the linear graph remains the only fitting parameter. When plotting cos θ a ( 1 / r a ) and cos θ c ( 1 / r c ) independently, each of them must be approximated by a straight line starting from an unknown point cos θ ∞ ( 0 . ) This point must be the same for both axisymmetric and cylindrical droplets. In principle, this should be taken into account to reduce the number of fitting parameters from 4 to 3. In contrast to experimental measurements, the computational expense of simulating droplet sizes with macroscopic contact angles can be significant. However, the macroscopic contact angle θ ∞ can be found, e. g., by employing the classical Young equation (11) and calculating σ αβ ∞ and ∆ σ γ ∞ in three simulations of planar interfaces [38].

The use of dual-geometry techniques to determine the line tension may offer another potential advantage, as suggested by experiments with nanosized droplets. For nanosized droplets, the graph cos θ a vs 1 / r a can be fairly linear; however, such a linear fit applicable to small droplets may yield a false value of cos θ ∞ at extrapolation down to 1 / r a → 0 (see Fig. 5 in [50]) and the slope that does not correspond to asymptotic Eq. (31). In many simulations, the droplets are also nanosized, and thus the same issue may arise with simulated droplets, both axisymmetric and cylindrical, resulting in miscalculation of κ a and κ c. The application of Eq. (35) would ultimately yield an incorrect value of κ ∞ . However, not only would the slope be incorrect in such a situation, but the extrapolated values of cos θ ∞ would also be different for cos θ a vs 1 / r a and cos θ c vs 1 / r c, indicating an inconsistency. The technique we proposed here would reveal such a situation as a set of points on the graph of quantity (38) vs 1 / r a that do not fit a straight line originating at the point ( 0 0 . Both these indicators imply , ) that the droplets are not sufficiently large for the asymptotic equations, Eqs. (31), (32), and (37), to be applicable, thereby precluding the determination of the contact angle θ ∞ and the line tension κ ∞ of macroscopic droplets.

## IV. OVERVIEW OF RESULTS AND DISCUSSION

The thermodynamic analysis performed demonstrates that droplet-size corrections to the contact-angle cosine for axisymmetric sessile droplets typically attributed to line tension alone, are generally determined by line tension κ a, its dependence on droplet size (directly on the droplet radius r a on the surface and via the chemical potential), and also by dependences of all three surface tensions on droplet size (via the chemical potential). The same applies to the corrections to the contact-angle cosine for cylindrical sessile droplets, except for the line-tension term κ c / r c that is absent in the generalized Young equation (10) for cylindrical droplets.

In the first order in 1 / r a (31), this gives the slope of cos θ a vs 1 / r a to be proportional not to the line tension κ ∞ of macroscopic droplets (as it is widely believed according to Eq. (3)) but rather to the so-called 'apparent line tension' κ a, a sum of this line tension κ ∞ and a combination of adsorptions at the three interfaces, the number density of molecules in the liquid phase, and the macroscopic contact angle (33).

For cylindrical droplets, the contact angle is found to de- pend on the droplet half-width r c on the substrate already in the first order in 1 / r c (32), thereby yielding the apparent line tension κ c. The latter is determined by the same combination of the adsorptions, the number density of molecules in the liquid phase, and the macroscopic contact angle (34) but lacks the numerical coefficient 2 compared to the axisymmetric case (33). This gives a possibility for determining the line tension κ ∞ from κ a and κ c found in dual-geometry measurements (for droplets simulated at the same conditions, i. e., under the same temperature with the same interaction potentials, etc.) as (35). This was illustrated with calculations of κ ∞ from κ a and κ c of MD-simulated droplets [20, 22] (see Fig. 2).

We have also proposed a technique for dual-geometry measurements based on the calculation of the slope of the difference cos θ c -cos θ a at the same value of the chemical potential vs 1 / r a. It can be implemented as a calculation of this difference (38) at R c = R a / 2 or r c = r a / 2. The slope will equal κ ∞ / σ αβ ∞ , with σ αβ ∞ the surface tension of the planar liquidgas interface. The advantages of this technique are discussed in Sec. III.

The Laplace and the generalized Young equations in both geometries, as well as the Gibbs adsorption equation (16) for the liquid-gas interface, were obtained and considered for the surface of tension as the liquid-gas dividing surface. At an arbitrary choice of this dividing surface, all these equations become more complicated [42, 51]. Even though most researchers use these equations in their simple forms, the real choice of the dividing surface is often different or even not specified explicitly. In simulations, this can be done, e. g., by employing a boundary detection algorithm, and different algorithms can yield [52] visibly different dependences of cos θ a vs 1 / r a. As demonstrated in [39], the dependence of cos θ c on R c is visibly different for the surface of tension and the equimolecular liquid-gas dividing surface in the same simulated argon cylindrical nanodroplets. Thus, the choice of the liquid-gas dividing surface can be considered as an additional source of errors that need to be estimated.

As for the solid-fluid dividing surfaces, the analog of the Gibbs adsorption equation (17) has been employed to express the difference ∆ σ γ ≡ σ αγ -σ βγ , implying certain conditions that are not necessarily met at an arbitrary choice of these dividing surfaces. However, if the same choice is made in both geometries, the quantities δ ∆ σ γ ≡ ∆ σ γ -∆ σ γ ∞ will cancel out, both with use of Eq. (35) and in the proposed technique where Eq. (36) will remain valid. In contrast, the line tensions κ a and κ c (and their macroscopic value κ ∞ ) do generally depend on the choice of the solid-fluid dividing surfaces. For any specific choice, the dual-geometry approach must work, giving the value of κ ∞ corresponding to the chosen dividing surfaces.

The idea of employing dual-geometry measurements to determine the line tension has in fact been exploited by Kanduˇ c et al. [22]. However, there are substantial differences from our findings in details and results. In that work, the line tension is supposed to depend on the contact angle θ but not the radius/half-width r of the droplet on the substrate. We consider the contact angle θ to be an inappropriate variable to characterize the size dependence of the line tension (see appendix A for further details). Another difference is that only

the liquid-gas surface tension σ αβ is considered to depend on the droplet size (specifically, R ) in [20, 22], while, in our consideration, all three surface tensions are generally dependent on the droplet size. The Tolman length δ T ≈ Γ∞ αβ / n α ∞ was obtained independently in [20] through modeling of a free cylindrical droplet. In our results, all necessary information to account for the effects of adsorptions (the droplet-size dependence of the surface tensions) is contained in the measured dependence cos θ c ( 1 / r c ) .

The equilibrium equations for the contact angles derived in [20, 22] are also different from ours, even at constant σ αγ , σ βγ , and κ . This discrepancy can be attributed primarily to the omission of the bulk term in the free-energy-based work of the droplet formation, for which then a stationary point was found under condition of constant droplet volume [20, 22]. The dual-geometry consideration allowed the authors to find the 'line-tension-stiffness' ∂κ / ∂θ correction, estimated as O 1 ( / r ) , from the Tolman length and the apparent line tension of cylindrical droplets [20]. This correction was subsequently adopted with axisymmetric droplets [22], implicitly assuming ∂κ a / ∂θ a = ∂κ c / ∂θ c . In the present work, we do not assume ∂κ a / ∂ r a = ∂κ c / ∂ r c; however, these derivatives are shown (30) to be O ( 1 / r 2 ) , i. e., negligible.

Finally, let us discuss if the proposed technique can be easily adapted to determine the size-dependent line tension κ a and/or κ c (not only their limit κ ∞ at r → ∞ ). Equation (36) does not contain solid-fluid surface tensions, but it contains the surface tensions σ αβ a,c . They can be calculated separately for free spherical/cylindrical droplets of the same radius R a,c of the liquid-gas surface as for the sessile droplet of the same type of geometry (thus, at the same value of chemical potential). However, since κ a ( µ , r a ) and κ c ( µ , r c ) are different functions, the expression ∂κ a / ∂ r a -∂κ c / ∂ r c on the right-hand side of Eq. (36) is still O ( ( δµ ) 2 ) = O 1 ( / r 2 a ) and therefore cannot be neglected even in the second order in 1 / r a. Thus, employing equation (36) is not a viable method for obtaining more detailed information on the line tension from dual-geometry measurements of cos θ ( 1 / r ) .

## ACKNOWLEDGMENTS

The authors thank N. A. Volkov and M.S. Polovinkin for valuable discussions and comments on the manuscript. The research has been carried out with financial support of the Russian Science Foundation, grant No. 22-13-00151, https://rscf.ru/en/project/22-13-00151/

## Appendix A: Droplet size variable(s) the line tension depends on

In Sec. II A, the equilibrium conditions (7) and (8) (for axisymmetric droplets), and (9) and (10) (for cylindrical droplets) were derived, supposing the line tension κ a,c to depend on the droplet size via its radius/half-width r a,c on the substrate. However, in the way these conditions were derived, this choice of the size variable(s) for the line tension is neither unique nor obvious.

What if we suppose the line tension to depend on the radius R of the liquid-gas interface and/or the contact angle θ ? If the liquid-gas interface is still assumed to be spherical/cylindrical, derivation considered in Sec. II A should be modified. For κ a,c = κ a,c ( R a,c , θ a,c ) , the equilibrium conditions (6b) applied to the functions Ω a,c ( R a,c , θ a,c ) yield, for axisymmetric droplets,

$$\underset { \text{erived} } { \text{assured} } \quad p ^ { \alpha } - p ^ { \beta } = \frac { 2 \sigma _ { a } ^ { \alpha \beta } } { R _ { a } } + \left [ \frac { \partial \kappa _ { a } } { \partial \theta _ { a } } \frac { \sin \theta _ { a } } { R _ { a } } - \frac { \partial \kappa _ { a } } { \partial R _ { a } } \cos \theta _ { a } \right ] \frac { 2 \sin \theta _ { a } } { ( 1 - \cos \theta _ { a } ) ^ { 2 } R _ { a } },$$

$$\overset { \bullet } { \text{irly to} } \quad, \quad \sigma _ { a } ^ { \alpha \beta } \cos \theta _ { a } = - \Delta \sigma ^ { \gamma } - \frac { \kappa _ { a } } { r _ { a } } + \frac { \partial \kappa _ { a } } { \partial \theta _ { a } } \frac { 2 + \cos \theta _ { a } } { R _ { a } } - \frac { \partial \kappa _ { a } } { \partial R _ { a } } \frac { ( 1 + \cos \theta ) ^ { 2 } } { \sin \theta _ { a } },$$

$$& \cdot + \left \lfloor \frac { \partial \kappa _ { a } } { \partial \theta _ { a } } \frac { \sin \theta _ { a } } { R _ { a } } - \frac { \partial \kappa _ { a } } { \partial R _ { a } } \cos \theta _ { a } \right \rfloor \frac { 2 \sin \theta _ { a } } { ( 1 - \cos \theta _ { a } ) ^ { 2 } R _ { a } }, \\ & \quad \omega \sigma ^ { \gamma } - \frac { \kappa _ { a } } { r _ { a } } + \frac { \partial \kappa _ { a } } { \partial \theta _ { a } } \frac { 2 + \cos \theta _ { a } } { R _ { a } } - \frac { \partial \kappa _ { a } } { \partial R _ { a } } \frac { ( 1 + \cos \theta ) ^ { 2 } } { \sin \theta _ { a } }, \\ & \cdot \cdot \cdot \cdot \cdot$$

and for cylindrical droplets,

$$p ^ { \alpha } - p ^ { \beta } = \frac { \sigma _ { c } ^ { \alpha \beta } } { R _ { c } } + \frac { \frac { \partial \kappa _ { c } } { \partial \theta _ { c } } \frac { \sin \theta _ { c } } { R _ { c } } - \frac { \partial \kappa _ { c } } { \partial R _ { c } } \cos \theta _ { c } } { ( \sin \theta _ { c } - \theta _ { c } \cos \theta _ { c } ) R _ { c } },$$

$$\underset { \underset { \underset { \underset { \text{acc}} { c } } } { \ u m e } } { \ u m e } \, \quad \sigma _ { c } ^ { \alpha \beta } \cos \theta _ { c } = - \Delta \sigma ^ { \gamma } + \frac { \frac { \partial \kappa _ { c } } { \partial \theta _ { c } } \frac { \theta _ { c } - \sin \theta _ { c } \cos \theta _ { c } } { R _ { c } } - \frac { \partial \kappa _ { c } } { \partial R _ { c } } \sin ^ { 2 } \theta _ { c } } { \sin \theta _ { c } - \theta _ { c } \cos \theta _ { c } }.$$

As can be seen, the equilibrium conditions for θ (the generalized Young equations) and for R (the Laplace equations) are subject to change and mix. The latter will only regain their original form at ( ∂κ / ∂ R ) cos θ = R -1 ( ∂κ / ∂θ ) sin θ , a condition that corresponds to the case κ = κ ( r ) . In all other cases, the effective force that tends to change R and/or θ due to dependence κ ( R , θ ) modifies the Laplace equation in both geometries. This nonlocal effect of the line tension has no reasonable explanation. The liquid-gas interface is not something rigid, and a force, tending to change it locally at the contact line, cannot affect the whole interface. The above equations are thus consequences of the assumption of a spherical/cylindrical shape for the interface.

Let us consider a variational derivation of the equilibrium conditions [7, 8, 21, 24], wherein the droplet is not initially assumed to be a spherical/cylindrical segment, and then discuss, within this more general context, why r is the only dropletsize variable on which the line tension can be assumed to depend.

We assume that the axisymmetric droplet has an arbitrary radial profile, m a ( x ) , where x is the radial coordinate in cylindrical coordinates. Similarly, we assume that the cylindrical droplet (in this case, a segment of a curvilinear cylinder) has an arbitrary transversal profile, m c ( x ) , where x is the transversal Cartesian coordinate (i. e., the coordinate along the substrate surface and perpendicular to the droplet's longitudinal axis). Refer to Figs. 3(a) and 3(b) for illustration. The profiles m a,c ( x ) are assumed to be twice differentiable. The grand potential Ω a,c of the system in both geometries can be expressed as a functional of m a,c ( x ) :

$$\text{to de-} \\ \text{in the} \\ \text{:rived}, \quad & \Omega _ { a, c } [ m _ { a, c } ( x ) ] = \int _ { 0 } ^ { r _ { a, c } } \omega _ { a, c } ( x, m _ { a, c } ( x ), m _ { a, c } ^ { \prime } ( x ) ) \, d x \\ \text{either} \quad & \quad + L _ { a, c } \kappa _ { a, c } \left ( r _ { a, c }, m _ { a, c } ^ { \prime } ( r _ { a, c } ) \right ) + \Omega ^ { \gamma } - p ^ { \beta } V _ { t } + \sigma ^ { \beta \gamma } A _ { t }$$

Figure 3. Axisymmetric (a) and cylindrical (b) shapes of sessile droplets in variational consideration. (c) A local alteration in the [axisymmetric or cylindrical] profile of a droplet near its three-phase contact line, which can significantly change the contact-angle-dependent line tension, making an arbitrarily small contribution to the integral part of the functional (A2) of the variation δ Ω . The altered profile with the slope m r ′ ( ) = -tan θ s (where θ s : ∂κ / ∂θ θ | = θ s = 0) corresponds to an extremum of the nonintegral line term, which makes the δ m r ′ ( ) term equal to zero, while the δ r term remains zero since the position of the contact line does not alter.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00547v2-figure-34b3e791-3ffc-4158-886c-11384d2019e7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00547v2-figure-098e228a-5898-4dc0-a372-7febec78ea3a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00547v2-figure-9e1f6999-e034-4e44-a095-c05c52734f23.jpg)

with with δ ( L a,c κ a,c ) for the axisymmetric/cylindrical profile as

$$\omega _ { a, c } \equiv \begin{cases} 2 \pi x \left [ \Delta \sigma ^ { \gamma } + \sigma _ { a } ^ { \alpha \beta } \sqrt { 1 + ( m _ { a } ^ { \prime } ( x ) ) ^ { 2 } } - ( p ^ { \alpha } - p ^ { \beta } ) m _ { a } ( x ) \right ], \\ L _ { c } \left [ \Delta \sigma ^ { \gamma } + \sigma _ { c } ^ { \alpha \beta } \sqrt { 1 + ( m _ { c } ^ { \prime } ( x ) ) ^ { 2 } } - ( p ^ { \alpha } - p ^ { \beta } ) m _ { c } ( x ) \right ], \quad & \delta ( L _ { a, c } \end{cases}$$

where the symbol ' ' indicates a derivative with respect to the ′ argument of the function, V t = V α a,c + V β a,c = const, A t = A αγ a,c + A βγ a,c = const. Here we take into account that

$$V _ { a } ^ { \alpha } = 2 \pi \int _ { 0 } ^ { r _ { a } } x m _ { a } ( x ) \, \mathrm d x, \ \ A _ { a } ^ { \alpha \beta } = 2 \pi \int _ { 0 } ^ { r _ { a } } x \, \sqrt { 1 + ( m _ { c } ^ { \prime } ( x ) ) ^ { 2 } } \, \mathrm d x, \quad \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

$$V _ { \mathrm a } ^ { \alpha } & = 2 \pi \int _ { 0 } ^ { r _ { \mathrm a } } x m _ { \mathrm a } ( x ) \, d x, \quad A _ { \mathrm a } ^ { \alpha \beta } = 2 \pi \int _ { 0 } ^ { r _ { \mathrm a } } x \, \sqrt { 1 + ( m _ { \mathrm c } ^ { \prime } ( x ) ) ^ { 2 } } \, d x, \quad \text{sidered} \\ V _ { \mathrm c } ^ { \alpha } & = L _ { \mathrm c } \int _ { 0 } ^ { r _ { \mathrm c } } m _ { \mathrm c } ( x ) \, d x, \quad A _ { \mathrm c } ^ { \alpha \beta } = L _ { \mathrm c } \int _ { 0 } ^ { r _ { \mathrm c } } \sqrt { 1 + ( m _ { \mathrm c } ^ { \prime } ( x ) ) ^ { 2 } } \, d x, \quad \text{point, } \cdot \\ A _ { \mathrm a } ^ { \alpha \gamma } & = \pi r _ { \mathrm a } ^ { 2 } = 2 \pi \int _ { 0 } ^ { r _ { \mathrm a } } x \, d x, \quad A _ { \mathrm c } ^ { \alpha \gamma } = L _ { \mathrm c } r _ { \mathrm c } = L _ { \mathrm c } \int _ { 0 } ^ { r _ { \mathrm c } } d x, \quad L _ { \mathrm a } = 2 \pi r _ { \mathrm a }, \quad \text{endpoi} \\ \dots \, \mathrm c \, \mathrm c \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \, \mathrm c \, \dots \}$$

and L c does not depend on r c.

In the functional (A1), not only the droplet profile m a,c ( x ) is variable, but also the position r a,c of the right endpoint of the interval, as well as the nonintegral term L a,c κ a,c with a variable L a = 2 π r a for an axisymmetric droplet and a nonvariable L c for a cylindrical droplet. In the context of our analysis, we suppose the line tension κ a,c to depend both on the radius/halfwidth r a,c of the droplet on the substrate and the contact angle θ a,c , i. e., the local slope m ′ a,c of the profile at x = r a,c , since m ′ a,c ( r a,c ) = -tan θ a,c .

To find the equilibrium equations, we need to write down the first variation of the functional (A1) containing both integral and nonintegral contributions, with the right endpoint that can freely move along the x axis. The boundary conditions for m a,c ( x ) are m ′ a,c ( 0 ) = 0 and m a,c ( r a,c ) = 0. Taking into account these boundary conditions, the first variation of the functional (A1) can be expressed as [53]

$$& \delta \Omega _ { a, c } \simeq \delta ( L _ { a, c } \kappa _ { a, c } ) + \int _ { 0 } ^ { r _ { a, c } } \left ( \frac { \partial \omega _ { a, c } } { \partial m _ { a, c } } \delta m + \frac { \partial \omega _ { a, c } } { \partial m ^ { \prime } _ { a, c } } \delta m ^ { \prime } \right ) \mathrm d x \\ & + \omega _ { a, c } \left ( r _ { a, c }, m _ { a, c } ( r _ { a, c } ), m ^ { \prime } _ { a, c } ( r _ { a, c } ) \right ) \delta r _ { a, c } = \delta ( L _ { a, c } \kappa _ { a, c } ) \quad \quad ( \mathrm A 2 )$$

$$+ \int _ { 0 } ^ { r _ { a, c } } \left ( \frac { \partial \omega _ { a, c } } { \partial m _ { a, c } } - \frac { d } { d x } \frac { \partial \omega _ { a, c } } { \partial m _ { a, c } ^ { \prime } } \right ) \delta m \, d x + \left ( \omega _ { a, c } - m _ { a, c } ^ { \prime } \frac { \partial \omega _ { a, c } } { \partial m _ { a, c } ^ { \prime } } \right ) \Big | _ { x = r _ { a, c } } \frac { \sqrt { ( \text{sphe} } ) } { \text{fixed} }.$$

$$& + \omega _ { a, c } \left ( r _ { a, c }, m _ { a, c } ( r _ { a, c } ), m _ { a, c } ^ { \prime } ( r _ { a, c } ) \right ) \delta r _ { a, c } = \delta ( L _ { a, c } \kappa _ { a, c } ) \quad ( A 2 ) \quad \text{Cons} \quad \text{and} \, E c \\ & + \int _ { 0 } ^ { r _ { a, c } } \left ( \frac { \partial \omega _ { a, c } } { \partial m _ { a, c } } - \frac { d } { d x } \frac { \partial \omega _ { a, c } } { \partial m _ { a, c } ^ { \prime } } \right ) \delta m \, d x + \left ( \omega _ { a, c } - m _ { a, c } ^ { \prime } \frac { \partial \omega _ { a, c } } { \partial m _ { a, c } ^ { \prime } } \right ) \right | \, \delta r _ { a, c } \quad \text{(spher} \quad \text{II)} \,,$$

$$\mathfrak { a } ( x ) \Big { ] }, \quad \ \ \ w i n { \nu } { \mathfrak { L } _ { a, c } } { \mathfrak { a } _ { a, c } } / \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { m } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { [ } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { p } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n  { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak{ n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak  { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { \n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { t } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak {  \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } _ { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n. } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { s } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { c } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { n } { \mathfrak { \mathfrak { n } {$$

Here ∂κ a,c / ∂ m ′ a,c ( r a,c ) = ( ∂κ a,c / ∂θ a,c )( d m ′ a,c ( r a,c ) / d θ a,c ) -1 = -( ∂κ a,c / ∂θ a,c ) cos 2 θ a,c is used.

The variation of the droplet profile, δ m x ( ) , can also be considered at a fixed position of the endpoint r , i. e., at δ r = 0. It is less clear that the variation of the profile slope at the endpoint, δ m r ′ ( ) , should also be considered as independent. A slight alteration in the profile m x ( ) within the vicinity of the endpoint can significantly modify δ m r ′ ( ) , while maintaining a negligible change in the integral term in (A2) and keeping δ r = 0 by fixing the position of the contact line, as depicted in Fig. 3(c). Thus, the stationary conditions for the functional (A1), δ Ω = 0, decompose into three equations [7, 8]:

I) The Euler-Lagrange equation ∂ω ∂ / m = ( d ∂ω ∂ / m ′ ) / d : x

$$\left ( \frac { \mathrm d } { \mathrm d x } + \frac { 1 } { x } \right ) \frac { m _ { \mathrm a } ^ { \prime } } { \sqrt { 1 + ( m _ { \mathrm a } ^ { \prime } ) ^ { 2 } } } = - \frac { p ^ { \alpha } - p ^ { \beta } } { \sigma _ { \mathrm a } ^ { \alpha \beta } }, \quad \ ( A 3 )$$

$$\frac { \mathrm d } { \mathrm d x } \frac { m _ { c } ^ { \prime } } { \sqrt { 1 + ( m _ { c } ^ { \prime } ) ^ { 2 } } } = - \frac { p ^ { \alpha } - p ^ { \beta } } { \sigma _ { c } ^ { \alpha \beta } }. \quad \quad ( \mathrm A 4 )$$

The left-hand side of each equation represents the sum of the principal curvatures of the droplet's free surface, which is constant and equal to -( p α -p β ) / σ αβ a,c . Since m c ( x ) is the transversal profile and the principal curvature along the longitudinal coordinate is zero, the sum of the principal curvatures in Eq. (A4) equals -1 / R c. Therefore, this equation is the Laplace equation (9) for a circular-cylindrical droplet. For an axisymmetric droplet, the condition of the constant sum of the principal curvatures, in conjunction with the boundary condition m ′ a ( 0 ) = 0, indicates that the surface is spherical. Consequently, the sum of its principal curvatures is -2 / R a, and Eq. (A3) is the Laplace equation (7) for the axisymmetric (spherical) droplet.

II) The transversality condition (taking into account the freedom of the endpoint to move along the x axis according

to the boundary condition m r ( ) = 0) written as a condition of the zero coefficient at δ r in the expression (A2):

$$\Delta \sigma ^ { \gamma } + \sigma _ { a } ^ { \alpha \beta } \sqrt { 1 + ( m _ { a } ^ { \prime } ( r _ { a } ) ) ^ { 2 } } + \frac { 1 } { r _ { a } } \frac { \partial ( \kappa _ { a } r _ { a } ) } { \partial r _ { a } } = \frac { \sigma _ { a } ^ { \alpha \beta } ( m _ { a } ^ { \prime } ( r _ { a } ) ) ^ { 2 } } { \sqrt { 1 + ( m _ { a } ^ { \prime } ( r _ { a } ) ) ^ { 2 } } }, \quad \underset { \text{size} } { \text{size} } }$$

$$\Delta \sigma ^ { \gamma } + \sigma _ { c } ^ { \alpha \beta } \sqrt { 1 + ( m _ { c } ^ { \prime } ( r _ { c } ) ) ^ { 2 } } + \frac { \partial \kappa _ { c } } { \partial r _ { c } } = \frac { \sigma _ { c } ^ { \alpha \beta } ( m _ { c } ^ { \prime } ( r _ { c } ) ) ^ { 2 } } { \sqrt { 1 + ( m _ { c } ^ { \prime } ( r _ { c } ) ) ^ { 2 } } }.$$

$$& \frac { ^ { 3 } ( m _ { a } ^ { \prime } ( r _ { a } ) ) ^ { 2 } } { \overline { + } ( m _ { a } ^ { \prime } ( r _ { a } ) ) ^ { 2 } }, \\ & \frac { \iota _ { c } ^ { \prime } ( r _ { c } ) ) ^ { 2 } } { m _ { c } ^ { \prime } ( r _ { c } ) ) ^ { 2 } }. \\ & \quad. \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \$$

Since we have already derived the Laplace equations (7) and (9) meaning that the droplets' free surfaces are spherical and cylindrical segments, respectively, m ′ a,c ( r a,c ) = -tan θ a,c ; thus, it can be demonstrated that the equations above yield the generalized Young equations (8) and (10).

III) The condition of the zero coefficient at δ m r ′ ( ) yields ∂κ / ∂ m r ′ ( ) = 0, i. e., ∂κ / ∂θ = 0. Given that the preceding two equations determine a spherical/cylindrical segment and the values of its two parameters, this additional equation makes the entire set of equations generally inconsistent.

Let us discuss again the question of the independence of the variation δ m r ′ ( ) from the variation of the integral in (A2). This presumed independence led to conditions (I) and (III) in their particular form and made them independent (in contrast to the approach used before, including the beginning of this section). However, this assumption apparently contradicts the 'regularity' of the function m x ′ ( ) at x = r . The dependence of the line tension κ on the contact angle θ leads to the common functional Ω [ m x ( )] being dependent on the derivative of the function m x ( ) at a single point x = r (i. e., the functional has a singular part similar to the derivative of the Dirac delta function). This generally makes the variational problem ill posed, in the sense that it lacks a solution within the domain of continuously differentiable functions. It is quite simple to construct a sequence of profiles { mi ( x ) } ∞ i = 1 , only slightly altering the regular one (given by the solution of the Laplace

- [1] T. Young, III. An essay on the cohesion of fluids, Philos. Trans. R. Soc. London 95 , 65 (1805).
- [2] A. I. Rusanov, Classification of line tension, Colloids Surf. A 156 , 315 (1999).
- [3] A. I. Rusanov, Surface thermodynamics revisited, Surf. Sci. Rep. 58 , 111 (2005).
- [4] V. S. Veselovskii and V. N. Pertsov, Prilipanie puzyrkov k tverdym poverhnostyam [Adhesion of bubbles to solid surfaces], Zh. Fiz. Khim. 8 , 245 (1936).
- [5] L. M. Shcherbakov and P. P. Ryazantsev, The effect of the energy of the wetting perimeter on the boundary conditions, in Research in Surface Forces , Vol. 2. Three-Dimensional Aspects of Surface Forces, edited by B. V. Deryagin (Consultants Bureau, New York, 1966) pp. 33-35.
- [6] B. A. Pethica, Contact angles, Rep. Prog. Appl. Chem. 46 , 14 (1961); The contact angle equilibrium, J. Colloid Interface Sci. 62 , 567 (1977).
- [7] A. I. Rusanov, On the wetting theory of elastically deformable bodies. 5. Reducing the deformation effects to linear tension, Colloid J. USSR 39 , 618 (1977).
- [8] A. I. Rusanov, Thermodynamics of solid surfaces, Surf. Sci.

and the generalized Young equations) in a very small vicinity of x = r in such a way that mi ′ ( r ) -- → -i → ∞ tan θ s , where θ s : ∂κ / ∂θ θ | = θ s = 0 (see Fig. 3(c)). However, such a sequence does not generally converge to a regular function, instead producing a singular contribution, clearly nonphysical in nature, to the equilibrium profile m x ( ) . A reasonable regularization procedure would involve the elimination of the singular part of the function, which corresponds to ignoring the condition (III) based on the contact-angle dependence of the line tension. Thus, the only case in which the problem is both physically meaningful and mathematically well posed is when ∂κ / ∂θ ≡ 0.

In regard to the possible dependence of the line tension κ a,c on the curvature radius R a,c of the liquid-gas interface, its consideration should lead to similar conclusions. Assuming such a dependence on the local value of the interface curvature/radius at the contact line, a similar construction with a local alteration of the droplet profile can be made, leading to the same conclusion. However, if the line tension is assumed to depend on the [nonlocal] mean curvature of the liquidgas interface, the rigorous variational problem statement remains unclear. This is because it is generally not known in advance that the solution will be a profile m a,c ( x ) of constant curvature. If this curvature is assumed to be constant and the line tension depends on its value, this would only affect the variation δ ( L a,c κ a,c ) in (A2), and specifically, the value of the line tension, but not the generalized Young equation for the contact angle. The curvature is still determined by the value of the chemical potential via the Gibbs-Duhem and the Laplace equations (which will be a consequence of the EulerLagrange equation). Therefore, the variation δ ( L a,c κ a,c ) is to be considered at a constant liquid-gas interface curvature, which yields the same outcome as if the curvature dependence of the line tension were incorporated directly into its chemical potential dependence.

## Rep. 29 , 173 (1996).

- [9] L. Boruvka and A. W. Neumann, Generalization of the classical theory of capillarity, J. Chem. Phys. 66 , 5464 (1977).
- [10] B. M. Law, S. P. McBride, J. Y. Wang, H. S. Wi, G. Paneru, S. Betelu, B. Ushijima, Y. Takata, B. Flanders, F. Bresme, H. Matsubara, T. Takiue, and M. Aratono, Line tension and its influence on droplets and particles at surfaces, Progr. Surf. Sci. 92 , 1 (2017).
- [11] F. Bresme and N. Quirke, Computer simulation study of the wetting behavior and line tensions of nanometer size particulates at a liquid-vapor interface, Phys. Rev. Lett. 80 , 3791 (1998).
- [12] A. I. Milchev and A. A. Milchev, Wetting behavior of nanodroplets: The limits of Young's rule validity, Europhys. Lett. 56 , 695 (2001).
- [13] H. K. Guo and H.-P. Fang, Drop size dependence of the contact angle of nanodroplets, Chin. Phys. Lett. 22 , 787 (2005).
- [14] G. Scocchi, D. Sergi, C. D'Angelo, and A. Ortona, Wetting and contact-line effects for spherical and cylindrical droplets on graphene layers: A comparative molecular-dynamics investigation, Phys. Rev. E 84 , 061602 (2011).

- [15] J. H. Weijs, A. Marchand, B. Andreotti, D. Lohse, and J. H. Snoeijer, Origin of line tension for a Lennard-Jones nanodroplet, Phys. Fluids 23 , 022001 (2011).
- [16] M. Barisik and A. Beskok, Wetting characterisation of silicon (1,0,0) surface, Mol. Simul 39 , 700 (2013).
- [17] H. Peng, G. R. Birkett, and A. V. Nguyen, The impact of line tension on the contact angle of nanodroplets, Mol. Simul 40 , 934 (2014).
- [18] M. Isaiev, S. Burian, L. Bulavin, M. Gradeck, F. Lemoine, and K. Termentzidis, Efficient tuning of potential parameters for liquid-solid interactions, Mol. Simul. 42 , 910 (2015).
- [19] M. Khalkhali, N. Kazemi, H. Zhang, and Q. Liu, Wetting at the nanoscale: A molecular dynamics study, J. Chem. Phys. 146 , 114704 (2017).
- [20] M. Kanduˇ c, Going beyond the standard line tension: Sizedependent contact angles of water nanodroplets, J. Chem. Phys. 147 , 174701 (2017).
- [21] D. V. Tatyanenko and A. K. Shchekin, Comparable effects of adsorption and line tension on contact angle of a nucleated droplet on a partially wettable substrate, Interfacial Phenom. Heat Transfer 5 , 113 (2017).
- [22] M. Kanduˇ c, L. Eixeres, S. Liese, and R. R. Netz, Generalized line tension of water nanodroplets, Phys. Rev. E 98 , 032804 (2018).
- [23] S. K. Das, S. A. Egorov, P. Virnau, D. Winter, and K. Binder, Do the contact angle and line tension of surface-attached droplets depend on the radius of curvature?, J. Phys.: Condens. Matter 30 , 255001 (2018).
- [24] D. V. Tatyanenko and A. K. Shchekin, Thermodynamic analysis of adsorption and line-tension contributions to contact angles of small sessile droplets, Colloid J. 81 , 455 (2019).
- [25] C. A. Ward and J. Wu, Effect of contact line curvature on solidfluid surface tensions without line tension, Phys. Rev. Lett. 100 , 256103 (2008).
- [26] C. A. Ward and K. Sefiane, Adsorption at the solid-liquid interface as the source of contact angle dependence on the curvature of the three-phase line, Adv. Colloid Interface Sci. 161 , 171 (2010).
- [27] J. Zhang, P. Wang, M. K. Borg, J. M. Reese, and D. Wen, A critical assessment of the line tension determined by the modified Young's equation, Phys. Fluids 30 , 082003 (2018).
- [28] M. Iwamatsu, A generalized Young's equation to bridge a gap between the experimentally measured and the theoretically calculated line tensions, J. Adhesion Sci. Technol. 32 , 2305 (2018).
- [29] B. Zhao, S. Luo, E. Bonaccurso, G. K. Auernhammer, X. Deng, Z. Li, and L. Chen, Resolving the apparent line tension of sessile droplets and understanding its sign change at a critical wetting angle, Phys. Rev. Lett. 123 , 094501 (2019).
- [30] N. Kubochkin and T. Gambaryan-Roisman, Wetting at nanoscale: Effect of surface forces and droplet size, Phys. Rev. Fluids 6 , 093603 (2021).
- [31] W. Klauser, F. T. von Kleist-Retzow, and S. Fatikow, Line tension and drop size dependence of contact angle at the nanoscale, Nanomaterials 12 , 369 (2022).
- [32] L. Schimmele and S. Dietrich, Line tension and the shape of nanodroplets, Eur. Phys. J. E 30 , 427 (2009).
- [33] S. M. Dammer and D. Lohse, Gas enrichment at liquid-wall interfaces, Phys. Rev. Lett. 96 , 206101 (2006).
- [34] F. Ould-Kaddour and D. Levesque, Molecular simulation of fluid-solid interfaces at nanoscale, J. Chem. Phys. 135 , 224705 (2011).
- [35] H. Jiang, F. Müller-Plathe, and A. Z. Panagiotopoulos, Contact angles from Young's equation in molecular dynamics simula-
- tions, J. Chem. Phys. 147 , 084708 (2017).
- [36] M. Kanduˇ c and R. R. Netz, Atomistic simulations of wetting properties and water films on hydrophilic surfaces, J. Chem. Phys. 146 , 164705 (2017).
- [37] M. Isaiev, S. Burian, L. Bulavin, W. Chaze, M. Gradeck, G. Castanet, S. Merabia, P. Keblinski, and K. Termentzidis, Gibbs adsorption impact on a nanodroplet shape: Modification of Young-Laplace equation, J. Phys. Chem. B 122 , 3176 (2018).
- [38] S. R. Carlson, O. Schullian, M. R. Becker, and R. R. Netz, Modeling water interactions with graphene and graphite via force fields consistent with experimental contact angles, J. Phys. Chem. Lett. 15 , 6325 (2024).
- [39] S. Burian, Y. Shportun, A. Yaroshchuk, L. Bulavin, D. Lacroix, and M. Isaiev, Size-dependent wetting contact angles at the nanoscale defined by equimolar surfaces and surfaces of tension, Sci. Rep. 14 , 31340 (2024).
- [40] S. Mechkov, G. Oshanin, M. Rauscher, M. Brinkmann, A. M. Cazabat, and S. Dietrich, Contact line stability of ridges and drops, EPL 80 , 66002 (2007).
- [41] A. I. Rusanov, D. V. Tatyanenko, and A. K. Shchekin, Comment on 'Size dependence of bubble wetting on surfaces: breakdown of contact angle match between small sized bubbles and droplets' by H. Zhang and X. Zhang, Nanoscale , 2019, 11, 2823, Nanoscale 13 , 4308 (2021).
- [42] A. I. Rusanov, A. K. Shchekin, and D. V. Tatyanenko, The line tension and the generalized Young equation: the choice of dividing surface, Colloids Surf. A 250 , 263 (2004).
- [43] P. Chen, The line adsorption equation: the one-dimensional counterpart of the Gibbs adsorption equation, Colloids Surf. A 161 , 23 (2000).
- [44] M. N. Popescu, G. Oshanin, S. Dietrich, and A.-M. Cazabat, Precursor films in wetting phenomena, J. Phys.: Condens. Matter 24 , 243102 (2012).
- [45] B. V. Derjaguin and Z. M. Zorin, Issledovaniya poverhnostnoy kondensatsii i adsorbtsii parov vblizi nasyscheniya opticheskim mikropolyarizatsionnym metodom. II. [Optical micropolarization studies of surface condensation and adsorption of vapors near saturation. II.], Zh. Fiz. Khim. 29 , 1755 (1955); Optical study of the adsorption and surface condensation of vapours in the vicinity of saturation on a smooth surface, Prog. Surf. Sci. 40 , 83 (1992).
- [46] A. C. Hall, Optical studies of thin films on surfaces of fused quartz, J. Phys. Chem. 74 , 2742 (1970).
- [47] A. Checco, P. Guenoun, and J. Daillant, Nonlinear dependence of the contact angle of nanodroplets on contact line curvature, Phys. Rev. Lett. 91 , 186101 (2003).
- [48] R. C. Tolman, The effect of droplet size on surface tension, J. Chem. Phys. 17 , 333 (1949).
- [49] T. Williams and C. Kelley, gnuplot 6.0. An interactive plotting program, http://gnuplot.info/docs\_6.0/Gnuplot\_ 6.pdf (2024).
- [50] L.-O. Heim and E. Bonaccurso, Measurement of line tension on droplets in the submicrometer range, Langmuir 29 , 14147 (2013).
- [51] L. Schimmele, M. Napiórkowski, and S. Dietrich, Conceptual aspects of line tensions, J. Chem. Phys. 127 , 164715 (2007).
- [52] M. S. Polovinkin, N. A. Volkov, D. V. Tatyanenko, and A. K. Shchekin, Contact angle calculations for argon and water sessile droplets on planar lyophilic and lyophobic surfaces within molecular dynamics modeling, Colloids Surf. A 702 , 134932 (2024).
- [53] I. M. Gelfand and S. V. Fomin, Calculus of Variations (PrenticeHall, Englewood Cliffs, NJ, 1963).
Dmitry V. Tatyanenko Konstantin D. Apitsin
Line tension from dual-geometry sessile droplet measurements: Combining contact-angle size-dependence data for axisymmetric and cylindrical droplets to determine the line tension
We perform a thermodynamic analysis of various contributions to the size
dependence of the contact-angle cosine for both axisymmetric and cylindrical
sessile droplets. This shows that a commonly used method to determine the line
tension from the slope of the contact-angle cosine dependence on the
three-phase contact-line curvature (for axisymmetric droplets) provides a
certain combination of the line tension, the adsorptions at the three
interfaces, and the macroscopic contact angle. To extract the contribution
related to the contact line in the leading order and determine the line
tension, we propose a simple technique using droplet-size dependences of the
contact angle for axisymmetric and cylindrical droplets under the same
conditions.