1
|
Crundwell FK. The impact of surface charge on the ionic dissociation of common salt (NaCl). Chem Eng Sci 2019. [DOI: 10.1016/j.ces.2019.04.050] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
2
|
Jung H, Yethiraj A. A simulation method for the phase diagram of complex fluid mixtures. J Chem Phys 2018; 148:244903. [PMID: 29960369 DOI: 10.1063/1.5033958] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The phase behavior of complex fluid mixtures is of continuing interest, but obtaining the phase diagram from computer simulations can be challenging. In the Gibbs ensemble method, for example, each of the coexisting phases is simulated in a different cell, and ensuring the equality of chemical potentials of all components requires the transfer of molecules from one cell to the other. For complex fluids such as polymers, successful insertions are rare. An alternative method is to simulate both coexisting phases in a single simulation cell, with an interface between them. The challenge here is that the interface position moves during the simulation, making it difficult to determine the concentration profile and coexisting concentrations. In this work, we propose a new method for single cell simulations that uses a spatial concentration autocorrelation function to (spatially) align instantaneous concentration profiles from different snapshots. This allows one to obtain average concentration profiles and hence the coexisting concentrations. We test the method by calculating the phase diagrams of two systems: the Widom-Rowlinson model and the symmetric blends of freely jointed polymer molecules for which phase diagrams from conventional methods are available. Excellent agreement is found, except in the neighborhood of the critical point where the interface is broad and finite size effects are important. The method is easy to implement and readily applied to any mixture of complex fluids.
Collapse
Affiliation(s)
- Hyuntae Jung
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| | - Arun Yethiraj
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| |
Collapse
|
3
|
Abstract
Desolvation barriers are present for solute-solvent exchange events, such as ligand binding to an enzyme active site, during protein folding, and at battery electrodes. For solution-grown crystals, desolvation at kink sites can be the rate-limiting step for growth. However, desolvation and the associated kinetic barriers are poorly understood. In this work, we use rare-event simulation techniques to investigate attachment/detachment events at kink sites of a NaCl crystal in water. We elucidate the desolvation mechanism and present an optimized reaction coordinate, which involves one solute collective variable and one solvent collective variable. The attachment/detachment pathways for Na+ and Cl- are qualitatively similar, with quantitative differences that we attribute to different ion sizes and solvent coordination. The attachment barriers primarily result from kink site desolvation, while detachment barriers largely result from breaking ion-crystal bonds. We compute ion detachment rates from kink sites and compare with results from an independent study. We anticipate that the reaction coordinate and desolvation mechanism identified in this work may be applicable to other alkali halides.
Collapse
|
4
|
Lovrić J, Duflot D, Monnerville M, Toubin C, Briquez S. Water-Induced Organization of Palmitic Acid at the Surface of a Model Sea Salt Particle: A Molecular Dynamics Study. J Phys Chem A 2016; 120:10141-10149. [DOI: 10.1021/acs.jpca.6b07792] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Josip Lovrić
- Laboratoire de
Physique des Lasers, Atomes et Molécules (PhLAM) CNRS, UMR
8523, Univ. Lille, F-59000 Lille, France
| | - Denis Duflot
- Laboratoire de
Physique des Lasers, Atomes et Molécules (PhLAM) CNRS, UMR
8523, Univ. Lille, F-59000 Lille, France
| | - Maurice Monnerville
- Laboratoire de
Physique des Lasers, Atomes et Molécules (PhLAM) CNRS, UMR
8523, Univ. Lille, F-59000 Lille, France
| | - Céline Toubin
- Laboratoire de
Physique des Lasers, Atomes et Molécules (PhLAM) CNRS, UMR
8523, Univ. Lille, F-59000 Lille, France
| | - Stéphane Briquez
- Laboratoire de
Physique des Lasers, Atomes et Molécules (PhLAM) CNRS, UMR
8523, Univ. Lille, F-59000 Lille, France
| |
Collapse
|
5
|
Akbarzadeh H, Abbaspour M, Salemi S, Masoumi A, Shamkhali AN. A modified thermodynamic insight to deliquescence of a void-containing nanocrystal confirmed by MD simulation. AIChE J 2016. [DOI: 10.1002/aic.15306] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Hamed Akbarzadeh
- Dept. of Chemistry, Faculty of Basic Sciences; Hakim Sabzevari University; 96179-76487 Sabzevar Iran
| | - Mohsen Abbaspour
- Dept. of Chemistry, Faculty of Basic Sciences; Hakim Sabzevari University; 96179-76487 Sabzevar Iran
| | - Sirous Salemi
- Dept. of Chemistry, Faculty of Basic Sciences; Hakim Sabzevari University; 96179-76487 Sabzevar Iran
| | - Azizeh Masoumi
- Dept. of Chemistry, Faculty of Basic Sciences; Hakim Sabzevari University; 96179-76487 Sabzevar Iran
| | - Amir Nasser Shamkhali
- Dept. of Chemistry, Faculty of Basic Science; University of Mohaghegh Ardabili; 56199-11367 Ardabil Iran
| |
Collapse
|
6
|
Affiliation(s)
- Gabriele Lanaro
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1
| | - G. N. Patey
- Department of Chemistry, University of British Columbia, Vancouver, British Columbia, Canada V6T 1Z1
| |
Collapse
|
7
|
Holmberg N, Chen JC, Foster AS, Laasonen K. Dissolution of NaCl nanocrystals: an ab initio molecular dynamics study. Phys Chem Chem Phys 2014; 16:17437-46. [DOI: 10.1039/c4cp00635f] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
NaCl nanocrystal dissolution was investigated in atomistic detail revealing a difference in the solvation of two different ionic species.
Collapse
Affiliation(s)
- Nico Holmberg
- Department of Chemistry
- Aalto University
- FI-00076 Aalto, Finland
| | - Jian-Cheng Chen
- Department of Applied Physics
- Aalto University
- FI-00076 Aalto, Finland
- COMP Centre of Excellence in Computational Nanoscience
- Aalto University
| | - Adam S. Foster
- Department of Applied Physics
- Aalto University
- FI-00076 Aalto, Finland
- COMP Centre of Excellence in Computational Nanoscience
- Aalto University
| | - Kari Laasonen
- Department of Chemistry
- Aalto University
- FI-00076 Aalto, Finland
- COMP Centre of Excellence in Computational Nanoscience
- Aalto University
| |
Collapse
|
8
|
Klimeš J, Bowler DR, Michaelides A. Understanding the role of ions and water molecules in the NaCl dissolution process. J Chem Phys 2013; 139:234702. [DOI: 10.1063/1.4840675] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
|
9
|
|
10
|
Drawing Out the Structural Information About the First Hydration Layer of the Isolated Cl− Anion Through the FTIR-ATR Difference Spectra. J SOLUTION CHEM 2013. [DOI: 10.1007/s10953-013-9970-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
|
11
|
Liu LM, Laio A, Michaelides A. Initial stages of salt crystal dissolution determined with ab initio molecular dynamics. Phys Chem Chem Phys 2011; 13:13162-6. [DOI: 10.1039/c1cp21077g] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
|
12
|
Yagasaki T, Saito S, Ohmine I. Effects of nonadditive interactions on ion solvation at the water/vapor interface: a molecular dynamics study. J Phys Chem A 2010; 114:12573-84. [PMID: 21077653 DOI: 10.1021/jp1084795] [Citation(s) in RCA: 22] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
The solvation of halide ions at the water/vapor interface is investigated by using molecular dynamics simulations with nonpolarizable molecular mechanical (MM), polarizable MM, and quantum mechanical (QM)/MM methods. The free energy profile of the ion solvation is decomposed into the energy and the entropic contributions along the ion displacement from inside to the surface of water. It is found that the surface affinity of the ion, relative to the bulk value, is determined by a subtle balance between the energetic destabilization and the entropic stabilization with the ion displacement. The amount of energetic destabilization is found to be reduced when nonadditive interactions are included, as in the polarizable MM and QM/MM models. The structure of water around the ion at the interface is also largely modified when the higher order effects are considered. For example, the induced dipole effect enhances the solvation structure around the ion at the interface significantly and thus reduces the amount of entropic stabilization at the interface, relative to in the bulk. It is found that this induced dipole effect causes the slowing in the ion-water hydrogen bond dynamics at the interface. On the other hand, the higher order induced multipole effects in the QM/MM method suppress both the excessive enhancement of the solvation structure and the slowing of the ion-water hydrogen bond dynamics at the interface. The present study demonstrates that not only the induced dipole moment but also the higher order induced multipole moments, which are neglected in standard empirical models, are essential for the correct description of the ion solvation at the water/vapor interface.
Collapse
Affiliation(s)
- Takuma Yagasaki
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, and The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | | | | |
Collapse
|
13
|
Smith PE. The effect of urea on the morphology of NaCl crystals: A combined theoretical and simulation study. FLUID PHASE EQUILIBRIA 2010; 290:36-42. [PMID: 20383314 PMCID: PMC2850219 DOI: 10.1016/j.fluid.2009.08.006] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2023]
Abstract
It has been known for over a century that the presence of cosolvents such as urea and formamide can alter the morphology of NaCl crystals grown from solution. To help understand this effect we have been developing a theoretical approach based on the Kirkwood-Buff (KB) theory of solutions, and have combined this with computer simulations of the interation of urea with different crystal faces of NaCl. In this way one can predict the effect of urea on the thermodynamic stability of different NaCl faces, with atomic level detail provided by the simulations. We observe that urea is preferentially excluded from 100 and 111 crystal faces, but is less excluded from 111 faces which present chloride ions at the surface. The results indicate that the 111 face is stabilized in urea solutions and promotes the formation of octahedral over cubic NaCl crystals. The approach is totally general and can be applied to understand a variety of interfacial properties. Furthermore, we apply KB theory to study several other issues regarding the simulation of crystal growth.
Collapse
Affiliation(s)
- Paul E Smith
- Department of Chemistry, 213 CBC Building, Kansas State University, Manhattan, KS 66506-0401, Tel: 785-532-5109, ,
| |
Collapse
|
14
|
Liu LM, Krack M, Michaelides A. Interfacial water: A first principles molecular dynamics study of a nanoscale water film on salt. J Chem Phys 2009; 130:234702. [DOI: 10.1063/1.3152845] [Citation(s) in RCA: 43] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
|
15
|
Lu PD, Wang F, Zhao LJ, Li WX, Li XH, Dong JL, Zhang YH, Lu GQ. Molecular events in deliquescence and efflorescence phase transitions of sodium nitrate particles studied by Fourier transform infrared attenuated total reflection spectroscopy. J Chem Phys 2008; 129:104509. [DOI: 10.1063/1.2973623] [Citation(s) in RCA: 27] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
|
16
|
Zhao LJ, Wang F, Zhang K, Zeng QX, Zhang YH. Deliquescence and Efflorescence Processes of Aerosol Particles Studied byin situFTIR and Raman Spectroscopy. CHINESE J CHEM PHYS 2008. [DOI: 10.1088/1674-0068/21/01/1-11] [Citation(s) in RCA: 15] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
17
|
Bahadur R, Russell LM. Water uptake coefficients and deliquescence of NaCl nanoparticles at atmospheric relative humidities from molecular dynamics simulations. J Chem Phys 2008; 129:094508. [DOI: 10.1063/1.2971040] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
|
18
|
Bahadur R, Russell LM, Alavi S. Surface Tensions in NaCl−Water−Air Systems from MD Simulations. J Phys Chem B 2007; 111:11989-96. [PMID: 17894485 DOI: 10.1021/jp075356c] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Surface tensions for liquid-vapor (lv), solid-liquid (sl), and solid-vapor (sv) interfaces are calculated from molecular dynamics simulations of the NaCl-water-air system. Three distinct calculation techniques based on thermodynamic properties are used to describe the multicomponent mixtures. Simulations of each bulk phase (including a liquid saturated solution) and various interfaces are carried out at both NPT and NVT conditions. The thermodynamic relation for energy difference between interface and bulk phases provides an upper bound to the surface tension, while the energy-integral and test area methods provide direct estimates. At 1 atm and 300 K, the best predictions for surface tensions are sigmasv (NaCl-air) of 114 mN m(-1), sigmasl (NaCl- soln) of 63 mN m(-1), sigmalv (soln-air) of 82 mN m(-1), and sigmalv (water-air) of 66 mN m(-1). The calculated surface tensions from simulations have uncertainties between 5 and 10%, which are higher than measurements for the liquid interfaces and lower than the measurement uncertainty for the solid interfaces. The calculated upper bounds for surface tensions of liquid interfaces compare well with experimental results but provide no improvement over existing measurements. However, the bounding values for solid interfaces lower uncertainty by as much as a factor of 10 as compared to the indirect experimental measurements currently available. The energy-integral and test area methods appear to underestimate the surface tension of water by 10%, which is consistent with previous studies using similar model potentials. The calculated upper bounds of surface tension show a weakly positive correlation with pressure in the 0.1-100 atm range for liquid-solid, liquid-vapor, and solid-vapor interfaces.
Collapse
Affiliation(s)
- Ranjit Bahadur
- Scripps Institution of Oceanography, University of California San Diego, La Jolla, CA 92093-0221, USA
| | | | | |
Collapse
|