A generalized scale of free energy of excess adsorption of solute and absolute composition of the interfacial phase

Moles of a surfactant (gamma2(1)) absorbed per unit area of the solid-liquid interface estimated analytically from the difference of the solute molality in the bulk phase before and after adsorption have been quantitatively related to the absolute compositions deltan1 and deltan2 of the solvent and solute forming the inhomogeneous surface phase in contact with the bulk phase of homogeneous composition. By use of isopiestic experiments, negative values of gamma2(1) for the adsorption of inorganic salts onto a solid-liquid interface have been calculated in the same manner. From the linear plot of gamma2(1) versus the ratio of the bulk mole fractions of the solute and solvent, values of deltan1 and deltan2 have been evaluated under a limited range of concentrations. For the adsorption of the surfactant and the inorganic salt respectively onto the fluid interface, gamma2(1) values have been evaluated from the surface tension concentration data using the Gibbs adsorption equation. Gamma2(1) based on the arbitrary placement of the Gibbs dividing plane near the fluid interface is quantitatively related to the composition of the inhomogeneous surface phase. Also, the Gibbs equation for multicomponent solutions has been appropriately expressed in terms of a suitably derived coefficient m. Integrating the Gibbs adsorption equation for a multicomponent system, the standard free energy change, deltaG degrees, per unit of surface area as a result of the maximum adsorption gamma2(m) of the surfactant at fluid interfaces due to the change of the activity alpha2 of the surfactant in the bulk from zero to unity have been calculated. A similar procedure has been followed for the calculation of deltaG degrees for the surfactant adsorption at solid-liquid interfaces using thermodynamically derived equations. deltaG degrees values for surfactant adsorption for all such systems are found to be negative. General expressions of deltaG degrees for negative adsorption of the salt on fluid and solid-liquid interfaces respectively have also been derived on thermodynamic grounds. deltaG degrees for all such systems are positive due to the excess spontaneous hydration of the interfacial phase in the presence of inorganic salt. Negative and positive values of deltaG degree for excess surfactant and salt adsorption respectively have been discussed in light of a generalized scale of free energy of adsorption.     

Interactions between hen egg-white lysozyme, PEG2,000, and PLA50 at the air-water interface

In this paper, we compared the efficiency of polymer films, made of a poly(ethylene glycol) (PEG2,000)/poly(d,l-lactide) (PLA50) mixture, or a PEG2,000-PLA50 copolymer, to prevent adsorption of a model protein, the hen egg-white lysozyme (HEWL), at the air-water interface. This was achieved by analyzing the surface pressure/surface area curves, and the X-ray reflectivity data of the polymer films spread on a Langmuir trough, obtained in absence or in presence of the protein. For both the mixture and the copolymer, the amount of protein adsorbed at the air-water interface decreases when the density of the polymer surface coverage increases. It was shown that even in a condensed state, the polymer film made by the mixture can not totally prevent HEWL molecules to adsorb and penetrate the polymer mixed film, but however, protein molecules would not be directly exposed to the more hydrophobic phase, i.e. the air phase. It was also shown that the configuration adopted by the copolymer at the interface in its condensed state would prevent adsorption of HEWL molecules for several hours; this would be due in particular to the presence of PEG segments in the interfacial film.     

The adsorption and unfolding kinetics determines the folding state of proteins at the air-water interface and thereby the equation of state

Unfolding of proteins has often been mentioned as an important factor during the adsorption process at air-water interfaces and in the increase of surface pressure at later stages of the adsorption process. This work focuses on the question whether the folding state of the adsorbed protein depends on the rate of adsorption to the interface, which can be controlled by bulk concentration. Therefore, the adsorption of proteins with varying structural stabilities at several protein concentrations was studied using ellipsometry and surface tensiometry. For beta-lactoglobulin the adsorbed amount (Gamma) needed to reach a certain surface pressure (Pi) decreased with decreasing bulk concentration. Ovalbumin showed no such dependence. To verify whether this difference in behavior is caused by the difference in structural stability, similar experiments were performed with cytochrome c and a destabilized variant of this protein. Both proteins showed identical Pi-Gamma, and no dependence on bulk concentration. From this work it was concluded that unfolding will only take place if the kinetics of adsorption is similar or slower than the kinetics of unfolding. The latter depends on the activation energy of unfolding (which is in the order of 100-300 kJ/mol), rather than the free energy of unfolding (typically 10-50 kJ/mol).     

Monte Carlo simulation study of the kinetics of brush formation by irreversible adsorption of telechelic polymers onto a solid substrate

The irreversible adsorption of telechelic polymer chains from solution and melts onto solid substrates has been studied using the bond fluctuation Monte Carlo model. Complex brush formation kinetics dominated by diffusion of chains to the substrate at short times (diffusion-limited regime or DLR) and by penetration of chains through the maturing brush at longer times (penetration-limited regime or PLR) were observed. During the entire adsorption process, the rate of chain adsorption decreases monotonically with time. In the DLR, characterized by a maximum in the concentration of singly bound chains and a rapidly increasing fraction of doubly bound chains (loops), this decrease is due primarily to the depletion of free chains near the substrate and the formation of concentration gradients of free (nonadsorbed) chains in the bulk solution. The DLR and PLR are separated by an intermediate regime during which the brush becomes dominated by doubly bound chains and both penetration of the maturing brush and diffusion of chains to the brush surface play a role in determining the kinetics of brush growth. The PLR is characterized by steep gradients of free chains within the growing brush and the disappearance of concentration gradients for free chains in the bulk solution. In the PLR, the concentration of singly bound chains is low and decreases slowly while surface coverage and the fraction of doubly bound chains increase slowly. The rates of adsorption of new chains and the formation of loops in the PLR slow dramatically with increasing surface coverage and increasing chain length and less dramatically with decreasing bulk concentration.     

Mechanisms of fibrinogen adsorption at solid substrates

Adsorption of fibrinogen, modeled as a linear chain of touching beads of various sizes, was theoretically studied using the random sequential adsorption (RSA) model. The adsorption process was assumed to consist of two steps: (i) formation of an irreversibly bound fibrinogen monolayer under the side-on orientation, which is independent of the bulk protein concentration and (ii) formation of the reversibly bound, end-on monolayer, whose coverage was dependent on the bulk concentration. Calculation based on the RSA model showed that the maximum surface concentration of the end-on (reversible) monolayer equals N(⊥∞) = 6.13 × 10(3) μm(-2) which is much larger than the previously found value for the side-on (irreversible) monolayer, equal to N(∞) = 2.27 × 10(3) μm(-2). Hence, the maximum surface concentration of fibrinogen in both orientations is determined to be 8.40 × 10(3) μm(-2) corresponding to the protein coverage of 5.70 mg m(-2) assuming 20% hydration. Additionally, the surface blocking function (ASF) was determined for the end-on fibrinogen adsorption, approximated for the entire range of coverage by the interpolating polynomial. For the coverage approaching the jamming limit, the surface blocking function (ASF) was shown to vanish proportionally to (θ(⊥∞) - θ(⊥))(2). These calculation allowed one to theoretically predict adsorption isotherms for the end-on regime of fibrinogen and adsorption kinetics under various transport conditions (diffusion and convection). Using these theoretical results, a quantitative interpretation of experimental data obtained by TIRF and ellipsometry was successfully performed. The equilibrium adsorption constant for the end-on adsorption regime was found to be 8.04 × 10(-3) m. On the basis of this value, the depth of the adsorption energy minimum, equal to -17.4 kT, was predicted, which corresponds to ΔG = -41.8 kJ mol(-1). This is in accordance with adsorption energy derived as the sum of the van der Waals and electrostatic interactions. Besides having significance for predicting fibrinogen adsorption, theoretical results derived in this work also have implications for basic science providing information on mechanisms of anisotropic protein molecule adsorption on heterogeneous surfaces.     

Lattice Monte Carlo simulation of polymer adsorption at an interface, 1. Monodisperse polymer

Adsorption of a monodisperse polymer at a solid-liquid interface is comprehensively studied by Monte Carlo simulation. The distributions of total segment density and different adsorption configurations including trains, loops and tails are obtained. Effects of reduced exchange interaction energies $ \tilde \varepsilon $, bulk concentrations ϕ^*, reduced adsorption energies $ \tilde \varepsilon_a $ and chain lengths r on those distributions are studied. Comparisons with predictions of the Scheutjens-Fleer (SF) theory are also provided. Generally, the chain molecules are more easily adsorbed at an interface in non-solvents than in good solvents. Longer chains are more likely to be adsorbed than shorter ones. The reduced adsorption energy and the bulk concentration have shown strong effects on the segment-density distributions. In addition, the thickness of the adsorption layer is mainly determined by the extension of tails into the bulk solution, which are in turn determined by the chain length. The trains, loops and tails are overwhelmingly short. On the other hand, the amounts of trains and loops are usually much greater than that of tails. Though not perfect, satisfactory agreement is found in comparison with the theoretical predictions of the SF theory.     

A comprehensive study of concepts and phenomena of the nonspecific adsorption of beta-lactoglobulin.

We investigate nonspecific protein adsorption processes by comparing experimentally measured adsorption kinetics of beta-lactoglobulin with mathematical models. The adsorption and desorption behavior of this protein on a hydrophilic glass surface in citrate buffer (pH 3.0), monitored for a large set of different bulk concentrations (0.5x10(-8) M-1.5x10(-6) M) using a supercritical angle fluorescence (SAF) biosensor, is reported. Increasing adsorption rates and overshootings in the beginning of the adsorption are observed as well as a transition to an almost irreversibly bound state of the protein in the long term. Furthermore, rinsing experiments prove that adsorbed proteins abruptly change their desorption behavior from irreversible to reversible when a critical surface coverage theta(crit) is reached. Based on all experimental observations, a mathematical model composed of three adsorbed states differing in their surface affinity is proposed. Terms to account for lateral interactions between surface-bound proteins are included, which yield an excellent fit of the measured kinetics. For the first time, several phenomena that have been discussed in theoretical studies are confirmed by comparing experimental data with a single model.     

Interfacial energetics of protein adsorption from aqueous buffer to surfaces with varying hydrophilicity

Adsorption isotherms constructed from time-and-concentration-dependent advancing contact angles thetaa show that the profound biochemical diversity among ten different blood proteins with molecular weight spanning 10-1000 kDa has little discernible effect on the amount adsorbed from aqueous phosphate-buffered saline (PBS) solution after 1 h contact with a particular test surface selected from the full range of observable water wettability (as quantified by PBS adhesion tension tauoa=gammaolv cos thetaoa; where gammaolv is the liquid-vapor interfacial tension and thetaoa is the advancing PBS contact angle). The maximum advancing spreading pressure, Pimaxa, determined from adsorption isotherms decreases systematically with tauoa for methyl-terminated self-assembled monolayers (CH3 SAM, tauo=-15 mN/m), polystyrene spun-coated onto electronic-grade SiOx wafers (PS, tauo=7.2 mN/m), aminopropyltriethoxysilane-treated SiOx surfaces (APTES, tauo = 42 mN/m), and fully water wettable SiOx (tauo=72 mN/m). Likewise, the apparent Gibbs' surface excess [Gammasl-Gammasv], which measures the difference in the amount of protein adsorbed Gamma (mol/cm2) at solid-vapor (SV) and solid-liquid (SL) interfaces, decreases with tauo from maximal values measured on the CH3 SAM surface through zero (no protein adsorption in excess of bulk solution concentration) near tauo=30 mN/m (thetaa=65 degrees). These latter results corroborate the conclusion drawn from independent studies that water is too strongly bound to surfaces with tauo>or=30 mN/m to be displaced by adsorbing protein and that, as a consequence, protein does not accumulate within the interfacial region of such surfaces at concentrations exceeding that of bulk solution ([Gammasl-Gammasv]=0 at tauo=30 mN/m). Results are collectively interpreted to mean that water controls protein adsorption to surfaces and that the mechanism of protein adsorption can be understood from this perspective for a diverse set of proteins with very different amino acid compositions.     

Amphiphilic derivatives of dextran: adsorption at air/water and oil/water interfaces

Ionic amphiphilic dextran derivatives were synthesized by the attachment of sodium sulfopropyl and phenoxy groups on the native polysaccharide. A family of dextran derivatives was thus obtained with varying hydrophobic content and charge density in the polymer chains. The surface-active properties of polymers were studied at the air-water and dodecane-water interfaces using dynamic surface/interfacial tension measurements. The adsorption was shown to begin in a diffusion-limited regime at low polymer concentrations, that is to say, with the diffusion of macromolecules in the bulk solution. In contrast, at long times the interfacial adsorption is limited by interfacial phenomena: adsorption kinetics or transfer into the adsorbed layer. A semiempirical equation developed by Filippov was shown to correctly fit the experimental curves over the whole time range. The presence of ionic groups in the chains strongly lowers the adsorption kinetics. This effect can be interpreted by electrostatic interactions between the free molecules and the already adsorbed ones. The adsorption kinetics at air-water and oil-water interfaces are compared.     

Nanoscale colloids in a freely adsorbing polymer solution: a Monte Carlo simulation study

A key issue in nanoscale materials and chemical processing is the need for thermodynamic and kinetic models covering colloid-polymer systems over the mesoscopic length scale (approximately 1-100 nm). We have applied Monte Carlo simulations to attractive nanoscale colloid-polymer mixtures toward developing a molecular basis for models of these complex systems. The expanded ensemble Monte Carlo simulation method is applied to calculate colloid chemical potentials (micro(c)) and polymer adsorption (gamma) in the presence of freely adsorbing Lennard-Jones (LJ) homopolymers (surface modifiers). gamma and micro(c) are studied as a function of nanoparticle diameter (sigma(c)), modifier chain length (n) and concentration, and colloid-polymer attractive strength over 0.3 < Rg/sigma(c) < 6 (Rg is the polymer radius of gyration). In the attractive regime, nanocolloid chemical potential decreases and adsorbed amount increases as sigma(c), or n is increased. The scaling of gamma with n from the simulations agrees with the theory of Aubouy and Raphael (Macromolecules 1998, 31, 4357) in the extreme limits of Rg/sigma(c). When Rg/sigma(c) is large, the "colloid" approaches a molecular size and interacts only locally with a few polymer segments and gamma approximately n. When Rg/sigma(c) is small, the system approaches the conventional colloid-polymer size regime where multiple chains interact with a single particle, and gamma approximately sigma(c)2, independent of n. In contrast, adsorption in the mesoscopic range of Rg/sigma(c) investigated here is represented well by a power law gamma approximately n(p), with 0 < p < 1 depending on concentration and LJ attractive strength. Likewise, the chemical potential from our results is fitted well with micro(c) approximately n(q)sigma(c)3, where the cubic term results from the sigma(c) dependence of particle surface area (approximately sigma(c)2) and LJ attractive magnitude (approximately sigma(c)). The q-exponent for micro(c) (micro(c) approximately n(q)) varies with composition and LJ attractive strength but is always very close to the power exponent for gamma (gamma approximately n(p)). This result leads to the conclusion that in attractive systems, polymer adsorption (and thus polymer-colloid attraction) dominates the micro(c) dependence on n, providing a molecular interpretation of the effect of adsorbed organic layers on nanoparticle stability and self-assembly.     

Recent progress in the determination of solid surface tensions from contact angles

Advancing contact angles of different liquids measured on the same solid surface fall very close to a smooth curve when plotted as a function of liquid surface tension, i.e., gamma(lv)costheta versus gamma(lv). Changing the solid surface, and hence gamma(sv), shifts the curve in a regular manner. These patterns suggest that gamma(lv)costheta depends only on gamma(lv) and gamma(sv). Thus, an "equation of state for the interfacial tensions" was developed to facilitate the determination of solid surface tensions from contact angles in conjunction with Young's equation. However, a close examination of the smooth curves showed that contact angles typically show a scatter of 1-3 degrees around the curves. The existence of the deviations introduces an element of uncertainty in the determination of solid surface tensions. Establishing that (i) contact angles are exclusively a material property of the coating polymer and do not depend on experimental procedures and that (ii) contact angle measurements with a sophisticated methodology, axisymmetric drop shape analysis (ADSA), are highly reproducible guarantees that the deviations are not experimental errors and must have physical causes. The contact angles of a large number of liquids on the films of four different fluoropolymers were studied to identify the causes of the deviations. Specific molecular interactions at solid-vapor and/or solid-liquid interfaces account for the minor contact angle deviations. Such interactions take place in different ways. Adsorption of vapor of the test liquid onto the solid surface is apparently the only process that influences the solid-vapor interfacial tension (gamma(sv)). The molecular interactions taking place at the solid-liquid interface are more diverse and complicated. Parallel alignment of liquid molecules at the solid surface, reorganization of liquid molecules at the solid-liquid interface, change in the configuration of polymer chains due to contact with certain probe liquids, and intermolecular interactions between solid and liquid molecules cause the solid-liquid interfacial (gamma(sl)) tension to be different from that predicted by the equation of state, i.e., gamma(sl) is not a precise function of gamma(lv) and gamma(sv). In other words, the experimental contact angles deviate from the "ideal" contact angle pattern. Specific criteria are proposed to identify probe liquids which eliminate specific molecular interactions. Octamethylcyclotetrasiloxane (OMCTS) and decamethylcyclopentasiloxane (DMCPS) are shown to meet those criteria, and therefore are the most suitable liquids to characterize surface tensions of low energy fluoropolymer films with an accuracy of +/-0.2 mJ/m2.     

反气相色谱法测定苯乙烯-氧乙烯-苯乙烯三嵌段聚合物的表面能

 采用反气相色谱法测定了苯乙烯 氧乙烯 苯乙烯三嵌段聚合物 (PS PEO PS)的色散成分的表面能 (γsd) ,研究探讨了温度及嵌段聚合物链段结构组成对γsd 的影响 ,并确定了γsd 与温度的数学关系式。研究结果表明 :在 70℃～ 12 0℃范围内 ,PS PEO PS的表面能较低 ;随着PS PEO PS表面组成中氧乙烯 (EO)成分的增加 ,色散成分的γsd 增大 ,且对温度的变化极其敏感 :随温度的升高 ,γsd 急剧地呈线性下降。     

Synthesis, characterization and protein adsorption behaviors of PLGA/PEG di-block co-polymer blend films

A series of poly(D,L-lactic-co-glycolic acid) (PLGA)/poly(ethyleneglycol) (PEG) di-block copolymers were synthesized by ring-opening polymerization of D,L-lactide and glycolide with different molecular weights of monomethoxy polyethyleneglycol (mPEG) 750, 2000 and 5000 as an initiator. The bulk properties of these co-polymers were characterized by using 1H NMR spectroscopy, gel permeation chromatography, differential scanning calorimetry (DSC). Electron spectroscopy for chemical analysis (ESCA) results, in which the blend films with the di-block copolymers showed increasing surface oxygen atomic percentage with increasing PEG chain length, indicate that PEG chain segment in the di-block copolymers is surface oriented and enriched onto the surface of the blend films. The extent of protein adsorption onto the surface of these blend films was studied, using iodine radio-labeled human serum albumin, gamma globulin and human growth hormone. The protein adsorption amount was reduced for the blend films prepared with PLGA/PEG 750 and 2000 di-block copolymers, but increased to a great extent for PLGA/PEG 5000 di-block copolymer. This is due to the increased water uptake capacity of the blend film, which absorbed more protein molecules into a swollen polymer matrix in addition to surface adsorption.     

Spreading of polymer molecules on solid surfaces

When a homopolymer adsorbs from dilute solution onto a solid surface it first attaches and then maximises its number of contacts with the substrate by means of a spreading process. Evidence for this spreading process can be obtained from experiments on the adsorption kinetics. We report on a case where the adsorption kinetics depend on the rate at which the polymer was supplied to the surface (a protein adsorbing onto silica). Also, we discuss competitive adsorption experiments in which one kind of chain molecule attempts to displace another one from the surface. In these experiments, the desorption rate of the displaced species reflects the spreading rate of the displacer. When this rate is slower than the supply of displacer molecules, oversaturated layers result that spontaneously eject polymer. We have measured the rate of these displacement-driven desorption processes in various cases and conclude that it depends strongly on the energy of the segment-surface bond. A model involving the diffusion of defects over the surface may account for this finding.     

Competitive adsorption of nonionic surfactant and nonionic polymer on silica

Competitive adsorption of the nonionic polymer poly(ethylene oxide) (PEO) and the nonionic surfactant of the type poly(ethylene oxide) alkyl ether from aqueous solutions on a silica surface is examined. From one-component solutions, both species readily adsorb onto silica and, in the bulk of mixed (two-component) solutions, polymer-surfactant complexes are not observed. Because both species bind by the same mechanism to silica, subtle differences in layer structure, or other species-specific parameters, determine whether one or both of the species will adsorb. It was found that various surfactants can displace PEO up to a certain critical molecular weight. Surfactants with a high aggregation number, in bulk and on the surface, can displace PEO with a higher molar mass than surfactants with a low aggregation number. As the molar mass of the polymer increases, the time a surfactant needs to completely displace the polymer increases. We can explain both the existence of the critical molar mass and the decrease in adsorption kinetics with a shift in the critical surface association concentration (CSAC).     

Dynamic adsorption of weakly interacting polymer/surfactant mixtures at the air/water interface

The dynamic adsorption of polymer/surfactant mixtures containing poly(ethylene oxide) (PEO) with either tetradecyltrimethylammonium bromide (C(14)TAB) or sodium dodecyl sulfate (SDS) has been studied at the expanding air/water interface created by an overflowing cylinder, which has a surface age of 0.1-1 s. The composition of the adsorption layer is obtained by a new approach that co-models data obtained from ellipsometry and only one isotopic contrast from neutron reflectometry (NR) without the need for any deuterated polymer. The precision and accuracy of the polymer surface excess obtained matches the levels achieved from NR measurements of different isotopic contrasts involving deuterated polymer, and requires much less neutron beamtime. The PEO concentration was fixed at 100 ppm and the electrolyte concentration at 0.1 M while the surfactant concentration was varied over three orders of magnitude. For both systems, at low bulk surfactant concentrations, adsorption of the polymer is diffusion-controlled while surfactant adsorption is under mixed kinetic/diffusion control. Adsorption of PEO is inhibited once the surfactant coverage exceeds 2 μmol m(-2). For PEO/C(14)TAB, polymer adsorption drops abruptly to zero over a narrow range of surfactant concentration. For PEO/SDS, inhibition of polymer adsorption is much more gradual, and a small amount remains adsorbed even at bulk surfactant concentrations above the cmc. The difference in behavior of the two mixtures is ascribed to favorable interactions between the PEO and SDS in the bulk solution and at the surface.     

Protein adsorption: A quest for a universal mechanism

Recent theoretical and experimental results pertinent to protein adsorption kinetics obtained for well-defined systems using direct experimental techniques are discussed. Attention is focused on albumins and fibrinogen, whose structure and physicochemical characteristic are well-known. It is confirmed that the experimental data obtained by AFM imaging, QCM, OWLS, XPS and electrokinetic techniques (streaming potential) are prone to a quantitative interpretation in terms of the coarse-grained and molecular dynamics modeling. This allows to derive reliable data concerning the mass transfer rates, hydration functions, maximum coverages and adsorption/desorption kinetic constants. These results confirm that the protein adsorption mechanism is governed by electrostatic interactions among heterogeneously distributed charges. The protein substrate interactions promote the molecule transfer through the surface layer, control the free energy and in consequence the residence time of the molecule on substrate surfaces. On the other hand, the interactions among adsorbed molecules control the maximum coverage and the formation of bilayer structures. As a result of this complex electrostatics, one often observes in protein adsorption studies the formation of irreversibly bound fraction of molecules that contact the substrate and a reversibly adsorbed fraction otherwise. This leads to the appearance of anomalous isotherms, characterized by considerable adsorption for negligible bulk protein concentration, which deviate from the Langmuir model.     

Multistage adsorption of diffusing macromolecules and viruses

We derive the equations that describe adsorption of diffusing particles onto a surface followed by additional surface kinetic steps before being transported across the interface. Multistage surface kinetics occurs during membrane protein insertion, cell signaling, and the infection of cells by virus particles. For example, viral entry into healthy cells is possible only after a series of receptor and coreceptor binding events occurs at the cellular surface. We couple the diffusion of particles in the bulk phase with the multistage surface kinetics and derive an effective, integrodifferential boundary condition that contains a memory kernel embodying the delay induced by the surface reactions. This boundary condition takes the form of a singular perturbation problem in the limit where particle-surface interactions are short ranged. Moreover, depending on the surface kinetics, the delay kernel induces a nonmonotonic, transient replenishment of the bulk particle concentration near the interface. The approach generalizes that of Ward and Tordai [J. Chem. Phys. 14, 453 (1946)] and Diamant and Andelman [Colloids Surf. A 183-185, 259 (2001)] to include surface kinetics, giving rise to qualitatively new behaviors. Our analysis also suggests a simple scheme by which stochastic surface reactions may be coupled to deterministic bulk diffusion.     

Free Energy of Bare and Capped Gold Nanoparticles Permeating through a Lipid Bilayer

Herein, we study the permeation free energy of bare and octane‐thiol‐capped gold nanoparticles (AuNPs) translocating through a lipid membrane. To investigate this, we have pulled the bare and capped AuNPs from bulk water to the membrane interior and estimated the free energy cost. The adsorption of the bare AuNP on the bilayer surface is energetically favorable but further loading inside it requires energy. However, the estimated free‐energy barrier for loading the capped AuNP into the lipid membrane is much higher compared to bare AuNP. We also demonstrate the details of the permeation process of bare and capped AuNPs. Bare AuNP induces the curvature in the lipid membrane whereas capped AuNP creates an opening in the interacting monolayer and get inserted into the membrane. The insertion of capped AuNP induces a partial unzipping of the lipid bilayer, which results in the ordering of the local lipids interacting with the nanoparticle. However, bare AuNP disrupts the lipid membrane by pushing the lipid molecules inside the membrane. We also analyze pore formation due to the insertion of capped AuNP into the membrane, which results in water molecules penetrating the hydrophobic region.