Molecular Dynamics Simulation of the Water/2-Heptanone Liquid-Liquid Interface

Water/organic liquid interface properties with amine, carboxyl, thiol, and methyl terminal groups as seen from MD simulations

Molecular dynamics simulations were performed to study structural and dynamic properties of polar butanamine/water/butanamine, pentanoic acid/water/pentanoic acid, butanethiol/water/butanethiol, and nonpolar pentane/water/pentane systems. The mass density profiles along the interface normal to the organic liquid/water system, the difference in the local structure of H2 O molecules in bulk and in the vicinity of interface, as well as the diffusion behavior of water molecules at the interface with above-mentioned organic liquids have been investigated. Our MD simulation has shown that the diffusion of water molecules across the water/organic liquid interface is influenced by the hydrogen bondsn HB between water molecules and the terminal groups of organic liquids. It was found that the loss of the hydrogen bondsn HB in the nonpolar organic liquid leads to a decrease in the value of the normal component of the diffusion coefficientD z , while the tangential diffusion coefficients, bothD x andD y , increase.

Molecular-Level Insight of CP52/NBR Damping Composites through a Combination of Molecular Dynamics Simulation and Experimental Method

To enhance the damping properties of nitrile butadiene rubber (NBR), the elastomer used was blended with chlorinated paraffin 52 (CP52) to prepare NBR/CP52 composites. The results showed that CP52 could significantly enhance the damping properties of NBR and shift the glass transition temperature (T_g) to lower temperatures. Molecular dynamics models of the CP52/NBR system were established, and the damping properties of the CP52-reinforced NBR were investigated using molecular dynamics (MD) simulations. Through the combination of MD simulations and the experimental results, the essential mechanism of the enhanced damping properties of the NBR was methodically expatiated and was ascribed to the Cl-CP-H····NC-NBR (type I) and CP-Cl····H-NBR-CN (type II) analogous hydrogen bonds formed between NBR and CP52. The higher the CP52 content, the higher the analogous hydrogen bond concentration, and the better the damping properties of the CP52/NBR composites. The experimental results were very consistent with the MD simulation results, meaning that the combination method can provide a new means to optimize the design of damping materials and broaden the application range of small polar molecules in the damping modification of polar rubber materials.

Liquid/liquid interface in periodic boundary condition

We study how surface phenomena can change the interface geometry in liquid-liquid two-phase systems with periodic boundary conditions. Without any curvature effect on surface tension, planar (slab), cylindrical, and spherical structures are successively obtained as a function of the total composition and elongation of the box, in accordance with molecular dynamics simulations for a water/heptane system. The curvature effects described by Tolman relationship desymmetrize the phase diagram by stabilizing a concavity but it leads to inconsistencies with high curvature. Helfrich model partially resolves this and predicts the possible presence of shells reflecting a frustrated system.

Molecular dynamics simulation study used in systems with supercritical water

Supercritical water (SCW) is a green solvent. The supercritical fluids have been increasingly concerned and studied in many areas such as SCW gasification, biofuel production, SCW hydrothermal conversion, organic wastes treatment and utilization, nanotechnology, etc. Because of the severe circumstances and rapid reactions in supercritical water, it is difficult for experimental researchers to disentangle various fundamental reaction steps from the intermediate and product distributions. From this perspective, molecular dynamics (MD) simulation based on quantum chemistry is an efficient tool for studying and exploring complex molecular systems. In recent years, molecular simulations and quantum chemical calculations have become powerful for illustrating the possible internal mechanism of a complex system. However, now there is no literature about the overview of MD simulation study of the system with SCW. Therefore, in this paper, an overview of MD simulation investigation applied in various systems with SCW is presented. In the current review we explore diverse research areas. Namely, the applications of MD simulation on investigating the properties of SCW, pyrolysis/gasification systems with SCW, dissolution systems and oxidation systems with SCW were summarized. And the corresponding problems in diverse systems were discussed. Furthermore, the advances and problems in MD simulation study were also discussed. Finally, possible directions for future research were outlined. This work is expected to be one reference for the further theoretical and molecular simulation investigations of systems involving SCW.

Microelectrochemical Techniques for Probing Kinetics at Liquid/Liquid Interfaces

We provide an overview of recent advances in microelectrochemical approaches to investigate the kinetics of various physicochemical processes that occur at the interface between two immiscible electrolyte solutions (ITIES). To place the advances in context, background material on the structure of the ITIES, derived from both experimental studies and computer simulation, is also provided. The main focus of the article is micro-ITIES techniques, single droplet measurements, microelectrochemical measurements at expanding droplets (MEMED) and scanning electrochemical microscopy (SECM). Recent developments in a combined SECM-Langmuir trough technique for probing diffusion processes across Langmuir monolayers at the water/air (W/A) interface are also highlighted, by considering an organic monolayer at a water surface as a special case of a liquid/liquid interface.

Quantitative Predictions of the Interfacial Tensions of Liquid-Liquid Interfaces through Atomistic and Coarse Grained Models

We report molecular simulations of oil-water liquid-liquid interfaces by using atomistic and coarse grained (CG) MARTINI force fields. We also apply the electronic continuum (EC) model to the MARTINI force field for the calculation of the interfacial tension of oil/water-salt systems. In a first step, we propose to calculate the interfacial tensions using thermodynamic and mechanical definitions of hydrocarbon-water interfacial systems modified by the addition of salts and alcohol. We also establish here the order of magnitude of the long-range corrections to the interfacial tension in fluid-fluid interfaces. Whereas the atomistic models are able to reproduce quantitatively the interfacial tension and the coexisting densities of oil-water systems, the coarse-description shows some deviations in the prediction of the interfacial tensions. Nevertheless, the physical features of these liquid-liquid interfaces are well-captured by this CG description. The CG force field offers then a very challenging alternative that will require however a more developed calibration of the parameters on the basis of liquid-liquid properties.

Calculation of the intrinsic solvation free energy profile of an ionic penetrant across a liquid-liquid interface with computer simulations

We introduce the novel concept of an intrinsic free energy profile, allowing one to remove the artificial smearing caused by thermal capillary waves, which renders difficulties for the calculation of free energy profiles across fluid interfaces in computer simulations. We apply this concept to the problem of a chloride ion crossing the interface between water and 1,2-dichloroethane and show that the present approach is able to reveal several important features of the free energy profile which are not detected with the usual, nonintrinsic calculations. Thus, in contrast to the nonintrinsic profile, a free energy barrier is found at the aqueous side of the (intrinsic) interface, which is attributed to the formation of a water "finger" the ion pulls with itself upon approaching the organic phase. Further, by the presence of a nonsampled region, the intrinsic free energy profile clearly indicates the coextraction of the first hydration shell water molecules of the ion when entering the organic phase.

Simple Ion Transfer at Liquid|Liquid Interfaces

The main aspects related to the charge transfer reactions occurring at the interface between two immiscible electrolyte solutions (ITIES) are described. The particular topics to be discussed involve simple ion transfer. Focus is given on theoretical approaches, numerical simulations, and experimental methodologies. Concerning the theoretical procedures, different computational simulations related to simple ion transfer are reviewed. The main conclusions drawn from the most accepted models are described and analyzed in regard to their relevance for explaining different aspects of ion transfer. We describe numerical simulations implementing different approaches for solving the differential equations associated with the mass transport and charge transfer. These numerical simulations are correlated with selected experimental results; their usefulness in designing new experiments is summarized. Finally, many practical applications can be envisaged regarding the determination of physicochemical properties, electroanalysis, drug lipophilicity, and phase-transfer catalysis.

Analysis of the subcritical carbon dioxide-water interface

We follow the evolution of the H(2)O/CO(2) interface at 300 K from the low pressure limit to near-critical pressures in molecular dynamics simulations using the SPC water and EPM2 carbon dioxide models. The intrinsic structure of the interface is elucidated by accumulating density profiles relative to the fluctuating capillary wave surface. Our main finding is that a carbon dioxide film of increasing density and thickness grows in two stages at the interface while the structure of the water surface barely changes. At low density, the entire film density profile grows linearly with the bulk CO(2) density. This regime continues up to a bulk CO(2) density of roughly 0.00095 Å(-3). At pressures above this point, we observe a distinct second peak in the CO(2) density, along with a tail of excess density that decays exponentially with distance from the interface. The decay length of the exponential tail diverges with increasing CO(2) pressure according to an inverse power law decay. Over the entire range of pressures, the CO(2) film had no detectable effect on the orientational order of the water surface. As expected, when the film of excess CO(2) at the interface grows, we find that the surface tension drops with increasing pressure. This is in qualitative accord with existing measurements, although the rate at which the surface tension falls with increasing pressure according to the SPC and EPM2 models is too small, indicating that the surface excess of CO(2) is underestimated by these models.

Molecular level properties of the free water surface and different organic liquid/water interfaces, as seen from ITIM analysis of computer simulation results

Molecular dynamics simulations of the interface of water with four different apolar phases, namely water vapour, liquid carbon tetrachloride, liquid dichloromethane (DCM) and liquid dichloroethane (DCE) are performed on the canonical ensemble at 298 K. The resulting configurations are analysed using the novel method of identification of the truly interfacial molecules (ITIM). Properties of the first three molecular layers of the liquid phases (e.g. width, spacing, roughness, extent of the in-layer hydrogen bonding network) as well as of the molecules constituting these layers (e.g., dynamics, orientation) are investigated in detail. In the analyses, particular attention is paid to the effect of the polarity of the non-aqueous phase and to the length scale of the effect of the vicinity of the interface on the various properties of the molecules. The obtained results show that increasing polarity of the non-aqueous phase leads to the narrowing of the interface, in spite of the fact that, at the same time, the truly interfacial layer of water gets somewhat broader. The influence of the nearby interface is found to extend only to the first molecular layer in many respects. This result is attributed to the larger space available for the truly interfacial than for the non-interfacial molecules (as the shapes of the two liquid surfaces are largely independent of each other, resulting in the presence of voids between the two phases), and to the fact that the hydrogen bonding interaction of the truly interfacial water molecules with other waters is hindered in the direction of the interface.

Water/hydrocarbon interfaces: effect of hydrocarbon branching on single-molecule relaxation

Water/hydrocarbon interfaces are studied using molecular dynamics simulations in order to understand the effect of hydrocarbon branching on the dynamics of the system at and away from the interface. A recently proposed procedure for studying the intrinsic structure of the interface in such systems is utilized, and dynamics are probed in the usual laboratory frame as well as the intrinsic frame. The use of these two frames of reference leads to insight into the effect of capillary waves at the interface on dynamics. The systems were partitioned into zones with a width of 5 A, and a number of quantities of dynamical relevance, namely, the residence times, mean squared displacements, the velocity auto correlation functions, and orientational time correlations for molecules of both phases, were calculated in the laboratory and intrinsic frames at and away from the interface. For the aqueous phase, translational motion is found to be (a) diffusive at long times and not anomalous as in proteins or micelles, (b) faster at the interface than in the bulk, and (c) faster upon reduction of the effect of capillary waves. The rotational motion of water is (a) more anisotropic at the interface than in the bulk and (b) dependent on the orientation of the covalent O-H bond with respect to the plane of the interface. The effect of hydrocarbon branching on aqueous dynamics was found to be small, a result similar to the effect on the interfacial water structure. The hydrocarbon phase shows a larger variation for all dynamical probes, a trend consistent with their interfacial structure.

Molecular dynamics simulation of nanoparticle self-assembly at a liquid-liquid interface

We have used molecular dynamics simulations to investigate the in situ self-assembly of modified hydrocarbon nanoparticles (mean diameter of 1.2 nm) at a water-trichloroethylene (TCE) interface. The nanoparticles were first distributed randomly in the water phase. The MD simulation shows the in situ formation of nanoparticle clusters and the migration of both single particles and clusters from the water phase to the trichloroethylene phase, possibly due to the hydrophobic nature of the nanoparticles. Eventually, the single nanoparticles or clusters equilibrate at the water-TCE interface, and the surrounding liquid molecules pack randomly when in contact with the nanoparticle surfaces. In addition, the simulations show that the water-TCE interfacial thickness analyzed from density profiles is influenced by the presence of nanoparticles either near or in contact with the interface but is independent of the number of nanoparticles present. The nanoparticles, water molecules, and TCE molecules all exhibit diffusion anisotropy.

Improved syntheses of bis(beta-cyclodextrin) derivatives, new carriers for gadolinium complexes

In the last decade a number of reports have been published on the synthesis and characterization of bridged cyclodextrin dimers (bis-CDs) connected with linkers of different lengths and structures. These dimers, having two hydrophobic cavities in close proximity, display much higher binding affinities and molecular selectivities than parent CDs, forming stable supramolecular adducts. We describe new synthetic protocols for the preparation of bis(beta-CDs) bearing 2-2', 3-3' and 6-6' bridges. Some of the critical steps were carried out either under high-intensity ultrasound (US) or microwave (MW) irradiation. Bis(beta-CDs) containing 6-6' ureido- and thioureido-bridges were prepared in high yields by a MW-promoted aza-Wittig reaction using polymer-bound triphenylphosphine, while those containing 2,2' and 3,3' bridges were prepared from mono-alkenyl beta-CDs by the cross-metathesis reaction (homodimerization) in the presence of 2(nd)-generation Grubbs catalyst under sonochemical conditions. By these improved protocols CD dimers could be obtained in gram amounts to prepare stable adducts of bis-CDs with contrast agents (CAs) containing gadolinium(iii) chelates. In the case of Gd(iii) chleate "G-1" the inclusion complexes were found to be 2 to 3 orders of magnitude more stable than that formed by beta-CD (K(ass) = 4.3 x 10(4) M(-1)vs 8.0 x 10(2) M(-1)). Relaxivity increased as well by factors of 3 and 4, viz. from 9.1 mM(-1) s(-1) (beta-CD) to 27.7 and 35 mM(-1) s(-1).

Computer Simulation Investigation of the Water−Benzene Interface in a Broad Range of Thermodynamic States from Ambient to Supercritical Conditions

The dependence of the properties of the water-benzene system on the thermodynamic conditions in a broad range of temperatures and pressures has been investigated by computer simulation methods. For this purpose, Monte Carlo simulations have been performed at 23 different thermodynamic states, ranging from ambient to supercritical conditions. The density profiles of the water and benzene molecules have been determined at each of the thermodynamic states investigated. Information on the dependence of the mutual solubility of the two components in each other as well as of the width of the interface on the temperature and pressure has been extracted from these profiles. The width of the interface has been found to increase with increasing temperature up to a certain point, where it diverges. The temperature of this divergence corresponds to the mixing of the two phases. The determination of the critical mixing temperature at various pressures allowed us to estimate the upper critical curve, separating the two-phase and one-phase liquid systems, of the phase diagram of the simulated water-benzene system. In analyzing the preferential orientation of the interfacial molecules relative to the interface, it has been found that the main orientational preference of the benzene molecules is to lie parallel with the plane of the interface, and the water molecules penetrated deepest into the benzene phase prefer to stay perpendicular to the interface, pointing by one of their O-H bonds almost straight toward the benzene phase, whereas the waters located at the aqueous side of the interface are preferentially aligned parallel with the interfacial plane. Although the strength of the observed orientational preferences decreases rapidly with increasing temperature, the preferred orientations themselves are found to be independent of the thermodynamic conditions. Remains of the orientational preferences of the molecules are found to be present up to temperatures as high as 650 K. The analysis of the relative orientation of the neighboring water-benzene pairs has revealed that the radius of the first hydration shell of the benzene molecules is independent of the thermodynamic conditions, even if the system consists of one single phase. It has been found that the nearest water neighbors of the benzene molecules are preferentially located above and below the benzene ring, whereas more distant water neighbors, belonging still to the first hydration shell, prefer to stay within the plane of the benzene molecule. In the two-phase systems the dipole vector of the nearest waters has been found to be preferentially perpendicular to the vector pointing from the center of the benzene molecule to the water O atom.

Orientational order of the water molecules at the vicinity of the water–benzene interface in a broad range of thermodynamic states, as seen from Monte Carlo simulations

Monte Carlo simulation of the water/benzene liquid-liquid interfacial system has been performed at six different thermodynamic state points, ranging from ambient conditions up to the vicinity of the critical point of water. The system has been found to consist of two immiscible liquid phases at every state point studied. The orientational preferences of the interfacial water molecules have been analysed in detail using the simulated configurations. The results obtained at ambient conditions are in agreement with previous results on various different water/apolar interfaces. Thus, interfacial water molecules have been found to have dual orientational preferences: the molecules located nearest to the organic phase prefer to stay perpendicular to the interface, pointing flatly toward the apolar phase by their dipole vectors, whereas the waters located somewhat farther from the organic phase prefer the parallel alignment with the interface. The observed orientational preferences are found to be rather stable with changing thermodynamic conditions: although the increase of the temperature has led, due to the increasing thermal motion of the molecules, to a gradual weakening of the orientational preferences, both preferences are found to exist up to at least 450 K, and found to be completely washed out at 575 K only. The pressure has not been found to influence the orientation of the water molecules noticeably.