Dynamics of multicellular swirling on micropatterned substrates

Chirality plays a crucial role in biology, as it is highly conserved and fundamentally important in the developmental process. To better understand the relationship between the chirality of individual cells and that of tissues and organisms, we develop a generalized mechanics model of chiral polarized particles to investigate the swirling dynamics of cell populations on substrates. Our analysis reveals that cells with the same chirality can form distinct chiral patterns on ring-shaped or rectangular substrates. Interestingly, our studies indicate that an excessively strong or weak individual cellular chirality hinders the formation of such chiral patterns. Our studies also indicate that there exists the influence distance of substrate boundaries in chiral patterns. Smaller influence distances are observed when cell-cell interactions are weaker. Conversely, when cell-cell interactions are too strong, multiple cells tend to be stacked together, preventing the formation of chiral patterns on substrates in our analysis. Additionally, we demonstrate that the interaction between cells and substrate boundaries effectively controls the chiral distribution of cellular orientations on ring-shaped substrates. This research highlights the significance of coordinating boundary features, individual cellular chirality, and cell-cell interactions in governing the chiral movement of cell populations and provides valuable mechanics insights into comprehending the intricate connection between the chirality of single cells and that of tissues and organisms.

Modeling Water Interactions with Graphene and Graphite via Force Fields Consistent with Experimental Contact Angles

Accurate simulation models for water interactions with graphene and graphite are important for nanofluidic applications, but existing force fields produce widely varying contact angles. Our extensive review of the experimental literature reveals extreme variation among reported values of graphene-water contact angles and a clustering of graphite-water contact angles into groups of freshly exfoliated (60° ± 13°) and not-freshly exfoliated graphite surfaces. The carbon-oxygen dispersion energy for a classical force field is optimized with respect to this 60° graphite-water contact angle in the infinite-force-cutoff limit, which in turn yields a contact angle for unsupported graphene of 80°, in agreement with the mean of the experimental results. Interaction force fields for finite cutoffs are also derived. A method for calculating contact angles from pressure tensors of planar equilibrium simulations that is ideally suited to graphite and graphene surfaces is introduced. Our methodology is widely applicable to any liquid-surface combination.

Multiscale Modeling of Aqueous Electric Double Layers

From the stability of colloidal suspensions to the charging of electrodes, electric double layers play a pivotal role in aqueous systems. The interactions between interfaces, water molecules, ions and other solutes making up the electrical double layer span length scales from Ångströms to micrometers and are notoriously complex. Therefore, explaining experimental observations in terms of the double layer's molecular structure has been a long-standing challenge in physical chemistry, yet recent advances in simulations techniques and computational power have led to tremendous progress. In particular, the past decades have seen the development of a multiscale theoretical framework based on the combination of quantum density functional theory, force-field based simulations and continuum theory. In this Review, we discuss these theoretical developments and make quantitative comparisons to experimental results from, among other techniques, sum-frequency generation, atomic-force microscopy, and electrokinetics. Starting from the vapor/water interface, we treat a range of qualitatively different types of surfaces, varying from soft to solid, from hydrophilic to hydrophobic, and from charged to uncharged.

Insights into the Rearrangement of the Molecular Assembly Structure of 6-Aminonicotinic Acid in Its Hydrated Environment

Water, as a ubiquitous and essential component of life, is known to have a significant impact on the structure and function of organic molecules. In this study, we investigate the role of water in the phase transition of organic molecular assembly structures by scanning tunneling microscopy at room temperature. The results show that the -O-H···O hydrogen induced by water molecules can lead to a significant structural transition in the molecular assembly, specifically through selective weakening of -C-H···O between 6-aminonicotinic acid and the formation of new -O-H···O bonds between 6-aminonicotinic acid and water molecules. Subsequent thermal treatment of these molecular assembly structures reveals that the formation of -N-H···O hydrogen bonds induced by water molecules has created a different pathway for the phase transition of the molecular assembly structure. This knowledge has important implications for the design of organic molecules with specific nanostructures that can be controlled through hydration.

Hydrophobicity of Self-Assembled Monolayers of Alkanes: Fluorination, Density, Roughness, and Lennard-Jones Cutoffs

The interplay of fluorination and structure of alkane self-assembled monolayers and how these affect hydrophobicity are explored via molecular dynamics simulations, contact angle goniometry, and surface-enhanced infrared absorption spectroscopy. Wetting coefficients are found to grow linearly in the monolayer density for both alkane and perfluoroalkane monolayers. The larger contact angles of monolayers of perfluorinated alkanes are shown to be primarily caused by their larger molecular volume, which leads to a larger nearest-neighbor grafting distance and smaller tilt angle. Increasing the Lennard-Jones force cutoff in simulations is found to increase hydrophilicity. Specifically, wetting coefficients scale like the inverse square of the cutoff, and when extrapolated to the infinite cutoff limit, they yield contact angles that compare favorably to experimental values. Nanoscale roughness is also found to reliably increase monolayer hydrophobicity, mostly via the reduction of the entropic part of the work of adhesion. Analysis of depletion lengths shows that droplets on nanorough surfaces partially penetrate the surface, intermediate between Wenzel and Cassie-Baxter states.

Experimental study and modelling of asphaltene deposition on metal surfaces with superhydrophobic and low sliding angle inner coatings

Inner coatings have emerged as a novel technique to prevent the deposition of paraffin, wax, scale, and corrosion of pipelines during oil production and transport. Few studies addressed this technique for preventing asphaltene deposition. In this study, two superhydrophobic inner coatings, including polytetrafluoroethylene (PTFE) coating and nanosilica coating, were fabricated on metal surfaces and the asphaltene deposition on these coated surfaces was examined. A model oil solution was prepared using asphaltene and heptol and the effect of static and dynamic flow states on the amount of asphaltene deposition on uncoated electrodes, PTFE coated electrodes, and nanosilica coated electrodes were investigated. The results showed that the PTFE coating is more effective in reducing asphaltene deposition than nanosilica coating. The PTFE coating could reduce 56% of the deposition in a static state and more than 70% in a dynamic state at an asphaltene concentration of 2000 ppm. For PTFE coating in a dynamic state, the deposition rate is negligible in long times. In addition, it was found that the type of flow state affects the asphaltene deposition kinetics. The results demonstrate that, in the static state, the nth-order kinetics model, and in the dynamic state, the double exponential models are in best agreement with the experimental data.

Carbon dioxide as a line active agent: Its impact on line tension and nucleation rate

By considering a water capillary bridge confined between two flat surfaces, we investigate the thermodynamics of the triple line delimiting this solid-liquid-vapor system when supplemented in carbon dioxide. In more detail, by means of atom-scale simulations, we show that carbon dioxide accumulates at the solid walls and, preferably, at the triple lines where it plays the role of a line active agent. The line tension of the triple line, which is quantitatively assessed using an original mechanical route, is shown to be driven by the line excess concentrations of the solute (carbon dioxide) and solvent (water). Solute accumulation at the lines decreases the negative line tension (i.e., more negative) while solvent depletion from the lines has the opposite effect. Such an unprecedented quantitative assessment of gas-induced line tension modifications shows that the absolute value of the negative line tension increases upon increasing the carbon dioxide partial pressure. As a striking example, for hydrophilic surfaces, the line tension is found to increase by more than an order of magnitude when the carbon dioxide pressure exceeds 3 MPa. By considering the coupling between line and surface effects induced by gaseous adsorption, we hypothesize from the observed gas concentration-dependent line tension a nontrivial impact on heterogeneous nucleation of nanometric critical nuclei.

Water-Repellent Characteristics of Beech Wood Coated with Parylene-N

In recent years, awareness regarding sustainability and the responsible usage of natural resources has become more important in our modern society. As a result, wood as a building material experiences a renaissance. However, depending on the use case, protective measures may be necessary to increase wood's durability and prolong its service life. The chemical vapor deposition (CVD) of parylene-N layers offers an interesting alternative to solvent-based and potentially environmentally harmful coating processes. The CVD process utilized in this study generated transparent, uniform barrier layers and can be applied on an extensive range of substrates without the involvement of any solvents. In this study, European beech wood samples (Fagus sylvatica L.) were coated with parylene-N using the CVD process, with paracyclophane as a precursor. The aim of the study was to analyze the water absorption of beech wood, in relation to the different layer thicknesses of parylene-N. Therefore, four different coating thicknesses from 0.5 to 40 μm were deposited, depending on the initial amount of precursor used. The deposited layers were analyzed by reflection interference spectroscopy and scanning electron microscopy, and their chemical structures and compositions were investigated by X-ray photoelectron spectroscopy and Fourier-transform infrared spectroscopy. Due to the chemical structure of parylene-N, the deposited layers led to a significantly increased water contact angle and reduced the water uptake by 25-34% compared to the uncoated reference samples. A linear correlation between layer thickness and water absorption was observed. The coating of wood with parylene-N provides a promising water barrier, even with thin layers.

How the hydroxylation state of the (110)-rutile TiO2 surface governs its electric double layer properties

The hydroxylation state of an oxide surface is a central property of its solid/liquid interface and its corresponding electrical double layer. This study integrated both a reactive force field (ReaxFF) and a non-reactive potential into a hierarchical framework within molecular dynamics (MD) simulations to reveal how the hydroxylation state of the (110)-rutile TiO2 surface affects the electrical double layer properties. The simulation results obtained in the ReaxFF framework have shown that, while water dissociation occurs only at the under-coordinated Ti5c sites on the pristine TiO2 surface, the presence of point defects on the surface facilitates water dissociation at the oxygen vacancy sites, leading to two protonated oxygen bridge atoms for each vacancy site. As a consequence of enhanced water dissociation at the vacancy sites, water dissociation is quenched at the under-coordinated Ti5c sites resulting in two competitive hydroxylation mechanisms on the (110)-TiO2 surface. Using non-reactive MD simulations with hydroxylation states derived from the ReaxFF analysis, we demonstrate that water dissociation at the vacancy sites is a central mechanism governing the structuring of water near the interface. While the structuring of water near the interface is the main contribution to the electric field, water dissociation at the vacancy site enhances the adsorption of the electrolyte ions at the interface. The adsorbed ions lead to an increase of the effective surface charge as well as surface (zeta) potentials which are in the range of experimental observations. Our work provides a hierarchical multiscale simulation approach, covering a series of results with in-depth discussion for atomic/molecular level understanding of water dissociation and its effect on electric double layer properties of TiO2 to advance water splitting.

Effect of Hydrogen Bonds on the Dielectric Properties of Interfacial Water

The dielectric constant for water is reduced under confinement. Although this phenomenon is well known, the underlying physical mechanism for the reduction is still in debate. In this work, we investigate the effect of the orientation of hydrogen bonds on the dielectric properties of confined water using molecular dynamics simulations. We find a reduced rotational diffusion coefficient for water molecules close to the solid surface. The reduced rotational diffusion arises due to the hindered rotation away from the plane parallel to the channel walls. The suppressed rotation in turn affects the orientational polarization of water, leading to a low value for the dielectric constant at the interface. We attribute the constrained out-of-plane rotation to originate from a higher density of planar hydrogen bonds formed by the interfacial water molecules.

Fluidity and phase transitions of water in hydrophobic and hydrophilic nanotubes

We put water flow under scrutiny to report radial distributions of water viscosity within hydrophobic and hydrophilic nanotubes as functions of the water-nanotube interactions ([Formula: see text]), surface wettability (θ), and nanotube size (R) using a proposed hybrid continuum-molecular mechanics. Based on the computed viscosity data, [Formula: see text] phase diagram of the phase transitions of confined water in nanotubes is developed. It is revealed that water exhibits different multiphase structures, and the formation of one of these structures depends on [Formula: see text] R parameters. A drag of water flow at the first water layer is revealed, which is conjugate to sharp increase in the viscosity and formation of an ice phase under severe confinement (R ≤ 3.5 nm) and strong water-nanotube interaction conditions. A vapor/vapor-liquid phase is observed at hydrophobic and hydrophilic interfaces. A state of confinement is revealed at which water exhibits different multiphase structures under the same flow rate. The derived viscosity functions are used to accurately determine factors of flow enhancement/inhibition of confined water.

On the equilibrium contact angle of sessile liquid drops from molecular dynamics simulations

We present a new methodology to estimate the contact angles of sessile drops from molecular simulations by using the Gaussian convolution method of Willard and Chandler [J. Phys. Chem. B 114, 1954-1958 (2010)] to calculate the coarse-grained density from atomic coordinates. The iso-density contour with average coarse-grained density value equal to half of the bulk liquid density is identified as the average liquid-vapor (LV) interface. Angles between the unit normal vectors to the average LV interface and unit normal vector to the solid surface, as a function of the distance normal to the solid surface, are calculated. The cosines of these angles are extrapolated to the three-phase contact line to estimate the sessile drop contact angle. The proposed methodology, which is relatively easy to implement, is systematically applied to three systems: (i) a Lennard-Jones (LJ) drop on a featureless LJ 9-3 surface; (ii) an SPC/E water drop on a featureless LJ 9-3 surface; and (iii) an SPC/E water drop on a graphite surface. The sessile drop contact angles estimated with our methodology for the first two systems are shown to be in good agreement with the angles predicted from Young's equation. The interfacial tensions required for this equation are computed by employing the test-area perturbation method for the corresponding planar interfaces. Our findings suggest that the widely adopted spherical-cap approximation should be used with caution, as it could take a long time for a sessile drop to relax to a spherical shape, of the order of 100 ns, especially for water molecules initiated in a lattice configuration on a solid surface. But even though a water drop can take a long time to reach the spherical shape, we find that the contact angle is well established much faster and the drop evolves toward the spherical shape following a constant-contact-angle relaxation dynamics. Making use of this observation, our methodology allows a good estimation of the sessile drop contact angle values even for moderate system sizes (with, e.g., 4000 molecules), without the need for long simulation times to reach the spherical shape.

Pressure increases the ice-like order of water at hydrophobic interfaces

Ice-like order of water at hydrophobic interfaces is quantified on different length scales based on molecular dynamics simulations.

Communication: Is a coarse-grained model for water sufficient to compute Kapitza conductance on non-polar surfaces?

Coarse-grained models have increasingly been used in large-scale particle-based simulations. However, due to their lack of degrees of freedom, it is a priori unlikely that they straightforwardly represent thermal properties with the same accuracy as their atomistic counterparts. We take a first step in addressing the impact of liquid coarse-graining on interfacial heat conduction by showing that an atomistic and a coarse-grained model of water may yield similar values of the Kapitza conductance on few-layer graphene with interactions ranging from hydrophobic to mildly hydrophilic. By design the water models employed yield similar liquid layer structures on the graphene surfaces. Moreover, they share common vibration properties close to the surfaces and thus couple with the vibrations of graphene in a similar way. These common properties explain why they yield similar Kapitza conductance values despite their bulk thermal conductivity differing by more than a factor of two.

The Hydrophobic Gap at High Hydrostatic Pressures

We have gained new insight into the so-called hydrophobic gap, a molecularly thin region of decreased electron density at the interface between water and a solid hydrophobic surface, by X-ray reflectivity experiments and molecular dynamics simulations at different hydrostatic pressures. Pressure variations show that the hydrophobic gap persists up to a pressure of 5 kbar. The electron depletion in the interfacial region strongly decreases with an increase in pressure, indicating that the interfacial region is compressed more strongly than bulk water. The decrease is most significant up to 2 kbar; beyond that, the pressure response of the depletion is less pronounced.

Water nanodroplet on a graphene surface-a new old system

The major subject of our study is the accuracy of contact angle calculations. Reporting new simulation data for graphene-water systems, we show that the majority of previously reported data should be treated with caution, since the proper contact angle can be recorded only after a sufficiently long simulation time. It has been proven that-if one wants to gain accuracy greater than 0.1°-long calculations (exceeding 50 ns) are required. Finally, we also show, using both a Groningen Machine for Chemical Simulations (GROMACS) package and our new molecular dynamics (MD) code, that the changes in the contact angle, caused by graphene bottom layer rotation, are within the range of calculation error. We also propose a novel definition of the bottom of the droplet as the height where the density is half the density of liquid water. This new definition is applied in the method of the contact angle calculation from the MD simulation data.

Dominance of Dispersion Interactions and Entropy over Electrostatics in Determining the Wettability and Friction of Two-Dimensional MoS2 Surfaces

The existence of partially ionic bonds in molybdenum disulfide (MoS₂), as opposed to covalent bonds in graphene, suggests that polar (electrostatic) interactions should influence the interfacial behavior of two-dimensional MoS₂ surfaces. In this work, using molecular dynamics simulations, we show that electrostatic interactions play a negligible role in determining not only the equilibrium contact angle on the MoS₂ basal plane, which depends solely on the total interaction energy between the surface and the liquid, but also the friction coefficient and the slip length, which depend on the spatial variations in the interaction energy. While the former is found to result from the exponential decay of the electric potential above the MoS₂ surface, the latter results from the trilayered sandwich structure of the MoS₂ monolayer, which causes the spatial variations in dispersion interactions in the lateral direction to dominate over those in electrostatic interactions in the lateral direction. Further, we show that the nonpolarity of MoS₂ is specific to the two-dimensional basal plane of MoS₂ and that other planes (e.g., the zigzag plane) in MoS₂ are polar with respect to interactions with water, thereby illustrating the role of edge effects, which could be important in systems involving vacancies or nanopores in MoS₂. Finally, we simulate the temperature dependence of the water contact angle on MoS₂ to show that the inclusion of entropy, which has been neglected in recent mean-field theories, is essential in determining the wettability of MoS₂. Our findings reveal that the basal planes in graphene and MoS₂ are unexpectedly similar in terms of their interfacial behavior.

A Relation for Nanodroplet Diffusion on Smooth Surfaces

In this work, we study the diffusion of nanodroplets on smooth surfaces through molecular dynamics (MD) simulations and theoretical analyses. Molecular dynamics simulations show that nanodroplet surface diffusion is different from that of single molecules and solid particles. The dependence of nanodroplet diffusion coefficient on temperature undergoes a transition from linear to nonlinear as the surface wettability is weakened due to the coupling of temperature and surface energy. We also develop a simple relation for the diffusion coefficient by using the contact angle and contact radius of the droplet. It works well for a wide range of surface wettabilities and different sized nanodroplets, as confirmed by MD simulations.

Fluorine-Terminated Diamond Surfaces as Dense Dipole Lattices: The Electrostatic Origin of Polar Hydrophobicity

Despite the pronounced polarity of C-F bonds, many fluorinated carbon compounds are hydrophobic: a controversial phenomenon known as "polar hydrophobicity". Here, its underlying microscopic mechanisms are explored by ab initio calculations of fluorinated and hydrogenated diamond (111) surfaces interacting with single water molecules. Gradient- and van der Waals-corrected density functional theory simulations reveal that "polar hydrophobicity" of the fully fluorinated surfaces is caused by a negligible surface/water electrostatic interaction. The densely packed C-F surface dipoles generate a short-range electric field that decays within the core repulsion zone of the surface and hence vanishes in regions accessible by adsorbates. As a result, water physisorption on fully F-terminated surfaces is weak (adsorption energies Ead < 0.1 eV) and dominated by van der Waals interactions. Conversely, the near-surface electric field generated by loosely packed dipoles on mixed F/H-terminated surfaces has a considerably longer range, resulting in a stronger water physisorption (Ead > 0.2 eV) that is dominated by electrostatic interactions. The suppression of electrostatic interactions also holds for perfluorinated molecular carbon compounds, thus explaining the prevalent hydrophobicity of fluorocarbons. In general, densely packed polar terminations do not always lead to short-range electric fields. For example, surfaces with substantial electron density spill-out give rise to electric fields with a much slower decay. However, electronic spill-out is limited in F/H-terminated carbon materials. Therefore, our ab initio results can be reproduced and rationalized by a simple classical point-charge model. Consequently, classical force fields can be used to study the wetting of F/H-terminated diamond, revealing a pronounced correlation between adsorption energies of single H2O molecules and water contact angles.

Bulk and interfacial liquid water as a transient network

The special macroscopic properties of liquid water stem from its structure as a complex network of molecules connected by hydrogen bonds. While the dynamics of single molecules within this network has been extensively investigated, only little attention has been paid to the closed loops (meshes) of hydrogen-bonded molecules which determine the network topology. Using molecular dynamics simulations we analyze the size, shape, geometrical arrangement, and dynamical stability of loops containing up to 10 hydrogen bonds. We find that six-membered loops in liquid water even at room temperature retain a striking similarity with the well-known structure of ice. Analyzing the network dynamics we find that rings of more than five hydrogen bonds are stabilized compared to a random collection containing the same number of single bonds. We finally show that in the vicinity of hydrophobic and hydrophilic interfaces loops arrange in a preferred orientation.

Dry-Surface Simulation Method for the Determination of the Work of Adhesion of Solid-Liquid Interfaces

We introduce a methodology, referred to as the dry-surface method, to calculate the work of adhesion of heterogeneous solid-liquid interfaces by molecular simulation. This method employs a straightforward thermodynamic integration approach to calculate the work of adhesion as the reversible work to turn off the attractive part of the actual solid-liquid interaction potential. It is formulated in such a way that it may be used either to evaluate the ability of force fields to reproduce reference values of the work of adhesion or to optimize force-field parameters with reference values of the work of adhesion as target quantities. The methodology is tested in the case of water on a generic model of nonpolar substrates with the structure of gold. It is validated through a quantitative comparison to phantom-wall calculations and against a previous characterization of the thermodynamics of the gold-water interface. It is found that the work of adhesion of water on nonpolar substrates is a nonlinear function of the microscopic solid-liquid interaction energy parameter. We also comment on the ability of mean-field approaches to predict the work of adhesion of water on nonpolar substrates. In addition, we discuss in detail the information on the solid-liquid interfacial thermodynamics delivered by the phantom-wall approach. We show that phantom-wall calculations yield the solid-liquid interfacial tension relative to the solid surface tension rather than the absolute solid-liquid interfacial tension as previously believed.

Progress on the Surface Nanobubble Story: What is in the bubble? Why does it exist?

Interfaces between aqueous solutions and hydrophobic solid surfaces are important in various areas of science and technology. Many researchers have found that forces between hydrophobic surfaces in aqueous solution are significantly different from the classical DLVO theory. Long-range attractive forces (non-DLVO forces) are thought to be affected by nanoscopic gaseous domains at the interfaces. This is a review of the latest research on nanobubbles at hydrophobic surfaces from experimental and simulation studies. The review focusses on non-intrusive optical view of surface nanobubbles and gas enrichment on solid surfaces by imaging and force mapping. By use of these recent experimental data in conjunction with molecular simulation work, all major theories on surface nanobubble formation and stability are critically reviewed. Even though the current body of research cannot comprehensively explain all properties of surface nanobubbles observed, the fundamental understanding has been significantly improved. Line tension has been shown to be incapable of explaining the contact angle of nanobubbles. Dense gas layer theory provides a new explanation on both large contact angle and long-time stability. The high density of gas in these domains may significantly affect the gas-water interface which is in line with some observation made on bulk nanobubbles. Along this line of inquiry, experimental and simulation effort should be focussed on measuring the density within surface nanobubbles and the properties of the gas water interface which may be the key to explaining the stability of these nanobubbles.

Detachment of semiflexible polymer chains from a substrate: a molecular dynamics investigation

Using Molecular Dynamics simulations, we study the force-induced detachment of a coarse-grained model polymer chain from an adhesive substrate. One of the chain ends is thereby pulled at constant speed off the attractive substrate and the resulting saw-tooth profile of the measured mean force ⟨f⟩ vs height D of the end-segment over the plane is analyzed for a broad variety of parameters. It is shown that the observed characteristic oscillations in the ⟨f⟩-D profile depend on the bending and not on the torsional stiffness of the detached chains. Allowing for the presence of hydrodynamic interactions (HI) in a setup with explicit solvent and dissipative particle dynamics-thermostat, rather than the case of Langevin thermostat, one finds that HI have little effect on the ⟨f⟩-D profile. Also the change of substrate affinity with respect to the solvent from solvophilic to solvophobic is found to play negligible role in the desorption process. In contrast, a changing ratio ε(s)(B)/ε(s)(A) of the binding energies of A- and B-segments in the detachment of an AB-copolymer from adhesive surface strongly changes the ⟨f⟩-D profile whereby the B-spikes vanish when ε(s)(B)/ε(s)(A)<0.15. Eventually, performing an atomistic simulation of (bio)-polymers, we demonstrate that the simulation results, derived from our coarse-grained model, comply favorably with those from the all-atom simulation.

Wettability of azobenzene self-assembled monolayers

The wettability properties of azobenzene self-assembled monolayers (SAMs), in the trans and cis forms, are investigated herein by classical Molecular Dynamics simulations of validated assembly structures described with a dedicated force field. The two different methodologies used for the calculation of the contact angle, one based on the Young's equation and the other on geometrical models, have provided a consistent description of the SAMs wettability in line with available experimental results. Furthermore, we provide an atomistic description of the first layers of water molecules at the solvent-SAM interface, which rationalizes the wettability difference between the cis- and trans-SAMs.

Influence of hydrogen bonding on the structure of the (001) corundum-water interface. Density functional theory calculations and Monte Carlo simulations

Density functional theory calculations and classical Monte Carlo simulations are applied to study the behavior of water in contact with a hydroxylated corundum (001) surface. Using DFT with periodic boundary conditions at T = 0 K, we systematically study the influence of the number of water molecules on the surface geometry and on the structure of the contact water layer. Only little effect of the thickness of the water layer on the geometry of the surface hydroxyl groups is observed. On the other hand, the molecules in the second layer have strong influence on the arrangement of water molecules in direct contact with the solid surface. In order to mimic macroscopic systems at room temperature, we perform inhomogeneous MC simulations of model corundum surface in contact with the water phase modeled by SPC/E model. The water molecules are classified according to their hydrogen-bonded partners into several groups. It is found that the preferential orientation of interfacial water molecules is primarily determined by the type of their hydrogen bonding. The hydroxyl groups at the corundum surface can serve as hydrogen bond donor or acceptor, depending on their orientation. No surface hydroxyls are found to coordinate two water molecules at the same time. On the other hand, water molecules coordinated by two different surface groups appear in MC simulations.

Automated optimization of water-water interaction parameters for a coarse-grained model

We have developed an automated parameter optimization software framework (ParOpt) that implements the Nelder-Mead simplex algorithm and applied it to a coarse-grained polarizable water model. The model employs a tabulated, modified Morse potential with decoupled short- and long-range interactions incorporating four water molecules per interaction site. Polarizability is introduced by the addition of a harmonic angle term defined among three charged points within each bead. The target function for parameter optimization was based on the experimental density, surface tension, electric field permittivity, and diffusion coefficient. The model was validated by comparison of statistical quantities with experimental observation. We found very good performance of the optimization procedure and good agreement of the model with experiment.

Influence of hydrophobic effects on streaming potential

We study the influence of hydrophobic effects on streaming potential mediated flow through a narrow confinement. In a clear departure from the approach used in prior works, we use a phase-field model to capture the hydrophobicity-induced depletion in the near wall region, and express the variation of viscosity and permittivity across the interfacial layer in terms of the phase-field variable. We then use these in the determination of the flow velocity, and highlight the sensitive interplay between the intrinsic length scale of the electrical double layer and that of the depletion in terms of the variations of an effective normalized viscosity that captures the electroviscous effect. We expect that this work will be an important step forward in the realistic continuum modeling of interfacial physics in the particular context of streaming potential mediated flows.

Solvation thermodynamics and heat capacity of polar and charged solutes in water

The solvation thermodynamics and in particular the solvation heat capacity of polar and charged solutes in water is studied using atomistic molecular dynamics simulations. As ionic solutes we consider a F(-) and a Na(+) ion, as an example for a polar molecule with vanishing net charge we take a SPC/E water molecule. The partial charges of all three solutes are varied in a wide range by a scaling factor. Using a recently introduced method for the accurate determination of the solvation free energy of polar solutes, we determine the free energy, entropy, enthalpy, and heat capacity of the three different solutes as a function of temperature and partial solute charge. We find that the sum of the solvation heat capacities of the Na(+) and F(-) ions is negative, in agreement with experimental observations, but our results uncover a pronounced difference in the heat capacity between positively and negatively charged groups. While the solvation heat capacity ΔC(p) stays positive and even increases slightly upon charging the Na(+) ion, it decreases upon charging the F(-) ion and becomes negative beyond an ion charge of q = -0.3e. On the other hand, the heat capacity of the overall charge-neutral polar solute derived from a SPC/E water molecule is positive for all charge scaling factors considered by us. This means that the heat capacity of a wide class of polar solutes with vanishing net charge is positive. The common ascription of negative heat capacities to polar chemical groups might arise from the neglect of non-additive interaction effects between polar and apolar groups. The reason behind this non-additivity is suggested to be related to the second solvation shell that significantly affects the solvation thermodynamics and due to its large spatial extent induces quite long-ranged interactions between solvated molecular parts and groups.

Interfacial entropy of water on rigid hydrophobic surfaces

A simple theoretical model is proposed for computing the interfacial entropy of water at rigid hydrophobic surfaces. The interfacial entropy, which is not considered in mean field models of static wettability, is evaluated from the fluctuations of the water-surface dispersion energy at the single particle level and represents the configurational bias imposed on the fluid molecules by the attractive external potential of a solid wall. A comparison with results obtained from molecular dynamics simulations shows that the model quantitatively describes the entropy loss of water when a water-vapor interface turns to water in contact with hydrophobic surfaces such as graphene, graphite, and diamond, while it overestimates this quantity on hydrophilic surfaces.

Molecular surface chemistry in marine bioadhesion

This review covers the in situ molecular physicochemical characterisation of bioadhesives at solid/liquid interfaces, with the aim of elucidating the adhesion strategies that lie at the root of marine biofouling. It focuses on three major foulers: mussels, algae and barnacles. The dispersal of these organisms, their colonisation of surfaces, and ultimately their survival rely critically on the ability of the organisms' larvae or spores to locate a favourable settlement site and undergo metamorphosis, thus initiating their sessile existence. Differences in the composition of adhesive secretions and the strategies employed for their temporary or permanent implementation exists between the larval and adult life stages. To date, only a few adhesive secretions from marine fouling organisms have been adequately described in terms of their chemical composition, and a survey revealed the presence of certain recurrent functional groups, specifically catechol, carboxylate, monoester-sulphate and -phosphate. This review will describe the binding modes of such functionalities to wet mineral/metal oxides surfaces. Such functionalities will be ranked based on their ability to bind to hydrophilic surfaces replacing surface-bound water (Langmuir adsorption constant) as well as other adsorbates (competitive adsorption). A plausible explanation for the propensity of the reviewed adhesive functionalities to bind to hydrous metal oxide surfaces will be given on the basis of the Hard and Soft Acids and Bases principle, Hofmeister effects and entropic considerations. From the in situ analysis of marine organism bioadhesives and adsorption studies of functionalities relevant to the bioadhesion process, insights can be gleaned for a knowledge-based innovation of antifouling strategies and the synthesis of strong, durable adhesive materials, which are suitable for implementation in wet environments.

Anionic and cationic Hofmeister effects on hydrophobic and hydrophilic surfaces

Using a two-step modeling approach, we address the full spectrum of direct, reversed, and altered ionic sequences as the charge of the ion, the charge of the surface, and the surface polarity are varied. From solvent-explicit molecular dynamics simulations, we extract single-ion surface interaction potentials for halide and alkali ions at hydrophilic and hydrophobic surfaces. These are used within Poisson-Boltzmann theory to calculate ion density and electrostatic potential distributions at mixed polar/unpolar surfaces for varying surface charge. The resulting interfacial tension increments agree quantitatively with experimental data and capture the Hofmeister series, especially the anomaly of lithium, which is difficult to obtain using continuum theory. Phase diagrams that feature different Hofmeister series as a function of surface charge, salt concentration, and surface polarity are constructed from the long-range force between two surfaces interacting across electrolyte solutions. Large anions such as iodide have a high hydrophobic surface affinity and increase the effective charge magnitude on negatively charged unpolar surfaces. Large cations such as cesium also have a large hydrophobic surface affinity and thereby compensate an external negative charge surface charge most efficiently, which explains the well-known asymmetry between cations and anions. On the hydrophilic surface, the size-dependence of the ion surface affinity is reversed, explaining the Hofmeister series reversal when comparing hydrophobic with hydrophilic surfaces.

Site-dipole field and vortices in confined water

We use molecular dynamics simulations to examine how the spatial patterns formed by the diffusive motion of water molecules are influenced by confinement between hydrophilic and hydrophobic surfaces. For bulk liquid water, Higo et al. found that the long-range orientational order of spatial dipole moments can form vortex-like spatial patterns for as long as 300 ps [Higo et al., Proc. Natl. Acad. Sci. USA 98, 5961 (2001)]. We perform a similar analysis for confined water and we find that the existence of vortices in these systems is dependent on both the surface separation and the surface hydrophilicity. Vortices perpendicular to the surface normal disappear when the surface separation is reduced to a thickness where the system is comprised of mostly interfacial water molecules. Vortices exist at slightly smaller separations for hydrophobic surfaces than for hydrophilic surfaces because the dipoles are less aligned at the hydrophobic surfaces. The dipole alignment that is induced by the hydrophilic surfaces is counter to the direction required by vortices.

Consistent description of electrohydrodynamics in narrow fluidic confinements in the presence of hydrophobic interactions

Electrohydrodynamics in the presence of hydrophobic interactions in narrow confinements is traditionally represented from a continuum viewpoint by a Navier slip-based conceptual paradigm, in which the slip length carries the sole burden of incorporating the effects of substrate wettability on interfacial electromechanics, precluding any explicit dependence of the interfacial potential distribution on the substrate wettability. Here we show that this traditional way of treating electrokinetics-wettability coupling may lead to serious discrepancies in predicting the resultant transport characteristics as manifested through an effective zeta potential. We suggest that an alternative consistent description of the underlying physics through a free-energy-based formalism, in conjunction with considerations of hydrodynamic and electrical property variations consistent with the pertinent phase-field description, may represent the underlying consequences in a more rational manner, as compared to the traditional slip-based model coupled with a two-layer description. Our studies further reveal that the above discrepancies may not occur solely due to the slip-based route of representing the interfacial wettability, but may be additionally attributed to the act of "discretizing" the interfacial phase fraction distribution through an artificial two-layer route.

Oxide/water interfaces: how the surface chemistry modifies interfacial water properties

The organization of water at the interface with silica and alumina oxides is analysed using density functional theory-based molecular dynamics simulation (DFT-MD). The interfacial hydrogen bonding is investigated in detail and related to the chemistry of the oxide surfaces by computing the surface charge density and acidity. We find that water molecules hydrogen-bonded to the surface have different orientations depending on the strength of the hydrogen bonds and use this observation to explain the features in the surface vibrational spectra measured by sum frequency generation spectroscopy. In particular, 'ice-like' and 'liquid-like' features in these spectra are interpreted as the result of hydrogen bonds of different strengths between surface silanols/aluminols and water.