Kenji Yasuoka

Coexistence and transition between Cassie and Wenzel state on pillared hydrophobic surface

Water droplets on rugged hydrophobic surfaces typically exhibit one of the following two states: (i) the Wenzel state [Wenzel RN (1936) Ind Eng Chem 28:988–994] in which water droplets are in full contact with the rugged surface (referred as the wetted contact) or (ii) the Cassie state [Cassie, ABD, Baxter S (1944) Trans Faraday Soc 40:546–551] in which water droplets are in contact with peaks of the rugged surface as well as the “air pockets” trapped between surface grooves (the composite contact). Here, we show large-scale molecular dynamics simulation of transition between Wenzel state and Cassie state of water droplets on a periodic nanopillared hydrophobic surface. Physical conditions that can strongly affect the transition include the height of nanopillars, the spacing between pillars, the intrinsic contact angle, and the impinging velocity of water nanodroplet (“raining” simulation). There exists a critical pillar height beyond which water droplets on the pillared surface can be either in the Wenzel state or in the Cassie state, depending on their initial location. The free-energy barrier separating the Wenzel and Cassie state was computed on the basis of a statistical-mechanics method and kinetic raining simulation. The barrier ranges from a few tenths of kBT0 (where kB is the Boltzmann constant, and T0 is the ambient temperature) for a rugged surface at the critical pillar height to ≈8 kBT0 for the surface with pillar height greater than the length scale of water droplets. For a highly rugged surface, the barrier from the Wenzel-to-Cassie state is much higher than from Cassie-to-Wenzel state. Hence, once a droplet is trapped deeply inside the grooves, it would be much harder to relocate on top of high pillars.

42 TFlops hierarchical N -body simulations on GPUs with applications in both astrophysics and turbulence

As an entry for the 2009 Gordon Bell price/performance prize, we present the results of two different hierarchical N-body simulations on a cluster of 256 graphics processing units (GPUs). Unlike many previous N-body simulations on GPUs that scale as O(N2), the present method calculates the O(N log N) treecode and O(N) fast multipole method (FMM) on the GPUs with unprecedented efficiency. We demonstrate the performance of our method by choosing one standard application -a gravitational N-body simulation- and one non-standard application -simulation of turbulence using vortex particles. The gravitational simulation using the treecode with 1,608,044,129 particles showed a sustained performance of 42.15 TFlops. The vortex particle simulation of homogeneous isotropic turbulence using the periodic FMM with 16,777,216 particles showed a sustained performance of 20.2 TFlops. The overall cost of the hardware was 228,912 dollars. The maximum corrected performance is 28.1TFlops for the gravitational simulation, which results in a cost performance of 124 MFlops/

Measurement of contact-angle hysteresis for droplets on nanopillared surface and in the Cassie and Wenzel states: A molecular dynamics simulation study

. This correction is performed by counting the Flops based on the most efficient CPU algorithm. Any extra Flops that arise from the GPU implementation and parameter differences are not included in the 124 MFlops/

Self-assembly of surfactants and polymorphic transition in nanotubes.

.

Molecular Mechanisms of How Mercury Inhibits Water Permeation through Aquaporin-1: Understanding by Molecular Dynamics Simulation

We perform large-scale molecular dynamics simulations to measure the contact-angle hysteresis for a nanodroplet of water placed on a nanopillared surface. The water droplet can be in either the Cassie state (droplet being on top of the nanopillared surface) or the Wenzel state (droplet being in contact with the bottom of nanopillar grooves). To measure the contact-angle hysteresis in a quantitative fashion, the molecular dynamics simulation is designed such that the number of water molecules in the droplets can be systematically varied, but the number of base nanopillars that are in direct contact with the droplets is fixed. We find that the contact-angle hysteresis for the droplet in the Cassie state is weaker than that in the Wenzel state. This conclusion is consistent with the experimental observation. We also test a different definition of the contact-angle hysteresis, which can be extended to estimate hysteresis between the Cassie and Wenzel state. The idea is motivated from the appearance of the hysteresis loop typically seen in computer simulation of the first-order phase transition, which stems from the metastability of a system in different thermodynamic states. Since the initial shape of the droplet can be controlled arbitrarily in the computer simulation, the number of base nanopillars that are in contact with the droplet can be controlled as well. We show that the measured contact-angle hysteresis according to the second definition is indeed very sensitive to the initial shape of the droplet. Nevertheless, the contact-angle hystereses measured based on the conventional and new definition seem converging in the large droplet limit.

Water proton configurations in structures I, II, and H clathrate hydrate unit cells

We study self-assembly and polymorphic transitions of surfactant molecules in water within a nanotube and the effect of water-nanotube interactions on the self-assembly morphologies. We present a simulation evidence of a cornucopia of polymorphic structures of surfactant assemblies--many of which have not been observed in bulk solutions--through adjusting the water-nanotube chemical interactions which range from hydrophilic to hydroneutral and to hydrophobic. The ability to control the morphologies of surfactant assemblies within nanoscale confinement can be used for patterning the interior surface of nanochannels for application in nanofluidics and nanomedical devices.

Fast multipole methods on a cluster of GPUs for the meshless simulation of turbulence

Aquaporin (AQP) functions as a water-conducting pore. Mercury inhibits the water permeation through AQP. Although site-directed mutagenesis has shown that mercury binds to Cys189 during the inhibition process, it is not fully understood how this inhibits the water permeation through AQP1. We carried out 40 ns molecular dynamics simulations of bovine AQP1 tetramer with mercury (Hg-AQP1) or without mercury (Free AQP1). In Hg-AQP1, Cys191 (Cys189 in human AQP1) is converted to Cys-SHg+ in each monomer. During each last 10 ns, we observed water permeation events occurred 23 times in Free AQP1 and never in Hg-AQP1. Mercury binding did not influence the whole structure, but did induce a collapse in the orientation of several residues at the ar/R region. In Free AQP1, backbone oxygen atoms of Gly190, Cys191, and Gly192 lined, and were oriented to, the surface of the water pore channel. In Hg-AQP1, however, the SHg+ of Cys191-SHg+ was oriented toward the outside of the water pore. As a result, the backbone oxygen atoms of Gly190, Cys191, and Gly192 became disorganized and the ar/R region collapsed, thereby obstructing the permeation of water. We suggest that mercury disrupts the water pore of AQP1 through local conformational changes in the ar/R region.

Origin of subdiffusion of water molecules on cell membrane surfaces

Position and orientation of water protons need to be specified when the molecular simulation studies are performed for clathrate hydrates. Positions of oxygen atoms in water are experimentally determined by X-ray diffraction analysis of clathrate hydrate structures, but positions of water hydrogen atoms in the lattice are disordered. This study reports a determination of the water proton coordinates in unit cell of structure I (sI), II (sII), and H (sH) clathrate hydrates that satisfy the ice rules, have the lowest potential energy configuration for the protons, and give a net zero dipole moment. Possible proton coordinates in the unit cell were chosen by analyzing the symmetry of protons on the hexagonal or pentagonal faces in the hydrate cages and generating all possible proton distributions which satisfy the ice rules. We found that in the sI and sII unit cells, proton distributions with small net dipole moments have fairly narrow potential energy spreads of about 1 kJ∕mol. The total Coulomb potential on a test unit charge placed in the cage center for the minimum energy∕minimum dipole unit cell configurations was calculated. In the sI small cages, the Coulomb potential energy spread in each class of cage is less than 0.1 kJ∕mol, while the potential energy spread increases to values up to 6 kJ∕mol in sH and 15 kJ∕mol in the sII cages. The guest environments inside the cages can therefore be substantially different in the sII case. Cartesian coordinates for oxygen and hydrogen atoms in the sI, sII, and sH unit cells are reported for reference.

DS-CUDA: A Middleware to Use Many GPUs in the Cloud Environment

Recent advances in the parallelizability of fast N-body algorithms, and the programmability of graphics processing units (GPUs) have opened a new path for particle based simulations. For the simulation of turbulence, vortex methods can now be considered as an interesting alternative to finite difference and spectral methods. The present study focuses on the efficient implementation of the fast multipole method and pseudo-particle method on a cluster of NVIDIA GeForce 8800 GT GPUs, and applies this to a vortex method calculation of homogeneous isotropic turbulence. The results of the present vortex method agree quantitatively with that of the reference calculation using a spectral method. We achieved a maximum speed of 7.48 TFlops using 64 GPUs, and the cost performance was near

Cutoff radius effect of the isotropic periodic sum and Wolf method in liquid-vapor interfaces of water.

9.4/GFlops. The calculation of the present vortex method on 64 GPUs took 4120 s, while the spectral method on 32 CPUs took 4910 s.