1
|
Liu HY, Mei KJ, Borrelli WR, Schwartz BJ. Simulating the Competitive Ion Pairing of Hydrated Electrons with Chaotropic Cations. J Phys Chem B 2024. [PMID: 39178349 DOI: 10.1021/acs.jpcb.4c04290] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/25/2024]
Abstract
Experiments show that the absorption spectrum of the hydrated electron (ehyd-) blue-shifts in electrolyte solutions compared with what is seen in pure water. This shift has been assigned to the ehyd-'s competitive ion-pairing interactions with the salt cation relative to the salt anion based on the ions' positions on the Hofmeister series. Remarkably, little work has been done investigating the ehyd-'s behavior when the salts have chaotropic cations, which should greatly change the ion-pairing interactions given that the ehyd- is a champion chaotrope. In this work, we remedy this by using mixed quantum/classical simulations to analyze the behavior of two different models of the ehyd- in aqueous RbF and RbI electrolyte solutions as a function of salt concentration. We find that the magnitude of the salt-induced spectral blue-shift is determined by a combination of the number of chaotropic Rb+ cations near the ehyd- and the number of salt anions near those cations so that the spectrum of the ehyd- directly reflects its local environment. We also find that the use of a soft-cavity ehyd- model predicts stronger competitive interactions with Rb+ relative to I- than a more traditional hard cavity model, leading to different predicted spectral shifts that should provide a way to distinguish between the two models experimentally. Our simulations predict that at the same concentration, salts with chaotropic cations should produce larger spectral blue-shifts than salts with kosmotropic cations. We also found that at high salt concentrations with chaotropic cations, the predicted blue-shift is greater when the salt anion is kosmotropic instead of chaotropic. Our goal is for this work to inspire experimentalists to make such measurements, which will help provide a spectroscopic means to distinguish between simulations models that predict different hydration structures for the ehyd-.
Collapse
Affiliation(s)
- Hannah Y Liu
- Department of Chemistry & Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Kenneth J Mei
- Department of Chemistry & Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - William R Borrelli
- Department of Chemistry & Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Benjamin J Schwartz
- Department of Chemistry & Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| |
Collapse
|
2
|
Markmann V, Pan J, Hansen BL, Haubro ML, Nimmrich A, Lenzen P, Levantino M, Katayama T, Adachi SI, Gorski-Bilke S, Temps F, Dohn AO, Møller KB, Nielsen MM, Haldrup K. Real-time structural dynamics of the ultrafast solvation process around photo-excited aqueous halides. Chem Sci 2024; 15:11391-11401. [PMID: 39055005 PMCID: PMC11268492 DOI: 10.1039/d4sc01912a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Accepted: 06/11/2024] [Indexed: 07/27/2024] Open
Abstract
This work investigates and describes the structural dynamics taking place following charge-transfer-to-solvent photo-abstraction of electrons from I- and Br- ions in aqueous solution following single- and 2-photon excitation at 202 nm and 400 nm, respectively. A Time-Resolved X-ray Solution Scattering (TR-XSS) approach with direct sensitivity to the structure of the surrounding solvent as the water molecules adopt a new equilibrium configuration following the electron-abstraction process is utilized to investigate the structural dynamics of solvent shell expansion and restructuring in real-time. The structural sensitivity of the scattering data enables a quantitative evaluation of competing models for the interaction between the nascent neutral species and surrounding water molecules. Taking the I0-O distance as the reaction coordinate, we find that the structural reorganization is delayed by 0.1 ps with respect to the photoexcitation and completes on a time scale of 0.5-1 ps. On longer time scales we determine from the evolution of the TR-XSS difference signal that I0: e- recombination takes place on two distinct time scales of ∼20 ps and 100 s of picoseconds. These dynamics are well captured by a simple model of diffusive evolution of the initial photo-abstracted electron population where the charge-transfer-to-solvent process gives rise to a broad distribution of electron ejection distances, a significant fraction of which are in the close vicinity of the nascent halogen atoms and recombine on short time scales.
Collapse
Affiliation(s)
- Verena Markmann
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Jaysree Pan
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Bianca L Hansen
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Morten L Haubro
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Amke Nimmrich
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
- Department of Chemistry and Molecular Biology, University of Gothenburg Gothenburg Sweden
| | - Philipp Lenzen
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Matteo Levantino
- European Synchrotron Radiation Facility CS40220 Grenoble 38043 Cedex 9 France
| | - Tetsuo Katayama
- Japan Synchrotron Radiation Research Institute Kouto 1-1-1, Sayo Hyogo 679-5198 Japan
- RIKEN SPring-8 Center 1-1-1 Kouto, Sayo Hyogo 679-5148 Japan
| | - Shin-Ichi Adachi
- Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK) 1-1 Oho, Tsukuba Ibaraki 305-0801 Japan
- Department of Materials Structure Science, School of High Energy Accelerator Science 1-1 Oho, Tsukuba Ibaraki 305-0801 Japan
| | | | - Friedrich Temps
- Christian-Albrechts-University Kiel Olshausenstr. 40 24098 Kiel Germany
| | - Asmus O Dohn
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
- Science Institute, University of Iceland 107 Reykjavík Iceland
| | - Klaus B Møller
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Martin M Nielsen
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Kristoffer Haldrup
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| |
Collapse
|
3
|
Singh N, Kashyap HK. Are NaTFSI and NaFSI Salt-Based Water-in-Salt Electrolytes Structurally Similar or Different? J Phys Chem B 2024. [PMID: 39045859 DOI: 10.1021/acs.jpcb.4c02863] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2024]
Abstract
Water-in-salt electrolytes (WiSEs) are a promising class of electrolytes due to their wide electrochemical stability window and nonflammability. In this study, we explore the structural organization of sodium bis(trifluoromethylsulfonyl)imide (NaTFSI) and sodium bis(fluorosulfonyl)imide (NaFSI) salt-based aqueous electrolytes, covering dilute to highly concentrated regions, by employing an all-atom molecular dynamics simulation. For the NaTFSI-based electrolyte, we observe that Na+ ions are mostly surrounded by water molecules at all the salt concentrations due to the very strong interaction between them. While TFSI anions weakly coordinate with Na+ ions and other TFSI anions, they also mostly prefer to be surrounded by water molecules. These interactions were found to have moderate dependence on the concentration of the NaTFSI salt. For the NaFSI-based electrolyte, while the Na+-water interaction is stronger at lower salt concentrations, the number of nearest neighbor FSI anions is found to be more than that of water at higher concentrations (≥20 m). This is because the increase in the salt concentration leads to expulsion of water molecules from the solvation shell of Na+ ions and enhances the interaction between Na+ ions and oxygen atoms of FSI. At the highest salt concentration (solubility limit), the bulk-like water structure is completely disrupted and dominated by an anionic network in the FSI-based electrolyte. In contrast, water-water hydrogen bonding network is still present even in the highly concentrated TFSI-based electrolyte. The simulated X-ray scattering pattern displays a low-q peak, revealing the presence of an intermediate range ordering due to alternating anion-rich and water/Na+-rich regions in both the electrolytes. However, the characteristic length scale corresponding to the low-q peak decreases with increasing the salt content in both the electrolytes.
Collapse
Affiliation(s)
- Navneet Singh
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Hemant K Kashyap
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|
4
|
Qin T, Yang H, Wang L, Xue W, Yao N, Li Q, Chen X, Yang X, Yu X, Zhang Q, Li H. Molecule Design for Non-Aqueous Wide-Temperature Electrolytes via the Intelligentized Screening Method. Angew Chem Int Ed Engl 2024:e202408902. [PMID: 38934230 DOI: 10.1002/anie.202408902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2024] [Revised: 06/06/2024] [Accepted: 06/26/2024] [Indexed: 06/28/2024]
Abstract
Operating a lithium-ion battery (LIB) in a wide temperature range is essential for ensuring a stable electricity supply amidst fluctuating temperatures caused by climate or terrain changes. Electrolyte plays a pivotal role in determining the temperature durability of batteries. However, specialized electrolytes designed for either low or high temperatures typically possess distinct features. Therefore, wide-temperature electrolytes (WTEs) are necessary as they encompass a combination of diverse properties, which complicates the clear instruction of WTE design. Here we represent an artificial intelligence (Al)-assisted workflow of WTE design through stepwise parameterizations and calculations. Linear mono-nitriles are identified as ideal wide-liquidus-range solvents that can "softly" solvate lithium ions by weak interactions. In addition, the explainable modules revealed the halogenoid similarity of cyanide as fluorine on the electrolyte properties (e.g. boiling point and dielectric constant). With the further introduction of an ether bond, 3-methoxypropionitrile (MPN) has been eventually determined as a main electrolyte solvent, enabling the battery operation from -60 to 120 °C. Particularly, a LiCoO2/Li cell using the proposed WTE can realize stable cycling with capacity retention reaching 72.3 % after 50 cycles under a high temperature of 100 °C.
Collapse
Affiliation(s)
- Tian Qin
- Beijing Frontier Research Center on Clean Energy, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
- Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Haoyi Yang
- Beijing Frontier Research Center on Clean Energy, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
| | - Lei Wang
- Beijing Frontier Research Center on Clean Energy, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
| | - Weiran Xue
- Beijing Frontier Research Center on Clean Energy, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
- Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Nan Yao
- Beijing Key Laboratory of Green Chemical Reaction Engineering and Technology, Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China
| | - Quan Li
- Beijing Frontier Research Center on Clean Energy, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
| | - Xiang Chen
- Beijing Key Laboratory of Green Chemical Reaction Engineering and Technology, Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China
| | - Xiukang Yang
- Key Laboratory of Environmentally Friendly Chemistry and Applications of Ministry of Education, National Base for International Science & Technology Cooperation, Hunan Province Key Laboratory for Electrochemical Energy Storage and Conversion, School of Chemistry, Xiangtan University, Hunan, Xiangtan, 411105, China
| | - Xiqian Yu
- Beijing Frontier Research Center on Clean Energy, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
- Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Qiang Zhang
- Beijing Key Laboratory of Green Chemical Reaction Engineering and Technology, Department of Chemical Engineering, Tsinghua University, Beijing, 100084, China
| | - Hong Li
- Beijing Frontier Research Center on Clean Energy, Institute of Physics, Chinese Academy of Sciences, Beijing, 100190, China
- Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing, 100049, China
| |
Collapse
|
5
|
Khammari A, Jeon J, Cho M. Solvation Structures and Transport Mechanisms of Cations in Water-in-Bisalt Electrolytes for High-Concentration Lithium-Ion Batteries Revealed by Molecular Dynamics Simulations. J Phys Chem B 2024; 128:5735-5745. [PMID: 38820050 DOI: 10.1021/acs.jpcb.4c02001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/02/2024]
Abstract
The development of safe and cost-effective electrolytes for rechargeable batteries is currently underway. While water-based electrolytes hold promise, their restricted electrochemical stability window poses a challenge. Combining multiple ionic species emerges as a promising strategy to broaden this stability window and optimize Li-ion battery performance. This study focuses on dual-cation electrolytes, which blend lithium and potassium acetates to enhance the electrochemical characteristics of the solution at high concentrations. We investigated the solvation structure of each ion and its interactions on a molecular level. Our analysis reveals that ion clusters and aggregates are formed through shared acetate and water molecules at high salt concentrations. Furthermore, the residence time analyses of atom pairs indicate that cations diffuse in vehicular mode at low concentrations. In contrast, they switch to a structural mode at high concentrations due to diminishing water content. This study offers a comprehensive model for exploring diverse solvation structures of cations and gaining insights into their diffusion mechanisms within water-in-bisalt electrolytes for aqueous Li-ion batteries.
Collapse
Affiliation(s)
- Anahita Khammari
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, Republic of Korea
- Korea University, Seoul 02841, Republic of Korea
| | - Jonggu Jeon
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, Republic of Korea
- Korea University, Seoul 02841, Republic of Korea
| | - Minhaeng Cho
- Center for Molecular Spectroscopy and Dynamics, Institute for Basic Science (IBS), Seoul 02841, Republic of Korea
- Department of Chemistry, Korea University, Seoul 02841, Republic of Korea
| |
Collapse
|
6
|
Ishraaq R, Das S. All-atom molecular dynamics simulations of polymer and polyelectrolyte brushes. Chem Commun (Camb) 2024; 60:6093-6129. [PMID: 38819435 DOI: 10.1039/d4cc01557f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Densely grafted polymer and polyelectrolyte (PE) brushes, owing to their significant abilities to functionalize surfaces for a plethora of applications in sensing, diagnostics, current rectification, surface wettability modification, drug delivery, and oil recovery, have attracted significant attention over the past several decades. Unfortunately, most of the attention has primarily focused on understanding the properties of the grafted polymer and the PE chains with little attention devoted to studying the behavior of the brush-supported ions (counterions needed to screen the PE chains) and water molecules. Over the past few years, our group has been at the forefront of addressing this gap: we have employed all-atom molecular dynamics (MD) simulations for studying a wide variety of polymer and PE brush systems with specific attention to unraveling the properties and behavior of the brush-supported water molecules and ions. Our findings have revealed some of the most fascinating properties of such brush-supported ions and water molecules, including the most remarkable control of nanofluidic transport afforded by the specific ion and water responses induced by the PE brushes grafted on the inner walls of the nanochannel. This feature article aims to summarize some of our key contributions associated with such atomistic simulations of polymer and PE brushes and brush-supported water molecules and counterions.
Collapse
Affiliation(s)
- Raashiq Ishraaq
- Department of Mechanical Engineering, University of Maryland, College Park, MD 20742, USA.
| | - Siddhartha Das
- Department of Mechanical Engineering, University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
7
|
Zhu Y, Lao Z, Zhang M, Hou T, Xiao X, Piao Z, Lu G, Han Z, Gao R, Nie L, Wu X, Song Y, Ji C, Wang J, Zhou G. A locally solvent-tethered polymer electrolyte for long-life lithium metal batteries. Nat Commun 2024; 15:3914. [PMID: 38724546 PMCID: PMC11082227 DOI: 10.1038/s41467-024-48078-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2023] [Accepted: 04/18/2024] [Indexed: 05/12/2024] Open
Abstract
Solid polymer electrolytes exhibit enhanced Li+ conductivity when plasticized with highly dielectric solvents such as N,N-dimethylformamide (DMF). However, the application of DMF-containing electrolytes in solid-state batteries is hindered by poor cycle life caused by continuous DMF degradation at the anode surface and the resulting unstable solid-electrolyte interphase. Here we report a composite polymer electrolyte with a rationally designed Hofmann-DMF coordination complex to address this issue. DMF is engineered on Hofmann frameworks as tethered ligands to construct a locally DMF-rich interface which promotes Li+ conduction through a ligand-assisted transport mechanism. A high ionic conductivity of 6.5 × 10-4 S cm-1 is achieved at room temperature. We demonstrate that the composite electrolyte effectively reduces the free shuttling and subsequent decomposition of DMF. The locally solvent-tethered electrolyte cycles stably for over 6000 h at 0.1 mA cm-2 in Li | |Li symmetric cell. When paired with sulfurized polyacrylonitrile cathodes, the full cell exhibits a prolonged cycle life of 1000 cycles at 1 C. This work will facilitate the development of practical polymer-based electrolytes with high ionic conductivity and long cycle life.
Collapse
Affiliation(s)
- Yanfei Zhu
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Zhoujie Lao
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Mengtian Zhang
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Tingzheng Hou
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China.
| | - Xiao Xiao
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Zhihong Piao
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Gongxun Lu
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Zhiyuan Han
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Runhua Gao
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Lu Nie
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Xinru Wu
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Yanze Song
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Chaoyuan Ji
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China
| | - Jian Wang
- Canadian Light Source, Saskatoon, S7N 2V3, Canada
| | - Guangmin Zhou
- Tsinghua-Berkeley Shenzhen Institute & Tsinghua Shenzhen International Graduate School, Shenzhen, 518055, PR China.
| |
Collapse
|
8
|
Thorat A, Chauhan R, Sartape R, Singh MR, Shah JK. Effect of K + Force Fields on Ionic Conductivity and Charge Dynamics of KOH in Ethylene Glycol. J Phys Chem B 2024; 128:3707-3719. [PMID: 38572661 PMCID: PMC11033864 DOI: 10.1021/acs.jpcb.3c08480] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2023] [Revised: 03/08/2024] [Accepted: 03/18/2024] [Indexed: 04/05/2024]
Abstract
Predicting ionic conductivity is crucial for developing efficient electrolytes for energy storage and conversion and other electrochemical applications. An accurate estimate of ionic conductivity requires understanding complex ion-ion and ion-solvent interactions governing the charge transport at the molecular level. Molecular simulations can provide key insights into the spatial and temporal behavior of electrolyte constituents. However, such insights depend on the ability of force fields to describe the underlying phenomena. In this work, molecular dynamics simulations were leveraged to delineate the impact of force field parameters on ionic conductivity predictions of potassium hydroxide (KOH) in ethylene glycol (EG). Four different force fields were used to represent the K+ ion. Diffusion-based Nernst-Einstein and correlation-based Einstein approaches were implemented to estimate the ionic conductivity, and the predicted values were compared with experimental measurements. The physical aspects, including ion-aggregation, charge distribution, cluster correlation, and cluster dynamics, were also examined. A force field was identified that provides reasonably accurate Einstein conductivity values and a physically coherent representation of the electrolyte at the molecular level.
Collapse
Affiliation(s)
- Amey Thorat
- School
of Chemical Engineering, Oklahoma State
University, Stillwater, Oklahoma 74078, United States
| | - Rohit Chauhan
- Department
of Chemical Engineering, University of Illinois
at Chicago, Chicago, Illinois 60608, United States
| | - Rohan Sartape
- Department
of Chemical Engineering, University of Illinois
at Chicago, Chicago, Illinois 60608, United States
| | - Meenesh R. Singh
- Department
of Chemical Engineering, University of Illinois
at Chicago, Chicago, Illinois 60608, United States
| | - Jindal K. Shah
- School
of Chemical Engineering, Oklahoma State
University, Stillwater, Oklahoma 74078, United States
| |
Collapse
|
9
|
Chio CC, Tse YLS. Reparameterization of Polarizable Force Fields for Studying Ion Transfer across Liquid-Liquid Interfaces. J Phys Chem B 2024; 128:1987-1999. [PMID: 38356148 DOI: 10.1021/acs.jpcb.3c07055] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2024]
Abstract
We have developed a general scheme for refining classical polarizable molecular dynamics (MD) force fields that can accurately describe the molecular interactions in systems with liquid-liquid interfaces. While ab initio MD (AIMD) simulations can naturally describe molecular interactions, they are often so computationally expensive that simulating large system sizes and/or long time scales is usually infeasible. To resolve this, we parameterized efficient and accurate classical polarizable force fields that use AIMD reference data by minimizing both the relative entropy and the root mean squared deviation in atomic forces. We utilized our new multiscale models to study chloride ion transfer across the water-dichloromethane (DCM) interface with and without the tetraethylammonium cation as the phase-transfer catalyst. Our calculated free-energy barrier for the water-DCM interface is consistent with the other reported simulation results. We further analyzed the ion-transfer process by studying the hydration shell structures around the chloride ion and the ion-pair formation to better understand the mechanism. We observed that electronic polarizability is an important factor for the studied phase-transfer catalyst to lower the free-energy barrier of the ion transfer. Using the water-benzene interface system as an additional example, we show that our parameterization scheme provides a general route for modeling different liquid-liquid interface systems even when the experimental data or force field parameters are not readily available.
Collapse
Affiliation(s)
- Chung Chi Chio
- Department of Chemistry, The Chinese University of Hong Kong, Sha Tin, New Territories, Hong Kong, China
| | - Ying-Lung Steve Tse
- Department of Chemistry, The Chinese University of Hong Kong, Sha Tin, New Territories, Hong Kong, China
| |
Collapse
|
10
|
Mondal J, Maji D, Biswas R. Temperature-dependent dielectric relaxation measurements of (acetamide + K/Na SCN) deep eutectic solvents: Decoding the impact of cation identity via computer simulations. J Chem Phys 2024; 160:084506. [PMID: 38421071 DOI: 10.1063/5.0193512] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2023] [Accepted: 02/07/2024] [Indexed: 03/02/2024] Open
Abstract
The impact of successive replacement of K+ by Na+ on the megahertz-gigahertz polarization response of 0.25[fKSCN + (1 - f)NaSCN] + 0.75CH3CONH2 deep eutectic solvents (DESs) was explored via temperature-dependent (303 ≤ T/K ≤ 343) dielectric relaxation (DR) measurements and computer simulations. Both the DR measurements (0.2 ≤ ν/GHz ≤ 50) and the simulations revealed multi-Debye relaxations accompanied by a decrease in the solution static dielectric constant (ɛs) upon the replacement of K+ by Na+. Accurate measurements of the DR response of DESs below 100 MHz were limited by the well-known one-over-frequency divergence for conducting solutions. This problem was tackled in simulations by removing the zero frequency contributions arising from the ion current to the total simulated DR response. The temperature-dependent measurements revealed a much stronger viscosity decoupling of DR times for Na+-containing DES than for the corresponding K+ system. The differential scanning calorimetry measurements indicated a higher glass transition temperature for Na+-DES (∼220 K) than K+-DES (∼200 K), implying more fragility and cooperativity for the former (Na+-DES) than the latter. The computer simulations revealed a gradual decrease in the average number of H bonds (⟨nHB⟩) per acetamide molecule and increased frustrations in the average orientational order upon the replacement of K+ by Na+. Both the measured and simulated ɛs values were found to decrease linearly with ⟨nHB⟩. Decompositions of the simulated DR spectra revealed that the cation-dependent cross interaction (dipole-ion) term contributes negligibly to ɛs and appears in the terahertz regime. Finally, the simulated collective single-particle reorientational relaxations and the structural H-bond fluctuation dynamics revealed the microscopic origin of the cation identity dependence shown by the measured DR relaxation times.
Collapse
Affiliation(s)
- Jayanta Mondal
- Department of Chemical and Biological Sciences, S. N. Bose National Centre for Basic Sciences, JD Block, Sector III, Salt Lake, Kolkata 700106, India
| | - Dhrubajyoti Maji
- Department of Chemical and Biological Sciences, S. N. Bose National Centre for Basic Sciences, JD Block, Sector III, Salt Lake, Kolkata 700106, India
| | - Ranjit Biswas
- Department of Chemical and Biological Sciences, S. N. Bose National Centre for Basic Sciences, JD Block, Sector III, Salt Lake, Kolkata 700106, India
| |
Collapse
|
11
|
Gullbrekken Ø, Gunnarshaug AF, Lervik A, Kjelstrup S, Schnell SK. Effect of the Ion, Solvent, and Thermal Interaction Coefficients on Battery Voltage. J Am Chem Soc 2024; 146:4592-4604. [PMID: 38340142 PMCID: PMC10885156 DOI: 10.1021/jacs.3c11589] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/12/2024]
Abstract
In order to increase the adoption of batteries for sustainable transport and energy storage, improved charging and discharging capabilities of lithium-ion batteries are necessary. To achieve this, accurate data that describe the internal state of the cells are essential. Several models have been derived, and transport coefficients have been reported for use in these models. We report for the first time a complete set of transport coefficients to model the concentration and temperature polarization in a lithium-ion battery ternary electrolyte, allowing us to test common assumptions. We include effects due to gradients in chemical potentials and in temperature. We find that the voltage contributions due to salt and solvent polarization are of the same order of magnitude as the ohmic loss and must be taken into account for more accurate modeling and understanding of battery performance. We report new Soret and Seebeck coefficients and find thermal polarization to be significant in cases relevant to battery research. The analysis is suitable for electrochemical systems, in general.
Collapse
Affiliation(s)
- Øystein Gullbrekken
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, N-7491 Trondheim, Norway
| | - Astrid Fagertun Gunnarshaug
- PoreLab, Department of Chemistry, Norwegian University of Science and Technology, NTNU, N-7491 Trondheim, Norway
| | - Anders Lervik
- PoreLab, Department of Chemistry, Norwegian University of Science and Technology, NTNU, N-7491 Trondheim, Norway
| | - Signe Kjelstrup
- PoreLab, Department of Chemistry, Norwegian University of Science and Technology, NTNU, N-7491 Trondheim, Norway
| | - Sondre Kvalvåg Schnell
- Department of Materials Science and Engineering, Norwegian University of Science and Technology, NTNU, N-7491 Trondheim, Norway
| |
Collapse
|
12
|
Wróbel P, Kubisiak P, Eilmes A. Ab Initio Molecular Dynamics Simulations of Aqueous LiTFSI Solutions─Structure, Hydrogen Bonding, and IR Spectra. J Phys Chem B 2024; 128:1001-1011. [PMID: 38235720 PMCID: PMC10839825 DOI: 10.1021/acs.jpcb.3c06633] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2023] [Revised: 12/30/2023] [Accepted: 01/05/2024] [Indexed: 01/19/2024]
Abstract
Simulations of Density Functional Theory-based ab initio molecular dynamics (AIMD) have been performed for a series of aqueous lithium bis(trifluoromethylsulfonyl)imide (LiTFSI) solutions with concentrations ranging from salt-in-water to water-in-salt systems. Analysis of the structure of electrolytes has revealed a preference of Li+ cations to interact with water molecules. In concentrated LiTFSI solutions, water molecules form small associates. The total number of hydrogen bonds (HBs) in the system decreases with salt concentrations, with bonds to water acceptors being only partially replaced by interactions with TFSI anions. Infrared (IR) spectra in the region of the O-H stretching frequency calculated from AIMD trajectories are in good agreement with experimental data. Statistics of oscillations of individual O-H bonds have shown correlations between vibrational frequencies and the structure of HBs formed by water. The changes in the IR spectrum have been related to the varying contributions of different local environments of the water molecules. The abundances of the three spectral components calculated from the simulations agree well with the decomposition of the experimental IR spectra reported in the literature.
Collapse
Affiliation(s)
- Piotr Wróbel
- Faculty of Chemistry, Jagiellonian
University, Gronostajowa 2, 30-387 Kraków, Poland
| | - Piotr Kubisiak
- Faculty of Chemistry, Jagiellonian
University, Gronostajowa 2, 30-387 Kraków, Poland
| | - Andrzej Eilmes
- Faculty of Chemistry, Jagiellonian
University, Gronostajowa 2, 30-387 Kraków, Poland
| |
Collapse
|
13
|
Blazquez S, Bourg IC, Vega C. Madrid-2019 force field: An extension to divalent cations Sr2+ and Ba2. J Chem Phys 2024; 160:046101. [PMID: 38270237 DOI: 10.1063/5.0186233] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2023] [Accepted: 12/27/2023] [Indexed: 01/26/2024] Open
Abstract
In this work, we present a parameterization of Sr2+ and Ba2+ cations, which expands the alkali earth set of cations of the Madrid-2019 force field. We have tested the model against the experimental densities of eight different salts, namely, SrCl2, SrBr2, SrI2, Sr(NO3)2, BaCl2, BaBr2, BaI2, and Ba(NO3)2. The force field is able to reproduce the experimental densities of all these salts up to their solubility limit. Furthermore, we have computed the viscosities for two selected salts, finding that the experimental values are overestimated, but the predictions are still reasonable. Finally, the structural properties for all the salts have been calculated with this model and align remarkably well with experimental observations.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Ian C Bourg
- Department of Civil and Environmental Engineering, Princeton University, Princeton, New Jersey 08544, USA
- High Meadows Environmental Institute, Princeton University, Princeton, New Jersey 08544, USA
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
14
|
Kwan V, Consta S, Malek SMA. Variation of Surface Propensity of Halides with Droplet Size and Temperature: The Planar Interface Limit. J Phys Chem B 2024; 128:193-207. [PMID: 38127582 DOI: 10.1021/acs.jpcb.3c05701] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2023]
Abstract
The radial number density profiles of halide and alkali ions in aqueous clusters with equimolar radius ≲1.4 nm, which correspond to ≲255 H2O molecules, have been extensively studied by computations. However, the surface abundance of Cl-, Br-, and I- relative to the bulk interior in these smaller clusters may not be representative of the larger systems. Indeed, here we show that the larger the cluster is, the lower the relative surface abundance of chaotropic halides is. In droplets with an equimolar radius of ≈2.45 nm, which corresponds to ≈2000 H2O molecules, the polarizable halides show a clear number density maximum in the droplet's bulk-like interior. A similar pattern is observed in simulations of the aqueous planar interface with halide salts at room temperature. At elevated temperature the surface propensity of Cl- decreases gradually, while that of I- is partially preserved. The change in the chaotropic halide location at higher temperatures than the room temperature may considerably affect photochemical reactivity in atmospheric aerosols, vapor-liquid nucleation and growth mechanisms, and salt crystallization via solvent evaporation. We argue that the commonly used approach of nullifying parameters in a force field in order to find the factors that determine the ion location does not provide transferable insight into other force fields.
Collapse
Affiliation(s)
- Victor Kwan
- Department of Chemistry, The University of Western Ontario, London, ON, Canada N6A 5B7
| | - Styliani Consta
- Department of Chemistry, The University of Western Ontario, London, ON, Canada N6A 5B7
| | - Shahrazad M A Malek
- Department of Physics and Physical Oceanography, Memorial University of Newfoundland, St. John's, NL, Canada A1B 3X7
| |
Collapse
|
15
|
Liu Y, Lin Y, Yang Z, Lin C, Zhang X, Chen S, Hu G, Sa B, Chen Y, Zhang Y. Stable Harsh-Temperature Lithium Metal Batteries Enabled by Tailoring Solvation Structure in Ether Electrolytes. ACS NANO 2023; 17:19625-19639. [PMID: 37819135 DOI: 10.1021/acsnano.3c01895] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/13/2023]
Abstract
For lithium metal batteries (LMBs), the elevated operating temperature results in severe capacity fading and safety issues due to unstable electrode-electrolyte interphases and electrolyte solvation structures. Therefore, it is crucial to construct advanced electrolytes capable of tolerating harsh environments to ensure stable LMBs. Here, we proposed a stable localized high-concentration electrolyte (LHCE) by introducing the highly solvating power solvent diethylene glycol dimethyl ether (DGDME). Computational and experimental evidence discloses that the original DGDME-LHCE shows favorable features for high-temperature LMBs, including high Li+-binding stability, electro-oxidation resistance, thermal stability, and nonflammability. The tailored solvated sheath structure achieves the preferred decomposition of anions, inducing the stable (cathode and Li anode)/interphases simultaneously, which enables a homogeneous Li plating-stripping behavior on the anode side and a high-voltage tolerance on the cathode side. For the Li||Li cells coupled with DGDME-LHCE, they showcase outstanding reversibility (a long lifespan of exceeding 1900 h). We demonstrate exceptional cyclic stability (∼95.59%, 250 cycles), high Coulombic efficiency (>99.88%), and impressive high-voltage (4.5 V) and high-temperature (60 °C) performances in Li||NCM523 cells using DGDME-LHCE. Our advances shed light on an encouraging ether electrolyte tactic for the Li-metal batteries confronted with stringent high-temperature challenges.
Collapse
Affiliation(s)
- Yongchuan Liu
- Key Laboratory of Optoelectronic Materials Chemistry and Physics, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou, Fujian 350002, People's Republic of China
| | - Yuansheng Lin
- Key Laboratory of Optoelectronic Materials Chemistry and Physics, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou, Fujian 350002, People's Republic of China
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, People's Republic of China
| | - Zhanlin Yang
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, People's Republic of China
| | - Changxin Lin
- Key Laboratory of Optoelectronic Materials Chemistry and Physics, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou, Fujian 350002, People's Republic of China
- University of Chinese Academy of Sciences, Beijing 100049, People's Republic of China
| | - Xiangxin Zhang
- Key Laboratory of Optoelectronic Materials Chemistry and Physics, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou, Fujian 350002, People's Republic of China
| | - Sujing Chen
- Key Laboratory of Optoelectronic Materials Chemistry and Physics, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou, Fujian 350002, People's Republic of China
| | - Guolin Hu
- Key Laboratory of Optoelectronic Materials Chemistry and Physics, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou, Fujian 350002, People's Republic of China
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, People's Republic of China
| | - Baisheng Sa
- Key Laboratory of Eco-materials Advanced Technology, College of Materials Science and Engineering, Fuzhou University, Fuzhou 350108, People's Republic of China
| | - Yuanqiang Chen
- Key Laboratory of Optoelectronic Materials Chemistry and Physics, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou, Fujian 350002, People's Republic of China
| | - Yining Zhang
- Key Laboratory of Optoelectronic Materials Chemistry and Physics, Fujian Institute of Research on the Structure of Matter, Chinese Academy of Sciences, 155 Yangqiao Road West, Fuzhou, Fujian 350002, People's Republic of China
| |
Collapse
|
16
|
Taylor PA, Stevens MJ. Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2023; 46:97. [PMID: 37831216 DOI: 10.1140/epje/s10189-023-00355-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Accepted: 09/26/2023] [Indexed: 10/14/2023]
Abstract
Strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties.
Collapse
Affiliation(s)
- Phillip A Taylor
- Sandia National Laboratories, Center for Integrated Nanotechnologies, Albuquerque, NM, 87123, USA
| | - Mark J Stevens
- Sandia National Laboratories, Center for Integrated Nanotechnologies, Albuquerque, NM, 87123, USA.
| |
Collapse
|
17
|
Vigil DL, Stevens MJ, Frischknecht AL. Accurate Calculation of Solvation Properties of Lithium Ions in Nonaqueous Solutions. J Phys Chem B 2023; 127:8002-8008. [PMID: 37676921 DOI: 10.1021/acs.jpcb.3c05591] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/09/2023]
Abstract
We perform all-atom molecular dynamics simulations of lithium triflate in 1,2-dimethoxyethane using six different literature force fields. This system is representative of many experimental studies of lithium salts in solvents and polymers. We show that multiple historically common force fields for lithium ions give qualitatively incorrect results when compared with those from experiments and quantum chemistry calculations. We illustrate the importance of correctly selecting force field parameters and give recommendations on the force field choice for lithium electrolyte applications.
Collapse
Affiliation(s)
- Daniel L Vigil
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Mark J Stevens
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Amalie L Frischknecht
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| |
Collapse
|
18
|
Blazquez S, Abascal JLF, Lagerweij J, Habibi P, Dey P, Vlugt TJH, Moultos OA, Vega C. Computation of Electrical Conductivities of Aqueous Electrolyte Solutions: Two Surfaces, One Property. J Chem Theory Comput 2023; 19:5380-5393. [PMID: 37506381 PMCID: PMC10448725 DOI: 10.1021/acs.jctc.3c00562] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2023] [Indexed: 07/30/2023]
Abstract
In this work, we computed electrical conductivities under ambient conditions of aqueous NaCl and KCl solutions by using the Einstein-Helfand equation. Common force fields (charge q = ±1 e) do not reproduce the experimental values of electrical conductivities, viscosities, and diffusion coefficients. Recently, we proposed the idea of using different charges to describe the potential energy surface (PES) and the dipole moment surface (DMS). In this work, we implement this concept. The equilibrium trajectories required to evaluate electrical conductivities (within linear response theory) were obtained by using scaled charges (with the value q = ±0.75 e) to describe the PES. The potential parameters were those of the Madrid-Transport force field, which accurately describe viscosities and diffusion coefficients of these ionic solutions. However, integer charges were used to compute the conductivities (thus describing the DMS). The basic idea is that although the scaled charge describes the ion-water interaction better, the integer charge reflects the value of the charge that is transported due to the electric field. The agreement obtained with experiments is excellent, as for the first time electrical conductivities (and the other transport properties) of NaCl and KCl electrolyte solutions are described with high accuracy for the whole concentration range up to their solubility limit. Finally, we propose an easy way to obtain a rough estimate of the actual electrical conductivity of the potential model under consideration using the approximate Nernst-Einstein equation, which neglects correlations between different ions.
Collapse
Affiliation(s)
- Samuel Blazquez
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jose L. F. Abascal
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Jelle Lagerweij
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Parsa Habibi
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD Delft, The Netherlands
| | - Poulumi Dey
- Department
of Materials Science and Engineering, Faculty of Mechanical, Maritime
and Materials Engineering, Delft University
of Technology, Mekelweg
2, 2628CD Delft, The Netherlands
| | - Thijs J. H. Vlugt
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Othonas A. Moultos
- Engineering
Thermodynamics, Process and Energy Department, Faculty of Mechanical,
Maritime and Materials Engineering, Delft
University of Technology, Leeghwaterstraat 39, 2628CB Delft, The Netherlands
| | - Carlos Vega
- Dpto.
Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
19
|
Bae J, Zhu Z, Yan J, Kim DM, Ko Y, Jain A, Helms BA. Closed-loop cathode recycling in solid-state batteries enabled by supramolecular electrolytes. SCIENCE ADVANCES 2023; 9:eadh9020. [PMID: 37566660 PMCID: PMC10421023 DOI: 10.1126/sciadv.adh9020] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/21/2023] [Accepted: 07/12/2023] [Indexed: 08/13/2023]
Abstract
Deconstructing solid-state batteries (SSBs) to physically separated cathode and solid-electrolyte particles remains intensive, as does the remanufacturing of cathodes and separators from the recovered materials. To address this challenge, we designed supramolecular organo-ionic (ORION) electrolytes that are viscoelastic solids at battery operating temperatures (-40° to 45°C) yet are viscoelastic liquids above 100°C, which enables both the fabrication of high-quality SSBs and the recycling of their cathodes at end of life. SSBs implementing ORION electrolytes alongside Li metal anodes and either LFP or NMC cathodes were operated for hundreds of cycles at 45°C with less than 20% capacity fade. Using a low-temperature solvent process, we isolated the cathode from the electrolyte and demonstrated that refurbished cells recover 90% of their initial capacity and sustain it for an additional 100 cycles with 84% capacity retention in their second life.
Collapse
Affiliation(s)
- Jiwoong Bae
- The Molecular Foundry, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
| | - Zhuoying Zhu
- Energy Storage and Distributed Resources Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
| | - Jiajun Yan
- Materials Sciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
| | - Dong-Min Kim
- The Molecular Foundry, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
- Joint Center for Energy Storage Research, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
| | - Youngmin Ko
- The Molecular Foundry, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
| | - Anubhav Jain
- Energy Storage and Distributed Resources Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
| | - Brett A. Helms
- The Molecular Foundry, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
- Materials Sciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA
| |
Collapse
|
20
|
Amărandi RM, Al-Matarneh MC, Popovici L, Ciobanu CI, Neamțu A, Mangalagiu II, Danac R. Exploring Pyrrolo-Fused Heterocycles as Promising Anticancer Agents: An Integrated Synthetic, Biological, and Computational Approach. Pharmaceuticals (Basel) 2023; 16:865. [PMID: 37375812 DOI: 10.3390/ph16060865] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Revised: 05/17/2023] [Accepted: 06/09/2023] [Indexed: 06/29/2023] Open
Abstract
Five new series of pyrrolo-fused heterocycles were designed through a scaffold hybridization strategy as analogs of the well-known microtubule inhibitor phenstatin. Compounds were synthesized using the 1,3-dipolar cycloaddition of cycloimmonium N-ylides to ethyl propiolate as a key step. Selected compounds were then evaluated for anticancer activity and ability to inhibit tubulin polymerization in vitro. Notably, pyrrolo[1,2-a]quinoline 10a was active on most tested cell lines, performing better than control phenstatin in several cases, most notably on renal cancer cell line A498 (GI50 27 nM), while inhibiting tubulin polymerization in vitro. In addition, this compound was predicted to have a promising ADMET profile. The molecular details of the interaction between compound 10a and tubulin were investigated through in silico docking experiments, followed by molecular dynamics simulations and configurational entropy calculations. Of note, we found that some of the initially predicted interactions from docking experiments were not stable during molecular dynamics simulations, but that configurational entropy loss was similar in all three cases. Our results suggest that for compound 10a, docking experiments alone are not sufficient for the adequate description of interaction details in terms of target binding, which makes subsequent scaffold optimization more difficult and ultimately hinders drug design. Taken together, these results could help shape novel potent antiproliferative compounds with pyrrolo-fused heterocyclic cores, especially from an in silico methodological perspective.
Collapse
Affiliation(s)
- Roxana-Maria Amărandi
- TRANSCEND Research Center, Regional Institute of Oncology Iasi, 2-4 General Henri Mathias Berthelot Street, 700483 Iasi, Romania
| | - Maria-Cristina Al-Matarneh
- "Petru Poni" Institute of Macromolecular Chemistry of Romanian Academy, 41A Grigore Ghica Voda Alley, 700487 Iasi, Romania
- Faculty of Chemistry, Alexandru Ioan Cuza University of Iasi, 11 Carol I, 700506 Iasi, Romania
| | - Lăcrămioara Popovici
- Faculty of Chemistry, Alexandru Ioan Cuza University of Iasi, 11 Carol I, 700506 Iasi, Romania
| | - Catalina Ionica Ciobanu
- Institute of Interdisciplinary Research-CERNESIM Centre, Alexandru Ioan Cuza University of Iasi, 11 Carol I, 700506 Iasi, Romania
| | - Andrei Neamțu
- TRANSCEND Research Center, Regional Institute of Oncology Iasi, 2-4 General Henri Mathias Berthelot Street, 700483 Iasi, Romania
| | - Ionel I Mangalagiu
- Faculty of Chemistry, Alexandru Ioan Cuza University of Iasi, 11 Carol I, 700506 Iasi, Romania
| | - Ramona Danac
- Faculty of Chemistry, Alexandru Ioan Cuza University of Iasi, 11 Carol I, 700506 Iasi, Romania
| |
Collapse
|
21
|
Maji D, Biswas R. Dielectric relaxation and dielectric decrement in ionic acetamide deep eutectic solvents: Spectral decomposition and comparison with experiments. J Chem Phys 2023; 158:2888209. [PMID: 37139998 DOI: 10.1063/5.0147378] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2023] [Accepted: 04/17/2023] [Indexed: 05/05/2023] Open
Abstract
Frequency-dependent dielectric relaxation in three deep eutectic solvents (DESs), (acetamide+LiClO4/NO3/Br), was investigated in the temperature range, 329 ≤ T/K ≤ 358, via molecular dynamics simulations. Subsequently, decomposition of the real and the imaginary components of the simulated dielectric spectra was carried out to separate the rotational (dipole-dipole), translational (ion-ion), and ro-translational (dipole-ion) contributions. The dipolar contribution, as expected, was found to dominate all the frequency-dependent dielectric spectra over the entire frequency regime, while the other two components together made tiny contributions only. The translational (ion-ion) and the cross ro-translational contributions appeared in the THz regime in contrast to the viscosity-dependent dipolar relaxations that dominated the MHz-GHz frequency window. Our simulations predicted, in agreement with experiments, anion-dependent decrement of the static dielectric constant (ɛs ∼ 20 to 30) for acetamide (ɛs ∼ 66) in these ionic DESs. Simulated dipole-correlations (Kirkwood g factor) indicated significant orientational frustrations. The frustrated orientational structure was found to be associated with the anion-dependent damage of the acetamide H-bond network. Single dipole reorientation time distributions suggested slowed down acetamide rotations but did not indicate presence of any "rotationally frozen" molecule. The dielectric decrement is, therefore, largely static in origin. This provides a new insight into the ion dependence of the dielectric behavior of these ionic DESs. A good agreement between the simulated and the experimental timescales was also noticed.
Collapse
Affiliation(s)
- Dhrubajyoti Maji
- Department of Chemical and Biological Sciences, S. N. Bose National Centre for Basic Sciences, Block-JD, Sector III, Salt Lake, Kolkata, West Bengal 700106, India
| | - Ranjit Biswas
- Department of Chemical and Biological Sciences, S. N. Bose National Centre for Basic Sciences, Block-JD, Sector III, Salt Lake, Kolkata, West Bengal 700106, India
| |
Collapse
|
22
|
Singh N, Jha MK, Dhattarwal HS, Kashyap HK. How NaFTA salt affects the structural landscape and transport properties of Pyrr 1,3FTA ionic liquid. J Chem Phys 2023; 158:104502. [PMID: 36922141 DOI: 10.1063/5.0133966] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/12/2023] Open
Abstract
Recently, it has been demonstrated that ionic liquids (ILs) with an asymmetric anion render a wider operational temperature range and can be used as a solvent in sodium ion batteries. In the present study, we examine the microscopic structure and dynamics of pure 1-methyl-1-propylpyrrolidinium fluorosulfonyl(trifluoromethylsulfonyl)amide (Pyrr1,3FTA) IL using atomistic molecular dynamics simulations. How the addition of the sodium salt (NaFTA) having the same anion changes the structural landscape and transport properties of the pure IL has also been explored. The simulated x-ray scattering structure functions reveal that the gradual addition of NaFTA salt (up to 1.2 molal) suppresses the charge alternating feature of the pure IL because of the replacement of the Pyrr+ cations with the Na+ ions. The Na+ ions are majorly found near the oxygen atoms of the anions, but the probability of finding the Na+ ions near these atoms slightly decreases with increasing salt concentration. As expected, the Na+ ions stay away from the Pyrr+ cations. However, the probability of finding the anions around anions increases with increasing salt concentration. The simulated self-diffusion coefficients of the ions in the pure IL reveal slightly faster diffusion of the Pyrr+ cations as compared to the FTA- anions. Interestingly, in the salt solution, despite having smaller size, the diffusion of the Na+ ions is found to be lesser than the Pyrr+ cations and the FTA- anions. The analysis of the ionic conductivity and transport numbers reveals that the fractional contribution of the FTA- anion to the overall conductivity remains nearly constant with increasing salt concentration, but the contribution of Pyrr+ cation decreases and Na+ ion increases.
Collapse
Affiliation(s)
- Navneet Singh
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Mrityunjay K Jha
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Harender S Dhattarwal
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Hemant K Kashyap
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|
23
|
Olsen R, Kvamme B. Impact of ethylene glycol on ions influencing corrosion in pores between iron oxide and calcium carbonate. MOLECULAR SIMULATION 2023. [DOI: 10.1080/08927022.2023.2184298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/11/2023]
Affiliation(s)
- Richard Olsen
- Department of Physics and Technology, University of Bergen, Bergen, Norway
| | | |
Collapse
|
24
|
Blazquez S, Conde MM, Vega C. Scaled charges for ions: An improvement but not the final word for modeling electrolytes in water. J Chem Phys 2023; 158:054505. [PMID: 36754806 DOI: 10.1063/5.0136498] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
In this work, we discuss the use of scaled charges when developing force fields for NaCl in water. We shall develop force fields for Na+ and Cl- using the following values for the scaled charge (in electron units): ±0.75, ±0.80, ±0.85, and ±0.92 along with the TIP4P/2005 model of water (for which previous force fields were proposed for q = ±0.85 and q = ±1). The properties considered in this work are densities, structural properties, transport properties, surface tension, freezing point depression, and maximum in density. All the developed models were able to describe quite well the experimental values of the densities. Structural properties were well described by models with charges equal to or larger than ±0.85, surface tension by the charge ±0.92, maximum in density by the charge ±0.85, and transport properties by the charge ±0.75. The use of a scaled charge of ±0.75 is able to reproduce with high accuracy the viscosities and diffusion coefficients of NaCl solutions for the first time. We have also considered the case of KCl in water, and the results obtained were fully consistent with those of NaCl. There is no value of the scaled charge able to reproduce all the properties considered in this work. Although certainly scaled charges are not the final word in the development of force fields for electrolytes in water, its use may have some practical advantages. Certain values of the scaled charge could be the best option when the interest is to describe certain experimental properties.
Collapse
Affiliation(s)
- S Blazquez
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M M Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - C Vega
- Dpto. Química Física I, Fac. Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
25
|
Banerjee S, Ghorai PK, Maji D, Biswas R. Difference in "Supercooling" Affinity between (Acetamide + Na/KSCN) Deep Eutectics: Reflections in the Simulated Anomalous Motions of the Constituents and Solution Microheterogeneity Features. J Phys Chem B 2022; 126:10146-10155. [PMID: 36414001 DOI: 10.1021/acs.jpcb.2c04994] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
Deep depression of freezing points of ionic amide deep eutectic solvents (DESs) is known to exhibit a significant dependence on the identity of ions present in those systems and the nature of the functional group attached to the host amide. This deep depression of the freezing point is sometimes termed as "supercooling". For (acetamide + electrolyte) DESs, experiments have revealed signatures of ion-dependent spatiotemporal heterogeneity features. The focus of this work is to provide microscopic explanations of these experimentally observed macroscopic system properties in terms of particle jumps and insights about the origin of the cation dependence. For this purpose, extensive molecular dynamics simulations have been performed employing (acetamide + Na/KSCN) deep eutectics as representative ionic systems at 303, 318, 333, and 348 K. The individual translational motions of acetamide and the ions are followed, and their connections to solution heterogeneity are explored. The center-of-mass motion for Na+ has been found to be more anomalous than that for K+. This difference corroborates well with experimental reports on heterogeneous relaxations in these systems. Simulated viscosity coefficients and dynamic heterogeneity features also reflect this difference. Moreover, simulated reorientational relaxations of acetamide molecules in these ionic DESs suggest that a Na+-containing DES is more heterogeneous than the corresponding K+-containing system. Estimated void and neck distributions for acetamide molecules differ as the alkali metal ions differ. In brief, this study provides a detailed microscopic view of the cation dependence of the microheterogeneous relaxation dynamics of these DESs reported repeatedly by different experiments.
Collapse
Affiliation(s)
- Swarup Banerjee
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata, Mohanpur 741246, India
| | - Pradip Kr Ghorai
- Department of Chemical Sciences, Indian Institute of Science Education and Research (IISER) Kolkata, Mohanpur 741246, India
| | - Dhrubajyoti Maji
- Department of Chemical, Biological and Macromolecular Sciences, S.N. Bose National Centre for Basic Sciences, Kolkata 700106, India
| | - Ranjit Biswas
- Department of Chemical, Biological and Macromolecular Sciences, S.N. Bose National Centre for Basic Sciences, Kolkata 700106, India
| |
Collapse
|
26
|
Effect of alkali metal ions on water structure: insights into the pressure-like effect. Struct Chem 2022. [DOI: 10.1007/s11224-022-02101-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
27
|
Tonti L, Floris FM. Hydrophilic Versus Hydrophobic Coupling in the Pressure Dependence of the Chemical Potential of Alkali Metal and Halide Ions in Water. J Phys Chem B 2022; 126:9325-9338. [PMID: 36326490 PMCID: PMC9677433 DOI: 10.1021/acs.jpcb.2c02373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
We computed the chemical potential for some alkali metal ions (K+, Rb+, and Cs+) and two halide ions (Br- and I-) in aqueous solution at ambient T and various pressures in the range 1-8000 atm. Results were obtained from classic Monte Carlo simulations in the NPT ensemble by means of the free energy perturbation method. Here, the chemical potential is computed as the sum of a term relative to a Lennard-Jones solute and a term relative to the process in which this solute is transformed into the ion. Hydrophobic and hydrophilic features of these two components of the chemical potential show opposite behaviors under isothermal compression. The increase in pressure determines an increase in the hydrophobic component, which becomes more positive with a stronger effect for larger ions. Correspondingly, the values of the hydrophilic component become more negative for alkali ions, whereas they are only slightly affected by compression for halide ions. Hydrophobic-hydrophilic quasi-compensation in the slopes is observed for Rb+. For a smaller ion, such as K+, the dependence on pressure of the hydrophilic component is slightly dominant. For a larger ion, as observed in the cases of Cs+, Br-, and I-, the hydrophobic component assumes the determinant role. Pressure dependence of the chemical potential is little affected by corrections introduced for molecular potential truncation. This view can change for possible boundary artifacts that could have affected the static electrostatic potential. Some inference is obtained from comparison with experimental data at 1 atm on the free energy of hydration. Discrepancies show the characteristic asymmetry between cations and anions. The further addition of a correction based on the static potential significantly reduces these discrepancies with important error cancellation on the sum of chemical potentials of ions of opposite charge. The correction is applied also at higher pressures, and results are compared with those obtained by adding an alternative correction that is based on the water number density. Regardless of the ion, changes of the chemical potential induced by an increase in pressure appear to be dominated by the hydrophobic component, in particular when using the alternative correction. For bromide and iodide electrolytes, the two corrections give chemical potentials in good agreement.
Collapse
Affiliation(s)
- Luca Tonti
- Department
of Chemical Engineering, The University
of Manchester, M13 9PLManchester, U.K.,
| | - Franca Maria Floris
- Dipartimento
di Chimica e Chimica Industriale, Università
di Pisa, Via Giuseppe
Moruzzi 13, 56124Pisa, Italy,
| |
Collapse
|
28
|
Barbosa GD, Turner CH. Molecular-level solvation and selectivity behavior of Na+, K+, and Li+ within glycerol-derived solvents. Chem Eng Sci 2022. [DOI: 10.1016/j.ces.2022.117992] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
29
|
Xi C, Zheng F, Gao G, Song Z, Zhang B, Dong C, Du XW, Wang LW. Ion Solvation Free Energy Calculation Based on Ab Initio Molecular Dynamics Using a Hybrid Solvent Model. J Chem Theory Comput 2022; 18:6878-6891. [PMID: 36253911 DOI: 10.1021/acs.jctc.1c01298] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Free energy calculation of small molecules or ion species in aqueous solvent is one of the most important problems in electrochemistry study. Although there are many previous approaches to calculate such free energies, they are based on ab initio methods and suffer from various limitations and approximations. In the current work, we developed a hybrid approach based on ab initio molecular dynamics (AIMD) simulations to calculate the ion solvation energy. In this approach, a small water cluster surrounding the central ion is used, and implicit solvent model is applied outside the water cluster. A dynamic potential well is used during AIMD to keep the water cluster together. Quasi-harmonic approximation is used to calculate the entropy contribution, while the total energy average is used to calculate the enthalpy term. The obtained solvation voltages of the bulk metal agree with experiments within 0.3 eV, and the simulation results for the solvation energies of gaseous ions are close to the experimental observations. Besides the free energies, radial pair distribution functions and coordination numbers of hydrated cations are also obtained. The remaining challenges of this method are also discussed.
Collapse
Affiliation(s)
- Cong Xi
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States.,Institute of New Energy Materials, School of Materials Science and Engineering, Tianjin University, Tianjin30072, People's Republic of China
| | - Fan Zheng
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Guoping Gao
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Zhigang Song
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Buyu Zhang
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Cunku Dong
- Institute of New Energy Materials, School of Materials Science and Engineering, Tianjin University, Tianjin30072, People's Republic of China
| | - Xi-Wen Du
- Institute of New Energy Materials, School of Materials Science and Engineering, Tianjin University, Tianjin30072, People's Republic of China
| | - Lin-Wang Wang
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| |
Collapse
|
30
|
Peng L, Wu X, Jia M, Qian W, Zhang X, Zhou N, Zhang L, Jian C, Zhang S. Solvating power regulation enabled low concentration electrolyte for lithium batteries. Sci Bull (Beijing) 2022; 67:2235-2244. [DOI: 10.1016/j.scib.2022.10.008] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2022] [Revised: 09/22/2022] [Accepted: 10/08/2022] [Indexed: 11/06/2022]
|
31
|
Pethes I. Towards the correct microscopic structure of aqueous CsCl solutions with a comparison of classical interatomic potential models. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.119660] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
32
|
Song YW, Shen L, Yao N, Li XY, Bi CX, Li Z, Zhou MY, Zhang XQ, Chen X, Li BQ, Huang JQ, Zhang Q. Cationic lithium polysulfides in lithium–sulfur batteries. Chem 2022. [DOI: 10.1016/j.chempr.2022.07.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
33
|
Dhattarwal HS, Kashyap HK. Heterogeneity and Nanostructure of Superconcentrated LiTFSI-EmimTFSI Hybrid Aqueous Electrolytes: Beyond the 21 m Limit of Water-in-Salt Electrolyte. J Phys Chem B 2022; 126:5291-5304. [PMID: 35819799 DOI: 10.1021/acs.jpcb.2c02822] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ionic liquids such as EmimTFSI (1-ethyl-3-methylimidazolium bis(trifluoromethanesulfonyl)imide) have been found to improve the solubility of LiTFSI salt in water-in-salt electrolyte (WiSE) from 21 to 60 m. However, the molecular origin of such enhancement in the solubility is still unknown. In the present work, we elucidate the microscopic structures of LiTFSI-EmimTFSI-based hybrid aqueous electrolytes and compare them with the structure of LiTFSI-based WiSE using molecular dynamics simulations. Our analysis reveals the presence of alternating water-rich clusters and TFSI-rich extended domains in the WiSE. In these clusters and domains, the Li+ ions reside such that the total number of oxygen atoms around them is conserved to four, where water contributes about three oxygen atoms. The addition of EmimTFSI in the WiSE results in removal of water from the nearest-neighbor solvation shell of TFSI- ions but not from the Li+ ions. Significant structural changes are observed when LiTFSI salt is further added to LiTFSI-EmimTFSI aqueous solution. In both the hybrid electrolytes, water and Emim+ cations are found to avoid each other. The simulated X-ray scattering structure factor reveals the presence of larger length-scale heterogeneity in the most concentrated solution of the hybrid aqueous electrolyte. We observe that this nanoscale heterogeneity originates from a water-TFSI-Emim-TFSI-water-TFSI-Emim-TFSI-like arrangement in which Li+ ions are dispersed such that the coordination number of oxygen atoms around them is enhanced to five, wherein the major contribution comes from the TFSI- ions. We envision that the enhanced LiTFSI solubility originates from the replacement of water molecules with TFSI- ions in the first solvation shell of Li+ ions.
Collapse
Affiliation(s)
- Harender S Dhattarwal
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Hemant K Kashyap
- Department of Chemistry, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| |
Collapse
|
34
|
Barbosa GD, Bara JE, Turner CH. Molecular simulation of glycerol-derived triether podands for lithium ion solvation. Phys Chem Chem Phys 2022; 24:9459-9466. [PMID: 35388849 DOI: 10.1039/d2cp00646d] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Solvate ionic liquids (ILs) are promising candidates for several applications due to their stability, high coulombic efficiency, and low volatility. In this work, we investigate the solvation of lithium-bistriflimide by different glycerol-derived triether solvents, using molecular dynamics simulations. Very strong interactions between Li+ and the solvent oxygen sites are found, leading to significant conformational changes in the solvent. By comparing the conformation of the neat solvents with their IL mixtures at different concentrations and temperatures, we find that the presence of Li+ induces a distinct crown-like structure in the solvent molecules. The Li+ cations and the surrounding solvent form a podand complex, which is stable even at elevated temperatures. These glycerol-derived solvents exhibit distinct interactions with Li+ cations which may be exploited in electrolytic applications or lithium recovery processes.
Collapse
Affiliation(s)
- Gabriel D Barbosa
- Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, AL 35487, USA.
| | - Jason E Bara
- Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, AL 35487, USA.
| | - C Heath Turner
- Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, AL 35487, USA.
| |
Collapse
|
35
|
Lin C, Wei H, Li H, Duan X. Structures of cationic and anionic polyelectrolytes in aqueous solutions: the sign effect. SOFT MATTER 2022; 18:1603-1616. [PMID: 35080232 DOI: 10.1039/d1sm01700d] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
In this study, we use molecular dynamics simulation to explore the structures of anionic and cationic polyelectrolytes in aqueous solutions. We first confirm the significantly stronger solvation effects of single anions compared to cations in water at the fixed ion radii, due to the reversal orientations of asymmetric dipolar H2O molecules around the ions. Based on this, we demonstrate that the solvation discrepancy of cations/anions and electrostatic correlations of ionic species can synergistically cause the nontrivial structural difference between single anionic and cationic polyelectrolytes. The cationic polyelectrolyte shows an extended structure whereas the anionic polyelectrolyte exhibits a collapsed structure, and their structural differences decline with increasing the counterion size. Furthermore, we corroborate that multiple cationic polyelectrolytes or multiple anionic polyelectrolytes can exhibit largely differential molecular architectures in aqueous solutions. In the solvation dominant regime, the polyelectrolyte solutions exhibit uniform structures; whereas, in the electrostatic correlation dominant regime, the polyelectrolyte solutions exhibit heterogeneous structures, in which the likely charged chains microscopically aggregate through counterion condensations. Increasing the intrinsic chain rigidity causes polyelectrolyte extension and hence moderately weakens the inter-chain clustering. Our work highlights the various, unique structures and molecular architectures of polyelectrolytes in solutions caused by the multi-body correlations between polyelectrolytes, counterions and asymmetric dipolar solvent molecules, which provides insights into the fundamental understanding of ion-containing polymers.
Collapse
Affiliation(s)
- Chengjiang Lin
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin 130022, China.
- School of Applied Chemistry and Engineering, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Hao Wei
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin 130022, China.
| | - Hongfei Li
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin 130022, China.
- School of Applied Chemistry and Engineering, University of Science and Technology of China, Hefei 230026, P. R. China
| | - Xiaozheng Duan
- State Key Laboratory of Polymer Physics and Chemistry, Changchun Institute of Applied Chemistry, Chinese Academy of Sciences, Changchun, Jilin 130022, China.
| |
Collapse
|
36
|
Thurston BA, Grest GS, Stevens MJ. Overlap Concentration of Sodium Polystyrene Sulfonate in Solution. ACS Macro Lett 2022; 11:217-222. [PMID: 35574772 DOI: 10.1021/acsmacrolett.1c00649] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The overlap concentration c* of sodium polystyrene sulfonate in water is calculated using multichain atomistic and coarse grained (CG) simulations for a range of chain lengths. Fully atomistic molecular dynamics simulations are carried out for N = 32-192 monomers. The CG model was parameterized to match the end-to-end distance from the atomistic simulations at small N and allows us to simulate a much larger N. Treating the hydrophobic backbone by inclusion of attraction between monomers is an essential addition to the CG model. The simulation c* are in agreement with experimental data, yet at c*, the chains are not fully stretched, even for N as large as 1200. This implies that none of the experimental systems are in the scaling regime and to reach the scaling regime for NaPSS chains much longer than N = 1200 are required.
Collapse
Affiliation(s)
- Bryce A. Thurston
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Gary S. Grest
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| | - Mark J. Stevens
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, New Mexico 87185, United States
| |
Collapse
|
37
|
Zhang W, Yamaguchi T, Fang C, Yoshida K, Zhou Y, Zhu F, Machida S, Hattori T, Li W. Structure of an aqueous RbCl solution in the gigapascal pressure range by neutron diffraction combined with empirical potential structure refinement modeling. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2021.118080] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
38
|
Atwi R, Chen Y, Han KS, Mueller KT, Murugesan V, Rajput NN. An automated framework for high-throughput predictions of NMR chemical shifts within liquid solutions. NATURE COMPUTATIONAL SCIENCE 2022; 2:112-122. [PMID: 38177518 DOI: 10.1038/s43588-022-00200-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/10/2021] [Accepted: 01/20/2022] [Indexed: 01/06/2024]
Abstract
Identifying stable speciation in multi-component liquid solutions is fundamentally important to areas from electrochemistry to organic chemistry and biomolecular systems. Here we introduce a fully automated, high-throughput computational framework for the accurate prediction of stable species in liquid solutions by computing the nuclear magnetic resonance (NMR) chemical shifts. The framework automatically extracts and categorizes hundreds of thousands of atomic clusters from classical molecular dynamics simulations, identifies the most stable species in solution and calculates their NMR chemical shifts via density functional theory calculations. Additionally, the framework creates a database of computed chemical shifts for liquid solutions across a wide chemical and parameter space. We compare our computational results to experimental measurements for magnesium bis(trifluoromethanesulfonyl)imide Mg(TFSI)2 salt in dimethoxyethane solvent. Our analysis of the Mg2+ solvation structural evolutions reveals key factors that influence the accuracy of NMR chemical shift predictions in liquid solutions. Furthermore, we show how the framework reduces the performance of over 300 13C and 600 1H density functional theory chemical shift predictions to a single submission procedure.
Collapse
Affiliation(s)
- Rasha Atwi
- Department of Materials Science and Chemical Engineering, Stony Brook University, Stony Brook, NY, USA
| | - Ying Chen
- The Joint Center for Energy Storage Research (JCESR), Pacific Northwest National Laboratory, Richland, WA, USA
| | - Kee Sung Han
- The Joint Center for Energy Storage Research (JCESR), Pacific Northwest National Laboratory, Richland, WA, USA
| | - Karl T Mueller
- The Joint Center for Energy Storage Research (JCESR), Pacific Northwest National Laboratory, Richland, WA, USA
| | - Vijayakumar Murugesan
- The Joint Center for Energy Storage Research (JCESR), Pacific Northwest National Laboratory, Richland, WA, USA
| | - Nav Nidhi Rajput
- Department of Materials Science and Chemical Engineering, Stony Brook University, Stony Brook, NY, USA.
| |
Collapse
|
39
|
Chen W, Gu J, Liu Q, Yang M, Zhan C, Zang X, Pham TA, Liu G, Zhang W, Zhang D, Dunn B, Morris Wang Y. Two-dimensional quantum-sheet films with sub-1.2 nm channels for ultrahigh-rate electrochemical capacitance. NATURE NANOTECHNOLOGY 2022; 17:153-158. [PMID: 34795438 DOI: 10.1038/s41565-021-01020-0] [Citation(s) in RCA: 23] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/09/2021] [Accepted: 09/29/2021] [Indexed: 06/13/2023]
Abstract
Dense, thick, but fast-ion-conductive electrodes are critical yet challenging components of ultrafast electrochemical capacitors with high volumetric power/energy densities1-4. Here we report an exfoliation-fragmentation-restacking strategy towards thickness-adjustable (1.5‒24.0 μm) dense electrode films of restacked two-dimensional 1T-MoS2 quantum sheets. These films bear the unique architecture of an exceptionally high density of narrow (sub-1.2 nm) and ultrashort (~6.1 nm) hydrophobic nanochannels for confinement ion transport. Among them, 14-μm-thick films tested at 2,000 mV s-1 can deliver not only a high areal capacitance of 0.63 F cm-2 but also a volumetric capacitance of 437 F cm-3 that is one order of magnitude higher than that of other electrodes. Density functional theory and ab initio molecular dynamics simulations suggest that both hydration and nanoscale channels play crucial roles in enabling ultrafast ion transport and enhanced charge storage. This work provides a versatile strategy for generating rapid ion transport channels in thick but dense films for energy storage and filtration applications.
Collapse
Affiliation(s)
- Wenshu Chen
- State Key Laboratory of Metal Matrix Composites, Shanghai Jiao Tong University, Shanghai, China
- School of Environmental Science and Nanjing Key Laboratory of Advanced Functional Materials, Nanjing Xiaozhuang University, Nanjing, China
| | - Jiajun Gu
- State Key Laboratory of Metal Matrix Composites, Shanghai Jiao Tong University, Shanghai, China.
| | - Qinglei Liu
- State Key Laboratory of Metal Matrix Composites, Shanghai Jiao Tong University, Shanghai, China.
| | - Mengzhao Yang
- State Key Laboratory of Metal Matrix Composites, Shanghai Jiao Tong University, Shanghai, China
| | - Cheng Zhan
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Xining Zang
- Department of Materials Science and Engineering, Massachusetts Institute of Technology, Cambridge, MA, USA
| | - Tuan Anh Pham
- Physical and Life Sciences Directorate, Lawrence Livermore National Laboratory, Livermore, CA, USA
| | - Guangxiang Liu
- School of Environmental Science and Nanjing Key Laboratory of Advanced Functional Materials, Nanjing Xiaozhuang University, Nanjing, China
| | - Wang Zhang
- State Key Laboratory of Metal Matrix Composites, Shanghai Jiao Tong University, Shanghai, China
| | - Di Zhang
- State Key Laboratory of Metal Matrix Composites, Shanghai Jiao Tong University, Shanghai, China.
| | - Bruce Dunn
- Department of Materials Science and Engineering, University of California, Los Angeles, Los Angeles, CA, USA
| | - Y Morris Wang
- Department of Materials Science and Engineering, University of California, Los Angeles, Los Angeles, CA, USA.
| |
Collapse
|
40
|
Blazquez S, Conde MM, Abascal JLF, Vega C. The Madrid-2019 force field for electrolytes in water using TIP4P/2005 and scaled charges: Extension to the ions F−, Br−, I−, Rb+, and Cs+. J Chem Phys 2022; 156:044505. [DOI: 10.1063/5.0077716] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- S. Blazquez
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - M. M. Conde
- Departamento de Ingeniería Química Industrial y Medio Ambiente, Escuela Técnica Superior de Ingenieros Industriales, Universidad Politécnica de Madrid, 28006 Madrid, Spain
| | - J. L. F. Abascal
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - C. Vega
- Departamento Química Física I, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| |
Collapse
|
41
|
Tolmachev D, Lukasheva N, Ramazanov R, Nazarychev V, Borzdun N, Volgin I, Andreeva M, Glova A, Melnikova S, Dobrovskiy A, Silber SA, Larin S, de Souza RM, Ribeiro MCC, Lyulin S, Karttunen M. Computer Simulations of Deep Eutectic Solvents: Challenges, Solutions, and Perspectives. Int J Mol Sci 2022; 23:645. [PMID: 35054840 PMCID: PMC8775846 DOI: 10.3390/ijms23020645] [Citation(s) in RCA: 28] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/04/2021] [Revised: 01/02/2022] [Accepted: 01/04/2022] [Indexed: 12/13/2022] Open
Abstract
Deep eutectic solvents (DESs) are one of the most rapidly evolving types of solvents, appearing in a broad range of applications, such as nanotechnology, electrochemistry, biomass transformation, pharmaceuticals, membrane technology, biocomposite development, modern 3D-printing, and many others. The range of their applicability continues to expand, which demands the development of new DESs with improved properties. To do so requires an understanding of the fundamental relationship between the structure and properties of DESs. Computer simulation and machine learning techniques provide a fruitful approach as they can predict and reveal physical mechanisms and readily be linked to experiments. This review is devoted to the computational research of DESs and describes technical features of DES simulations and the corresponding perspectives on various DES applications. The aim is to demonstrate the current frontiers of computational research of DESs and discuss future perspectives.
Collapse
Affiliation(s)
- Dmitry Tolmachev
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Natalia Lukasheva
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Ruslan Ramazanov
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Victor Nazarychev
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Natalia Borzdun
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Igor Volgin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Maria Andreeva
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Artyom Glova
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Sofia Melnikova
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Alexey Dobrovskiy
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Steven A. Silber
- Department of Physics and Astronomy, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada;
- The Centre of Advanced Materials and Biomaterials Research, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada
| | - Sergey Larin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Rafael Maglia de Souza
- Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, Avenida Professor Lineu Prestes 748, São Paulo 05508-070, Brazil; (R.M.d.S.); (M.C.C.R.)
| | - Mauro Carlos Costa Ribeiro
- Departamento de Química Fundamental, Instituto de Química, Universidade de São Paulo, Avenida Professor Lineu Prestes 748, São Paulo 05508-070, Brazil; (R.M.d.S.); (M.C.C.R.)
| | - Sergey Lyulin
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
| | - Mikko Karttunen
- Institute of Macromolecular Compounds, Russian Academy of Sciences, Bolshoy pr. 31, 199004 St. Petersburg, Russia; (N.L.); (R.R.); (V.N.); (N.B.); (I.V.); (M.A.); (A.G.); (S.M.); (A.D.); (S.L.); (S.L.)
- Department of Physics and Astronomy, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada;
- The Centre of Advanced Materials and Biomaterials Research, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada
- Department of Chemistry, The University of Western Ontario, 1151 Richmond Street, London, ON N6A 5B7, Canada
| |
Collapse
|
42
|
Hou LP, Zhang XQ, Yao N, Chen X, Li BQ, Shi P, Jin CB, Huang JQ, Zhang Q. An encapsulating lithium-polysulfide electrolyte for practical lithium–sulfur batteries. Chem 2022. [DOI: 10.1016/j.chempr.2021.12.023] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
43
|
Qiu Y, Jiang Y, Zhang Y, Zhang H. Rational Design of Nonbonded Point Charge Models for Monovalent Ions with Lennard-Jones 12-6 Potential. J Phys Chem B 2021; 125:13502-13518. [PMID: 34860517 DOI: 10.1021/acs.jpcb.1c09103] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ions are of central importance in nature, and a variety of potential models was proposed to model ions in different phases for an in-depth exploration of ion-related systems. Here, we developed point charge models of 14 monovalent ions with the traditional 12-6 Lennard-Jones (LJ) potential for use in conjunction with 11 water models of TIP3P, OPC3, SPC/E, SPC/Eb, TIP3P-FB, a99SB-disp, TIP4P-Ew, OPC, TIP4P/2005, TIP4P-D, and TIP4P-FB. The designed models reproduced the real hydration free energy (HFE) of ions and the ion-oxygen distance (IOD) in the first hydration shell accurately and simultaneously, a performance similar to the previously reported 12-6-4 LJ-type ion models (12-6 LJ plus an attractive C4 term for cations or a repulsive one for anions). This work, along with our previous work on di-, tri-, and tetravalent metal cations (J. Chem. Inf. Model. 2021, 61, 4031-4044; J. Chem. Inf. Model. 2021, 61, 4613-4629), demonstrates the feasibility of the simple 12-6 LJ potential in ion modeling. In order for the 12-6 LJ potential to reproduce both the HFE and IOD, the LJ R parameters need to be close to Shannon's ionic radii for the highly charged cations and to the Stokes's van der Waals (vdW) radii for the monovalent ions. With an additional C4 term, the R parameters of 12-6-4 LJ ion models agree well with the Stokes's vdW radii and have a more physical meaning. It appears that the C4 term can be merged into the 12-6 LJ potential by a rational tuning of R and the LJ well depth. Simulations of the osmotic coefficients of alkali chloride solutions and the properties of gaseous and solid alkali halides indicate the necessity of further optimizing ion-ion interactions via, for instance, targeting more properties or using a more physical (polarizable) model.
Collapse
Affiliation(s)
- Yejie Qiu
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park, Pennsylvania 16802, United States
| | - Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, 100083 Beijing, China
| |
Collapse
|
44
|
Mukherjee K, Das S, Rajbangshi J, Tarif E, Barman A, Biswas R. Temperature-Dependent Dielectric Relaxation in Ionic Acetamide Deep Eutectics: Partial Viscosity Decoupling and Explanations from the Simulated Single-Particle Reorientation Dynamics and Hydrogen-Bond Fluctuations. J Phys Chem B 2021; 125:12552-12567. [PMID: 34752087 DOI: 10.1021/acs.jpcb.1c07299] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report here temperature-dependent (293 ≤ T (K) ≤ 336) dielectric relaxation (DR) measurements of (acetamide + LiBr/NO3-/ClO4-) deep eutectic solvents (DESs) in the frequency window of 0.2 ≤ ν (GHz) ≤ 50 and explore, via molecular dynamics simulations, the relative roles for the collective single-particle reorientational relaxations and the H-bond dynamics of acetamide in the measured DR response. In addition, DR measurements of neat molten acetamide were performed. Recorded DR spectra of these DESs require multi-Debye fits and produce well-separated DR time scales that are spread over several picoseconds to ∼1 ns. Simulations suggest DR time scales derive contributions from both the collective reorientational (Cl(t)) relaxation and structural H-bond (CHB(t)) dynamics of acetamide. A good correlation between the measured and simulated activation energies further reveals a strong connection between the measured DR and the simulated Cl(t) and CHB(t). Average DR times exhibit a strong fractional viscosity dependence, suggesting substantial microheterogeneity in these media. Simulations of Cl(t) and CHB(t) reveal strong stretched exponential relaxations with a stretching exponent, 0.4 ≤ β ≤ 0.7. The ratio between the average reorientational correlation times of first and second ranks, ⟨τ⟩l=1/⟨τ⟩l=2, deviates appreciably from Debye's l(l+1) law for homogeneous media. Importantly, a pronounced translation-rotation decoupling between the simulated reorientation and center-of-mass diffusion times was observed.
Collapse
|
45
|
Hou T, Fong KD, Wang J, Persson KA. The solvation structure, transport properties and reduction behavior of carbonate-based electrolytes of lithium-ion batteries. Chem Sci 2021; 12:14740-14751. [PMID: 34820089 PMCID: PMC8597828 DOI: 10.1039/d1sc04265c] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/03/2021] [Accepted: 09/16/2021] [Indexed: 11/21/2022] Open
Abstract
Despite the extensive employment of binary/ternary mixed-carbonate electrolytes (MCEs) for Li-ion batteries, the role of each ingredient with regards to the solvation structure, transport properties, and reduction behavior is not fully understood. Herein, we report the atomistic modeling and transport property measurements of the Gen2 (1.2 M LiPF6 in ethylene carbonate (EC) and ethyl methyl carbonate (EMC)) and EC-base (1.2 M LiPF6 in EC) electrolytes, as well as their mixtures with 10 mol% fluoroethylene carbonate (FEC). Due to the mixing of cyclic and linear carbonates, the Gen2 electrolyte is found to have a 60% lower ion dissociation rate and a 44% faster Li+ self-diffusion rate than the EC-base electrolyte, while the total ionic conductivities are similar. Moreover, we propose for the first time the anion–solvent exchange mechanism in MCEs with identified energetic and electrostatic origins. For electrolytes with additive, up to 25% FEC coordinates with Li+, which exhibits a preferential reduction that helps passivate the anode and facilitates an improved solid electrolyte interphase. The work provides a coherent computational framework for evaluating mixed electrolyte systems. The different roles of the anion, cyclic and linear carbonates, and additive in mixed-carbonate electrolytes are revealed. The anion–solvent exchange mechanism and factors influencing the solid electrolyte interphase (SEI) formation are deciphered.![]()
Collapse
Affiliation(s)
- Tingzheng Hou
- Department of Materials Science and Engineering, University of California Berkeley 210 Hearst Mining Building Berkeley California 94720 USA.,Energy Technologies Area, Lawrence Berkeley National Laboratory Berkeley California 94720 USA
| | - Kara D Fong
- Energy Technologies Area, Lawrence Berkeley National Laboratory Berkeley California 94720 USA.,Department of Chemical and Biomolecular Engineering, University of California Berkeley CA 94720 USA
| | - Jingyang Wang
- Department of Materials Science and Engineering, University of California Berkeley 210 Hearst Mining Building Berkeley California 94720 USA.,Materials Sciences Division, Lawrence Berkeley National Laboratory Berkeley California 94720 USA
| | - Kristin A Persson
- Department of Materials Science and Engineering, University of California Berkeley 210 Hearst Mining Building Berkeley California 94720 USA.,The Molecular Foundry, Lawrence Berkeley National Laboratory Berkeley California 94720 USA
| |
Collapse
|
46
|
Kubisiak P, Wróbel P, Eilmes A. How Temperature, Pressure, and Salt Concentration Affect Correlations in LiTFSI/EMIM-TFSI Electrolytes: A Molecular Dynamics Study. J Phys Chem B 2021; 125:12292-12302. [PMID: 34706539 PMCID: PMC8591607 DOI: 10.1021/acs.jpcb.1c07782] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2021] [Revised: 10/13/2021] [Indexed: 11/29/2022]
Abstract
Classical polarizable molecular dynamics simulations have been performed for LiTFSI solutions in the EMIM-TFSI ionic liquid. Different temperature or pressure values and salt concentrations have been examined. The structure and dynamics of the solvation shell of Li+ cations, diffusion coefficients of ions, conductivities of the electrolytes, and correlations between motions of ions have been analyzed. The results indicated that regardless of the conditions, significant correlations are present in all systems. The degree of correlations depends mainly on the salt fraction in the electrolyte and is much less affected by temperature and pressure changes. A positive correlation between motions of Li+ cations and TFSI anions, leading to the occurrence of negative Li+ transference numbers, exists for all conditions, although temperature and pressure changes affect the speed of anion exchange in Li+ solvation shells.
Collapse
Affiliation(s)
- Piotr Kubisiak
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Kraków, Poland
| | - Piotr Wróbel
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Kraków, Poland
| | - Andrzej Eilmes
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Kraków, Poland
| |
Collapse
|
47
|
Man VH, He X, Gao J, Wang J. Effects of All-Atom Molecular Mechanics Force Fields on Amyloid Peptide Assembly: The Case of PHF6 Peptide of Tau Protein. J Chem Theory Comput 2021; 17:6458-6471. [PMID: 34491058 DOI: 10.1021/acs.jctc.1c00028] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Molecular dynamics (MD) simulations play a vital role in revealing the mechanism of amyloid aggregation that is crucial to the therapeutic agent development for Alzheimer's Disease. However, the accuracy of MD simulation results strongly depends on the force field employed. In our previous benchmark for 17 all-atom force fields on modeling of amyloid aggregation using the Aβ16-22 dimer, we showed that AMBER14SB and CHARMM36m are suitable force fields for amyloid aggregation simulation, while GROMOS54a7 and OPLSAA are not good for the task. In this work, we continue assessing the applicability of atomistic force fields on amyloid aggregation using the VQIVYK (PHF6) peptide which is essential for tau-protein aggregation. Although, both Aβ16-22 and PHF6 peptides formed fibrils in vitro, the PHF6 fibrils are parallel β-sheets, while the Aβ16-22 fibrils are antiparallel β-sheets. We performed an all-atom large-scale MD simulation in explicit water on the PHF6 dimer and octa-peptides systems using five mainstream force fields, including AMBER99SB-disp, AMBER14SB, CHARMM36m, GROMOS54a7, and OPLSAA. The accumulated simulation time is 0.2 ms. Our result showed that the β-sheet structures of PHF6 peptides sampled by AMBER99SB-disp, AMBER14SB, GROMOS54a7, and OPLSAA are in favor of the antiparallel β-sheets, while the dominant type of β-sheet structures is parallel β-sheet by using CHARMM36m. Among the five force fields, CHARMM36m provides the strongest CH-π interaction that was observed in an NMR study. The comparison between our results and experimental observation indicates that CHARMM36m achieved the best performance on modeling the aggregation of PHF6 peptides. In summary, CHARMM36m is currently the most suitable force field for studying the aggregation of both amyloid-β and Tau through MD simulations.
Collapse
Affiliation(s)
- Viet Hoang Man
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Xibing He
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Jie Gao
- Department of Neuroscience, The Ohio State University Wexner Medical Center, Columbus, Ohio 43210, United States
| | - Junmei Wang
- Department of Pharmaceutical Sciences and Computational Chemical Genomics Screening Center, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| |
Collapse
|
48
|
Credidio B, Pugini M, Malerz S, Trinter F, Hergenhahn U, Wilkinson I, Thürmer S, Winter B. Quantitative electronic structure and work-function changes of liquid water induced by solute. Phys Chem Chem Phys 2021; 24:1310-1325. [PMID: 34604895 PMCID: PMC8768487 DOI: 10.1039/d1cp03165a] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Recent advancement in quantitative liquid-jet photoelectron spectroscopy enables the accurate determination of the absolute-scale electronic energetics of liquids and species in solution. The major objective of the present work is the determination of the absolute lowest-ionization energy of liquid water, corresponding to the 1b1 orbital electron liberation, which is found to vary upon solute addition, and depends on the solute concentration. We discuss two prototypical aqueous salt solutions, NaI(aq) and tetrabutylammonium iodide, TBAI(aq), with the latter being a strong surfactant. Our results reveal considerably different behavior of the liquid water 1b1 binding energy in each case. In the NaI(aq) solutions, the 1b1 energy increases by about 0.3 eV upon increasing the salt concentration from very dilute to near-saturation concentrations, whereas for TBAI the energy decreases by about 0.7 eV upon formation of a TBAI surface layer. The photoelectron spectra also allow us to quantify the solute-induced effects on the solute binding energies, as inferred from concentration-dependent energy shifts of the I− 5p binding energy. For NaI(aq), an almost identical I− 5p shift is found as for the water 1b1 binding energy, with a larger shift occurring in the opposite direction for the TBAI(aq) solution. We show that the evolution of the water 1b1 energy in the NaI(aq) solutions can be primarily assigned to a change of water's electronic structure in the solution bulk. In contrast, apparent changes of the 1b1 energy for TBAI(aq) solutions can be related to changes of the solution work function which could arise from surface molecular dipoles. Furthermore, for both of the solutions studied here, the measured water 1b1 binding energies can be correlated with the extensive solution molecular structure changes occurring at high salt concentrations, where in the case of NaI(aq), too few water molecules exist to hydrate individual ions and the solution adopts a crystalline-like phase. We also comment on the concentration-dependent shape of the second, 3a1 orbital liquid water ionization feature which is a sensitive signature of water–water hydrogen bond interactions. Significant differences are observed in liquid-water's lowest electron binding energy with increasing solute concentration in archetypal aqueous solutions. For NaI(aq) and TBAI(aq), the energy changes extend to +0.3 eV and −0.7 eV, respectively.![]()
Collapse
Affiliation(s)
- Bruno Credidio
- Molecular Physics Department, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany. .,Institute for Chemical Sciences and Engineering (ISIC), École Polytechnique Fédérale de Lausanne (EPFL), 1015 Lausanne, Switzerland
| | - Michele Pugini
- Molecular Physics Department, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany.
| | - Sebastian Malerz
- Molecular Physics Department, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany.
| | - Florian Trinter
- Molecular Physics Department, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany. .,Institut für Kernphysik, Goethe-Universität, Max-von-Laue-Straße 1, 60438 Frankfurt am Main, Germany
| | - Uwe Hergenhahn
- Molecular Physics Department, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany.
| | - Iain Wilkinson
- Department of Locally-Sensitive & Time-Resolved Spectroscopy, Helmholtz-Zentrum Berlin für Materialien und Energie, Hahn-Meitner-Platz 1, 14109 Berlin, Germany
| | - Stephan Thürmer
- Department of Chemistry, Graduate School of Science, Kyoto University, Kitashirakawa-Oiwakecho, Sakyo-Ku, Kyoto 606-8502, Japan.
| | - Bernd Winter
- Molecular Physics Department, Fritz-Haber-Institut der Max-Planck-Gesellschaft, Faradayweg 4-6, 14195 Berlin, Germany.
| |
Collapse
|
49
|
Wróbel P, Kubisiak P, Eilmes A. NaFSI and NaTFSI Solutions in Ether Solvents from Monoglyme to Poly(ethylene oxide)-A Molecular Dynamics Study. J Phys Chem B 2021; 125:10293-10303. [PMID: 34494439 PMCID: PMC8450900 DOI: 10.1021/acs.jpcb.1c05793] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Classical molecular dynamics simulations have been performed for a series of electrolytes based on sodium bis(fluorosulfonyl)imide or sodium bis(trifluoromethylsulfonyl)imide salts and monoglyme, tetraglyme, and poly(ethylene oxide) as solvents. Structural properties have been assessed through the analysis of coordination numbers and binding patterns. Residence times for Na-O interactions have been used to investigate the stability of solvation shells. Diffusion coefficients of ions and electrical conductivity of the electrolytes have been estimated from molecular dynamics trajectories. Contributions to the total conductivity have been analyzed in order to investigate the role of ion-ion correlations. It has been found that the anion-cation interactions are more probable in the systems with NaTFSI salts. Accordingly, the degree of correlations between ion motions is larger in NaTFSI-based electrolytes.
Collapse
Affiliation(s)
- Piotr Wróbel
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Kraków, Poland
| | - Piotr Kubisiak
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Kraków, Poland
| | - Andrzej Eilmes
- Faculty of Chemistry, Jagiellonian University, Gronostajowa 2, 30-387 Kraków, Poland
| |
Collapse
|
50
|
Shpigel N, Chakraborty A, Malchik F, Bergman G, Nimkar A, Gavriel B, Turgeman M, Hong CN, Lukatskaya MR, Levi MD, Gogotsi Y, Major DT, Aurbach D. Can Anions Be Inserted into MXene? J Am Chem Soc 2021; 143:12552-12559. [PMID: 34357752 DOI: 10.1021/jacs.1c03840] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Despite the continuous progress in the research and development of Ti3C2Tx (MXene) electrodes for high-power batteries and supercapacitor applications, the role of the anions in the electrochemical energy storage and their ability to intercalate between the MXene sheets upon application of positive voltage have not been clarified. A decade after the discovery of MXenes, the information about the possibility of anion insertion into the restacked MXene electrode is still being questioned. Since the positive potential stability range in diluted aqueous electrolytes is severely limited by anodic oxidation of the Ti, the possibility of anion insertion was evaluated in concentrated aqueous electrolyte solutions and aprotic electrolytes as well. To address this issue, we have conducted in situ gravimetric electrochemical quartz crystal microbalance with dissipation monitoring (EQCM-D) measurements in highly concentrated LiCl and LiBr electrolytes, which enable a significant extension of the operation range of the MXene electrodes toward positive potentials. Also, halogens are among the smallest anions and should be easier to intercalate between MXene layers, in comparison to multiatomic anions. On the basis of mass change variations in the positive voltage range and complementary density functional theory calculations, it was demonstrated that insertion of anionic species into MXene, within the range of potentials of interest for capacitive energy storage, is not likely to occur. This can be explained by the strong negative charge on Ti3C2Tx sheets terminated by functional groups.
Collapse
Affiliation(s)
- Netanel Shpigel
- Department of Chemistry and BINA-BIU Center for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Arup Chakraborty
- Department of Chemistry and BINA-BIU Center for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Fyodor Malchik
- Center for Physical and Chemical Methods of Research and Analysis, al-Farabi Kazakh National University, 050040 Almaty, Kazakhstan
| | - Gil Bergman
- Department of Chemistry and BINA-BIU Center for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Amey Nimkar
- Department of Chemistry and BINA-BIU Center for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Bar Gavriel
- Department of Chemistry and BINA-BIU Center for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Meital Turgeman
- Department of Chemistry and BINA-BIU Center for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Chulgi Nathan Hong
- Department of Mechanical and Process Engineering, ETH Zurich, Sonneggstrasse 3, Zürich 8092 Switzerland
| | - Maria R Lukatskaya
- Department of Mechanical and Process Engineering, ETH Zurich, Sonneggstrasse 3, Zürich 8092 Switzerland
| | - Mikhael D Levi
- Department of Chemistry and BINA-BIU Center for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Yury Gogotsi
- Department of Materials Science and Engineering, and A.J. Drexel Nanomaterials Institute, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | - Dan T Major
- Department of Chemistry and BINA-BIU Center for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| | - Doron Aurbach
- Department of Chemistry and BINA-BIU Center for Nanotechnology and Advanced Materials, Bar-Ilan University, Ramat-Gan 5290002, Israel
| |
Collapse
|