Open AccessArticle

On the Calculation of Solid-Fluid Contact Angles from Molecular Dynamics

Erik E. Santiso

^1,2,

Carmelo Herdes

¹ and

Erich A. Müller

^1,*

Department of Chemical Engineering, Imperial College London, South Kensington Campus, London, SW7 2AZ, UK

Department of Chemical and Biomolecular Engineering, North Carolina State University, Raleigh, NC 27695, USA

Author to whom correspondence should be addressed.

Entropy 2013, 15(9), 3734-3745; https://doi.org/10.3390/e15093734

Submission received: 21 July 2013 / Revised: 2 September 2013 / Accepted: 3 September 2013 / Published: 6 September 2013

(This article belongs to the Special Issue Molecular Dynamics Simulation)

Download

Browse Figures

Graphical abstract
"> Figure 1
Schematic of a liquid drop on a solid surface showing the contact angle. "> Figure 2
Frequency of contact angle values of water on graphite reported in literature; both from experimental results and numerical simulations [<a href="#B1-entropy-15-03734" class="html-bibr">1</a>]. "> Figure 3
Two-dimensional projections of a given configuration of a water droplet on a surface. The contact angles, measured using the auxiliary lines depicted in black, are: (a) 63.77° (b) 60.52° (c) 64.56° (d) 54.93°. "> Figure 4
(a) A snapshot from a MD simulation showing a droplet on top of a given surface (not shown). (b) The discretized density profile from the same system obtained using cubic subcells of width 3 nm. Density values are in molecules/Å [<a href="#B3-entropy-15-03734" class="html-bibr">3</a>]. "> Figure 5
(a) Water contact angle as a function of the water molecules on Wall05, inset shows the correspondent drop diameter, dashed lines are guide to the eye. (b) Snapshots of the drop interfaces for the smallest and the biggest system studied. "> Figure 6
(a) Water contact angle as a function of the fluid-substrate interactions, solid circles are simulation results, dashed red line marks the hydrophobic-hydrophilic threshold. (b) Corresponding equilibrium interface snapshots depicting interfacial beads. ">

Versions Notes

Abstract

A methodology for the determination of the solid-fluid contact angle, to be employed within molecular dynamics (MD) simulations, is developed and systematically applied. The calculation of the contact angle of a fluid drop on a given surface, averaged over an equilibrated MD trajectory, is divided in three main steps: (i) the determination of the fluid molecules that constitute the interface, (ii) the treatment of the interfacial molecules as a point cloud data set to define a geometric surface, using surface meshing techniques to compute the surface normals from the mesh, (iii) the collection and averaging of the interface normals collected from the post-processing of the MD trajectory. The average vector thus found is used to calculate the Cassie contact angle (i.e., the arccosine of the averaged normal z-component). As an example we explore the effect of the size of a drop of water on the observed solid-fluid contact angle. A single coarse-grained bead representing two water molecules and parameterized using the SAFT-γ Mie equation of state (EoS) is employed, meanwhile the solid surfaces are mimicked using integrated potentials. The contact angle is seen to be a strong function of the system size for small nano-droplets. The thermodynamic limit, corresponding to the infinite size (macroscopic) drop is only truly recovered when using an excess of half a million water coarse-grained beads and/or a drop radius of over 26 nm.

Keywords:

cloud data set; interfacial tension; coarse-graining; water; line tension; graphene

Graphical Abstract

1. Introduction

Wetting is the ability of a liquid to maintain contact with a solid surface, resulting from intermolecular interactions when the two are brought together. The presence of a liquid drop on a rigid surface is a reflection of the force balance between adhesive and cohesive forces and is commonly used to determine the wettability (the degree of wetting) of the solid-fluid system in terms of the solid-fluid contact angle, θ (see Figure 1). In this context, hydrophobicity is commonly referred to as the ability of a solid surface to repel water: if the water contact angle is smaller than 90°, the solid surface is considered hydrophilic and if the water contact angle is larger than 90°, the solid surface is considered hydrophobic.

Figure 1. Schematic of a liquid drop on a solid surface showing the contact angle.

Despite the fact of being such a well-defined problem the amount of conflicting (both theoretical and experimental) reported values for a given system is intriguing (see Figure 2, data from reference [1]).

Figure 2. Frequency of contact angle values of water on graphite reported in literature; both from experimental results and numerical simulations [1].

For instance, in the case of the graphite-water system, contact angles have been addressed extensively by experimental and theoretical approaches; however, a single general value has not been accepted [2,3,4,5,6,7,8,9].

A variety of causes for the discrepancies can be enumerated: heterogeneity and/or impurities at the surfaces or in the fluids, different methodological unstandardized calibration of equipment in experiments, possible system size effects and the distinct interaction potentials in simulations.

At the molecular scale, the main hurdle is that estimating contact angles for nanodroplets on surfaces is complicated by the fact that there are significant fluctuations in the shape of the droplet, and its geometry at a given step is often not axially symmetric. Furthermore, for very small nanoclusters, the fluid interfacial tension is a function of the curvature, and the planar limit is, in some cases, not recovered even after drop radii of 14 times the molecular diameter [10,11]. The change in the line tension with curvature (discussed in the latter part of this manuscript) is also an important factor affecting the result. It is seen that the contact angle will be, in the nanoscopic limit, a strong function of the system size [2]. In the analysis of simulations, contact angles [12,13] are commonly determined by using two-dimensional slices of the droplet and fitting its density profile to an empirical function, usually a circular section [14]. Such an approach, although appealing from the simplicity of the method, provides inconsistent results, particularly for small droplets [11]. Understandingly, different methods of increasing complexity have been devised for this purpose [15,16,17].

Figure 3 illustrates the difficulties associated with defining the contact angle using two-dimensional slices of molecular snapshots, especially for small droplets. Such droplets are often asymmetrical, and their shape changes substantially with time, due to capillary fluctuations. Using one or a few two-dimensional projections of a few molecular snapshots can thus potentially lead to large errors on the measured contact angle: for example, the values shown on Figure 3 vary over range of about 10°.

Figure 3. Two-dimensional projections of a given configuration of a water droplet on a surface. The contact angles, measured using the auxiliary lines depicted in black, are: (a) 63.77° (b) 60.52° (c) 64.56° (d) 54.93°.

In addition to the inherent fluctuations of the drop’s shape, there is a degree of arbitrariness in defining a single line separating the two phases in the two-dimensional projection. The interface layer is actually diffuse, as seen in Figure 3, spanning a width of several molecular diameters. This introduces another potential source of error, as the contact angle measured can change over a range of 5° or more depending on the choice of the contact line.

In this paper, we propose a method that does not make any a priori assumption about the shape of the drop, and uses the complete three-dimensional structure of the droplet near the surface to estimate the contact angle. For such analysis, we propose a geometrical estimation of the contact angle based on cloud point data breakdown of an equilibrated molecular dynamics (MD) trajectory of a drop on a given surface, this eliminates the effect of both spatial and temporal fluctuations on the estimated value of contact angle, yielding a value that better represents the average shape of the drop, as explained in detail in the next section.

2. Methodology

In order to estimate the average contact angle during an MD simulation, we first define a contact layer at each time step. This layer is defined as the set of molecules within the liquid-vapor interface that are close to (i.e., within a maximum distance z_max from) the solid surface. We then estimate the normal to the interface for each molecule in the contact layer by finding the plane that best fits the local shape of the interface at that point. The average of all the local normal vectors for all the time steps, <n>, can be used to calculate the average contact angle as θ = cos⁻¹(<n>·n_surf), where n_surf is the normal to the solid surface.

In practice, at each time step in the simulation, we carry out a two-step calculation: In the first step we use a discretized density profile to identify the molecules belonging to the liquid-vapor interface. In the second step, we estimate the local surface normal vectors for the molecules in the contact layer. The procedure is explained in detail below.

2.1. Identification of Interfacial Molecules

We identify the molecules belonging to the liquid-vapor interface using the following procedure:

(1): We divide the simulation box into subcells, and calculate the local number density in each subcell.
(2): We mark subcells as liquid if their number density is greater than a cutoff value ρ_c, otherwise we mark them as vapor cells. The cutoff density is chosen close to the average between the liquid and vapor number densities at the conditions of interest.
(3): Every subcell that is adjacent to at least one liquid cell and one vapor cell is marked as an interface cell.
(4): Finally, all molecules contained within the above-determined interface cells are marked as interfacial molecules.

The results obtained using the interface-sensing procedure are only weakly dependent on the choice of cutoff density, as long as this density is between the liquid and vapor densities. The number of subcells, on the other hand, should be chosen carefully: having more subcells identifies the interface with a higher resolution. However, if the number of subcells is too large, the density fluctuations within the droplet may cause the algorithm to incorrectly label too many cells as part of the interface. This can be easily detected by visualizing the interface molecules at one or a few simulation steps. Alternatively, one can use the average coordination number or the Minkowski-Bouligand dimension (obtained from a box-counting algorithm) [18,19] as a quantitative measure to select the appropriate subcell size. As the cells become smaller and the interface-sensing algorithm starts fitting the molecular-level fluctuations in the density rather than the actual interface, both the fractal dimension and the average coordination number will start deviating from the value expected for the larger subcells, indicating that the subcell size has become too small.

Figure 4. (a) A snapshot from a MD simulation showing a droplet on top of a given surface (not shown). (b) The discretized density profile from the same system obtained using cubic subcells of width 3 nm. Density values are in molecules/Å [3].

2.2. Estimation of Local Contact Angles

Having identified the molecules belonging to the interface, we choose a subset of the interface molecules as the interface contact layer. This layer contains the interface molecules that are within a given distance z_max from the solid surface. We then estimate the local contact angle at the position of each molecule i using the following procedure:

(a): Find all the interface molecules within a given cutoff radius, r_c, of molecule i.
(b): Find the average position r_avg of all the interface molecules found in step (a), including molecule i.
(c): Subtract the average position from the position of all the neighboring molecules, including i.
(d): Construct the covariance matrix of the centered positions, $r_{j}^{cent}$ found in step (c):

(1)

where denotes the outer (Kronecker) product.
(e): Find the eigenvector of Ω corresponding to its smallest eigenvalue. This is the normal to the plane that best fits the set of molecules in step (c). The sign of this normal is chosen to point away from the center of the droplet [20].

We average all the local normal vectors obtained at all the time steps to obtain the average normal. The component of this average normal perpendicular to the solid surface is the cosine of the average contact angle. The eigenvalues of the covariance matrix measure the variability of the positions along the three orthogonal directions defined by the eigenvectors. These can be used as a measure of how 1-, 2- or 3-dimensional the interface is. If one of the eigenvectors is close to zero, the interface is close to planar in that region.

The procedure described above to estimate the local contact angle requires arbitrarily setting a distance from the solid surface, z_max, to define the contact layer, and a cutoff radius r_c to estimate the tangent plane. In theory the contact angle would correspond to the limit when z_max→ 0. However, if z_max is chosen to be too small, the local density fluctuations within the contact layer cause the distribution of local contact angle values to be too noisy. Thus, one should choose the smallest value that gives a reasonable estimation error for the average contact angle. For the examples shown in this work, we have found z_max = 0.5 nm to give reasonable results.

For the cutoff radius, r_c, a compromise is also needed: if the value is too small, local density fluctuations cause the error in the estimated normal to be too large, whereas if the value is too large the orientation of the tangent plane will be affected by molecules far from the position of interest. The choice of r_c should be made by considerations similar to those used in choosing z_max. For the examples in this work, we have found r_c = 1 nm to be a reasonable value. A more detailed discussion of how to choose this cutoff can be found in reference [21]. Finally, we estimate the statistical error associated with the average contact angle computed by our method, we use the bootstrapping method [22].

2.3. Molecular Dynamics Details

Atomistic MD simulations are generally limited to nanometer-sized water droplets [2] Consequently, the apparent contact angle is usually drop-size dependent. To explore bigger systems, and aiming to find the optimal size for the MD contact angle calculation, we adopted a coarse-grained (CG) approach to describe the solid-water system in order to test the system size dependence beyond the atomistic limits, in a reasonable time [2].

We used the GROMACS simulation open source [23] suite to calculate the MD, which is well suited to implement Mie potentials [24]. Here, a single CG isotropic bead represents two water molecules [25]. Although several options are available for choosing the number of water molecules in a CG bead [26], this choice and the current parameterization produces sensible results, including a melting point, the surface tension, liquid densities and vapour pressures close to the experimental. The parameterization was carried out using SAFT-γ Mie approach, [27,28,29] where the water parameters were obtained by fitting to macroscopic properties, namely, the planar limit interfacial tension and liquid state density of water in a range from 0 °C to 40 °C. The SAFT EoS is a perturbation approach based on a well-defined Hamiltonian; here the CG beads are represented in the theory by a Mie potential, u:

u (r) = (\frac{λ_{r}}{λ_{r} - λ_{a}}) (\frac{λ_{r}}{λ_{a}}) ε [{(\frac{σ}{r})}^{λ_{r}} - {(\frac{σ}{r})}^{λ_{a}}] = A r^{- λ_{r}} - C r^{- λ_{a}}

(2)

where r is the intermolecular distance, and ε, and σ, are the adjustable parameters relating to the energy and distance scales, λ_a is the dispersion exponent and λ_r is the short-range repulsion. The potential is expressed in terms of two constants A and C for ease of tabulating in MD codes. Solid walls are modeled implicitly and described by an integrated potential of the same form in the z-dimension. Eight implicit solid surfaces of increasing ε _Wall#-W/k_B increments of □10 K, from 60 K to 130 K (labeled Wall01 to Wall08 respectively) are employed [30]. Table 1 summarizes the selected coarse-grained parameters.

Table 1. Coarse-grained parameters. W refers to a CG representation of water that accounts for two water molecules.

**Table 1.** Coarse-grained parameters. W refers to a CG representation of water that accounts for two water molecules.
Interaction	σ [nm]	ε/k_B [K]	λ_r, λ_a	C ×10² [kJ mol⁻¹ nm^λ^a]	A × 10⁴ [kJ mol⁻¹ nm^λ^r]
ε_Wall01-W/k_B	0.38716	60	10, 4	3.44107	1.15888
ε_Wall02-W/k_B		70		4.01458	1.35203
ε_Wall03-W/k_B		80		4.58810	1.54518
ε_Wall04-W/k_B		90		5.16161	1.73832
ε_Wall05-W/k_B		100		5.73512	1.93147
ε_Wall06-W/k_B		110		6.30863	2.12462
ε_Wall07-W/k_B		120		6.88214	2.31776
ε_Wall08-W/k_B		130		7.45565	2.51091
ε_W-W/k_B	0.37459	399.96	8, 6	8.71139	1.222380 × 10²

The systems are run under a canonical (NVT) ensemble, where the total volume, concentration and temperature are kept constant. Periodic boundary conditions were applied in xy dimensions, meanwhile an attractive wall was placed at z = 0 and a repulsive one (C = 5.73512 × 10⁻⁴ kJ mol⁻¹ nm⁴, A = 1.93147 × 10⁻⁶ kJ mol⁻¹ nm¹⁰ at the maximum height of the box. The number density of the atoms for each wall was set in 5 nm⁻² (c.f. in a Si crystal the number of atoms per nm² on the (100), (111) and (110) planes are 6.78, 7.83 and 9.59 nm⁻², respectively [31]). The simulations are thermostated to 298.15 K every 1 ps by a Nose-Hoover algorithm, all non-bonded interactions were truncated at 2.0 nm. The trajectories were recorded every 1000 time-steps (Δt = 0.01 ps) for at least 2000 ps after equilibrium.

3. Results and Discussions

We simulate 16,000, 32,000, 64,000, 128,000, 256,000, 512,000 and 1,024,000 water molecules (8,000, 16,000, 32,000, 64,000, 128,000 and 512,000 beads) on the eight solid surfaces. The simulation boxes dimensions can be seen in Table 2.

Table 2. Simulation box size for each system.

**Table 2.** Simulation box size for each system.
Water molecules	xy-dimensions [nm]	z-dimension [nm]
16,000	60 × 60	18
32,000
64,000
128,000
256,000	80 × 80	37.2
512,000	80 × 80	37.2
1,024,000	144 × 144	48

Although arbitrary, the reason guiding the choice of the simulation box size is to prevent the interaction between the sample and its periodic images. The meshing was done by dividing the simulation domain with a 1.2 × 1.2 × 1.2 nm³ subcell. The water droplet density contour is obtained from the cloud point data set analysis explained in the methodology section.

To test the size dependence of the water contact angle on a given surface we chose the intermediate Wall05 (see Table 1). The contact angles obtained using a moderately hydrophilic substrate (Wall05) as a function of the drop size can be seen in the Figure 5. The water contact angle over Wall05 exhibits a marked system-size dependence even up to 256,000 water molecules (128,000 beads) Larger drops exhibit a less pronounced but nevertheless noticeable size dependence. It is interesting to note that the effect of this scale up is an increase in the hydrophobicity of this surface in correspondence with previous simulations [32] and experimental studies at the macroscale. For this substrate, a system with 1,024,000 water molecules shows an apparent limiting contact angle of 74.30° [8]. The Wall05 parameters are chosen to loosely relate to graphene, a hydrophilic substrate.

Figure 5. (a) Water contact angle as a function of the water molecules on Wall05, inset shows the correspondent drop diameter, dashed lines are guide to the eye. (b) Snapshots of the drop interfaces for the smallest and the biggest system studied.

Figure 6. (a) Water contact angle as a function of the fluid-substrate interactions, solid circles are simulation results, dashed red line marks the hydrophobic-hydrophilic threshold. (b) Corresponding equilibrium interface snapshots depicting interfacial beads.

Following from the above, we use the system of 256,000 water molecules to study its contact angle as a function of the fluid-substrate interactions; the results can be seen in Figure 6. As expected, increasing the interaction energy between water molecules and the attractive wall diminishes the solid-fluid contact angle. For this particular coarse-grained surface model the functionality is linear, and a value of ε_Wall-W/k_B ~ 85 K can be taken as the boundary between hydrophilic and hydrophobic surface behavior.

An ancillary quantity of interest in the context of the drops on surfaces is the line tension. The line tension is the relation between the energy associated with the three phase contact line and the length of this line. This quantity is inherently dependent on the size of the drop [33] and is discussed as one of the reasons why the calculated tensions appear to be size-dependent. For large drops, one can obtain an estimate of the line tension using the approximation described by Weijs et al. [34]. For a typical case as shown above: the case of a drop of roughly 26 nm in diameter and considering the planar limit of the fluid-vapor tension (72 mN·m⁻¹) results in a contact tension line strength of O(10⁻⁸) J·m⁻¹. While this number seems to be different from those estimated from micrometer drop experiments [35] or from molecular simulation of atomistic water models [36,37]; it is appropriate to point out that values as low as ×10⁻¹¹ J·m⁻¹ and as high as ×10⁻⁵ J·m⁻¹ have been reported [33], which place our results in the correct context. Obviously, there is scope for much more detailed research into this topic.

4. Conclusions

We have proposed and validated a methodology for the unambiguous calculation of the solid-fluid contact angle from molecular dynamics simulations. We have tested model coarse-grained water-solid systems far beyond the limits commonly taken in atomistic simulation and showed, that for this particular model, more than 500,000 effective beads and/or drop diameters in excess of 50 nm would be required in order to obtain a result which is invariant of system size. So far the methodology has been applied over homogenous surfaces. However, is a well-known fact that surface roughness and energetic heterogeneities will have a profound effect on the contact angle calculations. The methodology presented is well suited to capture those effects.

Supplementary Materials

Supplementary File 1

Acknowledgments

This work was supported by the EPSRC through research grants (EP/I018212, EP/J014958 and EP/J010502).

Conflicts of Interest

The authors declare no conflict of interest.

References

Mattia, D. Templated growth and characterization of carbon nanotubes for nanofluidic applications. Ph.D. Thesis, Drexel University, Philadelphia, PA, USA, 2007. [Google Scholar]
Werder, T.; Walther, J.H.; Jaffe, R.L.; Halicioglu, T.; Koumoutsakos, P. On the water−carbon interaction for use in molecular dynamics simulations of graphite and carbon nanotubes. J. Phys. Chem. B 2003, 107, 1345–1352. [Google Scholar] [CrossRef]
Fowkes, F.M.; Harkins, W.D. The state of monolayers adsorbed at the interface solid-aqueous solution. J. Am. Chem. Soc. 1940, 62, 3377–3386. [Google Scholar] [CrossRef]
Schrader, M.E. Ultrahigh-vacuum techniques in the measurement of contact angles. LEED study of the effect of structure on the wettability of graphite. J. Phys. Chem. 1980, 84, 2774–2779. [Google Scholar] [CrossRef]
Hirvi, J.T.; Pakkanee, T.A. Molecular dynamics simulations of water droplets on polymer surfaces. J. Chem. Phys. 2006, 125, 144712. [Google Scholar] [CrossRef] [PubMed]
Hirvi, J.T.; Pakkanee, T.A. Enhanced hydrophobicity of rough polymer surfaces. J. Phys. Chem. B 2007, 111, 3336–3341. [Google Scholar] [CrossRef] [PubMed]
Hirvi, J.T.; Pakkanee, T.A. Wetting of nanogrooved polymer surfaces. Langmuir 2007, 23, 7724–7729. [Google Scholar] [CrossRef] [PubMed]
Ponter, A.B.; Yekta-Fard, M. The influence of environment on the drop size-contact angle relationship. Colloid Polym. Sci. 1985, 263, 673–681. [Google Scholar] [CrossRef]
Li, H.; Zeng, X.C. Wetting and interfacial properties of water nanodroplets in contact with graphene and monolayer boron-nitride sheets. ACS NANO 2012, 6, 2401–2409. [Google Scholar] [CrossRef] [PubMed]
Sampayo, J.G.; Malijevský, A.; Müller, E.A.; de Miguel, E.; Jackson, G. Communications: Evidence for the role of fluctuations in the thermodynamics of nanoscale drops and the implications in computations of the surface tension. J. Chem. Phys. 2010, 132, 141101. [Google Scholar] [CrossRef] [PubMed]
Malijevský, A.; Jackson, G. A perspective on the interfacial properties of nanoscopic liquid drops. J. Phys. Cond. Matt. 2012, 24, 464121. [Google Scholar] [CrossRef] [PubMed]
Ingebrigtsen, T.; Toxvaerd, S. Contact angles of Lennard-Jones liquids and droplets on planar surfaces. J. Phys. Chem. C 2007, 111, 8518–8523. [Google Scholar] [CrossRef]
Shi, B.; Dhir, V.K. Molecular dynamics simulation of the contact angle of liquids on solid surfaces. J. Chem. Phys. 2009, 130, 034705. [Google Scholar] [CrossRef] [PubMed]
De Ruijter, M.J.; Blake, T.D.; de Coninck, J. Dynamic wetting studied by molecular modeling simulations of droplet spreading. Langmuir 1999, 15, 7836–7847. [Google Scholar] [CrossRef]
Hautman, J.; Klein, M.L. Microscopic wetting phenomena. Phys Rev. Lett. 1991, 67, 1763–1766. [Google Scholar] [CrossRef] [PubMed]
Sergi, D.; Scocchi, G.; Ortona, A. Molecular dynamics simulations of the contact angle between water droplets and graphite surfaces. Fluid Phase Equilib. 2012, 332, 173–177. [Google Scholar] [CrossRef]
Das, S.K.; Binder, K. Does Young’s equation hold on the nanocale? A monte Carlo test for the binary Lennard-Jones fluid. Europhys. Lett. 2006, 92, 26006. [Google Scholar] [CrossRef]
Falconer, K. Fractal Geometry: Mathematical Foundations and Applications; Wiley: Chichester, UK, 1990. [Google Scholar]
Mandelbrot, B. The Fractal Geometry of Nature; Macmillan: New York, NY, USA, 1983. [Google Scholar]
Berkmann, J.; Caelli, T. Computation of surface geometry and segmentation using covariance techniques. IEEE Trans. Pattern Anal. 1994, 16, 1114–1116. [Google Scholar] [CrossRef]
Mitra, N.J.; Nguyen, A.; Guibas, L. Estimating surface normals in noisy point cloud data. Int. J. Comput. Geom. Appl. 2004, 14, 261–276. [Google Scholar] [CrossRef]
Efron, B.; Tibshirani, T. An Introduction to The Bootstrap; Chapman & Hall/CRC Monographs on Statistics and Applied Probability; Chapman and Hall/CRC: Boca Raton, FL, USA, 1993. [Google Scholar]
Van der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef] [PubMed]
Mie, G. Zur kinetischen Theorie der einatomigen Körper. Ann. Phys. 1903, 11, 657–697. (In German) [Google Scholar] [CrossRef]
Lobanova, O.; Avendaño, C.; Jackson, G.; Müller, E.A. SAFT-gamma force field for the simulation of molecular fluids: 5. A single-site coarse grained model of water. 2013. In preparation. [Google Scholar]
Hadley, K.R.; McCabe, C. Coarse-grained molecular models of water: A review. Mol. Simul. 2012, 38, 671–681. [Google Scholar] [CrossRef] [PubMed]
Avendaño, C.; Lafitte, T.; Galindo, A.; Adjiman, C.S.; Jackson, G.; Müller, E.A. SAFT-gamma force field for the simulation of molecular fluids: 1. A single-site coarse grained model of carbon dioxide. J. Phys. Chem. B 2011, 115, 11154–11169. [Google Scholar] [CrossRef] [PubMed]
Avendaño, C.; Lafitte, T.; Adjiman, C.S.; Galindo, A.; Müller, E.A.; Jackson, G. SAFT-γ force field for the simulation of molecular fluids: 2. Coarse-grained models of greenhouse gases, refrigerants, and long alkanes. J. Phys. Chem. B 2013, 117, 2717–2733. [Google Scholar] [CrossRef] [PubMed]
Lafitte, T.; Avendaño, C.; Papaioannou, V.; Galindo, A.; Adjiman, C.S.; Jackson, G.; Müller, E.A. SAFT-gamma force field for the simulation of molecular fluids: 3. Coarse-grained models of benzene and hetero-group models of n-decylbenzene. Mol. Phys. 2012, 110, 1189–1203. [Google Scholar] [CrossRef]
Israelachvili, J.N. Intermolecular and Surface Forces, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2011; Chapter 11; pp. 208–211. [Google Scholar]
Shakouri, A. EE145 Spring 2002. Available online: http://classes.soe.ucsc.edu/ee145/Spring02/EE145hmwk1Sol.pdf (accessed on 10 July 2013).
Scocchi, G.; Sergi, D.; D’Angelo, C.; Ortona, A. Wetting and contact-line effects for spherical and cylindrical droplets on graphene layers: A comparative molecular-dynamics investigation. Phys. Rev. E 2011, 84, 061602. [Google Scholar] [CrossRef]
Tadmor, R. Line energy, line tension and drop size. Surf. Sci. Lett. 2008, 602, L108–L111. [Google Scholar] [CrossRef]
Weijs, J.H.; Marchand, A.; Andreotti, B.; Lohse, D.; Snoeijer, J.H. Origin of line tension for a Lennard-Jones nanodroplet. Phys. Fluids 2011, 23, 022001. [Google Scholar] [CrossRef]
Pompe, T.; Herminghaus, S. Three-phase contact line energies from nanosclae surface topographies. Phys. Rev. Lett. 2000, 85, 1930–1933. [Google Scholar] [CrossRef] [PubMed]
Ritchie, J.A.; Yazdi, J.S.; Bratko, D.; Luzar, A. Metastable sessile nanodroplets on nanopatterned surfaces. J. Phys. Chem. C 2012, 116, 8634–8641. [Google Scholar] [CrossRef]
Daub, C.D.; Bartko, D.; Luzar, A. Electric control of wetting by salty nanodrops: Molecular dynamics simulations. J. Phys. Chem. C 2011, 115, 22393–22399. [Google Scholar] [CrossRef]

© 2013 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).

Share and Cite

MDPI and ACS Style

Santiso, E.E.; Herdes, C.; Müller, E.A. On the Calculation of Solid-Fluid Contact Angles from Molecular Dynamics. Entropy 2013, 15, 3734-3745. https://doi.org/10.3390/e15093734

AMA Style

Santiso EE, Herdes C, Müller EA. On the Calculation of Solid-Fluid Contact Angles from Molecular Dynamics. Entropy. 2013; 15(9):3734-3745. https://doi.org/10.3390/e15093734

Chicago/Turabian Style

Santiso, Erik E., Carmelo Herdes, and Erich A. Müller. 2013. "On the Calculation of Solid-Fluid Contact Angles from Molecular Dynamics" Entropy 15, no. 9: 3734-3745. https://doi.org/10.3390/e15093734

Article Menu

On the Calculation of Solid-Fluid Contact Angles from Molecular Dynamics

Abstract

1. Introduction

2. Methodology

2.1. Identification of Interfacial Molecules

2.2. Estimation of Local Contact Angles

2.3. Molecular Dynamics Details

3. Results and Discussions

4. Conclusions

Supplementary Materials

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI