Pore nucleation in mechanically stretched bilayer membranes

We report a computer-simulation study of the free-energy barrier for the nucleation of pores in the bilayer membrane under constant stretching lateral pressure. We find that incipient pores are hydrophobic but as the lateral size of the pore nucleus becomes comparable with the molecular length, the pore becomes hydrophilic. In agreement with previous investigations, we find that the dynamical process of growth and closure of hydrophilic pores is controlled by the competition between the surface tension of the membrane and the line tension associated with the rim of the pore. We estimate the line tension of a hydrophilic pore from the shape of the computed free-energy barriers. The line tension thus computed is in a good agreement with available experimental data. We also estimate the line tension of hydrophobic pores at both macroscopic and microscopic levels. The comparison of line tensions at these two different levels indicates that the "microscopic" line tension should be carefully distinguished from the "macroscopic" effective line tension used in the theoretical analysis of pore nucleation. The overall shape of the free-energy barrier for pore nucleation shows no indication for the existence of a metastable intermediate during pore nucleation.     

Nucleation free energy of pore formation in an amphiphilic bilayer studied by molecular dynamics simulations

The formation of a pore in a membrane requires a considerable rearrangement of the amphiphilic molecules about to form the bilayer edge surrounding the pore, and hence is accompanied by a steep increase of the free energy. Recent rupture and conductance experiments suggest that this reshuffling process is also responsible for a small energy barrier that stabilizes "prepores" with diameters of less than 1 nm, rendering both the opening and closing of pores an activated process. We use the potential of mean constraint force method to study this free energy profile, as a function of pore radius, in a coarse grained bilayer model. The calculations show that the free energy rises by (15-20) kT during pore opening, making it an extremely rare nucleation event. Although we do not observe a barrier to pore closure, the results do make the existence of such a barrier plausible. For larger pores we find a smooth transition to Litster's model, from which a line tension coefficient of about 3.7 x 10(-11) J m(-1) is deduced.     

Pores in bilayer membranes of amphiphilic molecules: coarse-grained molecular dynamics simulations compared with simple mesoscopic models

We investigate pores in fluid membranes by molecular dynamics simulations of an amphiphile-solvent mixture, using a molecular coarse-grained model. The amphiphilic membranes self-assemble into a lamellar stack of amphiphilic bilayers separated by solvent layers. We focus on the particular case of tensionless membranes, in which pores spontaneously appear because of thermal fluctuations. Their spatial distribution is similar to that of a random set of repulsive hard disks. The size and shape distribution of individual pores can be described satisfactorily by a simple mesoscopic model, which accounts only for a pore independent core energy and a line tension penalty at the pore edges. In particular, the pores are not circular: their shapes are fractal and have the same characteristics as those of two-dimensional ring polymers. Finally, we study the size-fluctuation dynamics of the pores, and compare the time evolution of their contour length to a random walk in a linear potential.     

Changes in density and surface tension of water in silica pores

 The density and surface tension of water in small pores of silicas have been investigated. These physical properties of water in the pores were calculated from a comparison of pore volumes and pore radii which were estimated from adsorption and desorption isotherms of nitrogen and water. Below a pore radius of about 5 nm both the density and the surface tension of water in the pores were smaller than those of the bulk liquid and decreased with a decrease in pore size. The density of water in the pores decreased with an increase in the concentration of surface hydroxyl groups. Similarly the surface tension of water in the pores is influenced by the surface hydroxyl groups. Anomalous changes in the density and surface tension of the water in the pores are attributed to the interaction of water molecules with surface hydroxyl groups and hydrogen-bond formation among water molecules. Received: 20 April 1999 Accepted in revised form: 17 November 1999     

Determination of the line tension of giant vesicles from pore-closing dynamics

Giant vesicles generated from synthetic and natural lipids such as phosphatidylcholines are useful models for understanding mechanical properties of cell membranes. Line tension is the one-dimensional force enabling the closing of transient pores on cell membranes. Transient pores were repeatedly and reproducibly formed on the membrane edge of giant vesicles generated from synthetic and natural phosphatidylcholines employing a nitrogen-pumped coumarin dye laser (440 nm). Line tension was determined at room temperature from closing of these pores that occurred over several seconds when the radius of the vesicle could be considered to be constant. The value of line tension depends on the nature of the lipid for single lipid systems, which, at room temperature, yielded a vesicle bilayer region in the gel, fluid, or mixed gel and fluid phases. The line tension for vesicles generated from phosphatidylcholines with saturated acyl chains of lengths of 12-18 carbon atoms ranges from 1 to 12 pN, exhibiting an increase with chain length. Vesicles generated from the natural Egg-PC, which is a mixture of lipids, are devoid of phase transition and exhibited the largest value of line tension (32 pN). This value is much larger than that estimated from the line tensions of vesicles obtained from lipids with homologous acyl chains. This study, to our knowledge, is the first to employ laser ablation to generate transient pores and determine line tension from the rate of pore closure and demonstrate a relationship between line tension and acyl chain length.     

Using wavelets to analyze AFM images of thin films: surface micelles and supported lipid bilayers

This paper presents micro- and nanoanalysis of thin films based on images obtained by atomic force microscopy (AFM). The analysis exploits the discrete wavelet transform and the resulting wavelet spectrum to study surface features. It is demonstrated that the wavelet technique can characterize micro- and nanosurface features and distinguish between similar surface structures. The use of a feature extraction method is shown. The method involves the separation of certain frequency content from the original AFM images and analyzing the data independently to gain quantitative information about the images. By using the feature extraction method, soft surfaces in water are analyzed and nanofeatures are measured. The packing of surface micelles of sodium dodecyl sulfate on a self-assembled monolayer is analyzed. The characteristics of pore formation, due to penetration of the antibacterial peptide protegrin, into a solid-supported lipid bilayer are quantified. The sizes of the pores are obtained, and it is observed that the line tension of the pores reduces the fluctuations of the lipid bilayer.     

SEM研究PET核孔膜的光接枝聚合

以PET核孔膜为基材 ,二苯甲酮为引发剂 ,采用光接枝方法实现了丙烯酸和丙烯酰胺在核孔膜上的接枝 ,用扫描电镜 (SEM)直接观察了接枝前后膜的表面形貌 ,考察了不同因素对于接枝位置和接枝效果的影响 .发现膜材料本身特性和接枝反应条件对接枝位置和接枝效果有较大影响 .通过光接枝能够实现膜孔的封盖、缩小、填堵等不同的效果 .采用正侧涂布法反应 ,标准直孔 ,特别是小孔径膜 (0 4 μm) ,不利于孔内的接枝 ,接枝主要在膜的表面 ,从而产生孔封盖效应 .双锥形的非标准直孔 ,由于孔壁的受光性好 ,容易发生孔壁上的接枝从而被填充 .大孔径的膜 (5 μm) ,需要加入交联剂才能在孔壁上形成厚的接枝层 .提出了一种新的反应方法 背侧吸附法 ,反应液依靠毛细作用由膜的底部吸入膜孔 ,膜的向光侧表面不存在反应液 ,接枝只发生在膜孔内 ,从而得到很好的填孔效果 .     

Molecular simulation of adsorption and intrusion in nanopores

This paper reports Monte Carlo simulations of the adsorption or intrusion in cylindrical silica nanopores. All the pores are opened at both ends towards an external bulk reservoir, so that they mimic real materials for which the confined fluid is always in contact with the external phase. This realistic model allows us to discuss the nature of the filling and emptying mechanisms. The adsorption corresponds to the metastable nucleation of the liquid phase, starting from a partially filled pore (a molecular thick film adsorbed at the pore surface). On the other hand, the desorption occurs through the displacement at equilibrium of a gas/liquid hemispherical interface (concave meniscus) along the pore axis. The intrusion of the non-wetting fluid proceeds through the invasion in the pore of the liquid/gas interface (convex meniscus), while the extrusion consists of the nucleation of the gas phase within the pore. In the case of adsorption, our simulation data are used to discuss the validity of the modified Kelvin equation (which is corrected for both the film adsorbed at the pore surface and the curvature effect on the gas/liquid surface tension).     

Modeling growth of organized nanoporous structures by anodic oxidation

Nanostructured porous oxides are produced by anodic dissolution of several metals. A scaling approach is introduced to explain pattern nucleation in an oxide layer, and a related microscopic model shows oxide growth with long nanopores. The scaling approach matches the time of ion transport across the thin oxide layer, which is related to metal corrosion, and the time of diffusion along the oxide/solution (OS) interface, which represents the extension of oxide dissolution. The selected pattern size is of order (dD(S)/v(O))(1/2), where d is the oxide thickness, v(O) is the migration velocity of oxygen ions across the oxide, and D(s) is the diffusion coefficient of H(+) ions along the oxide/solution interface. This result is consistent with available experimental data for those quantities, predicts the increase of pore size with the external voltage, and suggests the independence of pore size with the solution pH. Subsequently, we propose a microscopic model that expresses the main physicochemical processes as a set of characteristic lengths for diffusion and surface relaxation. It shows a randomly perturbed OS interface at short times, its evolution to pore nucleation and to stable growth of very long pores, in agreement with the mechanistic scenario suggested by two experimental groups. The decrease of the size of the walls between the pores with the interface tension is consistent with arguments for formation of titania nanotube arrays instead of nanopores. These models show that pattern nucleation and growth depend on matching a small number of physicochemical parameters, which is probably the reason for the production of nanostructured porous oxides from various materials under suitable electrochemical conditions.     

The effect of gas surface diffusion on the asymmetric permeability of two-layer porous membranes

The permeability is calculated for two-layer membranes composed of porous layers with nano- and microsized pores. It is found that, in nanosized pores, the gas transfer occurs in the free-molecular regime in the pore volume and via diffusion along the adsorption layer. The degree of adsorption layer filling is determined from the Langmuir isotherm. The dependence of the diffusion coefficient in the adsorption layer on the degree of surface coverage is taken into account. The transfer in the microsized pored is described in a hydrodynamic regime. The values of the membrane permeability are determined at different orientations with respect to the direction of a gas flow. It is shown that the difference between the permeability values may be as large as 60%.     

Nucleation in cylindrical capillaries

We use a local density functional theory in the square gradient approximation to explore the properties of critical nuclei for the liquid-vapor transition of van der Waals fluids in cylindrical capillaries. The proposed model allows us to investigate the effect of pore size, surface field, and supersaturation on the behavior of the system. Our calculations predict the existence of at least three different pathways for the nucleation of droplets and bubbles in these confined fluids: axisymmetric annular bumps and lenses, and asymmetric droplets. The morphological transition between these different structures is driven by the existence of states of zero compressibility in the capillary. We show that the classical capillarity theory provides surprisingly accurate predictions for the work of formation of critical nuclei in cylindrical pores when line tension contributions to the free energy are taken into account.     

Computational modelling of nematic phase ordering by film and droplet growth over heterogeneous substrates

This paper presents a computational study of defect nucleation associated with the kinetics of the isotropic-to-nematic phase ordering transition over heterogeneous substrates, as it occurs in new liquid crystal biosensor devices, based on the Landau-de Gennes model for rod-like thermotropic nematic liquid crystals. Two regimes are identified due to interfacial tension inequalities: (i) nematic surface film nucleation and growth normal to the heterogeneous substrate, and (ii) nematic surface droplet nucleation and growth. The former, known as wetting regime, leads to interfacial defect shedding at the moving nematic-isotropic interface. The latter droplet regime, involves a moving contact line, and exhibits two texturing mechanisms that also lead to interfacial defect shedding: (a) small and large contact angles of drops spreading over a heterogeneous substrate, and (b) small drops with large curvature growing over homogeneous patches of the substrate. The numerical results are consistent with qualitative defect nucleation models based on the kinematics of the isotropic-nematic interface and the substrate-nematic-isotropic contact line. The results extend current understanding of phase ordering over heterogeneous substrates by elucidating generic defect nucleation processes at moving interfaces and moving contact lines.     

Effects of radially dependent parameters on proton transport in polymer electrolyte membrane nanopores

A three-dimensional continuum model is explored to investigate the effects of radially dependent system parameters, such as relative permittivity and viscosity, on the transport of proton and water in nanoscale cylindrical pores of a fully hydrated polymer electrolyte membrane (PEM). The model employs Poisson, Nernst-Planck, and Stokes equations. Based on evidence from the literature for the presence of a stagnant water layer near the pore surface, we assume that a no-slip surface is located inside the pore, a few Angstroms from the pore wall. To solve the system numerically, the steady-state solution for the transport of protons and water is considered to be a perturbation around the equilibrium solution. Our results indicate that a radial variation of relative permittivity has the greatest influence on pore conductivity, reducing it by about 50% when compared to that of constant permittivity. On the other hand, viscosity plays the dominant role when the effective water drag within such pores is considered. We conclude that a continuum approach, including constant viscosity, is applicable in nanoscale models provided that the location of the no-slip surface is properly specified and the radial variation of the relative permittivity is taken into consideration.     

Modeling of pore formation in phase inversion processes: analysis of pore formation mechanism

The phase inversion process is the most important preparation process of porous polymer membranes. Recently, a numerical model based on first principles to predict pore structures has been proposed (Hopp-Hirschler and Nieken in J Membr Sci 564:820–831, 2018). This model enables a detailed investigation of the mechanism of pore formation in porous polymer membranes. This follow-up presents investigations of the mechanism of nucleation of pores during the phase inversion process in 1D. Pores originate due to accumulation of over-saturated mixtures inside a diffuse interface between homogeneous and demixed polymer solutions behind the precipitation front. This is caused by an expansion of the width of the diffuse interface and time-dependent concentration profiles which finally lead to a change of sign of total diffusive mass flux inside of the diffuse interface. As a result, oscillating compositions behind the precipitation front lead to formation of pores. It is concluded that large surface tension leads to initially small pore sizes. In the second part, a detailed discussion of directional diffusion behind the precipitation front is presented in 2D, which is responsible for different pore structures, e.g., finger or sponge pores. Depending on the dominant direction of diffusion finger pores, lamella structures or sponge pores are formed. This picture can straightforwardly be extended to 3D structures.

     

Nucleation in hydrophobic cylindrical pores: a lattice model

We consider the nucleation process associated with capillary condensation of a vapor in a hydrophobic cylindrical pore (capillary evaporation). The liquid-vapor transition is described within the framework of a simple lattice model. The phase properties are characterized both at the mean-field level and with Monte Carlo simulations. The nucleation process for the liquid to vapor transition is then specifically considered. Using umbrella sampling techniques, we show that nucleation occurs through the condensation of an asymmetric vapor bubble at the pore surface. Even for highly confined systems, good agreement is found with macroscopic considerations based on classical nucleation theory. The results are discussed in the context of recent experimental work on the extrusion of water in hydrophobic pores.     

Molecular mechanism for lipid flip-flops

Transmembrane lipid translocation (flip-flop) processes are involved in a variety of properties and functions of cell membranes, such as membrane asymmetry and programmed cell death. Yet, flip-flops are one of the least understood dynamical processes in membranes. In this work, we elucidate the molecular mechanism of pore-mediated transmembrane lipid translocation (flip-flop) acquired from extensive atomistic molecular dynamics simulations. On the basis of 50 successful flip-flop events resolved in atomic detail, we demonstrate that lipid flip-flops may spontaneously occur in protein-free phospholipid membranes under physiological conditions through transient water pores on a time scale of tens of nanoseconds. While the formation of a water pore is induced here by a transmembrane ion density gradient, the particular way by which the pore is formed is irrelevant for the reported flip-flop mechanism: the appearance of a transient pore (defect) in the membrane inevitably leads to diffusive translocation of lipids through the pore, which is driven by thermal fluctuations. Our findings strongly support the idea that the formation of membrane defects in terms of water pores is the rate-limiting step in the process of transmembrane lipid flip-flop, which, on average, requires several hours. The findings are consistent with available experimental and computational data and provide a view to interpret experimental observations. For example, the simulation results provide a molecular-level explanation in terms of pores for the experimentally observed fact that the exposure of lipid membranes to electric field pulses considerably reduces the time required for lipid flip-flops.     

Computational modelling of nematic phase ordering by film and droplet growth over heterogeneous substrates

This paper presents a computational study of defect nucleation associated with the kinetics of the isotropic‐to‐nematic phase ordering transition over heterogeneous substrates, as it occurs in new liquid crystal biosensor devices, based on the Landau–de Gennes model for rod‐like thermotropic nematic liquid crystals. Two regimes are identified due to interfacial tension inequalities: (i) nematic surface film nucleation and growth normal to the heterogeneous substrate, and (ii) nematic surface droplet nucleation and growth. The former, known as wetting regime, leads to interfacial defect shedding at the moving nematic‐isotropic interface. The latter droplet regime, involves a moving contact line, and exhibits two texturing mechanisms that also lead to interfacial defect shedding: (a) small and large contact angles of drops spreading over a heterogeneous substrate, and (b) small drops with large curvature growing over homogeneous patches of the substrate. The numerical results are consistent with qualitative defect nucleation models based on the kinematics of the isotropic–nematic interface and the substrate–nematic–isotropic contact line. The results extend current understanding of phase ordering over heterogeneous substrates by elucidating generic defect nucleation processes at moving interfaces and moving contact lines.     

耗散粒子动力学研究双层膜线张力与抗弯刚度的关系

利用耗散力粒子动力学模拟方法研究了拉伸状态下双层膜的抗弯刚度与膜孔线张力之间的关系.通过对双层膜在不同投影面积约束下的系统模拟,观察到3个区域:自由振动膜、伸展膜和穿孔膜.由前两个区域的应力张量计算得到的膜的面张力(σΣ)与拟合膜波动性质得到的面张力(σfluc)吻合的很好,除去在两区域的转变点附近σfluc略大于σΣ.当考虑在经典Helfrich弹性膜模型中被忽略的膜厚度时,线张力可以和抗弯刚度用一个简单的模型联系起来.通过对穿孔膜区域的数据进行分析,证明由本模型得到的抗弯刚度与拟合膜波动性质得到的抗弯刚度符合的很好.由此提出一种简便的测量方法,通过计算拉伸膜孔的面张力和统计膜厚度,拟合这个简单模型来测量膜的抗弯刚度.     

Thermodynamics of heterogeneous crystal nucleation in contact and immersion modes

One of the most intriguing problems of heterogeneous crystal nucleation in droplets is its strong enhancement in the contact mode (when the foreign particle is presumably in some kind of contact with the droplet surface) compared to the immersion mode (particle immersed in the droplet). Heterogeneous centers can have different nucleation thresholds when they act in contact or immersion modes. The underlying physical reasons for this enhancement have remained largely unclear. In this paper we present a model for the thermodynamic enhancement of heterogeneous crystal nucleation in the contact mode compared to the immersion one. To determine if and how the surface of a liquid droplet can thermodynamically stimulate its heterogeneous crystallization, we examine crystal nucleation in the immersion and contact modes by deriving and comparing with each other the reversible works of formation of crystal nuclei in these cases. The line tension of a three-phase contact gives rise to additional terms in the formation free energy of a crystal cluster and affects its Wulff (equilibrium) shape. As an illustration, the proposed model is applied to the heterogeneous nucleation of hexagonal ice crystals on generic macroscopic foreign particles in water droplets at T = 253 K. Our results show that the droplet surface does thermodynamically favor the contact mode over the immersion one. Surprisingly, the numerical evaluations suggest that the line tension contribution (from the contact of three water phases (vapor-liquid-crystal)) to this enhancement may be of the same order of magnitude as or even larger than the surface tension contribution.