Next Article in Journal
Polysaccharide from Okra (Abelmoschus esculentus (L.) Moench) Improves Antioxidant Capacity via PI3K/AKT Pathways and Nrf2 Translocation in a Type 2 Diabetes Model
Next Article in Special Issue
Management of Exciton for Highly-Efficient Hybrid White Organic Light-Emitting Diodes with a Non-Doped Blue Emissive Layer
Previous Article in Journal
Fabrication and Characterisation of Stimuli Responsive Piezoelectric PVDF and Hydroxyapatite-Filled PVDF Fibrous Membranes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Dynamics Simulation of Distribution and Diffusion Behaviour of Oil–Water Interfaces

1
Key Laboratory of Energy Thermal Conversion and Control of Ministry of Education, School of Energy and Environment, Southeast University, Nanjing 210096, China
2
School of Hydraulic, Energy and Power Engineering, Yangzhou University, Yangzhou 225127, China
3
Dalian National Laboratory for Clean Energy, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China
4
Jiangsu Key Laboratory of Micro and Nano Heat Fluid Flow Technology and Energy Application, School of Environmental Science and Engineering, Suzhou University of Science and Technology, Suzhou 215009, China
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(10), 1905; https://doi.org/10.3390/molecules24101905
Submission received: 24 April 2019 / Revised: 11 May 2019 / Accepted: 16 May 2019 / Published: 17 May 2019

Abstract

:
The distribution and diffusion behaviors of microscopic particles at fluorobenzene–water and pentanol–water interfaces are investigated using molecular dynamics simulation. The influences of Na+/Cl ions and the steric effects of organic molecules are examined. The concentration distributions of different species, the orientations of oil molecules at the interface, and oil–water interface morphology as well as the diffusion behaviors of water molecules are explored and analyzed. The results indicate that a few fluorobenzene molecules move into the water phase influenced by Na+/Cl ions, while the pentanol molecules at the interface prefer orientating their hydrophilic groups toward the water phase due to their large size. The water molecules more easily burst into the pentanol phase with larger molecular spaces. As the concentration of ions in the water phase increases, more water molecules enter into the pentanol molecules, leading to larger interface roughness and interface thickness. In addition, a lower diffusion coefficient for water molecules at the fluorobenzene–water interface are observed when introducing Na+/Cl ions in the water phase, while for the pentanol–water system, the mobility of interfacial water molecules are enhanced with less ions and inhibited with more ions.

1. Introduction

Oil–water interfaces exist widely in physical, chemical, and biological processes, such as interface adsorption [1], extraction [2], emulsification [3,4], phase transfer [5,6], and component transport [7,8]. From a microscopic point of view, the particle diffusion at the oil–water interface is simultaneously affected by the same and heterogeneous phases nearby in the interface, thus presenting distribution and diffusion characteristics that differ from the particles inside the substance. This directly affects the micromorphology of the oil–water interface and the diffusion efficiency at the interface. Therefore, the investigation of the distribution and diffusion of microscopic particles at the oil–water interface is not only of scientific significance for understanding the structure and physicochemical characteristics of the oil–water interface, but also to provide key technical support for the regulation and optimization of various interface processes in practical applications [9,10,11,12].
Since the 19th century, extensive experimental studies have been conducted on interfacial tension, solution adsorption, liquid membrane separation, and other interfacial properties and processes [13,14]. The rapid development of modern surface/interface analysis techniques (such as surface microscopy, surface spectroscopy, etc.) has pushed surface/interface science into the era of molecular/atomic analysis. Atkin et al. [15] used atomic force microscopy (AFM) and tunneling scanning microscopy (TSM) to find a unique multilayer solvation layer at the interface between gold and ionic liquids. Han et al. [16] studied the electrocatalytic performance of platinum–carbon nanoparticles by means of modern electrochemical analysis. It was indicated that the electrocatalytic activity of platinum–carbon nanoparticles was closely related to their size and irregular surface. However, these surface/interface analysis techniques are mainly used to analyze and characterize the microscopic characteristics of fixed surfaces/interfaces. It is difficult to study the distribution and diffusion behavior of microscopic particles at liquid–liquid interfaces, such as oil–water interfaces.
Fortunately, molecular simulation provides another important means for exploring the distribution and movement of microparticles in the free interface region, and to some extent, makes up for the deficiency of modern surface/interface analysis technology. In recent years, molecular dynamics simulation has been used to study the distribution and movement of microscopic particles at the oil–water interface. Patel et al. [17] studied the structure orientation and hydrogen bonding of water molecules at the octane–water interface. The results indicated that the average number of hydrogen bonds and dielectric constant of water molecules at the interface were less than those in the water phase. Through the study of the hydrogen bond interaction between water molecules at the interface of dodecane phase and water phase, Chen et al. [18] pointed out that the slow hydrogen bond exchange between water molecules delayed the rearrangement of water molecules at the interface when the water phase and oil phase contact each other. Vatente et al. [19] used nonequilibrium molecular dynamics simulation to find that when potassium ions transfer from the water phase to the oil phase, it is necessary to strip the hydration layer molecules and overcome a certain energy barrier. Kadota et al. [20] pointed out that the mobility of ions at oil–water interface and the formation of organic molecular clusters could be effectively regulated by changing the types of organic molecules in the oil phase.
In conclusion, it has been confirmed that there are differences in the distribution and diffusion characteristics of microscopic particles between the oil–water interface and the oil–water phase, and some progress has been made in the molecular dynamics simulation of oil–water interfaces. However, the microscopic behaviors of particle distribution and diffusion in oil–water systems remain unclear. Especially, the mechanism behind the microscopic morphology of oil–water interfaces and the diffusion characteristics of microscopic particles near the interface when adding electrolyte ions in the aqueous phase are still to be explored. Therefore, molecular dynamics simulation is performed to study the distribution and diffusion behavior of microscopic particles at the interface of two phases in two oil–water systems of fluorobenzene–water and 1-pentanol–water, which are commonly used in pharmaceutical manufacturing. The effect of adding electrolyte ions (Na+/Cl) to water on the interface morphology and water diffusion characteristics was analyzed in order to provide a clearer understanding of the distribution and diffusion of microscopic particles at the oil–water interface, which provides a theoretical basis for the design and control of physicochemical interface processes in engineering practice.

2. Molecular Dynamics Simulation

The structure of the oil–water system, the microstructure of water molecules, Na+/Cl ions, and oil phase molecules used in the molecular dynamics simulation study are depicted in Figure 1. As shown in the figure, the oil–water system consists of an upper water phase and a lower oil phase, forming an oil–water interface perpendicular to the z direction. The upper water phase contains 2000 water molecules (1.0 g/cm3), while the lower oil phase contains 300 fluorobenzene molecules (1.0 g/cm3) or 200 1-pentanol molecules (0.8 g/cm3). In the simulation, periodic boundary conditions are applied in the x, y, and z directions. The total thickness of the oil and water phases along the z direction is 9.49 nm, of which the oil phase thickness is about 5.00 nm, the water phase thickness is about 4.49 nm, and the area of the oil–water interface perpendicular to the z direction is 4.00 × 4.00 nm2. In this way, the oil and water phases in the z direction have enough thickness so that the interaction between them and the influence caused by the oil–water interface in the adjacent periodic cell can be neglected. Different quantities of electrolyte ions (Na+/Cl) were added to the simulated water phase to investigate the impact of electrolyte ions on the distribution and diffusion characteristics of the microscopic particles at the oil–water interface. The number of molecules and ions in each simulation system is listed in Table 1. As is known, the mass concentration of saturated NaCl solution is 26.47% at room temperature. In our work, the simulation systems were all dilute solutions of NaCl (from 0.00% to 6.10%).
The interactions of microscopic particles decide the distribution and diffusion characteristics of particles at the oil–water interface; thus, the choice of molecular force field is quite significant in the calculation these properties. COMPASS (Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies) force field [21] is an ab initio molecular force field, which can better deal with the interactions among inorganic and organic molecules and ions, and is suitable for the calculation of molecular structure and thermodynamic properties of condensed matter. Therefore, a COMPASS force field is used to calculate the interactions among water molecules, ions, and oil molecules. Bond length, bond angle, and dihedral angle are used to describe molecular structures. Lennard–Jones (LJ) potential and Coulomb potential are used to describe nonbonding interactions [22]. In order to eliminate unreasonable interaction between molecules and ions and optimize the initial particle distribution of the model system, steepest descent, conjugate gradient, and quasi-Newton methods were comprehensively used to optimize the initial configuration. On the basis of the initial configuration, we use canonical ensemble (NVT) to make the system balance via 1 ns of dynamic simulation, and then perform another 1 ns of molecular dynamics simulation. The simulation time step of which is 1 fs and the atomic coordinate information is saved every 10,000 steps. In the whole simulation process, the atom-based addition method is used to calculate Van der Waals forces, and the truncation radius is 1.2 nm. The Ewald summation method is applied in Coulomb force calculation [23] and the accuracy was 10−3 kcal/mol. A Nose heat bath [24] is adopted to control the temperature of the whole system during simulation at 298 K. All the sets and calculations stated above were conducted in the Materials Studio 6.0 packages [25].
In order to verify the simulation method in our work, we simulate the molecular distributions of an octane–water system as provided by Valente [19]. The concentration distributions are compared in Figure 2. The agreement of the interfacial distribution profiles of water and octane molecules verifies the rationality of the method used in our work.

3. Results and Discussion

3.1. Microscopic Particle Distribution and Oil–Water Interface Morphology

In order to establish an intuitive understanding of the distribution characteristics of various microscopic particles at the oil–water interface, the interfacial microstructures of the fluorobenzene–water system and the pentanol–water system at equilibrium state with different electrolyte ion contents are presented in Figure 3. As shown in Figure 3a–d, fluorobenzene molecules are closely packed in the oil phase in the fluorobenzene–water system due to their planar configuration and small volume. Owing to the attraction of hydrogen bonds of adjacent water molecules, the arrangement of fluorobenzene molecules at the interface is looser than that of fluorobenzene molecules in the oil phase. When sodium ions and chloride ions exist in the aqueous phase, a small amount of fluorobenzene molecules enter the aqueous phase (as shown in Figure 3b) because of the additional ion–fluorobenzene electrostatic interaction. However, with the further increase of electrolyte ions (Na+/Cl) in the aqueous phase, the density of aqueous phase layer increases, the electrostatic interaction between water molecules and anion and cation increases, thus the binding between water molecules and ions makes it difficult for fluorobenzene molecules to enter the aqueous phase (as shown in Figure 3d).
Compared with the fluorobenzene–water system, the pentanol molecule in the pentanol–water system has a larger volume and needs to overcome much larger steric hindrance in the process of movement and migration. Therefore, no pentanol molecule has been observed in the water phase, regardless of whether ions are added or not, as shown in Figure 3e–h. In addition, the amphiphilic nature of pentanol molecules makes the pentanol molecule at the interface point to the aqueous phase mainly with hydrophilic hydroxyl groups, as shown in Figure 3e. Moreover, owing to the large volume of pentanol molecules, loose accumulation, and large gaps between pentanol molecules at the interface, water molecules often enter the pentanol phase, as shown in Figure 3g.
Various particle concentration curves at the oil–water interface directly reflect the microscopic distribution characteristics of particles. The relative concentration distributions of various particles along the z direction is directedly shown in Figure 4 and Figure 5, respectively (Cr, the ratio of the concentration of microscopic particles at the interface to the average concentration in the aqueous phase). The concentration curves of fluorine atoms (F) and α carbon atoms (C1) characterize the concentration distributions of fluorobenzene molecules and the concentration curves of hydroxyl (-OH) and methyl (-CH3) groups characterize the concentration distribution of pentanol molecules, while the concentration curve of oxygen atom (O) represents the concentration distribution of water molecules.
It can be seen from the figures that there is a interfacial zone between the oil and water phase (defined as the area where the concentration of oil phase molecule or water phase molecule is 10%–90% of their respective bulk phase concentration), that is, the oil–water interface region [26]. It can be seen that the oil–water interface thickness of the fluorobenzene–water system decreases with the increase of electrolyte ions (Na+/Cl), while the oil–water interface thickness of the pentanol–water system increases with the increase of sodium ions and chloride ions, which is because the pentanol molecule has a longer alkyl chain than the fluorobenzene molecule, larger steric hindrance, and looser accumulation. In this way, when the ions in the water phase increase, the water molecules are squeezed, and the number of water molecules entering the interior of the pentanol phase increases, resulting in an increase of the volume of the oil–water interface. In addition, compared with the concentration curves of fluorine atoms (F) and α carbon atoms (C1), the concentration distribution curves of hydroxyl (-OH) and alkyl (-CH3) groups of pentanol molecules have concentration peaks near the oil–water interface, which are near the water and oil phases, respectively, indicating again that the amphiphilicity of pentanol molecule leads to an obvious orientation selectivity at the oil–water interface, that is, hydroxyl groups point towards the water phase and alkyl chains point towards the oil phase. Moreover, as can be seen from Figure 4 and Figure 5, the concentration curves of F and C1 at the interface basically coincide, which indicates that the fluorobenzene molecules at the interface do not show a specific orientation due to their small steric hindrance and easily achievable movement and rotation.
The rough morphology of the oil–water interface is also an important feature of the microscopic distribution of particles at the interface. Therefore, based on the spatial distribution of various particles in the oil–water two-phase systems, the oil–water interface boundary curves are generated by projecting the oil–water two-phase system to the x-y-z plane, as shown in Figure 6 and Figure 7. Here, the root mean square roughness Rms [27] is especially defined to quantitatively characterize the roughness of an oil–water interface:
R ms = 1 L 0 L [ z ( x ) z ¯ ( x ) ] 2 d x
Among them, L is the length in the x direction and z ( x ) z ¯ ( x ) is the difference between the z coordinates of each position on the interface boundary curve and the average of z. The z coordinates of different positions were obtained by Get Data software. Thus, for the fluorobenzene–water system (as shown in Figure 6), the roughness of the oil–water interface curve fluctuates slightly (Rms = 0.07–0.17), and the root mean square roughness reaches the minimum (Rms = 0.07) when the number of ions in the water phase reaches 40. It can be inferred that adding more ions in the water phase can reduce the roughness of the interface, which is due to the tight accumulation of molecules at the interface caused by the extrusion of water relative to the oil phase. Compared with the fluorobenzene–water system, the interfacial boundary curve of the pentanol–water system (as shown in Figure 7) fluctuates more obviously and its roughness is greater, which is mainly due to the large void of pentanol molecules, loose accumulation, and easy entry of water molecules into the void to form an irregular “erosion” morphology.

3.2. Diffusion Characteristics of Microscopic Particles at the Oil–water Interface

Diffusion and migration of microscopic particles at oil–water interfaces are important for the distribution and rough morphology of those particles. The mean squared displacement (MSD) of microscopic particles can be used to characterize the difficulty of particle diffusion. The MSD (formula 2) of water molecules in the water phase and at the interface of two oil–water systems is calculated in this paper, and the MSD curves are plotted in Figure 8 and Figure 9.
MSD = 1 N i = 1 N | r i ( t ) r i ( 0 ) | 2
On this basis, the diffusion coefficient (D) of water molecules at the interface is calculated as following [28]:
D = 1 6 lim t d M S D d t
In Formulas (2) and (3), the position coordinates of ri(0) and ri(t) are the position coordinates of the i-th water molecule at the initial time and the current time, respectively, and D is the diffusion coefficient of the water molecule. Accordingly, the corresponding diffusion coefficients are given in Table 2 (Di is the diffusion coefficients of water molecules at the interface and Db is the diffusion coefficients of water molecules in the water phase). From Formulas (2) and (3), it should be noted that for normal diffusion, the MSD-t plots should be linear and then the correct diffusion coefficient D can be obtained according to the hypothesis of Einstein function. Therefore, we extracted the data at the equilibrium stage.
By calculating the MSD curves and diffusion coefficients, it can be seen that adding electrolyte ions can reduce the diffusivity of water molecules in the fluorobenzene–water and pentanol–water systems (as shown in Figure 8, Figure 9 and Table 2), which is due to the fact that when electrolyte ions are added, the movement of water molecules not only overcomes the hydrogen bonding between them, but also overcomes the electrostatic interaction exerted by electrolyte ions. However, for water molecules at the oil–water interface, there are quite different situations in these two systems, that is, the diffusion ability of water molecules at the fluorobenzene–water interface decreases with the increase of the number of ions in the water phase, while the diffusion ability of water molecules at the pentanol–water interface increases first and then decreases with the increase of the number of ions in the water phase (as shown in Figure 8a,b and Table 2). In the fluorobenzene–water system, the extrusion of water molecules at the interface increases with the increase of the number of ions in the water phase. However, due to the tight accumulation of fluorobenzene molecules, this kind of squeezing action compresses the space of water molecules, resulting in a decrease in the diffusion capacity (as shown in Figure 3a–d and Figure 8). However, in the pentanol–water system, although the number of ions in the aqueous phase will still squeeze the water molecules at the interface, the water molecules can more easily enter the pentanol because of the loose accumulation of pentanol molecules at the interface. Therefore, this extrusion promotes the diffusion of water molecules into the voids of pentanol molecules, resulting in the enhancement of the diffusion ability of water molecules as the number of ions increases (as shown in Figure 3e–h and Figure 9). With the further increase of electrolyte ions in the aqueous phase, not only are the water molecules at the interface are squeezed, but so are the pentanol molecules. The gaps between pentanol molecules decrease, and the difficulty for water molecules to enter pentanol increases, resulting in a downward trend of diffusion ability of water molecules at the interface (Figure 9b).
In addition, the abnormal increase of the diffusion ability of water molecules at the pentanol–water interface with the increase of the number of ions in the water phase can be explained by radial distribution function (RDF) analysis. Figure 10a shows the g(r) of the hydrophilic and hydrophobic parts of fluorobenzene molecules around O (water) atoms for the fluorobenzene–water system with (without) Na+/Cl ions in the water phase. Larger peaks can be observed at ~0.3 nm in the g(r) of the hydrophilic groups and water molecules (F-O), especially for the system with Na+/Cl ions. This means the addition of NaCl leads to the compression of the interface and densely packing of the molecules. For the pentanol–water system (Figure 10b), more pronounced first peak was presented at the g(r) of hydrophilic part and water molecules (head-O). However, there is no obvious change in the g(r) profiles with and without NaCl in the water phase. It is suggested that the pentanol molecules at the interface orient their hydroxyl group toward the water phase. When adding NaCl in the water phase, some water molecules burst into the space of pentanol and the pentanol molecules retain regular configurations at the interface. Thus, the interfacial packing is not obviously changed.

4. Conclusions

In this paper, the distribution and diffusion behavior of microscopic particles at oil–water (fluorobenzene–water, 1-pentanol–water) two-phase interfaces were studied based on molecular dynamics simulation. The microscopic particle distribution, molecular orientation, interfacial roughness, and particle diffusion ability were analyzed. The effects of water molecular distribution and diffusion characteristics exerted on the types of organic molecules in the oil phase and the concentration of electrolyte ions (Na+/Cl) in water on the oil–water interface were examined. The conclusions are as follows:
(1) The addition of a small amount of electrolyte ions in the water phase will promote the movement of fluorobenzene molecules towards the water phase at the oil–water interface, while the pentanol molecules at the oil–water interface cannot easily enter the water phase because of their large molecular size. The pentanol molecules show obvious orientation selectivity at the oil–water interface, i.e., hydroxyl groups orientate towards the water phase and alkyl chains orientate towards the oil phase. For water molecules, the close accumulation of fluorobenzene makes it difficult for water molecules to enter the fluorobenzene phase, whereas the larger voids between pentanol molecules facilitate the entry of water molecules.
(2) There is a microscopic particle interfacial zone at the oil–water interface. With the increase of the number of electrolyte ions, the water molecules at the interface are squeezed, and the accumulation of fluorobenzene is denser, which is not conducive to the entry of water molecules, and the thickness of the interfacial zone decreases. In contrast, the increase of the number of electrolyte ions promotes more water molecules to enter the pentanol phase with larger voids, which leads to the increase of the thickness of the interfacial zone.
(3) The root mean square roughness of oil–water interface increases first and then decreases with increasing of ion concentration in the water phase, and the interface roughness of the oil–water interface is larger than that of pentanol molecules.
(4) Whether in the fluorobenzene–water system or in the pentanol–water system, the diffusion coefficient of water molecules in the water phase decreases with the increase of the number of ions; the diffusion coefficient of water molecules in the fluorobenzene–water system decreases with the increase of the number of ions at the oil–water interface, while the diffusion coefficient of water molecules increases first and then decreases with the increase of the number of ions at the pentanol–water interface.

Author Contributions

Investigation, C.Z., H.D. and P.L.; writing—original draft preparation, L.W. and B.Z.; supervision, C.Y.

Funding

This research was funded by National Natural Science Foundation of China (Grant No. 51776037 and 51806147), Natural Science Foundation of Jiangsu Province (Grant No. BK20170082), “Six talent peaks” Project of Jiangsu Province(Grant No. XNY-042), and “Zhishan Young Scholar” Program of Southeast University.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Naskar, B.; Diat, O.; Nardello-Rataj, V.; Bauduin, P. Nanometer-size polyoxometalate anions adsorb strongly on neutral soft surfaces. J. Phys. Chem. C 2015, 119, 20985–20992. [Google Scholar] [CrossRef]
  2. Vichapong, J.; Burakham, R.; Srijaranai, S. Vortex-assisted surfactant-enhanced-emulsification liquid-liquid microextraction with solidification of floating organic droplet combined with hplc for the determination of neonicotinoid pesticides. Talanta 2013, 117, 221–228. [Google Scholar] [CrossRef]
  3. Chaudhary, J.P.; Vadodariya, N.; Nataraj, S.K.; Meena, R. Chitosan-based aerogel membrane for robust oil-in-water emulsion separation. Acs Appl. Mater. Interfaces 2015, 7, 24957–24962. [Google Scholar] [CrossRef]
  4. Chen, Y.; Liu, X.; Shi, M. Hydrodynamics of double emulsion droplet in shear flow. Appl. Phys. Lett. 2013, 102, 051609. [Google Scholar] [CrossRef]
  5. Zhao, Q.Q.; Sun, J.; Liu, B.J.; He, J.X. Novel kinetics model for third-liquid phase-transfer catalysis system of the “complex” carbanion: Competitive role between catalytic cycles. Chem. Eng. J. 2015, 280, 782–795. [Google Scholar] [CrossRef]
  6. De-La-Torre, P.; Osorio, E.; Alzate-Morales, J.H.; Caballero, J.; Trilleras, J.; Astudillo-Saavedra, L.; Brito, I.; Cardenas, A.; Quiroga, J.; Gutierrez, M. Ultrasound-assisted phase-transfer catalysis method in an aqueous medium to promote the knoevenagel reaction: Advantages over the conventional and microwave-assisted solvent-free/catalyst-free method. Ultrason. Sonochem. 2014, 21, 1666–1674. [Google Scholar] [CrossRef]
  7. Rodriguez-Aller, M.; Guillarme, D.; Veuthey, J.L.; Gurny, R. Strategies for formulating and delivering poorly water-soluble drugs. J. Drug Deliv. Sci. Technol. 2015, 30, 342–351. [Google Scholar] [CrossRef]
  8. Chen, Y.P.; Deng, Z.L. Hydrodynamics of a droplet passing through a microfluidic t-junction. J. Fluid Mech. 2017, 819, 401–434. [Google Scholar] [CrossRef]
  9. Sharipova, A.; Aidarova, S.; Fainerman, V.B.; Stocco, A.; Cernoch, P.; Miller, R. Dynamics of adsorption of polyallylamine hydrochloride/sodium dodecyl sulphate at water/air and water/hexane interfaces. Colloid Surf. A 2011, 391, 112–118. [Google Scholar] [CrossRef]
  10. Dani, A.; Keiser, G.; Yeganeh, M.; Maldarelli, C. Hydrodynamics of particles at an oil–water interface. Langmuir 2015, 31, 13290–13302. [Google Scholar] [CrossRef] [PubMed]
  11. Park, B.J.; Furst, E.M. Fabrication of unusual asymmetric colloids at an oil–water interface. Langmuir 2010, 26, 10406–10410. [Google Scholar] [CrossRef] [PubMed]
  12. Abimannan, P.; Selvaraj, V.; Rajendran, V. Sonication effect on the reaction of 4-bromo-1-methylbenzene with sodium sulfide in liquid-liquid multi-site phase-transfer catalysis condition - kinetic study. Ultrason. Sonochem. 2015, 23, 156–164. [Google Scholar] [CrossRef]
  13. Wang, J.; Sun, L.; Zou, M.; Gao, W.; Liu, C.; Shang, L.; Gu, Z.; Zhao, Y. Bioinspired shape-memory graphene film with tunable wettability. Sci. Adv. 2017, 3, e1700004. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, Y.; Gao, W.; Zhang, C.; Zhao, Y. Three-dimensional splitting microfluidics. Lab Chip 2016, 16, 1332–1339. [Google Scholar] [CrossRef]
  15. Atkin, R.; El Abedin, S.Z.; Hayes, R.; Gasparotto, L.H.S.; Borisenko, N.; Endres, F. Afm and stm studies on the surface interaction of bmp tfsa and (emim tfsa ionic liquids with au(111). J. Phys. Chem. C 2009, 113, 13266–13272. [Google Scholar] [CrossRef]
  16. Han, K.S.; Moon, Y.S.; Han, O.H.; Hwang, K.J.; Kim, I.; Kim, H. Heat treatment and potential cycling effects on surface morphology, particle size, and catalytic activity of pt/c catalysts studied by c-13 nmr, tem, xrd and cv. Electrochem. Commun. 2007, 9, 317–324. [Google Scholar] [CrossRef]
  17. Patel, H.A.; Nauman, E.B.; Garde, S. Molecular structure and hydrophobic solvation thermodynamics at an octane-water interface. J. Chem. Phys. 2003, 119, 9199–9206. [Google Scholar] [CrossRef]
  18. Chen, M.; Lu, X.C.; Liu, X.D.; Hou, Q.F.; Zhu, Y.Y.; Zhou, H.Q. Retardation of water reorientation at the oil/water interface. J. Phys. Chem. C 2015, 119, 16639–16648. [Google Scholar] [CrossRef]
  19. Valente, M.; Sousa, S.F.; Magalhaes, A.L.; Freire, C. Transfer of the k+ cation across a water/dichloromethane interface: A steered molecular dynamics study with implications in cation extraction. J. Phys. Chem. B 2012, 116, 1843–1849. [Google Scholar] [CrossRef] [PubMed]
  20. Kadota, K.; Wake, T.; Gonda, K.; Kitayama, A.; Tozuka, Y.; Shimosaka, A.; Shirakawa, Y.; Hidaka, J. Effect of organic solvent on mutual diffusion and ionic behavior near liquid-liquid interface by molecular dynamics simulations. J. Mol. Liq. 2014, 197, 243–250. [Google Scholar] [CrossRef]
  21. Sun, H. Compass: An ab initio force-field optimized for condensed-phase applications - overview with details on alkane and benzene compounds. J. Phys. Chem. B 1998, 102, 7338–7364. [Google Scholar] [CrossRef]
  22. Rigby, D.; Sun, H.; Eichinger, B.E. Computer simulations of poly(ethylene oxide): Force field, pvt diagram and cyclization behaviour. Polym. Int. 1997, 44, 311–330. [Google Scholar] [CrossRef]
  23. Nymand, T.M.; Linse, P. Ewald summation and reaction field methods for potentials with atomic charges, dipoles, and polarizabilities. J. Chem. Phys. 2000, 112, 6152–6160. [Google Scholar] [CrossRef]
  24. Nose, S. A unified formulation of the constant temperature molecular-dynamics methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef]
  25. Materials Studio Accelrys, Release 6.0; Accelrys Software Inc.: San Diego, CA, USA, 2012.
  26. Choudhuri, J.R.; Chandra, A. Structure and dynamics of the liquid-liquid interface of an aqueous nacl solution with liquid carbon tetrachloride from first-principles simulations. J. Phys. Chem. C 2014, 118, 23083–23091. [Google Scholar] [CrossRef]
  27. Li, C.G.; Dong, S.; Zhang, G.X. Evaluation of the root-mean-square slope of 3d surface topography. Int. J. Mach. Tool Manu. 2000, 40, 445–454. [Google Scholar] [CrossRef]
  28. Abdi, N.; Abdi, Y.; Alemipour, Z.; NedaaeeOskoee, E. Chemical diffusion coefficient in dye sensitized solar cells as a function of porosity and surface roughness. Sol. Energy 2016, 135, 506–511. [Google Scholar] [CrossRef]
Sample Availability: Samples of the fluorobenzene–water and pentanol–water are not available from the authors.
Figure 1. Model and molecular structure of the oil–water interface.
Figure 1. Model and molecular structure of the oil–water interface.
Molecules 24 01905 g001
Figure 2. Comparison of concentration profiles of water and octane molecules at the oil–water interface between the present simulation and Valente’s work.
Figure 2. Comparison of concentration profiles of water and octane molecules at the oil–water interface between the present simulation and Valente’s work.
Molecules 24 01905 g002
Figure 3. Simulation of oil–water interface in different systems after 1 ns molecular dynamics simulation: (ad) Fluorobenzene–NaCl aqueous solution system containing different concentrations of ions; (eh) 1-Pentanol–NaCl aqueous solution system containing different concentrations of ions.
Figure 3. Simulation of oil–water interface in different systems after 1 ns molecular dynamics simulation: (ad) Fluorobenzene–NaCl aqueous solution system containing different concentrations of ions; (eh) 1-Pentanol–NaCl aqueous solution system containing different concentrations of ions.
Molecules 24 01905 g003
Figure 4. Concentration distribution curve of F, C, O, Na+, and Cl at the oil–water interface in the fluorobenzene–NaCl aqueous solution system: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4.
Figure 4. Concentration distribution curve of F, C, O, Na+, and Cl at the oil–water interface in the fluorobenzene–NaCl aqueous solution system: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4.
Molecules 24 01905 g004
Figure 5. Concentration distribution curve of F, C, O, Na+, and Cl at the oil–water interface in the 1-pentanol–NaCl aqueous solution system: (a) Case 5; (b) Case 6; (c) Case 7; (d) Case 8.
Figure 5. Concentration distribution curve of F, C, O, Na+, and Cl at the oil–water interface in the 1-pentanol–NaCl aqueous solution system: (a) Case 5; (b) Case 6; (c) Case 7; (d) Case 8.
Molecules 24 01905 g005
Figure 6. Interface morphology of the fluorobenzene–NaCl aqueous solution system with different contrations of ions: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4.
Figure 6. Interface morphology of the fluorobenzene–NaCl aqueous solution system with different contrations of ions: (a) Case 1; (b) Case 2; (c) Case 3; (d) Case 4.
Molecules 24 01905 g006
Figure 7. Interface morphology of the 1-pentanol–NaCl aqueous solution system with different concentrations of ions: (a) Case 5; (b) Case 6; (c) Case 7; (d) Case 8.
Figure 7. Interface morphology of the 1-pentanol–NaCl aqueous solution system with different concentrations of ions: (a) Case 5; (b) Case 6; (c) Case 7; (d) Case 8.
Molecules 24 01905 g007
Figure 8. Mean displacement curve of water molecules in: (a) Fluorobenzene–NaCl aqueous solution system; (b) 1-Pentanol–NaCl aqueous solution system.
Figure 8. Mean displacement curve of water molecules in: (a) Fluorobenzene–NaCl aqueous solution system; (b) 1-Pentanol–NaCl aqueous solution system.
Molecules 24 01905 g008
Figure 9. Mean displacement curve of water molecules at the interface: (a) Fluorobenzene–NaCl aqueous solution system; (b) 1-Pentanol–NaCl aqueous solution system.
Figure 9. Mean displacement curve of water molecules at the interface: (a) Fluorobenzene–NaCl aqueous solution system; (b) 1-Pentanol–NaCl aqueous solution system.
Molecules 24 01905 g009
Figure 10. The radial distribution functions g(r) of hydrophilic and hydrophobic parts of organic solvents and water molecules at different oil–water interfaces: (a) The g(r) of C1-O and F-O in the fluorobenzene–water system; (b) The g(r) of head-O and tail-O in 1-the pentanol–water system.
Figure 10. The radial distribution functions g(r) of hydrophilic and hydrophobic parts of organic solvents and water molecules at different oil–water interfaces: (a) The g(r) of C1-O and F-O in the fluorobenzene–water system; (b) The g(r) of head-O and tail-O in 1-the pentanol–water system.
Molecules 24 01905 g010
Table 1. The number of microscopic particles contained in the various oil–water systems.
Table 1. The number of microscopic particles contained in the various oil–water systems.
Fluorobenzene–NaCl Aqueous SolutionCase NumberCase 1Case 2Case 3Case 4
Fluorobenzene300300300300
Water Molecule2000200020002000
Sodium Ion0102040
Chloride0102040
Mass Concentration------1.60%3.15%6.10%
1-Pentanol–NaCl Aqueous SolutionCase NumberCase 5Case 6Case 7Case 8
1-Pentanol200200200200
Water Molecule2000200020002000
Sodium Ion0102040
Chloride0102040
Mass Concentration------1.60%3.15%6.10%
Table 2. Diffusion coefficients (in three-dimensional space) of water molecules in the water phase and at the interface with different oil–water interface systems (unit: 10−10 m2/s).
Table 2. Diffusion coefficients (in three-dimensional space) of water molecules in the water phase and at the interface with different oil–water interface systems (unit: 10−10 m2/s).
Diffusion CoefficientDiDb
Fluoro benzene–NaCl Aqueous SolutionCase 10.450.50
Case 20.370.39
Case 40.280.32
Diffusion CoefficientDiDb
1-Pentanol–NaCl Aqueous SolutionCase 50.330.43
Case 60.470.40
Case 80.380.35

Share and Cite

MDPI and ACS Style

Zhang, C.; Dai, H.; Lu, P.; Wu, L.; Zhou, B.; Yu, C. Molecular Dynamics Simulation of Distribution and Diffusion Behaviour of Oil–Water Interfaces. Molecules 2019, 24, 1905. https://doi.org/10.3390/molecules24101905

AMA Style

Zhang C, Dai H, Lu P, Wu L, Zhou B, Yu C. Molecular Dynamics Simulation of Distribution and Diffusion Behaviour of Oil–Water Interfaces. Molecules. 2019; 24(10):1905. https://doi.org/10.3390/molecules24101905

Chicago/Turabian Style

Zhang, Chengbin, Hanhui Dai, Pengfei Lu, Liangyu Wu, Bo Zhou, and Cheng Yu. 2019. "Molecular Dynamics Simulation of Distribution and Diffusion Behaviour of Oil–Water Interfaces" Molecules 24, no. 10: 1905. https://doi.org/10.3390/molecules24101905

APA Style

Zhang, C., Dai, H., Lu, P., Wu, L., Zhou, B., & Yu, C. (2019). Molecular Dynamics Simulation of Distribution and Diffusion Behaviour of Oil–Water Interfaces. Molecules, 24(10), 1905. https://doi.org/10.3390/molecules24101905

Article Metrics

Back to TopTop