Efficient tuning of potential parameters for liquid–solid interactions

Spherical and cylindrical water droplets on silicon surface are studied to tune the silicon–oxygen interaction. We use molecular dynamics simulations to estimate the contact angle of two different shaped droplets. We found that the cylindrical droplets are independent of the line tension as their three phases curvature is equal zero. Additionally, we compare an analytical model, taking into account or not the Tolman length and we show that for spherical small size droplets, this length is important to be included, in contrast to cylindrical droplets in which the influence of the Tolman length is negligible. We demonstrate that the usual convenient way to exclude linear tension in the general case can give wrong results. Here, we consider cylindrical droplets, since their contact angle does not depend on the droplet size in the range of few to 10ths of nanometres. The droplets are stabilised due to the periodic boundary conditions. This allows us to propose a new parameterisation for nanoscale droplets, which is independent the size of the droplets or its shape, minimising at the same time the calculation procedure. With the proposed methodology, we can extract the epsilon parameter of the interaction potential between a liquid and a solid from the nanoscaled molecular simulation with only as input the macrosized experimental wetting angle for a given temperature.     

Wetting characterisation of silicon (1,0,0) surface

Water droplets on bare silicon surfaces are studied to examine the wetting behaviour as a function of the surface energy and to parameterise water–silicon interactions in order to recover the hydrophobic behaviour measured by experiments. Two different wetting regimes characterised by a critical interaction strength value are observed. At a threshold value of the water–silicon interaction parameter, water molecules start penetrating into the first layer of silicon surface under thermally vibrating walls, resulting in two distinct wetting behaviours. Fixed (cold) silicon walls do not exhibit the two different wetting characteristics. Size effects are studied for nano-scale droplets, and line tension influence is observed depending on the surface wettability. Decrease in the droplet size increases the contact angle values for the low wetting cases, while contact angles decrease for smaller droplets on the high wetting surfaces. Considering the line tension effects and droplet size, ?_Si–O for water–silicon interactions to recover the hydrophobic behaviour of silicon surfaces is estimated to be 12.5% of the value predicted using the Lorentz–Berthelot mixing rule.     

Nanoscale droplet vaporisation by molecular dynamics

Nanoscale droplet vaporisation was studied by molecular dynamics, which allows the calculation of properties for droplets statistically without considering the discontinuous interface between a liquid droplet and surrounding gas. An argon droplet was created and immersed inside its vapour. After equilibration, the periphery of the system was heated by a carrier gas to vaporise the droplet. Replications were conducted to check the variation in the phenomenon. Thermodynamic properties such as the density, pressure and temperature profiles were sampled for each interval. The evolution of the surface tension of the droplet undergoing vaporisation was investigated. Moreover, the vaporisation rate of nanodroplets was compared with the kinetic theory-based Hertz–Knudsen–Langmuir equation and two diffusion-based models, which are the D² evaporation law and Kincaid and Longley model [Kincaid DC, Longley TS. A water droplet evaporation and temperature model. Trans ASAE. 1989; 32(2):457–463]. The kinetic model underestimates the vaporisation rate by one order of magnitude whereas the two diffusion-based models overestimate the rate by one order of magnitude.     

Wettability of polymers by mucin aqueous solutions

The wettability of poly(methyl methacrylate) and polyethylene by water and aqueous mucin solutions have been studied by sessile drop and under-water captive air bubble contact angles, respectively. From the sessile drop and octane under-water contact angles the polymer-water interfaces have been characterized in terms of works of adhesion and acid-base (polar) interactions. A large water-air contact angle hysteresis observed with poly(methyl methacrylate) surfaces has been attributed to side-chain beta relaxations of polymer ester methyl groups. The wettabilities of the polymers by mucin aqueous solutions have been studied as a function of protein concentration and related to the surface tensions. A positive slope of adhesion tension vs surface tension line, characteristic of polar surfaces, was found with poly(methyl methacrylate). By contrast, a change in the slope, explained as a change in mucin relative adsorption densities at solid/liquid and solid/vapour interfaces, was observed with polyethylene. This adhesion tension behavior appeared to be in agreement with previous data we have published concerning the quantity and state of mucin which are adsorbed to polymers characterized by different surface properties.     

Evaluation of the contact angle from molecular simulations

An alternative method for determination of the contact angle of droplets at a solid underlay from molecular simulations is proposed. The method is based on a recently developed general method of identification the surface molecules of a molecular system with the interface of an arbitrary shape and on a subsequent parametrisation of the surface of the droplet by a smooth function. The method has been verified first by considering two artificial systems with the exactly known contact angles and then by comparison with literature data for two realistic systems.     

Calculation of the surface tension of planar interfaces by molecular simulations: from Lennard-Jones fluids to binary mixtures

The relative longevity of the research in the field of the molecular simulations of the liquid–vapour interfaces of Lennard-Jones (LJ) particles can be explained by the dependence of the surface tension on many methodological factors. After a few illustrations on the parameters that can impact the results of surface tension on the LJ interfaces, we establish the ability of the current methodologies to quantitatively predict the surface tension of various liquid–vapour interfaces of pure components at different temperatures. We also show that the methods perform very well for the reproduction of the interfacial tension of binary mixtures in a wide range of pressures.     

Fouling of nanostructured insect cuticle: adhesion of natural and artificial contaminants

The adhesional properties of contaminating particles of scales of various lengths were investigated for a wide range of micro- and nanostructured insect wing cuticles. The contaminating particles consisted of artificial hydrophilic (silica) and spherical hydrophobic (C₁₈) particles, and natural pollen grains. Insect wing cuticle architectures with an open micro-/nanostructure framework demonstrated topographies for minimising solid–solid and solid–liquid contact areas. Such structuring of the wing membranes allows for a variety of removal mechanisms to contend with particle contact, such as wind and self-cleaning droplet interactions. Cuticles exhibiting high contact angles showed considerably lower particle adhesional forces than more hydrophilic insect surfaces. Values as low as 3 nN were recorded in air for silica of ～28 nm in diameter and <25 nN for silica particles 30 μm in diameter. A similar adhesional trend was also observed for contact with pollen particles.     

The role of nano-roughness in antifouling

Nano-engineered superhydrophobic surfaces have been investigated for potential fouling resistance properties. Integrating hydrophobic materials with nanoscale roughness generates surfaces with superhydrophobicity that have water contact angles (θ) >150° and concomitant low hysteresis (<10°). Three superhydrophobic coatings (SHCs) differing in their chemical composition and architecture were tested against major fouling species (Amphora sp., Ulva rigida, Polysiphonia sphaerocarpa, Bugula neritina, Amphibalanus amphitrite) in settlement assays. The SHC which had nanoscale roughness alone (SHC 3) deterred the settlement of all the tested fouling organisms, compared to selective settlement on the SHCs with nano- and micro-scale architectures. The presence of air incursions or nanobubbles at the interface of the SHCs when immersed was characterized using small angle X-ray scattering, a technique sensitive to local changes in electron density contrast resulting from partial or complete wetting of a rough interface. The coating with broad spectrum antifouling properties (SHC 3) had a noticeably larger amount of unwetted interface when immersed, likely due to the comparatively high work of adhesion (60.77 mJ m^?2 for SHC 3 compared to 5.78 mJ m^?2 for the other two SHCs) required for creating solid/liquid interface from the solid/vapour interface. This is the first example of a non-toxic, fouling resistant surface against a broad spectrum of fouling organisms ranging from plant cells and non-motile spores, to complex invertebrate larvae with highly selective sensory mechanisms. The only physical property differentiating the immersed surfaces is the nano-architectured roughness which supports longer standing air incursions providing a novel non-toxic broad spectrum mechanism for the prevention of biofouling.     

Towards automatic measurement of anteversion and neck–shaft angles in human femurs using CT images

Automatic assessment of human femur morphology may provide useful clinical information with regard to hip and knee surgery, prosthesis design and management of hip instability. To this end, neck–shaft and anteversion angles are usually used. We propose a full automatic method to estimate these angles in human femurs. Multislice CT images from 18 dried bones were analysed. The algorithm fits 3D cylinders to different regions of the bone to estimate the angles. A manual segmentation and a conventional angle assessment were used for validation. We found anteversion angle as 20 ± 7° and neck–shaft angle as 130 ± 9°. Mean distances from femur surface to cylinders were 5.5 ± 0.6, 3.5 ± 0.6 and 2.4 ± 0.4 mm for condyles, diaphysis and neck regions, respectively. Automatic and conventional angles were positively correlated (r²>0.85). Manual and automatic segmentations did not differ. The method was fast and 100% reproducible. A robust in vivo segmentation algorithm should be integrated to advance towards a clinically compliant methodology.     

网粒体超疏水性在茶小绿叶蝉防杀虫剂渗透中的屏障效应

【目的】茶小绿叶蝉Empoasca onukii体表覆盖的网粒体具有超疏水性，杀虫剂喷雾触碰虫体后药滴动态是否受网粒体影响尚未完全清楚。本研究旨在明确网粒体在茶小绿叶蝉成虫抵御杀虫剂雾滴渗透的屏障作用。【方法】以罗丹明B(RhB)作为指示剂添加到测试的杀虫剂(联苯菊酯和茚虫威)中，利用可拍照显微镜观察记录联苯菊酯(1.25 mg/L和0.05 mg/L)和茚虫威(0.006 mg/L和0.0009 mg/L)喷雾处理后24 h，茶小绿叶蝉成虫翅面药滴滚落、蒸发、被抖动扫除等行为动态，分析翅面药滴大小与蒸发后固化形态的关系；测定网粒体移除前后药滴与翅面的接触角，统计不同疏水性翅面上的网粒体分布密度；收集并利用扫描电镜分析叶蝉体表抖落的药滴及药剂颗粒是否含有网粒体，同时观察网粒体与翅面残留溶质接触的显微形态。【结果】药滴动态观察显示，圆球状药滴在茶小绿叶蝉成虫翅面不会自行滚落，72.0%成虫静止等待翅面药滴蒸发，蒸发后形成药剂颗粒或不规则药斑与药滴大小无关，而与虫体翅面的疏水类型有关，蒸发后24 h内翅面的药剂颗粒都被叶蝉抖动扫除；在叶蝉疏水性强翅面上，药滴的静态接触角为141.63±8.06°，药滴蒸发后形成药剂颗粒，网粒体分布密度为6.1±1.2粒/μm²，而疏水性弱的翅面上药滴蒸发形成药斑，网粒体分布密度为2.2±0.9粒/μm²；SEM图片显示被茶小绿叶蝉抖落的药滴和药剂颗粒表面均带有网粒体，药斑和药剂颗粒的显微结构显示网粒体出现聚集并与残留溶质相融合。【结论】超疏水性网粒体的均匀分布决定药滴触碰茶小绿叶蝉成虫翅面后形成圆球状，网粒体的亲油性及团聚性促使药滴蒸发后形成药剂颗粒，网粒体的脱落性使药剂颗粒可被茶小绿叶蝉成虫抖动扫除。     

Wetting control through topography and surface hydrophilic/hydrophobic property changes by coarse grained simulation

Abstract

The changes of wetting state of water droplet on the solid surface featuring pillared structures are quantitatively studied by Coarse Grained simulation. Our results demonstrate that wetting state changes with the different topography (surface roughness), and it depends on the intrinsic hydrophilic/hydrophobic property of surface as well. Only if the contact angle of water droplet on the smooth surface is larger than 93.13°, the wetting state translates from the Wenzel state to the Cassie state on the rough surface with certain pillar height and width, and the contact angle climb up to the highest point and then remain almost unchanged with the increasing of pillar height and the same pillar distance. However, the wetting state does not change if the contact angle on the smooth surface is 85.1° or less, no matter what pillar structure the surface has. Additionally, the contact angles will remain almost unchanged if the pillar height is higher than a certain value. Our simulation results provide a quantitative understanding about the wetting state of water droplet on solid rough surfaces, and the results show the wetting state can be controlled by combining rough structure design and hydrophilic/hydrophobic property change of surfaces.     

Multi-body optimization with subject-specific knee models: performance at high knee flexion angles

When estimating knee kinematics from skin markers and stereophotogrammetry, multi-body optimization (MBO) has provided promising results for reducing soft tissue artefacts (STA), but can still be improved. The goal of this study was to assess the performance of MBO with subject-specific knee models at high knee flexion angles (up to 110°) against knee joint kinematics measured by magnetic resonance imaging. Eight subjects were recruited. MBO with subject-specific knee models was more effective in compensating STA compared to no kinematic and spherical constraints, in particular for joint displacements. Moreover, it seems to be more reliable over large ranges of knee flexion angle. The ranges of root mean square errors for knee rotations/displacements were 3.0°–9.2°/1.3–3.5 mm for subject-specific knee models, 6.8°–8.7°/6.0–12.4 mm without kinematic constraint and 7.1°–9.8°/4.9–12.5 mm for spherical constraints.     

Fluctuations and mechanical strength of α-helices of polyglycine and poly(L-alanine)

By using the static correlations of fluctuations in the dihedral angles of the α-helices of polyglycine and poly(L -alanine) calculated previously, geometrical fluctuations of a section (consisting of up to 18 peptide units) of the α-helices of infinite length are calculated. These fluctuations are found to differ in some respects (e.g., the dependence of amplitudes on the length of section) from those of a circular rod made of homogeneous continuous material. However, the moduli of the mechanical strengths (tensile Young's modulus, bending Young's modulus, and the shear modulus) of a circular rod are calculated, whose geometrical fluctuations are approximately equal to the fluctuations of a section consisting of 18 peptide units. They are of the order of 10¹¹ dyn/cm². The tensile rigidity, flexural rigidity, and torsional rigidity are calculated to be 1.20 × 10^?3 dyn, 2.46 × 10^?19 dyn·cm² and 1.79 × 10^?19 dyn·cm² for polyglycine, and 1.96 × 10^?3 dyn, 4.05 × 10^?19 dyn·cm² and 3.28 × 10^?19 dyn·cm² for poly(L -alanine), respectively.     

Molecular dynamics simulation for shape change of water-in-oil droplets

We performed molecular dynamics (MD) simulations of water-in-oil droplet shape transformations induced by the addition of polymer chains. In a prior experiment, transformations of spherical droplets to rod-like, worm-like and network-like droplets were observed. In our previous study, we reproduced rod-like droplets via coarse-grained MD simulations, and the mechanism for the droplet shape change was elucidated by considering the contact area between the chains and the surfactant head groups. However, in that simulation model, we could not reproduce the worm-like and network-like droplets. In this study, we improved the simulation model. For a small number of chains, several spherical droplets were obtained. As the number of chains increased, the spherical droplets were transformed to rod-like, worm-like and network-like shapes by coalescence of the droplets. The calculated and experimental results agreed well, and we verified that the mechanism for the droplet shape transformations observed in the present simulations could be explained by the mechanism suggested in the previous study.     

Wetting transition of nanodroplets of water on textured surfaces: a molecular dynamics study

The crossover behaviour of water droplet's state from the Wenzel state to the Cassie state with varying pillar height and surface fraction is examined critically using molecular dynamics. We report the effect of the system size on the wetting behaviour of water droplets by examining the contact angle for both regimes. We observe that when the droplet size is comparable to the pillar dimension, the contact angle of droplets fluctuates with increasing droplet size because of the contact line pinning, which is more pronounced in the Wenzel regime. We further demonstrate the phantom-wall method to evaluate free energy of intermediate wetting states.     

Growth kinetics of microalgae in microfluidic static droplet arrays

We investigated growth kinetics of microalgae, Chlorella vulgaris, in immobilized arrays of nanoliter‐scale microfluidic drops. These static drop arrays enabled simultaneous monitoring of growth of single as well as multiple cells encapsulated in individual droplets. To monitor the growth, individual drop volumes were kept nearly intact for more than a month by controlling the permeation of water in and out of the microfluidic device. The kinetic growth parameters were quantified by counting the increase in the number of cells in each drop over time. In addition to determining the kinetic parameters, the cell‐size distribution of the microalgae was correlated with different stages of the growth. The single‐cell growth kinetics of C. vulgaris showed significant heterogeneity. The specific growth rate ranged from 0.55 to 1.52 day^?1 for different single cells grown in the same microfluidic device. In comparison, the specific growth rate in bulk‐scale experiment was 1.12 day^?1. It was found that the average cell size changes significantly at different stages of the cell growth. The mean cell‐size increased from 5.99 ± 1.08 to 7.33 ± 1.3 µm from exponential to stationary growth phase. In particular, when multiple cells are grown in individual drops, we find that in the stationary growth phase, the cell size increases with the age of cell suggesting enhanced accumulation of fatty acids in older cells. Biotechnol. Bioeng. 2012; 109: 2987–2996. © 2012 Wiley Periodicals, Inc.     

Surface tension measurements on filamentous broths

Summary The surface tension of microorganisms is an important parameter in biotechnological processes. In this work contact angles and surface tensions were determined for filamentous micro-organisms Streptomyces levoris and Aspergillus niger. The contact angles were 39.7° and 43.1°, and the surface tensions 59.0 mJ/m² and 57.3 mJ/m², respectively. The contact angles for the filamentous broths were measured with the axisymmetric drop shape analysis - contact diameter (ADSA-CD) method which is particularly well suited to rough surfaces. The contact angle results obtained were significantly larger than those earlier reported for bacteria.Visiting scientist from Helsinki University of Technology, Espoo, Finland.     

Interrelationship between water film thicknesses and contact angles and a model for CO2 adhesion

Thermodynamics of contact angle phenomena is strongly affected by the presence of thin liquid films. However, at present, studies for CO₂/brine/mineral systems only consider the films apart from contact angles. In this paper, molecular dynamics (MD) simulations have been performed to simultaneously investigate the interrelationship between water film thicknesses and water contact angles. Two types of contact angles were considered namely Young’s contact angle (no water film is present) and contact angle with film (a stable film is present). The results showed that as Young’s contact angle increased, film thickness decreased which leading to increasing of contact angle with film. The effects of CO₂-mineral pre-contact have also been investigated and it has been found that on mediate hydrophilic surfaces (Q³), water films were present when CO₂ droplets were placed above the surfaces, however, water films were absent when CO₂ droplets directly contact with the surfaces. This phenomenon implies that water films on mineral surfaces have a possibility to rupture and a film rupture mechanism for CO₂ adhesion on hydrated mineral surfaces was proposed. These results may provide new information on interactions among CO₂, water/brine and mineral to better understand the behaviour of CO₂ during geologic sequestration.     

Molecular dynamics simulation of the spreading of the nanosized droplet on a graphene-coated substrate: the effect of the contact line forces

Surface characteristics of graphene have an important impact on its performance. Substantial attention has been devoted to studying the static wetting behaviour of a graphene-coated substrate with little attention to the dynamic wetting behaviour. The impact of contact line forces (CLFs) on the droplet-spreading process has not been revealed completely yet. A series of molecular dynamics (MD) simulation is performed to investigate the spreading process of the water droplet on the graphene-coated substrate in this research. The increase of interaction potential parameter between substrate and water droplet makes the spreading radius getting bigger and the final static contact angle smaller. Apart from that, the higher hydrophilicity of underlying substrate can lead to the greater averaged forces of atoms near contact line. CLFs correlate well with the variation of kinetic energy of water molecules located in the contact line region. Surface tensions of water droplets on graphene-coated substrates are also examined. The liquid-vapour and solid-vapour surface tensions are constant. An increase in the surface tension of liquid-solid lead to the increase of balanced contact angles of water on the substrate. The results are useful for understanding the effect of CLFs on the dissipation of kinetic energy of water molecules.     

Determination of contact angle by molecular simulation using number and atomic density contours

The contact angles of Lennard-Jones fluid droplets on a structureless solid surface, simulated using Monte Carlo simulation, are calculated by fitting isochoric surfaces and making a number of assumptions about the droplet. The results show that there are significant uncertainties in the calculated contact angles due to the choice of these assumptions, such as the grid size used in tracking the isochoric density profile, the omission of isochoric data points near the surface and the function used to fit the isochoric profile. In this study, we propose a new method of calculating density contours based on atomic density instead of number density. This method results in a much smaller variation in contact angle when applying different assumptions than using number density for isochoric contours. The most consistent results, across a range of assumptions about the droplet and the contact angle, come from averaging the contact angle from several isochoric density profiles. In addition, this gives a measurement of the variation due to the choice of isochoric density.