Tuning the dynamics and molecular distribution of the self-spreading lipid bilayer

The self-spreading dynamics of lipid bilayers were investigated at controlled electrolyte concentrations. The self-spreading velocity increased when the concentration of NaCl was increased from 1 to 100 mM. Comparing the experimentally determined spreading energy with that estimated from theoretical models, we found that the self-spreading dynamics were well explained by considering the van der Waals interaction, double layer interaction and hydration interaction energies between the self-spreading bilayer and the substrate. The characteristic behavior at high concentration is attributable to the increase in the density of the lipid layer, originating from the effective shielding of the molecular charges by the electrolyte ions in solution. The distribution of doped dye-labeled molecule within the spreading bilayer was also controllable by tuning the electrolyte concentration. All of these findings were explained by systematic changes in bilayer-substrate or inter-molecular interactions depending on the electrolyte concentration.     

Controlling molecular diffusion in self-spreading lipid bilayer using periodic array of ultra-small metallic architecture on solid surface

Diffusion of target molecules incorporated in the self-spreading lipid bilayer was controlled by the introduction of periodic array of metallic architecture on solid surface. Retardation of the progress of target molecules became significant when the size of gap between small metal architectures was less than a few hundred nanometers. The self-spreading dynamics of the lipid bilayer depending on the size of the small gap were analyzed semiquantitatively. Estimated change in the driving force of the spreading layer suggests that highly localized compression of the spreading layer causes selective segregation of molecules.     

纳米结构表面浸润性质的分子动力学研究

采用分子动力学方法研究了氩纳米液滴在铂金属及其模型固体表面的浸润现象,获得了液滴在平滑表面和三角纳米结构阵列表面的接触角和展布特性.研究表明,液滴与壁面的势能作用较强时,液滴与纳米结构表面为均匀浸润,但是由于迟滞效应,接触角受表面纳米结构的影响不明显;势能作用较弱时,纳米结构间隙中存在类似蒸汽的低密度相,液滴与纳米结构表面为非均匀浸润,接触角受纳米结构的影响而增大;表面纳米结构可以使表面具有超疏水性.     

Molecular dynamics simulation of nanodroplet spreading enhanced by linear surfactants

We utilize molecular dynamics simulations to probe the surfactant-mediated spreading of a Lennard-Jones liquid droplet on a solid surface. The surfactants are linear hexamers that are insoluble in the liquid and reduce the surface tension of the liquid-vapor interface. We study how the interaction of the surfactant hexamers with the solid substrate influences spreading, as well as the dependence of spreading on surfactant concentration. We find that the spreading speed is strongly influenced by the attraction of the hydrophobic surfactant tail to the solid surface. When this attraction is sufficiently strong, surfactant molecules partition to the liquid-solid interface and facilitate spreading. This partitioning can lead to an inhomogeneous distribution of surfactant over the liquid-vapor interface, which could drive the Marangoni convection. We also observe that the surfactant molecules can assemble into micelles on the solid surface. The repulsion between micelles at the liquid-solid interface can lead to break-off and migration of the micelles from the liquid-solid to the gas-solid interface and spreading is facilitated in this way. Our model system contains features that are believed to underlie superspreading in experimental studies of droplet spreading.     

Properties of a water layer on hydrophilic and hydrophobic self-assembled monolayer surfaces: A molecular dynamics study

The microscopic behaviors of a water layer on different hydrophilic and hydrophobic surfaces of well ordered self-assembled monolayers (SAMs) are studied by molecular dynamics simulations. The SAMs consist of 18-carbon alkyl chains bound to a silicon(111) substrate, and the characteristic of its surface is tuned from hydrophobic to hydrophilic by using different terminal functional groups ( CH 3 , COOH). In the simulation, the properties of water membranes adjacent to the surfaces of SAMs were reported by comparing pure water in mobility, structure, and orientational ordering of water molecules. The results suggest that the mobility of water molecules adjacent to hydrophilic surface becomes weaker and the molecules have a better ordering. The distribution of hydrogen bonds indicates that the number of water-water hydrogen bonds per water molecule tends to be lower. However, the mobility of water molecules and distribution of hydrogen bonds of a water membrane in hydropho- bic system are nearly the same as those in pure water system. In addition, hydrogen bonds are mainly formed between the hydroxyl of the COOH group and water molecules in a hydrophilic system, which is helpful in understanding the structure of interfacial water.     

Molecular dynamics simulations of polyelectrolyte adsorption

We have performed molecular dynamics simulations of polyelectrolyte adsorption at oppositely charged surfaces from dilute polyelectrolyte solutions. In our simulations, polyelectrolytes were modeled by chains of charged Lennard-Jones particles with explicit counterions. We have studied the effects of the surface charge density, surface charge distribution, solvent quality for the polymer backbone, strength of the short-range interactions between polymers and substrates on the polymer surface coverage, and the thickness of the adsorbed layer. The polymer surface coverage monotonically increases with increasing surface charge density for almost all studied systems except for the system of hydrophilic polyelectrolytes adsorbing at hydrophilic surfaces. In this case the polymer surface coverage saturates at high surface charge densities. This is due to additional monomer-monomer repulsion between adsorbed polymer chains, which becomes important in dense polymeric layers. These interactions also preclude surface overcharging by hydrophilic polyelectrolytes at high surface charge densities. The thickness of the adsorbed layer shows monotonic dependence on the surface charge density for the systems of hydrophobic polyelectrolytes for both hydrophobic and hydrophilic surfaces. Thickness is a decreasing function of the surface charge density in the case of hydrophilic surfaces while it increases with the surface charge density for hydrophobic substrates. Qualitatively different behavior is observed for the thickness of the adsorbed layer of hydrophilic polyelectrolytes at hydrophilic surfaces. In this case, thickness first decreases with increasing surface charge density, then it begins to increase.     

Molecular modeling of interfaces between cellulose crystals and surrounding molecules: Effects of caprolactone surface grafting

A technical problem in cellulosic nanocomposite materials is the weak interaction between hydrophilic cellulose and hydrophobic polymer matrices. One approach to solve this difficulty is to chemically graft monomers of the matrix polymer onto the cellulose surface. An important question is to understand the effect such surface modification has on the interfacial properties. Semi-empirical approaches to estimate work of adhesion based on surface energies do not provide information on specific molecular interactions. Details about these interactions were obtained using molecular dynamics (MD) simulation. Cellulose interfaces with water and caprolactone medium were modeled with different amounts of grafted caprolactone. The modification lead to an increased work of adhesion between the surface and its surrounding medium. Furthermore, the MD simulations showed that the interaction between cellulose, both modified and non-modified, and surrounding medium is dominated by Coulomb interactions, predominantly as hydrogen bonds.     

An atomistic study of a poly(styrene sulfonate)/poly(diallyldimethylammonium) bilayer: the role of surface properties and charge reversal

In this work, we present the first molecular dynamics simulation on the formation of a polyelectrolyte bilayer resolved in atomistic detail, extending a previous study of the adsorption of poly(styrene sulfonate) [B. Qiao, J. J. Cerdà and C. Holm, Macromolecules, 2011, 44, 1707-1718.] to the formation of a poly(styrene sulfonate)/poly(diallyldimethylammonium) bilayer. Extensive molecular dynamics simulations of the adsorption process on different substrates (hydrophilic/hydrophobic, charged/neutral) were performed. Our results seem to indicate that a high enough surface charge density (0.164 C m(-2) here) may be required to achieve a multilayer linear growth in the framework of the electrostatic driven mechanism for PEM growth. We furthermore demonstrate that the use of stiff hydroxyl groups for creating a hydrophilic surface from a hydrophilic one can lead to severe simulation artifacts, and we discuss a simple remedy for this problem.     

Liquid–paper interactions during liquid drop impact and recoil on paper surfaces

The spreading and recoiling of water drops on several flat and macroscopically smooth model surfaces and on sized paper surfaces were studied over a range of drop impaction velocities using a high-speed CCD camera. The water drop spreading and recoiling results on several model hydrophobic and hydrophilic surfaces were found to be in agreement with observations reported in the literature. The maximum drop spreading diameter for those model surfaces at impact was found to be dependent upon the initial drop kinetic energy and the degree of hydrophobicity/hydrophilicity of the surface. The extent of the maximum drop recoiling was found to be much weaker for hydrophilic substrates than for hydrophobic substrates. Sized papers, however, showed an interesting switch of behaviour in the process of water drop impaction. They behave like a hydrophobic substrate when a water drop impacts on it, but like a hydrophilic substrate when water drop recoils. Although the contact angle between water and hydrophilic or hydrophobic non-porous surfaces changes from advancing to receding as reported in literature, the change of contact angle during water impact on paper surface is unique in that the level of sizing was found to have a smaller than expected influence on the degree of recoil. Atomic force microscopy (AFM) was used to probe fibres on a sized filter paper surface under water. The AFM data showed that water interacted strongly with the fibre even though the paper was heavily sized. Implications of this phenomenon were discussed in the context of inkjet print quality and of the surface conditions of sized papers. Results of this study are very useful in the understanding of inkjet ink droplet impaction on paper surfaces which sets the initial condition for ink penetration into paper after impaction.     

无规共聚物自组装球形胶束表面亲水性的耗散粒子动力学模拟

建立了含不同亲疏水粒子比的双亲性无规共聚物粗粒化模型. 采用耗散粒子动力学方法模拟了两亲性无规共聚物选择性溶剂自组装球形胶束表面的亲水性能. 模拟结果表明, 无规共聚物在选择性溶剂中自组装得到实心球形胶束, 球形胶束表面的亲水性与聚合物链亲水粒子含量、溶剂的选择性有关. 随着聚合物链所含亲水粒子增加, 球形胶束表面的亲水性增强. 球形胶束表面的亲水性随着疏水粒子与溶剂粒子间的排斥参数增大而增强, 模拟结果与实验结论一致. 该模拟方法给出的胶束微结构信息可以为双亲无规共聚物分子设计及自组装双亲胶束制备提供一定的理论指导.     

The hydrophobic effect: molecular dynamics simulations of water confined between extended hydrophobic and hydrophilic surfaces

Structural and dynamic properties of water confined between two parallel, extended, either hydrophobic or hydrophilic crystalline surfaces of n-alkane C(36)H(74) or n-alcohol C(35)H(71)OH, are studied by molecular dynamics simulations. Electron density profiles, directly compared with corresponding experimental data from x-ray reflectivity measurements, reveal a uniform weak de-wetting characteristic for the extended hydrophobic surface, while the hydrophilic surface is weakly wetted. These microscopic data are consistent with macroscopic contact angle measurements. Specific water orientation is present at both surfaces. The ordering is characteristically different between the surfaces and of longer range at the hydrophilic surface. Furthermore, the dynamic properties of water are different at the two surfaces and different from the bulk behavior. In particular, at the hydrophobic surface, time-correlation functions reveal that water molecules have characteristic diffusive behavior and orientational ordering due to the lack of hydrogen bonding interactions with the surface. These observations suggest that the altered dynamical properties of water in contact with extended hydrophobic surfaces together with a partial drying of the surfaces are more indicative of the hydrophobic effect than structural ordering, which we suggest to be independent of surface topology.     

Why are water-hydrophobic interfaces charged?

We report ab initio molecular dynamics simulations of hydroxide and hydronium ions near a hydrophobic interface, indicating that both ions behave like amphiphilic surfactants that stick to a hydrophobic hydrocarbon surface with their hydrophobic side. We show that this behavior originates from the asymmetry of the molecular charge distribution which makes one end of the ions strongly hydrophobic while the other end is even more hydrophilic than the regular water (H2O) molecules. The effect is more pronounced for the hydroxide than for the hydronium. Our results are consistent with several experimental observations and explain why hydrophobic surfaces in contact with water acquire a net negative charge, a phenomenon that has important implications for biology and polymer science.     

Kinetics of a polysoap collapse

We study the dynamics of collapse of a polysoap by means of large-scale molecular dynamics simulation and scaling arguments. A polysoap consists of a hydrophilic backbone and hydrophobic side chains attached at regular intervals along the backbone. In selective solvent conditions, the hydrophobic components aggregate, forcing the hydrophilic backbone to form loops anchored at the surface of the core, ultimately forming a micelle. The kinetics of polysoap collapse includes two major mechanisms: (1) early aggregation of the hydrophobic side chains controlled by first-order kinetics whose rate constant is given by a contact probability and (2) coalescence into larger clusters which requires activation to overcome energy barriers due to excluded volume repulsions between intermediate micelle coronas. In the late stage, the energy barrier is increasing as p(3/2), with p the number of aggregated side chains in an intermediate micelle. The corresponding late-stage rate constant decays exponentially as approximately exp(-p(3/2)).     

Hydrophobic-hydrophilic interactions of water with alkanethiolate chains from first-principles calculations.

The interaction of water with alkanethiolate chains is studied from first principles. A detailed analysis is performed by optimizing the structure of small water clusters, one-dimensional water chains, and ordered and disordered thin water layers adsorbed on hydroxyl(OH)- and methyl(CH3)-terminated alkanethiol monolayers. The hydrophilic/hydrophobic character of these two different substrates is investigated by means of an energetic analysis combined with hydrogen-bond counting. On the hydrophilic OH-terminated alkanethiol surface, zigzag, one-dimensional water chains and disordered thin water layers are the energetically favored structures. The ab initio results can be used to determine the optimal value of the empirical parameters characterizing a suitable force field to be used in classic molecular dynamics simulations.     

Origin of the Hydrophobic Behaviour of Hydrophilic CeO2

The nature of the hydrophobicity found in rare-earth oxides is intriguing. The CeO₂ (100) surface, despite its strongly hydrophilic nature, exhibits hydrophobic behaviour when immersed in water. In order to understand this puzzling and counter-intuitive effect we performed a detailed analysis of the confined water structure and dynamics. We report here an ab-initio molecular dynamics simulation (AIMD) study which demonstrates that the first adsorbed water layer, in immediate contact with the hydroxylated CeO₂ surface, generates a hydrophobic interface with respect to the rest of the liquid water. The hydrophobicity is manifested in several ways: a considerable diffusion enhancement of the confined liquid water as compared with bulk water at the same thermodynamic condition, a weak adhesion energy and few H-bonds above the hydrophobic water layer, which may also sustain a water droplet. These findings introduce a new concept in water/rare-earth oxide interfaces: hydrophobicity mediated by specific water patterns on a hydrophilic surface.     

Hydrophobic similarity between molecules: a MST-based hydrophobic similarity index

A similarity index based on the hydrophilic/hydrophobic properties of molecules is presented. Such an index is defined based on the fractional partition of the free energy of solvation developed within the framework of the self-consistent reaction field MST model, which divides the free energy of solvation or the free energy of transfer into contributions assigned to the surface elements defining the solute/solvent interface. These surface contributions can be integrated to derive atomic or group contributions. The suitability of the index to compute the molecular similarity based on hydrophobic/hydrophilic properties is examined by considering their application in a variety of test systems, including structure-activity relationships, absorption properties, and molecular recognition. The similarity index is expected to be a very powerful tool in molecular similarity studies for compounds of chemical, biochemical, and pharmaceutical interest.     

Transport of a liquid water and methanol mixture through carbon nanotubes under a chemical potential gradient

In this work, we report a dual-control-volume grand canonical molecular dynamics simulation study of the transport of a water and methanol mixture under a fixed concentration gradient through nanotubes of various diameters and surface chemistries. Methanol and water are selected as fluid molecules since water represents a strongly polar molecule while methanol is intermediate between nonpolar and strongly polar molecules. Carboxyl acid (-COOH) groups are anchored onto the inner wall of a carbon nanotube to alter the hydrophobic surface into a hydrophilic one. Results show that the transport of the mixture through hydrophilic tubes is faster than through hydrophobic nanotubes although the diffusion of the mixture is slower inside hydrophilic than hydrophobic pores due to a hydrogen network. Thus, the transport of the liquid mixture through the nanotubes is controlled by the pore entrance effect for which hydrogen bonding plays an important role.     

丙烯酸水溶液聚合法制备PTFE/PAA/Nafion膜及其性能

采用廉价的多孔聚四氟乙烯(PTFE)膜作为基底, 用少量的Nafion与PTFE膜复合可制备低成本的质子膜. 但疏水性的PTFE膜与亲水性的Nafion膜结合性不佳. 基于此, 本文对疏水性的PTFE膜材料表面进行设计, 先采用丙烯酸对疏水性的PTFE膜表面进行亲水性改性, 再喷涂亲水性Nafion膜, 完成低成本PTFE/PAA/Nafion膜的制备. 实验结果表明, 改性前的PTFE膜材料水接触角为150°, 改性后的膜接触角变为55.6°, 亲水性大幅上升, 膜的机械强度和尺寸稳定性(断裂强度为25.2 MPa, 80 ℃下的溶胀率为11.9%)均优于Nafion117膜, 而 Nafion用量则节省了60%. PTFE/PAA/Nafion膜具有高质子导通率(80 ℃下达到131.9 mS/cm), 接近于Nafion117膜, 最大功率密度可以达到404.2 mW/cm².     

Vibrational density of states of hydration water at biomolecular sites: hydrophobicity promotes low density amorphous ice behavior

Inelastic neutron scattering experiments and molecular dynamics simulations have been used to investigate the low frequency modes, in the region between 0 and 100 meV, of hydration water in selected hydrophilic and hydrophobic biomolecules. The results show changes in the plasticity of the hydrogen-bond network of hydration water molecules depending on the biomolecular site. At 200 K, the measured low frequency density of states of hydration water molecules of hydrophilic peptides is remarkably similar to that of high density amorphous ice, whereas, for hydrophobic biomolecules, it is comparable to that of low density amorphous ice behavior. In both hydrophilic and hydrophobic biomolecules, the high frequency modes show a blue shift of the libration mode as compared to the room temperature data. These results can be related to the density of water molecules around the biological interface, suggesting that the apparent local density of water is larger in a hydrophilic environment.