1
|
Dufrêche JF, Siboulet B, Duvail M. Chemical models for dense solutions. Faraday Discuss 2024. [PMID: 39041303 DOI: 10.1039/d4fd00084f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/24/2024]
Abstract
Here we examine the question of the chemical models widely used to describe dense solutions, particularly ionic solutions. First, a simple macroscopic analysis shows that, in the case of weak interactions, taking into account aggregated species amounts to modelling an effective attraction between solutes, although the stoichiometry used does not necessarily correspond to atomic reality. We then use a rigorous microscopic analysis to explain how, in the very general case, chemical models can be obtained from an atomic physical description. We show that there are no good or bad chemical models as long as we consider exact calculations. To obtain the simplest possible description, it is nevertheless advisable to take the speciation criterion that minimises the excess terms. Molecular simulations show that, very often, species can be defined simply by grouping ions which are in direct contact. In some cases, the appearance of macroscale clusters can be predicted.
Collapse
Affiliation(s)
- J-F Dufrêche
- ICSM, University of Montpellier, CEA, CNRS, ENSCM, Marcoule, France.
| | - B Siboulet
- ICSM, University of Montpellier, CEA, CNRS, ENSCM, Marcoule, France.
| | - M Duvail
- ICSM, University of Montpellier, CEA, CNRS, ENSCM, Marcoule, France.
| |
Collapse
|
2
|
Coste A, Slejko E, Zavadlav J, Praprotnik M. Developing an Implicit Solvation Machine Learning Model for Molecular Simulations of Ionic Media. J Chem Theory Comput 2024; 20:411-420. [PMID: 38118122 PMCID: PMC10782447 DOI: 10.1021/acs.jctc.3c00984] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2023] [Revised: 12/04/2023] [Accepted: 12/04/2023] [Indexed: 12/22/2023]
Abstract
Molecular dynamics (MD) simulations of biophysical systems require accurate modeling of their native environment, i.e., aqueous ionic solution, as it critically impacts the structure and function of biomolecules. On the other hand, the models should be computationally efficient to enable simulations of large spatiotemporal scales. Here, we present the deep implicit solvation model for sodium chloride solutions that satisfies both requirements. Owing to the use of the neural network potential, the model can capture the many-body potential of mean force, while the implicit water treatment renders the model inexpensive. We demonstrate our approach first for pure ionic solutions with concentrations ranging from physiological to 2 M. We then extend the model to capture the effective ion interactions in the vicinity and far away from a DNA molecule. In both cases, the structural properties are in good agreement with all-atom MD, showcasing a general methodology for the efficient and accurate modeling of ionic media.
Collapse
Affiliation(s)
- Amaury Coste
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Ljubljana SI-1001, Slovenia
| | - Ema Slejko
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Ljubljana SI-1001, Slovenia
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana SI-1000, Slovenia
| | - Julija Zavadlav
- Professorship
of Multiscale Modeling of Fluid Materials, TUM School of Engineering
and Design, Technical University of Munich, Garching Near Munich DE-85748, Germany
| | - Matej Praprotnik
- Laboratory
for Molecular Modeling, National Institute of Chemistry, Ljubljana SI-1001, Slovenia
- Department
of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana SI-1000, Slovenia
| |
Collapse
|
3
|
Gun'ko VM, Turov VV, Krupska TV, Borysenko MV. Surroundings effects on the interfacial and temperature behaviors of NaOH/water bound to hydrophilic and hydrophobic nanosilicas. J Colloid Interface Sci 2023; 634:93-109. [PMID: 36535173 DOI: 10.1016/j.jcis.2022.12.027] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2022] [Revised: 11/30/2022] [Accepted: 12/06/2022] [Indexed: 12/14/2022]
Abstract
HYPOTHESIS It has been assumed that the temperature and interfacial behaviors of concentrated alkali solutions under confined space effects may depend on adsorbent surface structure, hydrophilicity/hydrophobicity, porosity of solids, and dispersion media properties causing kosmotropic or chaotropic effects onto hydrogen bond network (HBN) in bound water and NaOH solution. EXPERIMENTS To analyze these effects, systems with NaOH/water (0.1 g/g/0.1 g/g) deposited onto compacted hydrophilic (A-300) and hydrophobic (AM1) nanosilicas were studied using 1H NMR spectroscopy (215-287 K). The materials were characterized using several experimental and theoretical methods. FINDINGS It has been shown that bound water and water/NaOH represent various clusters and domains whose characteristics depend strongly on nanosilica hydrophilicity/hydrophobicity, dispersion media (air, CDCl3, DMSO, CDCl3/DMSO), subsequent or simultaneous deposition of NaOH and water, and temperature. Water amount (0.1 g/g) was selected too small to completely dissolve NaOH (0.1 g/g) under confined space effects and low temperatures. Chaotropic hydrophobic AM1 and CDCl3 enhance water clusterization and HBN disorder (weakly associated water, WAW appears) in contrast to kosmotropic hydrophilic A-300, NaOH, and DMSO reducing the clusterization and HBN disorder in bound water (WAW disappears). Several aspects related to the interfacial and temperature behaviors of water and co-adsorbates bound to the nanosilicas were elucidated.
Collapse
Affiliation(s)
- Volodymyr M Gun'ko
- Chuiko Institute of Surface Chemistry, 17 General Naumov Street, 03164 Kyiv, Ukraine.
| | - Volodymyr V Turov
- Chuiko Institute of Surface Chemistry, 17 General Naumov Street, 03164 Kyiv, Ukraine
| | - Tetyana V Krupska
- Chuiko Institute of Surface Chemistry, 17 General Naumov Street, 03164 Kyiv, Ukraine
| | - Mykola V Borysenko
- Chuiko Institute of Surface Chemistry, 17 General Naumov Street, 03164 Kyiv, Ukraine
| |
Collapse
|
4
|
Habibi P, Rahbari A, Blazquez S, Vega C, Dey P, Vlugt TJH, Moultos OA. A New Force Field for OH – for Computing Thermodynamic and Transport Properties of H 2 and O 2 in Aqueous NaOH and KOH Solutions. J Phys Chem B 2022; 126:9376-9387. [DOI: 10.1021/acs.jpcb.2c06381] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Parsa Habibi
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
- Department of Materials Science and Engineering, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Mekelweg 2, 2628 CDDelft, The Netherlands
| | - Ahmadreza Rahbari
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
| | - Samuel Blazquez
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040Madrid, Spain
| | - Carlos Vega
- Depto. Química Física, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040Madrid, Spain
| | - Poulumi Dey
- Department of Materials Science and Engineering, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Mekelweg 2, 2628 CDDelft, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
| | - Othonas A. Moultos
- Engineering Thermodynamics, Process & Energy Department, Faculty of Mechanical, Maritime and Materials Engineering, Delft University of Technology, Leeghwaterstraat 39, 2628 CBDelft, The Netherlands
| |
Collapse
|
5
|
Le Crom S, Dourdain S, Pellet-Rostaing S, Duvail M. Long-Range Organization Study of Piperidinium-Based Ionic Liquids by Polarizable Molecular Dynamics Simulations. J Phys Chem B 2022; 126:3355-3365. [PMID: 35471118 DOI: 10.1021/acs.jpcb.2c01123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The nanoscale organization of some classes of ionic liquids is responsible for their singular properties. In this paper, we use polarizable molecular dynamics simulations and small-angle X-ray scattering to probe the structure of two piperidinium- and (trifluoromethylsulfonyl)imide-based ionic liquids ([EBPip+][NTf2-] and [EOPip+][NTf2-]) that differ in the alkyl chain length of their cation. The X-ray scattering intensities calculated numerically, from the radial distribution functions, are in excellent agreement with the experimental data. The analysis of the different contributions of the X-ray scattering data allowed us to highlight the correlations responsible for the low q peak observed for the long-chain alkyl cations. New angular analyses showed that anions were more likely to align with alkyl chains as their size increased, inducing angular correlation between anions at larger distances. They also showed that the long alkyl chains of the cations aligned more with each other than the short ones. These more aligned alkyl chains induce a smaller volume of the apolar microdomains compared to the well-studied imidazolium-based ionic liquids, leading to the smaller correlation distance for piperidinium-based ionic liquids.
Collapse
Affiliation(s)
| | | | | | - Magali Duvail
- ICSM, Univ Montpellier, CEA, CNRS, ENSCM, Bagnols-sur-Cèze, France
| |
Collapse
|
6
|
Chen Y, Fu X, Yu S, Quan K, Zhao C, Shao Z, Ye D, Qi H, Chen P. Parameterization of classical nonpolarizable force field for hydroxide toward the large‐scale molecular dynamics simulation of cellulose in pre‐cooled alkali/urea aqueous solution. J Appl Polym Sci 2021. [DOI: 10.1002/app.51477] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Yu Chen
- Beijing Engineering Research Centre of Cellulose and Its Derivatives School of Materials Science and Engineering, Beijing Institute of Technology Beijing China
| | - Xiaotong Fu
- Key Laboratory of Eco‐Textiles, Ministry of Education Jiangnan University Wuxi Jiangsu Province China
| | - Shuxian Yu
- Beijing Engineering Research Centre of Cellulose and Its Derivatives School of Materials Science and Engineering, Beijing Institute of Technology Beijing China
| | - Kun Quan
- China Institute of Marine Technology and Economy Beijing China
| | - Changjun Zhao
- Beijing Engineering Research Centre of Cellulose and Its Derivatives School of Materials Science and Engineering, Beijing Institute of Technology Beijing China
| | - Ziqiang Shao
- Beijing Engineering Research Centre of Cellulose and Its Derivatives School of Materials Science and Engineering, Beijing Institute of Technology Beijing China
| | - Dongdong Ye
- School of Textile Materials and Engineering Wuyi University Jiangmen Guangdong Province China
| | - Haisong Qi
- State Key Laboratory of Pulp and Paper Engineering South China University of Technology Guangzhou Guangdong Province China
| | - Pan Chen
- Beijing Engineering Research Centre of Cellulose and Its Derivatives School of Materials Science and Engineering, Beijing Institute of Technology Beijing China
| |
Collapse
|
7
|
Bertin M, Gomes Rodrigues D, Pierlot C, Albert-Mercier C, Davy C, Lambertin D, Nardello-Rataj V. Influence of cetyltrimethylammonium bromide and hydroxide ions on the interfacial tension and stability of emulsions of dodecane in aqueous silicate solutions. Colloids Surf A Physicochem Eng Asp 2021. [DOI: 10.1016/j.colsurfa.2021.127306] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
8
|
Moreno Martinez D, Acher E, Vatin M, Dourdain S, Guillaumont D, Guilbaud P. Aggregation of Bifunctional Extractants Used for Uranium(VI) Separation. J Phys Chem B 2021; 125:10759-10771. [PMID: 34320806 DOI: 10.1021/acs.jpcb.1c03529] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
DEHCNPB (butyl-N,N-di(2-ethylhexyl)carbamoyl-nonylphosphonate) is an amido-phosphonic acid that has remarkable properties for the separation of uranium from wet phosphoric acid. Despite previous studies, a detailed description of the DEHCNPB organic solutions at the supramolecular and molecular scales is missing. In the present work, we use classical Molecular Dynamics (MD) combined with SANS and SAXS experimental data in order to describe the aggregation of the bifunctional extractant DEHCNPB as well as the speciation of uranium(VI) in such systems. We provide a fine description of the molecular species in the organic solution and of the interactions within the aggregates formed, shedding light on solvent extraction mechanisms. Without uranium, the organic phase is highly composed of dimers and trimers H-bonded through phosphonate functions and without water molecules. With uranium, two to three extractant molecules coordinate directly the uranyl cation by their phosphonate groups. Uranyl is not fully dehydrated in this organic solution, and the amide groups of the extractants are found to form H-bonds with the water molecules bound to uranyl. These H-bond networks around the metallic cation stabilize the complexes and facilitate the extraction. These results underline the importance of considering weak interactions in the understanding of extraction processes and demonstrate how molecular simulations provide essential insights into such complex organic phase chemistry with a high number of species.
Collapse
Affiliation(s)
| | - Eléonor Acher
- CEA, DES, ISEC, DMRC, Univ Montpellier, Bagnols sur Ceze 30207, France
| | - Marin Vatin
- ICSM, Univ Montpellier, CEA, CNRS, ENSCM, Marcoule, France
| | | | | | - Philippe Guilbaud
- CEA, DES, ISEC, DMRC, Univ Montpellier, Bagnols sur Ceze 30207, France
| |
Collapse
|
9
|
Vatin M, Porro A, Sator N, Dufrêche JF, Berthoumieux H. Electrostatic interactions in water: a nonlocal electrostatic approach. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1825849] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- M. Vatin
- Laboratoire de Physique Théorique de la Matière Condensée (LPTMC), Sorbonne Université, CNRS, Paris, France
- ICSM/LMCT Site de Marcoule, Bagnols sur Céze Cedex, France
| | - A. Porro
- Laboratoire de Physique Théorique de la Matière Condensée (LPTMC), Sorbonne Université, CNRS, Paris, France
| | - N. Sator
- Laboratoire de Physique Théorique de la Matière Condensée (LPTMC), Sorbonne Université, CNRS, Paris, France
| | - J.-F. Dufrêche
- ICSM/LMCT Site de Marcoule, Bagnols sur Céze Cedex, France
| | - H. Berthoumieux
- Laboratoire de Physique Théorique de la Matière Condensée (LPTMC), Sorbonne Université, CNRS, Paris, France
| |
Collapse
|
10
|
Ramamoorthy RK, Levesque M, Belloni L, Carriere D. Structure Factor of EuCl 3 Aqueous Solutions via Coupled Molecular Dynamics Simulations and Integral Equations. J Phys Chem B 2020; 124:1787-1793. [PMID: 32026687 DOI: 10.1021/acs.jpcb.9b11537] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Identifying the structure of an aqueous solution is essential to rationalize various phenomena such as crystallization in solution, chemical reactivity, extraction of rare earth elements, and so forth. Despite this, the efforts to describe the structure of an aqueous solution have been hindered by the difficulty to retrieve structural data both from experiments and simulations. To overcome this, first, undersaturated EuCl3 aqueous solutions of concentrations varying from 0.15 to 1.8 mol/kg were studied using X-ray scattering. Second, for the first time, the theoretical X-ray signal of 1.8 mol/kg EuCl3 aqueous solution was simulated, with precise details for the complete range of scattering vectors using coupled molecular dynamics and hypernetted chain integral equations, and satisfactorily compared with the 1.8 mol/kg experimental X-ray scattering signal. The theoretical calculations demonstrate that the experimental structure factor is dominated by Eu3+-Eu3+ correlations.
Collapse
Affiliation(s)
| | - Maximilien Levesque
- PASTEUR, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Luc Belloni
- Université Paris-Saclay, CEA, CNRS, NIMBE, LIONS, 91191 Gif-sur-Yvette, France
| | - David Carriere
- Université Paris-Saclay, CEA, CNRS, NIMBE, LIONS, 91191 Gif-sur-Yvette, France
| |
Collapse
|
11
|
Dembowski M, Snyder MM, Delegard CH, Reynolds JG, Graham TR, Wang HW, Leavy II, Baum SR, Qafoku O, Fountain MS, Rosso KM, Clark SB, Pearce CI. Ion-ion interactions enhance aluminum solubility in alkaline suspensions of nano-gibbsite (α-Al(OH) 3) with sodium nitrite/nitrate. Phys Chem Chem Phys 2020; 22:4368-4378. [PMID: 31850442 DOI: 10.1039/c9cp05856g] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
Despite widespread industrial importance, predicting metal solubilities in highly concentrated, multicomponent aqueous solutions is difficult due to poorly understood ion-ion and ion-solvent interactions. Aluminum hydroxide solid phase solubility in concentrated sodium hydroxide (NaOH) solutions is one such case, with major implications for ore refining, as well as processing of radioactive waste stored at U.S. Department of Energy legacy sites, such as the Hanford Site, Washington State. The solubility of gibbsite (α-Al(OH)3) is often not well predicted because other ions affect the activity of hydroxide (OH-) and aluminate (Al(OH)4-) anions. In the present study, we systematically examined the influence of key anions, nitrite (NO2-) and nitrate (NO3-), as sodium salts on the solubility of α-Al(OH)3 in NaOH solutions taking care to establish equilibrium from both under- and oversaturation. Rapid equilibration was enabled by use of a highly pure and crystalline synthetic nano-gibbsite of well-defined particle size and shape. Measured dissolved aluminum concentrations were compared with those predicted by an α-Al(OH)3 solubility model derived for simple Al(OH)4-/OH- systems. Specific anion effects were expressed as an enhancement factor (Alenhc) conveying the excess of dissolved aluminum. At 45 °C, NaNO2 and NaNO3-containing systems exhibited Alenhc values of 2.70 and 1.88, respectively, indicating significant enhancement. The solutions were examined by Raman and high-field 27Al NMR spectroscopy, indicating specific interactions including Al(OH)4--Na+ contact ion pairing and Al(OH)4--NO2-/NO3- ion-ion interactions. Dynamic evolution of the α-Al(OH)3 particles including growth and agglomeration was observed revealing the importance of dissolution/reprecipitation in establishing equilibrium. These studies indicate that incomplete ion hydration, as a result of the low water activity in these concentrated electrolytes, results in: (i) enhanced reactivity of the hydroxide ion with respect to α-Al(OH)3; (ii) increased concentrations of Al(OH)4- in solution; and (iii) stronger ion-ion interactions that act to stabilize the supersaturated solutions. This information on the mechanisms by which α-Al(OH)3 becomes supersaturated is essential for more energy-efficient aluminum processing technologies, including the treatment of millions of gallons of Al(OH)4--rich high-level radioactive waste.
Collapse
Affiliation(s)
- Mateusz Dembowski
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA.
| | | | | | | | | | | | | | | | | | | | | | | | | |
Collapse
|