1
|
Lin M, Tian B, Huang R, Xiao C. Study on the Transport Properties of SO 2 and NO at the Interface of H 2O 2 Solutions Using Molecular Dynamics. J Phys Chem B 2024. [PMID: 38656112 DOI: 10.1021/acs.jpcb.4c00013] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2024]
Abstract
Gas-liquid interfaces have a unique structure different from the bulk phase due to the complex intermolecular interactions within them and are regarded as barriers that prevent gases from entering solution or as channels that affect gas reactions. In this study, the adsorption and mass-transfer mechanisms of sulfur dioxide and nitric oxide at the gas-liquid interface of a H2O2 solution were comprehensively analyzed using molecular dynamics (MD) simulations. The analysis on molecule angle showed that H2O molecules tended to align parallel to the solution surface on the surface of the H2O2 solution. Regardless of whether the gas was adsorbed on the surface of the solution or not, H2O2 molecules were always perpendicular to the interface of the solution. The analysis on molecule angle and radial distribution function revealed that the H atoms of H2O molecules had a corresponding turn, and SO2 molecules were greatly affected by the attraction of H2O2 molecules during the adsorption of gas molecules on the interface. Steered MD was utilized to investigate the mass-transfer process of SO2 and NO molecules across the gas-liquid interface. The S atoms of SO2 molecules were significantly influenced by H2O2 molecules, while the O atoms of NO molecules gradually transitioned from the gas phase to the liquid phase. The results provided information on how gas molecules interacted with the surface of the solution and the specific details of the molecular orientation at the solution surface.
Collapse
Affiliation(s)
- Mingqi Lin
- Department of Energy and Power Engineering, College of Electrical Engineering, Guizhou University, Huaxi District, Guiyang 550025, China
| | - Bobing Tian
- Department of Energy and Power Engineering, College of Electrical Engineering, Guizhou University, Huaxi District, Guiyang 550025, China
| | - Ren Huang
- Department of Energy and Power Engineering, College of Electrical Engineering, Guizhou University, Huaxi District, Guiyang 550025, China
| | - Chao Xiao
- Department of Energy and Power Engineering, College of Electrical Engineering, Guizhou University, Huaxi District, Guiyang 550025, China
| |
Collapse
|
2
|
Zhao X, Wang L, Liu S, Luo X, Zhang M, Fu F, Gao X, Yan T. Dissolution Behavior of Polycyclic Aromatic Hydrocarbons in Heavy Oil in the Presence of Supercritical Cyclohexane. ACS OMEGA 2024; 9:252-263. [PMID: 38222523 PMCID: PMC10785659 DOI: 10.1021/acsomega.3c04101] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/10/2023] [Revised: 11/20/2023] [Accepted: 11/22/2023] [Indexed: 01/16/2024]
Abstract
Supercritical cyclohexane (SC-cyclohexane) shows significant advantages in mild operating conditions and the modulation of product distribution. To gain insights into the upgrading process of heavy oil in SC-cyclohexane, the dissolution process of polycyclic aromatic hydrocarbons (PAHs) contained in heavy oil was simulated based on molecular dynamics with the use of naphthalene, benzopyrene, and mixtures of naphthalene and benzopyrene as the model compounds. As indicated by the radial distribution function results, in SC-cyclohexane exhibiting low density, cyclohexane formed a solvent shell around PAHs such that the local concentration was reduced and the aggregation of PAHs was inhibited. The results of the solvation free energy suggested that van der Waals forces between PAHs and cyclohexane were mainly dominant. As revealed by the dissolution process of the model compounds in SC-cyclohexane, a low density and a suitable temperature contributed to the solubilization of PAHs. An appropriate temperature and a low density can be selected for the upgrading reaction to limit coke formation.
Collapse
Affiliation(s)
- Xiangbo Zhao
- Shaanxi
Key Laboratory of Chemical Reaction Engineering, School of Chemistry
and Chemical Engineering, Yan’an
University, Yan’an 716000, Shaanxi, China
- Shaanxi
Gold Group Xi’an Qinjin Co., Ltd., Xi’an 710300, Shaanxi, China
| | - Litao Wang
- Project
Office of Heavy to Light Conversion, Petrochemical Research Institute, Beijing 100007, China
| | - Shuai Liu
- Shaanxi
Key Laboratory of Chemical Reaction Engineering, School of Chemistry
and Chemical Engineering, Yan’an
University, Yan’an 716000, Shaanxi, China
| | - Xuan Luo
- Shaanxi
Key Laboratory of Chemical Reaction Engineering, School of Chemistry
and Chemical Engineering, Yan’an
University, Yan’an 716000, Shaanxi, China
| | - Mengran Zhang
- Shaanxi
Key Laboratory of Chemical Reaction Engineering, School of Chemistry
and Chemical Engineering, Yan’an
University, Yan’an 716000, Shaanxi, China
| | - Feng Fu
- Shaanxi
Key Laboratory of Chemical Reaction Engineering, School of Chemistry
and Chemical Engineering, Yan’an
University, Yan’an 716000, Shaanxi, China
| | - Xiaoming Gao
- Shaanxi
Key Laboratory of Chemical Reaction Engineering, School of Chemistry
and Chemical Engineering, Yan’an
University, Yan’an 716000, Shaanxi, China
| | - Ting Yan
- Shaanxi
Key Laboratory of Chemical Reaction Engineering, School of Chemistry
and Chemical Engineering, Yan’an
University, Yan’an 716000, Shaanxi, China
| |
Collapse
|
3
|
Rezlerová E, Moučka F, Předota M, Lísal M. Structure and self-diffusivity of alkali-halide electrolytes in neutral and charged graphene nanochannels. Phys Chem Chem Phys 2023; 25:21579-21594. [PMID: 37548441 DOI: 10.1039/d3cp03027j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/08/2023]
Abstract
Understanding the microscopic behaviour of aqueous electrolyte solutions in graphene-based ultrathin nanochannels is important in nanofluidic applications such as water purification, fuel cells, and molecular sensing. Under extreme confinement (<2 nm), the properties of water and ions differ drastically from those in the bulk phase. We studied the structural and diffusion behaviour of prototypical aqueous solutions of electrolytes (LiCl, NaCl, and KCl) confined in both neutral and positively-, and negatively-charged graphene nanochannels. We performed molecular dynamics simulations of the solutions in the nanochannels with either one, two- or three-layer water structures using the effectively polarisable force field for graphene. We analysed the structure and intermolecular bond network of the confined solutions along with their relation to the self-diffusivity of water and ions. The simulations show that Na and K cations can more easily rearrange their solvation shells under the graphene nanoconfinement and adsorb on the graphene surfaces or dissolve in the confinement-induced layered water than the Li cation. The negative surface charge together with the presence of ions orient water molecules with hydrogens towards the graphene surfaces, which in turn weakens the intermolecular bond network. The one-layer nanochannels have the biggest effect on the water structure and intermolecular bonding as well as on the adsorption of ions with only co-ions entering these nanochannels. The self-diffusivity of confined water is strongly reduced with respect to the bulk water and decreases with diminishing nanochannel heights except for the negatively-charged one-layer nanochannel. The self-diffusivity of ions also decreases with the reducing the nanochannel heights except for the self-diffusivity of cations in the negatively-charged one-layer nanochannel, evidencing cooperative diffusion of confined water and ions. Due to the significant break-up of the intermolecular bond network in the negatively-charged one-layer nanochannel, self-diffusion coefficients of water and cations exceed those for the two- and three-layer nanochannels and become comparable to the bulk values.
Collapse
Affiliation(s)
- Eliška Rezlerová
- Research Group of Molecular and Mesoscopic Modelling, The Czech Academy of Sciences, Institute of Chemical Process Fundamentals, Prague, Czech Republic.
- Department of Physics, Faculty of Science, Jan Evangelista Purkyně University in Úst nad Labem, Ústín. Lab., Czech Republic
| | - Filip Moučka
- Research Group of Molecular and Mesoscopic Modelling, The Czech Academy of Sciences, Institute of Chemical Process Fundamentals, Prague, Czech Republic.
- Department of Physics, Faculty of Science, Jan Evangelista Purkyně University in Úst nad Labem, Ústín. Lab., Czech Republic
| | - Milan Předota
- Department of Physics, Faculty of Science, University of South Bohemia, České Budě jovice, Czech Republic
| | - Martin Lísal
- Research Group of Molecular and Mesoscopic Modelling, The Czech Academy of Sciences, Institute of Chemical Process Fundamentals, Prague, Czech Republic.
- Department of Physics, Faculty of Science, Jan Evangelista Purkyně University in Úst nad Labem, Ústín. Lab., Czech Republic
| |
Collapse
|
4
|
Sookhaki E, Zolghadr AR, Namazian M. The effect of various partial atomic charges on the bulk and liquid/vacuum interface properties of iodobenzene derivatives at their melting points. J Mol Graph Model 2023; 119:108400. [PMID: 36586348 DOI: 10.1016/j.jmgm.2022.108400] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 11/04/2022] [Accepted: 12/20/2022] [Indexed: 12/28/2022]
Abstract
In this work, various atomic charge schemes including natural bond orbital (NBO), electrostatic potential based (CHELPG), and σ-hole model charges were applied in the OPLS-AA force field to investigate their effects on the thermophysical and structural properties of iodobenzene and its derivatives. Molecular dynamics simulations presented in this work show that the studied structural and thermophysical properties are in good agreement with experiments when the CHELPG charge was coupled with the OPLS-AA force field. Also, the arrangement of iodobenzene derivatives in the liquid phase was investigated via combined radial/angular distribution functions (CDFs) analyses and halogen-bonding theory. The most probable orientation of iodobenzene derivatives at the liquid/vacuum interface was assigned by atom density profile and bivariate orientational distribution maps. For all studied iodobenzene derivatives, benzene rings are oriented such that the iodine atoms tend toward the vacuum phase.
Collapse
Affiliation(s)
- Elham Sookhaki
- Department of Chemistry, Yazd University, Yazd, Iran; Department of Chemistry, Shiraz University, Shiraz, 71946-84795, Iran
| | | | | |
Collapse
|
5
|
Honti B, Fábián B, Idrissi A, Jedlovszky P. Surface Properties of N,N-Dimethylformamide-Water Mixtures, As Seen from Computer Simulations. J Phys Chem B 2023; 127:1050-1062. [PMID: 36652674 DOI: 10.1021/acs.jpcb.2c07572] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
The liquid-vapor interface of N,N-dimethylformamide (DMF)-water mixtures, spanning the entire composition range, is investigated in detail at 298 K by molecular dynamics simulation and intrinsic surface analysis. DMF molecules are found to adsorb strongly at the liquid surface, but this adsorption extends only to the first molecular layer. Water and DMF molecules mix with each other on the molecular scale even in the surface layer; thus, no marked self-association of any of the components is seen at the liquid surface. The major surface component prefers such orientation in which the molecular dipole vector lays parallel with the macroscopic plane of the surface. On the other hand, the preferred orientation of the minor component is determined, at both ends of the composition range, by the possibility of H-bond formation with the major component. The lack of H-donating ability of DMF leads to a rapid breakup of the percolating H-bond network at the surface; due to the strong adsorption of DMF, this breakup occurs below the bulk phase DMF mole fraction of 0.03. The disruption of the surface H-bond network also accelerates the exchange of both species between the liquid surface and bulk liquid phase, although, for water, this effect becomes apparent only above a bulk phase DMF mole fraction of 0.4. H-bonds formed by a DMF and a water molecule live, on average, 25-60% longer than those formed by two water molecules at the liquid surface. A similar, but smaller (i.e., about 10-20%) difference is seen in the bulk liquid phase. The enhanced surface mobility of the molecules results in 2-6 times larger diffusion coefficient and 2-5 times shorter H-bond lifetime values at the liquid surface than in the bulk liquid phase. The diffusion of both molecules is slowed down in the presence of the other species; in the case of DMF, this effect is caused by the formation of water-DMF H-bonds, whereas for water, steric hindrances imposed by the bulky DMF neighbors are responsible for this slowing down.
Collapse
Affiliation(s)
- Barbara Honti
- Department of Organic Chemistry and Technology, Budapest University of Technology and Economics, Budapest 1111, Hungary
| | - Balázs Fábián
- Department of Theoretical Biophysics, Max Planck Institute of Biophysics, Frankfurt am Main 60438, Germany
| | - Abdenacer Idrissi
- University of Lille, CNRS UMR 8516 -LASIRe - Laboratoire Avancé de Spectroscopie pour les Interactions la Réactivité et l'environnement, 59000 Lille, France
| | - Pál Jedlovszky
- Department of Chemistry, Eszterházy Károly University, Leányka u. 6, H-3300 Eger, Hungary
| |
Collapse
|
6
|
González-Barramuño B, Cea-Klapp E, Piñeiro MM, Polishuk I, Quinteros-Lama H, Garrido JM. Surface anomalies in ethanol plus n-octane mixture: An effect of molecular orientations and hydrogen bonds. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119630] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
7
|
Balbisi M, Horváth RA, Szőri M, Jedlovszky P. Computer simulation investigation of the adsorption of acetamide on low density amorphous ice. An astrochemical perspective. J Chem Phys 2022; 156:184703. [PMID: 35568547 DOI: 10.1063/5.0093561] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The adsorption of acetamide on low density amorphous (LDA) ice is investigated by grand canonical Monte Carlo computer simulations at the temperatures 50, 100, and 200 K, characteristic of certain domains of the interstellar medium (ISM). We found that the relative importance of the acetamide-acetamide H-bonds with respect to the acetamide-water ones increases with decreasing temperature. Thus, with decreasing temperature, the existence of the stable monolayer, characterizing the adsorption at 200 K, is gradually replaced by the occurrence of marked multilayer adsorption, preceding even the saturation of the first layer at 50 K. While isolated acetamide molecules prefer to lay parallel to the ice surface to maximize their H-bonding with the surface water molecules, this orientational preference undergoes a marked change upon saturation of the first layer due to increasing competition of the adsorbed molecules for H-bonds with water and to the possibility of their H-bond formation with each other. As a result, molecules stay preferentially perpendicular to the ice surface in the saturated monolayer. The chemical potential value corresponding to the point of condensation is found to decrease linearly with increasing temperature. We provide, in analogy with the Clausius-Clapeyron equation, a thermodynamic explanation of this behavior and estimate the molar entropy of condensed phase acetamide to be 34.0 J/mol K. For the surface concentration of the saturated monolayer, we obtain the value 9.1 ± 0.8 µmol/m2, while the heat of adsorption at infinitely low surface coverage is estimated to be -67.8 ± 3.0 kJ/mol. Our results indicate that the interstellar formation of peptide chains through acetamide molecules, occurring at the surface of LDA ice, might well be a plausible process in the cold (i.e., below 50 K) domains of the ISM; however, it is a rather unlikely scenario in its higher temperature (i.e., 100-200 K) domains.
Collapse
Affiliation(s)
- Mirjam Balbisi
- Institute of Chemistry, Eötvös Loránd University, Pázmány P. stny. 1/A, H-1117 Budapest, Hungary
| | - Réka A Horváth
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budafoki út 8, H-1111 Budapest, Hungary
| | - Milán Szőri
- Institute of Chemistry, University of Miskolc, Egyetemváros A/2, H-3515 Miskolc, Hungary
| | - Pál Jedlovszky
- Department of Chemistry, Eszterházy Károly University, Leányka u. 6, H-3300 Eger, Hungary
| |
Collapse
|
8
|
Balbisi M, Horváth RA, Szőri M, Jedlovszky P. Adsorption of acetamide on crystalline and amorphous ice under atmospheric conditions. A grand canonical Monte Carlo simulation study. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.118870] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
9
|
Honti B, Szőri M, Jedlovszky P. Description of the Interfacial Behavior of Benzonitrile at Icy Surfaces by Grand Canonical Monte Carlo Simulations. J Phys Chem A 2022; 126:1221-1232. [PMID: 35168326 DOI: 10.1021/acs.jpca.1c10749] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The adsorption of benzonitrile at the surface of crystalline (Ih) and low-density amorphous (LDA) ice has been investigated by grand canonical Monte Carlo simulations at temperatures ranging from 50 to 200 K. It is found that, in spite of its rather large dipole moment of 4.5 D, benzonitrile molecules can only form a highly unsaturated monolayer on LDA ice, reaching not more than 50% of the surface concentration of the saturated monolayer even at the lowest temperature considered, and they practically do not adsorb on Ih ice. In spite of the observed weak ability of the benzonitrile molecules for being adsorbed, the estimated heat of adsorption at an infinitely low surface concentration of -66.8 ± 2.2 kJ/mol is rather large. This value includes the contribution of roughly -30 to -35 kJ/mol of a benzene ring, about -10 kJ/mol of a large molecular dipole moment, and about -20 to -25 kJ/mol of a benzonitrile-water H-bond, as estimated from comparisons with the heat of adsorption values of similar molecules. The surprisingly weak ability of benzonitrile for adsorption is thus attributed to the unusually strong cohesion between the molecules, considerably exceeding their adhesion to ice, as reflected in the 70-80 kJ/mol difference of the lateral and ice contributions to the binding energy of surface benzonitrile molecules in the presence of condensed benzonitrile.
Collapse
Affiliation(s)
- Barbara Honti
- Budapest University of Technology and Economics, Szt. Gellért tér 4, H-1111 Budapest, Hungary
| | - Milán Szőri
- Institute of Chemistry, University of Miskolc, Egyetemváros A/2, H-3515 Miskolc, Hungary
| | - Pál Jedlovszky
- Department of Chemistry, Eszterházy Károly University, Leányka u. 6, H-3300 Eger, Hungary
| |
Collapse
|
10
|
|
11
|
Voloshin VP, Medvedev NN. ORIENTATION OF WATER MOLECULES NEAR A GLOBULAR PROTEIN. J STRUCT CHEM+ 2021. [DOI: 10.1134/s002247662105005x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
12
|
A molecular simulation study on solvation free energy and structural properties of polycyclic aromatic hydrocarbons in supercritical water environment. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.114274] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
13
|
Wohlfahrt O, Dellago C, Sega M. Ab initio structure and thermodynamics of the RPBE-D3 water/vapor interface by neural-network molecular dynamics. J Chem Phys 2020; 153:144710. [DOI: 10.1063/5.0021852] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Oliver Wohlfahrt
- University of Vienna, Faculty of Physics, Boltzmanngasse 5, A-1090 Vienna, Austria
| | - Christoph Dellago
- University of Vienna, Faculty of Physics, Boltzmanngasse 5, A-1090 Vienna, Austria
| | - Marcello Sega
- University of Vienna, Faculty of Physics, Boltzmanngasse 5, A-1090 Vienna, Austria
- Forschungszentrum Jülich GmbH, Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Fürther Str. 248, D-90429 Nürnberg, Germany
| |
Collapse
|
14
|
Jin H, Ding W, Bai B, Cao C. Molecular dynamics simulation study used in systems with supercritical water. REV CHEM ENG 2020. [DOI: 10.1515/revce-2019-0068] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Abstract
Abstract
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.
Collapse
Affiliation(s)
- Hui Jin
- State Key Laboratory of Multiphase Flow in Power Engineering , Xi'an JiaoTong University , 710049, Shaanxi , China
| | - Weijing Ding
- State Key Laboratory of Multiphase Flow in Power Engineering , Xi'an JiaoTong University , 710049, Shaanxi , China
| | - Bin Bai
- State Key Laboratory of Multiphase Flow in Power Engineering , Xi'an JiaoTong University , 710049, Shaanxi , China
| | - Changqing Cao
- State Key Laboratory of Multiphase Flow in Power Engineering , Xi'an JiaoTong University , 710049, Shaanxi , China
| |
Collapse
|
15
|
A molecular dynamics simulation investigation on the solubility of polycyclic aromatic hydrocarbons in supercritical water. J Mol Liq 2020. [DOI: 10.1016/j.molliq.2020.112464] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
16
|
Heydari Dokoohaki M, Zolghadr AR, Klein A. Impact of the chemical structure on the distribution of neuroprotective N-alkyl-9H-carbazoles at octanol/water interfaces. NEW J CHEM 2020. [DOI: 10.1039/c9nj04251b] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
All-atom molecular dynamics (MD) simulations were performed on distribution and agglomeration dynamics of neuroprotective N-(3-anilinopropyl)-9H-carbazoles at octanol/water interfaces.
Collapse
Affiliation(s)
| | | | - Axel Klein
- Department of Chemistry
- Shiraz University
- Shiraz
- Iran
- Department für Chemie
| |
Collapse
|
17
|
Horváth RA, Fábián B, Szőri M, Jedlovszky P. Investigation of the liquid-vapour interface of aqueous methylamine solutions by computer simulation methods. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.110978] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
18
|
Kiss B, Picaud S, Szőri M, Jedlovszky P. Adsorption of Formamide at the Surface of Amorphous and Crystalline Ices under Interstellar and Tropospheric Conditions. A Grand Canonical Monte Carlo Simulation Study. J Phys Chem A 2019; 123:2935-2948. [PMID: 30839213 DOI: 10.1021/acs.jpca.9b00850] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The adsorption of formamide is studied both at the surface of crystalline (Ih) ice at 200 K and at the surface of low density amorphous (LDA) ice in the temperature range of 50-200 K by grand canonical Monte Carlo (GCMC) simulation. These systems are characteristic of the upper troposphere and of the interstellar medium (ISM), respectively. Our results reveal that while no considerable amount of formamide is dissolved in the bulk ice phase in any case, the adsorption of formamide at the ice surface under these conditions is a very strongly preferred process, which has to be taken into account when studying the chemical reactivity in these environments. The adsorption is found to lead to the formation of multimolecular adsorption layer, the occurrence of which somewhat precedes the saturation of the first molecular layer. Due to the strong lateral interaction acting between the adsorbed formamide molecules, the adsorption isotherm does not follow the Langmuir shape. Adsorption is found to be slightly stronger on LDA than Ih ice under identical thermodynamic conditions, due to the larger surface area exposed to the adsorption. Indeed, the monomolecular adsorption capacity of the LDA and Ih ice surfaces is found to be 10.5 ± 0.7 μmol/m2 and 9.4 μmol/m2, respectively. The first layer formamide molecules are very strongly bound to the ice surface, forming typically four hydrogen bonds with each other and the surface water molecules. The heat of adsorption at infinitely low surface coverage is found to be -105.6 kJ/mol on Ih ice at 200 K.
Collapse
Affiliation(s)
- Bálint Kiss
- Institute of Chemistry , University of Miskolc , Egyetemváros A/2 , H-3515 Miskolc , Hungary.,University of Lille, Faculty of Sciences and Technologies, LASIR (UMR CNRS 8516), 59655 Villeneuve d'Ascq , France
| | - Sylvain Picaud
- Institut UTINAM (CNRS UMR 6213), Université Bourgogne Franche-Comté, 16 Route de Gray , F-25030 Besançon , France
| | - Milán Szőri
- Institute of Chemistry , University of Miskolc , Egyetemváros A/2 , H-3515 Miskolc , Hungary
| | - Pál Jedlovszky
- Department of Chemistry , Eszterházy Károly University , Leányka u. 6 , H-3300 Eger , Hungary
| |
Collapse
|
19
|
Horváth RA, Hantal G, Picaud S, Szőri M, Jedlovszky P. Adsorption of Methylamine on Amorphous Ice under Interstellar Conditions. A Grand Canonical Monte Carlo Simulation Study. J Phys Chem A 2018. [PMID: 29537265 DOI: 10.1021/acs.jpca.8b01591] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The adsorption of methylamine at the surface of amorphous ice is studied at various temperatures, ranging from 20 to 200 K, by grand canonical Monte Carlo simulations under conditions that are characteristic to the interstellar medium (ISM). The results are also compared with those obtained earlier on crystalline ( Ih) ice. We found that methylamine has a strong ability of being adsorbed on amorphous ice, involving also multilayer adsorption. The decrease of the temperature leads to a substantial increase of this adsorption ability; thus, considerable adsorption is seen at 20-50 K even at bulk gas phase concentrations that are comparable with that of the ISM. Further, methylamine molecules can also be dissolved in the bulk amorphous ice phase. Both the adsorption capacity of amorphous ice and the strength of the adsorption on it are found to be clearly larger than those corresponding to crystalline ( Ih) ice, due to the molecular scale roughness of the amorphous ice surface as well as to the lack of clear orientational preferences of the water molecules at this surface. Thus, the surface density of the saturated adsorption monolayer is estimated to be 12.6 ± 0.4 μmol/m2, 20% larger than the value of 10.35 μmol/m2, obtained earlier for Ih ice, and at low enough surface coverages the adsorbed methylamine molecules are found to easily form up to three hydrogen bonds with the surface water molecules. The estimated heat of adsorption at infinitely low surface coverage is calculated to be -69 ± 5 kJ/mol, being rather close to the estimated heat of solvation in the bulk amorphous ice phase of -74 ± 7 kJ/mol, indicating that there are at least a few positions at the surface where the adsorbed methylamine molecules experience a bulk-like local environment.
Collapse
Affiliation(s)
- Réka A Horváth
- Apáczai Csere János School of the ELTE University , Papnövelde u. 4 , H-1053 Budapest , Hungary
| | - György Hantal
- Faculty of Physics , University of Vienna , Boltzmanngasse 5, A-1090 Vienna , Austria
| | - Sylvain Picaud
- Institut UTINAM (CNRS UMR 6213), Université Bourgogne Franche-Comté , 16 Route de Gray , F-25030 Besançon , France
| | - Milán Szőri
- Institute of Chemistry , University of Miskolc , Egyetemváros A/2 , H-3515 Miskolc , Hungary
| | - Pál Jedlovszky
- Department of Chemistry , Eszterházy Károly University , Leányka u. 6 , H-3300 Eger , Hungary
| |
Collapse
|
20
|
Dependence of the adsorption of halogenated methane derivatives at the ice surface on their chemical structure. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.05.110] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
21
|
Zavadlav J, Bevc S, Praprotnik M. Adaptive resolution simulations of biomolecular systems. EUROPEAN BIOPHYSICS JOURNAL: EBJ 2017; 46:821-835. [PMID: 28905203 DOI: 10.1007/s00249-017-1248-0] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2017] [Revised: 06/12/2017] [Accepted: 08/15/2017] [Indexed: 10/18/2022]
Abstract
In this review article, we discuss and analyze some recently developed hybrid atomistic-mesoscopic solvent models for multiscale biomolecular simulations. We focus on the biomolecular applications of the adaptive resolution scheme (AdResS), which allows solvent molecules to change their resolution back and forth between atomistic and coarse-grained representations according to their positions in the system. First, we discuss coupling of atomistic and coarse-grained models of salt solution using a 1-to-1 molecular mapping-i.e., one coarse-grained bead represents one water molecule-for development of a multiscale salt solution model. In order to make use of coarse-grained molecular models that are compatible with the MARTINI force field, one has to resort to a supramolecular mapping, in particular to a 4-to-1 mapping, where four water molecules are represented with one coarse-grained bead. To this end, bundled atomistic water models are employed, i.e., the relative movement of water molecules that are mapped to the same coarse-grained bead is restricted by employing harmonic springs. Supramolecular coupling has recently also been extended to polarizable coarse-grained water models with explicit charges. Since these coarse-grained models consist of several interaction sites, orientational degrees of freedom of the atomistic and coarse-grained representations are coupled via a harmonic energy penalty term. The latter aligns the dipole moments of both representations. The reviewed multiscale solvent models are ready to be used in biomolecular simulations, as illustrated in a few examples.
Collapse
Affiliation(s)
- Julija Zavadlav
- Department of Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, 1001, Ljubljana, Slovenia.,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000, Ljubljana, Slovenia.,Chair of Computational Science, ETH Zurich, Clausiusstrasse 33, 8092, Zurich, Switzerland
| | - Staš Bevc
- Department of Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, 1001, Ljubljana, Slovenia
| | - Matej Praprotnik
- Department of Molecular Modeling, National Institute of Chemistry, Hajdrihova 19, 1001, Ljubljana, Slovenia. .,Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Jadranska 19, 1000, Ljubljana, Slovenia.
| |
Collapse
|
22
|
Sedghamiz T, Bahrami M, Ghatee MH. Enantiospecific adsorption of propranolol enantiomers on naturally chiral copper surface: A molecular dynamics simulation investigation. Chem Phys 2017. [DOI: 10.1016/j.chemphys.2017.02.007] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
23
|
Zolghadr AR, Boroomand S. Spontaneous assembly of HSP90 inhibitors at water/octanol interface: A molecular dynamics simulation study. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2016.12.020] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
|
24
|
Zolghadr AR, Heydari Dokoohaki M. Self-assembly of neuroprotective carbazolium based small molecules at octane/water interface: A simulation investigation. Chem Phys 2016. [DOI: 10.1016/j.chemphys.2016.10.008] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
25
|
Debbarma R, Malani A. Comparative Study of Water Adsorption on a H(+) and K(+) Ion Exposed Mica Surface: Monte Carlo Simulation Study. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2016; 32:1034-1046. [PMID: 26727708 DOI: 10.1021/acs.langmuir.5b04131] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
Clay minerals are used in variety of applications ranging from composites to electronic devices. For their efficient use in such areas, understanding the effect of surface-active agents on interfacial properties is essential. We investigated the role of surface ions in the adsorption of water molecules by using a muscovite mica surface populated with two different, H(+) and K(+), surface ions. A series of grand canonical Monte Carlo (GCMC) simulations at various relative vapor pressures (p/p0) were performed to obtain the water structure and adsorption isotherm on the H(+)-exposed mica (H-mica) surface. The obtained results were compared to the recent simulation data of water adsorption on the K(+)-exposed mica (K-mica) surface reported by Malani and Ayyappa (Malani, A.; Ayappa, K. G. J. Phys. Chem. B 2009, 113, 1058-1067). Water molecules formed two prominent layers adjacent to the H-mica surface, whereas molecular layering was observed adjacent to the K-mica surface. The adsorption isotherm of water on the K-mica surface was characterized by three stages that corresponded to rapid adsorption in the initial regime below p/p0 = 0.1, followed by a linear development regime for p/p0 = 0.1-0.7 and rapid film thickening for p/p0 ≥ 0.7, whereas only latter two regimes were observed in the H-mica system. In addition, the film thickness of adsorbed water molecules for p/p0 < 0.7 was lower as compared to that for the K-mica surface and comparable beyond. The film thickness obtained from the MC simulations was in excellent agreement with the interferometry experimental data of Balmer et al. (Balmer, T. E.; Christenson, H. K.; Spencer, N. D.; Heuberger, M. Langmuir 2008, 24, 1566-1569). It was observed that the hydration behaviors of the two ions were completely different and depended on the size of their hydration shell and their ability to form hydrogen bonds. The behavior of water adsorption between these two cases was illustrated using the water density distribution, orientational distributions, hydrogen bonding, and isosteric heats of adsorption.
Collapse
Affiliation(s)
- Rousan Debbarma
- Department of Chemical Engineering, Indian Institute of Technology Bombay , Mumbai, India
| | - Ateeque Malani
- Department of Chemical Engineering, Indian Institute of Technology Bombay , Mumbai, India
| |
Collapse
|
26
|
McCutchen M, Chen LG, Bermudez H, Matysiak S. Interplay of Dynamical Properties between Ionic Liquids and Ionic Surfactants: Mechanism and Aggregation. J Phys Chem B 2015; 119:9925-32. [DOI: 10.1021/acs.jpcb.5b05151] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Michael McCutchen
- Fischell
Department of Bioengineering, University of Maryland, College Park, Maryland 20742, United States
| | - Lang G. Chen
- Department
of Polymer Science and Engineering, University of Massachusetts, Amherst, Massachusetts 01003, United States
| | - Harry Bermudez
- Department
of Polymer Science and Engineering, University of Massachusetts, Amherst, Massachusetts 01003, United States
| | - Silvina Matysiak
- Fischell
Department of Bioengineering, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
27
|
|
28
|
Idrissi A, Hantal G, Jedlovszky P. Properties of the liquid–vapor interface of acetone–methanol mixtures, as seen from computer simulation and ITIM surface analysis. Phys Chem Chem Phys 2015; 17:8913-26. [DOI: 10.1039/c4cp05974c] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The intrinsic surface of acetone–methanol mixtures is studied by computer simulation and ITIM analysis.
Collapse
Affiliation(s)
- Abdenacer Idrissi
- Laboratoire de Spectrochimie Infrarouge et Raman (UMR CNRS A8516)
- Université Lille 1
- Science et Technologies
- 59655 Villeneuve d’Ascq Cedex
- France
| | - György Hantal
- EKF Department of Chemistry
- H-3300 Eger
- Hungary
- Institut für Computergestützte Biologische Chemie
- University of Vienna
| | - Pál Jedlovszky
- EKF Department of Chemistry
- H-3300 Eger
- Hungary
- Laboratory of Interfaces and Nanosize Systems
- Institute of Chemistry
| |
Collapse
|
29
|
Jedlovszky P, Jójárt B, Horvai G. Properties of the intrinsic surface of liquid acetone, as seen from computer simulations. Mol Phys 2014. [DOI: 10.1080/00268976.2014.968227] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
30
|
Zavadlav J, Melo MN, Cunha AV, de Vries AH, Marrink SJ, Praprotnik M. Adaptive Resolution Simulation of MARTINI Solvents. J Chem Theory Comput 2014; 10:2591-8. [DOI: 10.1021/ct5001523] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Julija Zavadlav
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Hajdrihova
19, SI-1001 Ljubljana, Slovenia
| | - Manuel N. Melo
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
| | - Ana V. Cunha
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
| | - Alex H. de Vries
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
| | - Siewert J. Marrink
- Groningen
Biomolecular Sciences and Biotechnology Institute and Zernike Institute
for Advanced Materials, University of Groningen, Nijenborgh 7, 9747 AG Groningen, Netherlands
| | - Matej Praprotnik
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Hajdrihova
19, SI-1001 Ljubljana, Slovenia
| |
Collapse
|
31
|
Adsorption of H2O2 at the surface of Ih ice, as seen from grand canonical Monte Carlo simulations. Chem Phys Lett 2014. [DOI: 10.1016/j.cplett.2014.03.050] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
32
|
Reprint of “Role of the fluidity of a liquid phase in determining the surface properties of the opposite phase”. J Mol Liq 2014. [DOI: 10.1016/j.molliq.2013.06.009] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
33
|
Nagarajan A, Junghans C, Matysiak S. Multiscale Simulation of Liquid Water Using a Four-to-One Mapping for Coarse-Graining. J Chem Theory Comput 2013; 9:5168-75. [DOI: 10.1021/ct400566j] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Anu Nagarajan
- Fischell
Department of Bioengineering, University of Maryland, College Park, Maryland 20742, United States
| | - Christoph Junghans
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Silvina Matysiak
- Fischell
Department of Bioengineering, University of Maryland, College Park, Maryland 20742, United States
| |
Collapse
|
34
|
Role of the fluidity of a liquid phase in determining the surface properties of the opposite phase at the liquid–liquid interface. J Mol Liq 2013. [DOI: 10.1016/j.molliq.2013.04.012] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
35
|
Sega M, Kantorovich SS, Jedlovszky P, Jorge M. The generalized identification of truly interfacial molecules (ITIM) algorithm for nonplanar interfaces. J Chem Phys 2013; 138:044110. [DOI: 10.1063/1.4776196] [Citation(s) in RCA: 62] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
36
|
Ghatee MH, Zolghadr AR, Moosavi F, Ansari Y. Studies of structural, dynamical, and interfacial properties of 1-alkyl-3-methylimidazolium iodide ionic liquids by molecular dynamics simulation. J Chem Phys 2012; 136:124706. [DOI: 10.1063/1.3696004] [Citation(s) in RCA: 55] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023] Open
|
37
|
Kiss P, Darvas M, Baranyai A, Jedlovszky P. Surface properties of the polarizable Baranyai-Kiss water model. J Chem Phys 2012; 136:114706. [DOI: 10.1063/1.3692602] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023] Open
Affiliation(s)
- Péter Kiss
- Institute of Chemistry, Eötvös Loránd University, Pázmány Péter stny. 1/a, H-1117 Budapest, Hungary
| | | | | | | |
Collapse
|
38
|
Darvas M, Lasne J, Laffon C, Parent P, Picaud S, Jedlovszky P. Adsorption of acetaldehyde on ice as seen from computer simulation and infrared spectroscopy measurements. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2012; 28:4198-4207. [PMID: 22320190 DOI: 10.1021/la204472k] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
Detailed investigation of the adsorption of acetaldehyde on I(h) ice is performed under tropospheric conditions by means of grand canonical Monte Carlo computer simulations and compared to infrared spectroscopy measurements. The experimental and simulation results are in a clear accordance with each other. The simulations indicate that the adsorption process follows Langmuir behavior in the entire pressure range of the vapor phase of acetaldehyde. Further, it was found that the adsorption layer is strictly monomolecular, and the adsorbed acetaldehyde molecules are bound to the ice surface by only one hydrogen bond, typically formed with the dangling H atoms at the ice surface, in agreement with the experimental results. Besides this hydrogen bonding, at high surface coverages dipolar attraction between neighboring acetaldehyde molecules also contributes considerably to the energy gain of the adsorption. The acetaldehyde molecules adopt strongly tilted orientations relative to the ice surface, the tilt angle being scattered between 50° and 90° (i.e., perpendicular orientation). The range of the preferred tilt angles narrows, and the preference for perpendicular orientation becomes stronger upon saturation of the adsorption layer. The CH(3) group of the acetaldehyde molecules points as straight away from the ice surface within the constraint imposed by the tilt angle adopted by the molecule as possible. The heat of adsorption at infinitely low coverage is found to be -36 ± 2 kJ/mol from the infrared spectroscopy measurement, which is in excellent agreement with the computer simulation value of -34.1 kJ/mol.
Collapse
Affiliation(s)
- Mária Darvas
- Institut UTINAM-UMR CNRS 6213, Faculté des Sciences, Université de Franche-Comté, Besançon, France
| | | | | | | | | | | |
Collapse
|
39
|
Rodríguez-Ropero F, Fioroni M. Effect of Na+, Mg2+, and Zn2+ chlorides on the structural and thermodynamic properties of water/n-heptane interfaces. J Comput Chem 2011; 32:1876-86. [DOI: 10.1002/jcc.21770] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2010] [Revised: 01/18/2011] [Accepted: 01/18/2011] [Indexed: 01/08/2023]
|
40
|
Ghatee MH, Zolghadr AR, Moosavi F, Pakdel L. The extent of molecular orientation at liquid/vapor interface of pyridine and its alkyl derivatives by molecular dynamics simulation. J Chem Phys 2011; 134:074707. [DOI: 10.1063/1.3554361] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
McFearin CL, Richmond GL. The unique molecular behavior of water at the chloroform-water interface. APPLIED SPECTROSCOPY 2010; 64:986-994. [PMID: 20828435 DOI: 10.1366/000370210792434288] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
The molecular bonding and orientation of water at the chloroform-water interface has been examined in this study using vibrational sum-frequency spectroscopy (VSFS). The results provide a key puzzle piece towards our understanding of the systematic changes in the interfacial bonding and orientation of water that occur with variations in the polarity of the organic phase, especially when compared with previous studies of different liquid-liquid interfacial systems. In these VSFS studies the OH spectral responses of interfacial water molecules are used to characterize the interactions between water and the organic phase. The spectral analysis, aided by isotopic dilution studies, shows that the moderate polarity of the chloroform phase results in a mixed interfacial region with stronger organic-water bonding and fewer bonding interactions between adjacent water molecules than was previously found for studies of non-polar organic liquid-water interfaces. Even with the more mixed interfacial region and stronger organic-water interactions, interfacial water retains a significant amount of orientational ordering. These results are compared with recent predictions from molecular dynamics simulations about how molecules behave at the chloroform-water interface.
Collapse
Affiliation(s)
- Cathryn L McFearin
- Department of Chemistry, University of Oregon, 1253 University of Oregon, Eugene, Oregon 97403, USA
| | | |
Collapse
|
42
|
Hantal G, Darvas M, Pártay LB, Horvai G, Jedlovszky P. Molecular level properties of the free water surface and different organic liquid/water interfaces, as seen from ITIM analysis of computer simulation results. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2010; 22:284112. [PMID: 21399284 DOI: 10.1088/0953-8984/22/28/284112] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
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.
Collapse
Affiliation(s)
- György Hantal
- Laboratory of Interfaces and Nanosize Systems, Institute of Chemistry, Eötvös Loránd University, Pázmány P Stny 1/A, H-1117 Budapest, Hungary
| | | | | | | | | |
Collapse
|
43
|
Petitjean M, Hantal G, Chauvin C, Mirabel P, Le Calvé S, Hoang PNM, Picaud S, Jedlovszky P. Adsorption of benzaldehyde at the surface of ice, studied by experimental method and computer simulation. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2010; 26:9596-9606. [PMID: 20329716 DOI: 10.1021/la100169h] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
Adsorption study of benzaldehyde on ice surfaces is performed by combining experimental and theoretical approaches. The experiments are conducted over the temperature range 233-253 K using a coated wall flow tube coupled to a mass spectrometric detector. Besides the experimental way, the adsorption isotherm is also determined by performing a set of grand canonical Monte Carlo simulations at 233 K. The experimental and calculated adsorption isotherms show a very good agreement within the corresponding errors. Besides, both experimental and theoretical studies permit us to derive the enthalpy of adsorption of benzaldehyde on ice surfaces DeltaH(ads), which are in excellent agreement: DeltaH(ads) = -61.4 +/- 9.7 kJ/mol (experimental) and DeltaH(ads) = -59.4 +/- 5.1 kJ/mol (simulation). The obtained results indicate a much stronger ability of benzaldehyde of being adsorbed at the surface of ice than that of small aliphatic aldehydes, such as formaldehyde or acetaldehyde. At low surface coverages the adsorbed molecules exclusively lie parallel with the ice surface. With increasing surface coverage, however, the increasing competition of the adsorbed molecules for the surface area to be occupied leads to the appearance of two different perpendicular orientations relative to the surface. In the first orientation, the benzaldehyde molecule turns its aldehyde group toward the ice phase, and, similarly to the molecules in the lying orientation, forms a hydrogen bond with a surface water molecule. In the other perpendicular orientation the aldehyde group turns to the vapor phase, and its O atom interacts with the delocalized pi system of the benzene ring of a nearby lying benzaldehyde molecule of the second molecular layer. In accordance with this observed scenario, the saturated adsorption layer, being stable in a roughly 1 kJ/mol broad range of chemical potentials, contains, besides the first molecular layer, also traces of the second molecular layer of adsorbed benzaldehyde.
Collapse
Affiliation(s)
- Mélanie Petitjean
- Laboratoire des Matériaux, Surfaces et Procédés pour la Catalyse (LMSPC, UMR 7515 CNRS/UdS), 25 rue Becquerel, 67087 Strasbourg Cedex 02, France
| | | | | | | | | | | | | | | |
Collapse
|
44
|
Darvas M, Pojják K, Horvai G, Jedlovszky P. Molecular dynamics simulation and identification of the truly interfacial molecules (ITIM) analysis of the liquid-vapor interface of dimethyl sulfoxide. J Chem Phys 2010; 132:134701. [DOI: 10.1063/1.3368111] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023] Open
|
45
|
Darvas M, Pártay LB, Jedlovszky P, Horvai G. Computer simulation and ITIM analysis of the surface of water–methanol mixtures containing traces of water. J Mol Liq 2010. [DOI: 10.1016/j.molliq.2009.06.004] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
46
|
Abstract
We describe and illustrate a simple procedure for identifying a liquid interface from atomic coordinates. In particular, a coarse-grained density field is constructed, and the interface is defined as a constant density surface for this coarse-grained field. In applications to a molecular dynamics simulation of liquid water, it is shown that this procedure provides instructive and useful pictures of liquid-vapor interfaces and of liquid-protein interfaces.
Collapse
Affiliation(s)
- Adam P Willard
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | |
Collapse
|
47
|
Szőri M, Jedlovszky P, Roeselová M. Water adsorption on hydrophilic and hydrophobic self-assembled monolayers as proxies for atmospheric surfaces. A grand canonical Monte Carlo simulation study. Phys Chem Chem Phys 2010; 12:4604-16. [DOI: 10.1039/b923382b] [Citation(s) in RCA: 41] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
48
|
LIU S, CHAI J, YANG X. β-Diketones at Water/Supercritical CO2 Interface: A Molecular Dynamics Simulation. Chin J Chem Eng 2009. [DOI: 10.1016/s1004-9541(08)60307-9] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
49
|
Pártay LB, Jedlovszky P, Vincze Á, Horvai G. Properties of Free Surface of Water−Methanol Mixtures. Analysis of the Truly Interfacial Molecular Layer in Computer Simulation. J Phys Chem B 2008; 112:5428-38. [DOI: 10.1021/jp711547e] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lívia B. Pártay
- Laboratory of Interfaces and Nanosize Systems, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter stny. 1/a, H-1117 Budapest, Hungary, HAS Research Group of Technical Analytical Chemistry, Szt. Gellért tér 4, H-1111 Budapest, Hungary, Department of NBC and Environmental Security, Zrínyi Miklós National Defense University, Hungária krt. 9-11, H-1581 Budapest, Hungary, and Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Szt. Gellért tér 4, H
| | - Pál Jedlovszky
- Laboratory of Interfaces and Nanosize Systems, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter stny. 1/a, H-1117 Budapest, Hungary, HAS Research Group of Technical Analytical Chemistry, Szt. Gellért tér 4, H-1111 Budapest, Hungary, Department of NBC and Environmental Security, Zrínyi Miklós National Defense University, Hungária krt. 9-11, H-1581 Budapest, Hungary, and Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Szt. Gellért tér 4, H
| | - Árpád Vincze
- Laboratory of Interfaces and Nanosize Systems, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter stny. 1/a, H-1117 Budapest, Hungary, HAS Research Group of Technical Analytical Chemistry, Szt. Gellért tér 4, H-1111 Budapest, Hungary, Department of NBC and Environmental Security, Zrínyi Miklós National Defense University, Hungária krt. 9-11, H-1581 Budapest, Hungary, and Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Szt. Gellért tér 4, H
| | - George Horvai
- Laboratory of Interfaces and Nanosize Systems, Institute of Chemistry, Eötvös Loránd University, Pázmány Péter stny. 1/a, H-1117 Budapest, Hungary, HAS Research Group of Technical Analytical Chemistry, Szt. Gellért tér 4, H-1111 Budapest, Hungary, Department of NBC and Environmental Security, Zrínyi Miklós National Defense University, Hungária krt. 9-11, H-1581 Budapest, Hungary, and Department of Inorganic and Analytical Chemistry, Budapest University of Technology and Economics, Szt. Gellért tér 4, H
| |
Collapse
|
50
|
Chowdhary J, Ladanyi BM. Water/hydrocarbon interfaces: effect of hydrocarbon branching on single-molecule relaxation. J Phys Chem B 2008; 112:6259-73. [PMID: 18324803 DOI: 10.1021/jp0769025] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
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.
Collapse
Affiliation(s)
- Janamejaya Chowdhary
- Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523-1872, USA
| | | |
Collapse
|