Classical density functional theory of orientational order at interfaces: application to water

A classical density functional formalism has been developed to predict the position-orientation number density of structured fluids. It is applied to the liquid-vapor interface of pure water, where it consists of a classical term, a gradient correction, and an anisotropic term that yields order through density gradients. The model is calibrated to predict that water molecules have their dipole moments almost parallel to a planar interface, while the molecular plane is parallel to it on the liquid side and perpendicular to it on the vapor side. For a planar interface, the surface tension obtained is twice its experimental value, while the surface potential is in qualitative agreement with that calculated by others. The model is also used to predict the orientation of water molecules near the surface of droplets, as well as the dependence of equilibrium vapor pressure around them on their size.     

Aqueous Heterogeneity at the Air/Water Interface Revealed by 2D‐HD‐SFG Spectroscopy

Water molecules interact strongly with each other through hydrogen bonds. This efficient intermolecular coupling causes strong delocalization of molecular vibrations in bulk water. We study intermolecular coupling at the air/water interface and find intermolecular coupling 1) to be significantly reduced and 2) to vary strongly for different water molecules at the interface—whereas in bulk water the coupling is homogeneous. For strongly hydrogen‐bonded OH groups, coupling is roughly half of that of bulk water, due to the lower density in the near‐surface region. For weakly hydrogen‐bonded OH groups that absorb around 3500 cm^?1, which are assigned to the outermost, yet hydrogen‐bonded OH groups pointing towards the liquid, coupling is further reduced by an additional factor of 2. Remarkably, despite the reduced structural constraints imposed by the interfacial hydrogen‐bond environment, the structural relaxation is slow and the intermolecular coupling of these water molecules is weak.     

A nonadditive methanol force field: bulk liquid and liquid-vapor interfacial properties via molecular dynamics simulations using a fluctuating charge model

We study the bulk and interfacial properties of methanol via molecular dynamics simulations using a CHARMM (Chemistry at HARvard Molecular Mechanics) fluctuating charge force field. We discuss the parametrization of the electrostatic model as part of the ongoing CHARMM development for polarizable protein force fields. The bulk liquid properties are in agreement with available experimental data and competitive with existing fixed-charge and polarizable force fields. The liquid density and vaporization enthalpy are determined to be 0.809 g/cm3 and 8.9 kcal/mol compared to the experimental values of 0.787 g/cm3 and 8.94 kcal/mol, respectively. The liquid structure as indicated by radial distribution functions is in keeping with the most recent neutron diffraction results; the force field shows a slightly more ordered liquid, necessarily arising from the enhanced condensed phase electrostatics (as evidenced by an induced liquid phase dipole moment of 0.7 D), although the average coordination with two neighboring molecules is consistent with the experimental diffraction study as well as with recent density functional molecular dynamics calculations. The predicted surface tension of 19.66+/-1.03 dyn/cm is slightly lower than the experimental value of 22.6 dyn/cm, but still competitive with classical force fields. The interface demonstrates the preferential molecular orientation of molecules as observed via nonlinear optical spectroscopic methods. Finally, via canonical molecular dynamics simulations, we assess the model's ability to reproduce the vapor-liquid equilibrium from 298 to 423 K, the simulation data then used to obtain estimates of the model's critical temperature and density. The model predicts a critical temperature of 470.1 K and critical density of 0.312 g/cm3 compared to the experimental values of 512.65 K and 0.279 g/cm3, respectively. The model underestimates the critical temperature by 8% and overestimates the critical density by 10%, and in this sense is roughly equivalent to the underlying fixed-charge CHARMM22 force field.     

Long-range Lennard-Jones and electrostatic interactions in interfaces: application of the isotropic periodic sum method

Molecular dynamics (MD) simulations of heptane/vapor, hexadecane/vapor, water/vapor, hexadecane/water, and dipalmitoylphosphatidylcholine (DPPC) bilayers and monolayers are analyzed to determine the accuracy of treating long-range interactions in interfaces with the isotropic periodic sum (IPS) method. The method and cutoff (rc) dependences of surface tensions, density profiles, water dipole orientation, and electrostatic potential profiles are used as metrics. The water/vapor, heptane/vapor, and hexadecane/vapor interfaces are accurately and efficiently calculated with 2D IPS (rc=10 A). It is demonstrated that 3D IPS is not practical for any of the interfacial systems studied. However, the hybrid method PME/IPS [Particle Mesh Ewald for electrostatics and 3D IPS for Lennard-Jones (LJ) interactions] provides an efficient way to include both types of long-range forces in simulations of large liquid/vacuum and all liquid/liquid interfaces, including lipid monolayers and bilayers. A previously published pressure-based long-range LJ correction yields results similar to those of PME/IPS for liquid/liquid interfaces. The contributions to surface tension of LJ terms arising from interactions beyond 10 A range from 13 dyn/cm for the hexadecane/vapor interface to approximately 3 dyn/cm for hexadecane/water and DPPC bilayers and monolayers. Surface tensions of alkane/vapor, hexadecane/water, and DPPC monolayers based on the CHARMM lipid force fields agree very well with experiment, whereas surface tensions of the TIP3P and TIP4P-Ew water models underestimate experiment by 16 and 11 dyn/cm, respectively. Dipole potential drops (DeltaPsi) are less sensitive to long-range LJ interactions than surface tensions. However, DeltaPsi for the DPPC bilayer (845+/-3 mV proceeding from water to lipid) and water (547+/-2 mV for TIP4P-Ew and 521+/-3 mV for TIP3P) overestimate experiment by factors of 3 and 5, respectively, and represent expected deficiencies in nonpolarizable force fields.     

Structural characterization of interfacial n-octanol and 3-octanol using molecular dynamic simulations

Structurally isomeric octanol interfacial systems, water/vapor, 3-octanol/vapor, n-octanol/vapor, 3-octanol/water, and n-octanol/water are investigated at 298 K using molecular dynamics simulation techniques. The present study is intended to investigate strongly associated liquid/liquid interfaces and probe the atomistic structure of these interfaces. The octanol and water molecules were initially placed randomly into a box and were equilibrated using constant pressure techniques to minimize bias within the initial conditions as well as to fully sample the structural conformations of the interface. An interface formed via phase separation during equilibration and resulted in a slab geometry with a molecularly sharp interface. However, some water molecules remained within the octanol phase with a mole fraction of 0.12 after equilibration. The resulting "wet" octanol interfaces were analyzed using density profiles and orientational order parameters. Our results support the hypothesis of an ordered interface only 1 or 2 molecular layers deep before bulk properties are reached for both the 3-octanol and water systems. However, in contrast to most other interfacial systems studied by molecular dynamics simulations, the n-octanol interface extends for several molecular layers. The octanol hydroxyl groups form a hydrogen-bonding network with water which orders the surface molecules toward a preferred direction and produces a hydrophilic/hydrophobic layering. The ordered n-octanol produces an oscillating low-high density of oxygen atoms out of phase with a high-low density of carbon atoms, consistent with an oscillating dielectric. In contrast, the isomeric 3-octanol has only a single carbon-rich layer directly proximal to the interface, which is a result of the different molecular topology. Both 3-octanol and n-octanol roughen the water interface with respect to the water/vapor interface. The "wet" octanol phases, in the octanol/water systems reach bulk properties in a shorter distance than the "dry" octanol/vapor interfaces.     

The extent of molecular orientation at liquid/vapor interface of pyridine and its alkyl derivatives by molecular dynamics simulation

In this paper, molecular dynamics simulation was performed to investigate the liquid∕vapor interfacial structure of neat polar liquids. Large-scale ensembles of liquid pyridine and its alkyl derivatives, 4-methylpyridine and 4-ethylpyridine, were simulated by classical molecular dynamics at 298 K. For the liquid system of low polarity, the surface density profile of the atoms meet exactly at the middle of interfacial region, and atoms of hydrophobic nature can be hardly discriminated from hydrophilic ones in either vapor or liquid sides. For a liquid system of high polarity, the density profile of atoms with different nature is highly discriminated all over the interfacial region, and as the polarity increases, a dense region of atomic density is clearly developed in the subsurface region. The recognized bivariate method was also used to study the molecular orientational distribution quantitatively. Orientational analysis of the three liquid systems indicates that the pyridine ring plane in the outmost surface tends to be vertical. Its tendency in the innermost interfacial region is parallel. The orientational states available to 4-ethylpyridine and pyridine are discriminated by predicting the possibility of a bisector-wise tumbling for the ring plane in pyridine and a side-wise tumbling in 4-ethylpyridine. The orientational distribution maps explain the trend of experimental surface tension and surface entropy. As the dipole moment of these liquids increases with the alkyl chain length, the surface structural profile changes from a regular definite one to a surface of complex atomic structure involving a dense phase near the interface. The development of dense region in alkyl derivatives is the result of segregation of molecules due to the alkyl group, which is captured and discriminated by molecular dynamics simulation even when the length of a short alkyl chain is increased by one carbon atom.     

Molecular Structure and Dynamics of Water at the Water–Air Interface Studied with Surface‐Specific Vibrational Spectroscopy

Water interfaces provide the platform for many important biological, chemical, and physical processes. The water–air interface is the most common and simple aqueous interface and serves as a model system for water at a hydrophobic surface. Unveiling the microscopic (<1 nm) structure and dynamics of interfacial water at the water–vapor interface is essential for understanding the processes occurring on the water surface. At the water interface the network of very strong intermolecular interactions, hydrogen‐bonds, is interrupted and the density of water is reduced. A central question regarding water at interfaces is the extent to which the structure and dynamics of water molecules are influenced by the interruption of the hydrogen‐bonded network and thus differ from those of bulk water. Herein, we discuss recent advances in the study of interfacial water at the water–air interface using laser‐based surface‐specific vibrational spectroscopy.     

Molecular level structure of the liquid/liquid interface. Molecular dynamics simulation and ITIM analysis of the water-CCl4 system

The molecular level properties of the liquid/liquid interface between water and CCl(4) are analysed in detail on the basis of molecular dynamics computer simulation. This analysis requires a full list of the molecules that are right at the interface in both phases. Such a list can be provided by the novel method for identifying truly interfacial molecules (ITIM). The full list of the truly interfacial molecules various properties (e.g., width, molecular level roughness) of the interface can be meaningfully analysed. The residence time of the molecules at the interface, the percolation of the water molecules at the interfacial layer as well as in the second layer beneath the surface, the preferred orientations of the interfacial water molecules and the dependence of these orientational preferences on the local curvature of the interface are also analysed and discussed in detail.     

Structure of the liquid-vapor interface of water-methanol mixtures as seen from Monte Carlo simulations

Monte Carlo simulation of the vapor-liquid interface of water-methanol mixtures of different compositions, ranging from pure water to pure methanol, have been performed on the canonical (N, V, T) ensemble at 298 K. The analysis of the systems simulated has revealed that the interface is characterized by a double layer structure: methanol is strongly adsorbed at the vapor side of the interface, whereas this adsorption layer is followed at its liquid side by a depletion layer of methanol of lower concentration than in the bulk liquid phase of the system. The dominant feature of the interface has been found to be the adsorption layer in systems of methanol mole fractions below 0.2, and the depletion layer in systems of methanol mole fractions between 0.25 and 0.5. The orientation of the molecules located at the depletion layer is found to be already uncorrelated with the interface, whereas the methanol molecules of the adsorption layer prefer to align perpendicular to the interface, pointing straight toward the vapor phase by their methyl group. Although both the preference of the molecular plane for a perpendicular alignment with the interface and the preference of the methyl group for pointing straight to the vapor phase are found to be rather weak, the preference of the methyl group for pointing as straight toward the vapor phase as possible within the constraint imposed by the orientation of the molecular plane is found to be fairly strong. One of the two preferred orientations of the interfacial water molecules present in the neat system is found to disappear in the presence of methanol, because methanol molecules aligned in their preferred orientation can replace these water molecules in the hydrogen-bonding pattern of the interface.     

Viscosity of the aqueous liquid/vapor interfacial region: 2D electrochemical measurements with a piperidine nitroxy radical probe

Surface partitioning of 2,2,6,6-tetramethyl-1-piperidynyloxy radical (Tempo) to the air/water interface follows a Langmuir isotherm. The partition constant was obtained by the surface tension measurements in the concentration range of 1.0 x 10(-4)-2.4 x 10(-3) M yielding K = 640 +/- 99 M(-1). The lateral mobility of Tempo at the air/water interface was measured electrochemically in the surface concentration range of 2.0 x 10(-11)-1.4 x 10(-10) mol/cm2, corresponding to ca. 7.3-50% full monolayer coverage. The measurements employed cyclic voltammetry with line microelectrodes touching the air/water interface. The Tempo lateral diffusion constant of (1.5 +/- 0.7) x 10(-4) cm2/s is independent of surface concentration below 4.0 x 10(-11) mol/cm2. The extent of Tempo water interactions was assessed by the electronic structure calculations. These calculations showed that, at most, two water molecules can hydrogen bond with the oxygen atom of the nitroxyl group of Tempo, and that a single water molecule forms a hydrogen bond that is ca. 30% stronger than the H2O-H2O hydrogen bond. These calculations led to a postulate that Tempo diffuses along the interface largely unimmersed, and that it is coupled to the interfacial water via hydrogen bonding with H2O. In view of this postulate, the viscosity of the aqueous liquid/vapor interfacial region obtained by interpreting the Tempo diffusion constant in the low concentration region is as much as 4 times smaller than that of bulk liquid water.     

Polarization and experimental configuration analyses of sum frequency generation vibrational spectra, structure, and orientational motion of the air/water interface

Here we report a detailed study on spectroscopy, structure, and orientational distribution, as well as orientational motion, of water molecules at the air/water interface, investigated with sum frequency generation vibrational spectroscopy (SFG-VS). Quantitative polarization and experimental configuration analyses of the SFG data in different polarizations with four sets of experimental configurations can shed new light on our present understanding of the air/water interface. Firstly, we concluded that the orientational motion of the interfacial water molecules can only be in a limited angular range, instead of rapidly varying over a broad angular range in the vibrational relaxation time as suggested previously. Secondly, because different vibrational modes of different molecular species at the interface has different symmetry properties, polarization and symmetry analyses of the SFG-VS spectral features can help the assignment of the SFG-VS spectra peaks to different interfacial species. These analyses concluded that the narrow 3693 cm(-1) and broad 3550 cm(-1) peaks belong to C(infinityv) symmetry, while the broad 3250 and 3450 cm(-1) peaks belong to the symmetric stretching modes with C2v symmetry. Thus, the 3693 cm(-1) peak is assigned to the free OH, the 3550 cm(-1) peak is assigned to the singly hydrogen-bonded OH stretching mode, and the 3250 and 3450 cm(-1) peaks are assigned to interfacial water molecules as two hydrogen donors for hydrogen bonding (with C2v symmetry), respectively. Thirdly, analysis of the SFG-VS spectra concluded that the singly hydrogen-bonded water molecules at the air/water interface have their dipole vector directed almost parallel to the interface and is with a very narrow orientational distribution. The doubly hydrogen-bonded donor water molecules have their dipole vector pointing away from the liquid phase.     

Polarizable contributions to the surface tension of liquid water

Surface tension, gamma, strongly affects interfacial properties in fluids. The degree to which polarizability affects gamma in water is thus far not well established. To address this situation, we carry out molecular dynamics simulations to study the interfacial forces acting on a slab of liquid water surrounded by vacuum using the Gaussian charge polarizable (GCP) model at 298.15 K. The GCP model incorporates both a fixed dipole due to Gaussian distributed charges and a polarizable dipole. We find a well-defined bulklike region forms with a width of approximately 31 A. The average density of the bulklike region agrees with the experimental value of 0.997 g/cm3. However, we find that the orientation of the molecules in the bulklike region is strongly influenced by the interfaces, even at a distance five molecular diameters from the interface. Specifically, the orientations of both the permanent and induced dipoles show a preferred orientation parallel to the interface. Near the interface, the preferred orientation of the dipoles becomes more pronounced and the average magnitude of the induced dipoles decreases monotonically. To quantify the degree to which molecular orientation affects gamma, we calculate the contributions to gamma from permanent dipolar interactions, induced dipolar interactions, and dispersion forces. We find that the induced dipole interactions and the permanent dipole interactions, as well as the cross interactions, have positive contributions to gamma, and therefore contribute stability to the interface. The repulsive core interactions result in a negative contribution to gamma, which nearly cancels the positive contributions from the dipoles. The large negative core contributions to gamma are the result of small oxygen-oxygen separation between molecules. These small separations occur due to the strong attractions between hydrogen and oxygen atoms. The final predicted value for gamma (68.65 m/Nm) shows a deviation of approximately 4% of the experimental value of 71.972 m/Nm. The inclusion of polarization is critical for this model to produce an accurate value.     

Effect of flexibility on surface tension and coexisting densities of water

Molecular dynamics simulations of pure water at the liquid-vapor interface are performed using direct simulation of interfaces in a liquid slab geometry. The effect of intramolecular flexibility on coexisting densities and surface tension is analyzed. The dipole moment profile across the liquid-vapor interface shows different values for the liquid and vapor phases. The flexible model is a polarizable model. This effect is minor for liquid densities and is large for surface tension. The liquid densities increase from 2% at 300 K to 9% at 550 K when the force field is changed from a fully rigid simple point charge extended (SPCE) model to that of a fully flexible model with the same intermolecular interaction parameters. The increases in surface tension at both temperatures are around 11% and 36%, respectively. The calculated properties of the flexible models are closer to the experimental data than those of the rigid SPCE. The effect of the maximum number of reciprocal vectors (h(z) (max)) and the surface area on the calculated properties at 300 K is also analyzed. The coexiting densities are not sensitive to those variables. The surface tension fluctuates with h(z) (max) with an amplitude larger than 10 mN m(-1). The effect of using small interfacial areas is slightly larger than the error in the simulations.     

Molecular dynamics simulation of the water|nitrobenzene interface

The structure and dynamics of the neat water|nitrobenzene liquid|liquid interface are studied at 300 K using molecular dynamics computer simulations. The water is modeled using the flexible SPC potential, and the nitrobenzene is modeled using an empirically determined nitrobenzene potential energy function. Although nitrobenzene is a polar liquid with a large dielectric constant, the structure of the interface is similar to other water|non-polar organic liquid interfaces. Among the main structural features we describe are an enhancement of interfacial water hydrogen bonds, the specific orientation of water dipoles and nitrobenzene molecules, and a rough surface that is locally sharp. Surface roughness is also characterized dynamically. The dynamics of molecular reorientation are shown to be only mildly modified at the interface. The effect due to the polarizable many-body potential energy functions of both liquids is investigated and is found to affect only mildly the above results.     

Surface tension of water droplets upon homogeneous droplet nucleation in water vapor

A method has been proposed for determining interfacial free energy from the data of molecular dynamics simulation. The method is based on the thermodynamic integration procedure and is distinguished by applicability to both planar interfaces and those characterized by a high curvature. The workability of the method has been demonstrated by the example of determining the surface tension for critical nuclei of water droplets upon condensation of water vapor. The calculation has been performed at temperatures of 273–373 K and a pressure of 1 atm, thus making it possible to determine the temperature dependence of the surface tension for water droplets and compare the results obtained with experimental data and the simulation results for a “planar” vapor–liquid interface.     

Adsorption of 1-octanol at the free water surface as studied by Monte Carlo simulation

The adsorption of 1-octanol at the free water surface has been investigated by Monte Carlo computer simulation. Six different systems, built up by an aqueous and a vapor phase, the latter also containing various number of octanol molecules, have been simulated. The number of the octanol molecules has been chosen in such a way that the octanol surface density varies in a broad range, between 0.27 and 7.83 micromol/m(2) in the six systems simulated. For reference, the interfacial system containing bulk liquid octanol in the apolar phase has also been simulated. The results have shown that the formation of hydrogen bonds between the interfacial water and adsorbed octanol molecules is of key importance in determining the properties of the adsorbed layer. At low octanol surface concentration values all the octanol molecules are strongly (i.e., by hydrogen bonds) bound to the aqueous phase, whereas their interaction with each other is negligibly small. Hence, they are preferentially oriented in such a way that their own binding energy (and thus their own free energy) is minimized. In this preferred orientation the O-H bond of the octanol molecule points flatly toward the aqueous phase, declining by about 30 degrees from the interfacial plane, irrespectively from whether the octanol molecule is the H-donor or the H-acceptor partner in the hydrogen bond. Hence, in its preferred orientation the octanol molecule can form at least two low energy hydrogen bonds with water: one as a H-donor and another one as a H-acceptor. Moreover, the preferred orientation of the hydrogen bonded water partners is close to one of the two preferred interfacial water alignments, in which the plane of the water molecule is parallel with the interface. When increasing the octanol surface density, the water surface gets saturated with hydrogen bonded octanols, and hence any further octanol molecule can just simply condense to the layer of the adsorbed octanols. The surface density value at which this saturation occurs is estimated to be about 1.7 micromol/m(2). Above this surface density value the hydrogen bonded octanols and their water partners are oriented in such a way that the number of the water-octanol hydrogen bonds is maximized. Hence, the preferred alignment of the O...O axes of these hydrogen bonds is perpendicular to the interface. This orientation is far from the optimal alignment of the individual octanol molecules, which is also reflected in the observed fact that, unlike in the case of many other adsorbents, the average molecular binding energy of the adsorbed octanol molecules increases (i.e., becomes less negative) with increasing octanol surface density.     

A molecular beam study of the evaporation of water from a liquid jet

A method to maintain a clean surface of a liquid in a high vacuum is described. Using a very thin and fast liquid jet it is not only possible to prevent freezing of the liquid but also to reduce the number of collisions between evaporating molecules to negligibly small values. Thus many of the standard, vacuum dependent, particle probing techniques for solid surfaces can be used for studies of rapidly vaporizing, high vapor pressure liquids. In a first molecular beam investigation we have used time-of-flight analysis to measure the velocity distribution of H₂O molecules vaporizing from thin jets of pure liquid water. The experiments were carried out for liquid jet diameters between 50 and 5 µm. In this range the expanding vapor is observed to undergo the transition to the collision-free molecular flow regime. From the measured velocity distributions the local surface temperature is determined to be less than 210 K. This appears to be the lowest temperature ever reported for supercooled liquid water.     

Hydrogen-bonded clusters on the vapor/ethanol-aqueous-solution interface

The structure of the vapor/ethanol-aqueous-solution interface has been carefully investigated focusing on an intermolecular hydrogen bond (HB) and molecular clusters bound by HBs. This paper is a continuation of our previous molecular dynamics (MD) study (Langmuir 2005, 21, 10885), and all analysis was performed based on five independent adsorption-equilibrated configurations of a slab of ethanol solution at 298.15 K, where the ethanol mole fraction of the solution, chi(e), is 0.0052, 0.012, 0.024, 0.057, and 0.12, respectively. The geometrical definition of HB enabled the detection of the HB between ethanol-ethanol, ethanol-water, and water-water molecules. The variations of the density of HB and the coordination number of HB across the vapor/solution interface were analyzed. Analysis on the density of HB reveals that a monolayer of adsorbed ethanol can be classified into two parts where ethanol molecules prefer to form HBs with each other and ethanol molecules prefer to form HBs with water molecules. Despite chi(e), the coordination number of ethanol-ethanol HB monotonically increases toward the vapor region, while those of ethanol-water and water-water HBs monotonically decrease. In addition, the variation of the mean size of both ethanol one-component clusters and ethanol/water binary clusters across the interface were analyzed. The mean size of an ethanol one-component cluster and that of an ethanol/water binary cluster are expressed as a maximum at the interface. These behaviors are linked with the size distributions of both one-component and binary clusters. A relatively large system in this calculation also enables detailed discussion about the molar dependency of the bulk structural properties of an ethanol solution.     

Molecular dynamics simulation of hydrated DPPC monolayers using charge equilibration force fields

We present results of molecular dynamics simulations of a model DPPC-water monolayer using charge equilibration (CHEQ) force fields, which explicitly account for electronic polarization in a classical treatment of intermolecular interactions. The surface pressure, determined as the difference between the monolayer and pure water surface tensions at 323 K, is predicted to be 22.92 ±1.29 dyne/cm, just slightly below the broad range of experimental values reported for this system. The surface tension for the DPPC-water monolayer is predicted to be 42.35 ±1.16 dyne/cm, in close agreement with the experimentally determined value of 40.9 dyne/cm. This surface tension is also consistent with the value obtained from DPPC monolayer simulations using state-of-the-art nonpolarizable force fields. The current results of simulations predict a monolayer-water potential difference relative to the pure water-air interface of 0.64 ±0.02 Volts, an improved prediction compared to the fixed-charge CHARMM27 force field, yet still overestimating the experimental range of 0.3 to 0.45 Volts. As the charge equilibration model is a purely charge-based model for polarization, the current results suggest that explicitly modeled polarization effects can offer improvements in describing interfacial electrostatics in such systems.     

Molecular simulation of an amorphous poly(methyl methacrylate)–poly(tetrafluoroethylene) interface

Molecular dynamics calculations of an amorphous interfacial system of poly(methyl methacrylate) (PMMA) and poly(tetrafluoroethylene) (PTFE) containing about 10,000 interaction sites were performed for 15 ns under constant pressure and constant temperature conditions. The time evolutions of the thickness, density and number of atomic pairs in the interfaces suggested that the interfaces reached their equilibrium states with an interfacial thickness of about 2 nm at 500 K. The molecular motion in the interface and bulk was compared using mean square displacement and torsional autocorrelation function. The separation at a PMMA/PTFE interface was mimicked using non-equilibrium molecular dynamics calculations by applying the potential energy to the MD cell in a direction perpendicular to the interface. Initially, the PTFE layer close to the interface was deformed, and before complete separation, some segments of the PTFE molecules extended from the bulk to the surface of the PMMA layer, which were attached by the intermolecular interaction. The remaining PTFE molecules were entangled in the bulk, which probably prevented the transfer of the PTFE molecules to the surfaces of the PMMA layers. On the other hand, the PMMA layer was only slightly deformed. This separation behavior can be explained by taking into account the intermolecular interaction, the barrier to the conformational changes of the backbones and the entanglement of the PTFE molecules in the bulk.