Charge and electric field fluctuations in aqueous NaCl electrolytes

Water Dimer under Electric Fields: An Ab Initio Investigation up to Quantum Accuracy

It is well-established that strong electric fields (EFs) can align water dipoles, partially order the H-bond network of liquid water, and induce water splitting and proton transfers. To illuminate the fundamental behavior of water under external EFs, we present the first benchmark, to the best of our knowledge, of DFT calculations of the water dimer exposed to intense EFs against coupled cluster calculations. The analyses of the vibrational Stark effect and electron density provide a consistent picture of the intermolecular charge transfer effects driven along the H-bond by the increasing applied field at all theory levels. However, our findings prove that at extreme field regimes (∼1-2 V/Å) DFT calculations significantly exaggerate by ∼10-30% the field-induced strengthening of the H-bond, both within the GGA, hybrid GGA, and hybrid meta-GGA approximations. Notably, a linear correlation emerges between the vibrational Stark effect on OH stretching and H-bond strengthening: a 1 kcal mol-1 increase corresponds to an 80 cm-1 red-shift in OH stretching frequency.

Scaling perspectives of underscreening in concentrated electrolyte solutions

We present a scaling view of underscreening observed in salt solutions in the range of concentrations greater than about 1 M, in which the screening length increases with concentration. The system consists of hydrated clusters of positive and negative ions with a single unpaired ion as suggested by recent simulations. The environment of this ion is more hydrated than average which leads to a self-similar situation in which the size of this environment scales with the screening length. The prefactor involves the local dielectric constant and the cluster density. The scaling arguments as well as the cluster model lead to scaling of the screening length with the ion concentration, in agreement with observations.

Electrical noise in electrolytes: a theoretical perspective

Seemingly unrelated experiments such as electrolyte transport through nanotubes, nano-scale electrochemistry, NMR relaxometry and surface force balance measurements, all probe electrical fluctuations: of the electric current, the charge and polarization, the field gradient (for quadrupolar nuclei) and the coupled mass/charge densities. The fluctuations of such various observables arise from the same underlying microscopic dynamics of the ions and solvent molecules. In principle, the relevant length and time scales of these dynamics are encoded in the dynamic structure factors. However, modelling the latter for frequencies and wavevectors spanning many orders of magnitude remains a great challenge to interpret the experiments in terms of physical processes such as solvation dynamics, diffusion, electrostatic and hydrodynamic interactions between ions, interactions with solid surfaces, etc. Here, we highlight the central role of the charge-charge dynamic structure factor in the fluctuations of electrical observables in electrolytes and offer a unifying perspective over a variety of complementary experiments. We further analyze this quantity in the special case of an aqueous NaCl electrolyte, using simulations with explicit ions and an explicit or implicit solvent. We discuss the ability of the standard Poisson-Nernst-Planck theory to capture the simulation results, and how the predictions can be improved. We finally discuss the contributions of ions and water to the total charge fluctuations. This work illustrates an ongoing effort towards a comprehensive understanding of electrical fluctuations in bulk and confined electrolytes, in order to enable experimentalists to decipher the microscopic properties encoded in the measured electrical noise.

Molecular Simulations of Aqueous Electrolytes: Role of Explicit Inclusion of Charge Transfer into Force Fields

We describe here simulations of aqueous salt solutions that are performed using an explicit charge transfer force field. The emphasis of the discussion is on the calculation of a dynamical property of the solutions: self-diffusion of water. While force fields that are based on pairwise additive potentials or on potentials with explicit inclusion of polarization or with scaled charges can provide at best a qualitative agreement with experiments, force fields with explicit inclusion of charge transfer can produce quantitative agreement with experiment for NaCl and KCl solutions. We argue that a force field with explicit charge transfer contains new physics absent in the previously used force fields described in recent reviews of molecular simulations of aqueous electrolytes.

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion-solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. Here we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

Electric fields and potentials in condensed phases

The electric fields and potentials inside and at the interface of matter are relevant to many branches of physics, chemistry, and biology. Accurate quantification of these fields and/or potentials is essential to control and exploit chemical and physical transformations. Before we understand the response of matter to external fields, it is first important to understand the intrinsic interior and interfacial fields and potentials, both classically and quantum mechanically, as well as how they are probed experimentally. Here we compare and contrast, beginning with the hydrogen atom in vacuum and ending with concentrated aqueous NaCl electrolyte, both classical and quantum mechanical electric potentials and fields. We make contact with experimental vibrational Stark, electrochemical, X-ray, and electron spectroscopic probes of these potentials and fields, outline relevant conceptual difficulties, and underscore the advantage of electron holography as a basis to better understand electrostatics in matter.

Modeling of the Response of Hydrogen Bond Properties on an External Electric Field: Geometry, NMR Chemical Shift, Spin-Spin Scalar Coupling

The response of the geometric and NMR properties of molecular systems to an external electric field has been studied theoretically in a wide field range. It has been shown that this adduct under field approach can be used to model the geometric and spectral changes experienced by molecular systems in polar media if the system in question has one and only one bond, the polarizability of which significantly exceeds the polarizability of other bonds. If this requirement is met, then it becomes possible to model even extreme cases, for example, proton dissociation in hydrogen halides. This requirement is fulfilled for many complexes with one hydrogen bond. For such complexes, this approach can be used to facilitate a detailed analysis of spectral changes associated with geometric changes in the hydrogen bond. For example, in hydrogen-bonded complexes of isocyanide C≡¹⁵N-¹H⋯X, ¹J(¹⁵N¹H) depends exclusively on the N-H distance, while δ(¹⁵N) is also slightly influenced by the nature of X.

Electric Field and Temperature Effects on the Ab Initio Spectroscopy of Liquid Methanol

Although many H-bonded systems have been extensively investigated by means of infrared (IR) spectroscopy, the vibrational response to externally applied electric fields of polar liquids remains poorly investigated. However, local electric fields along with quantum-mechanical interactions rule the behavior of H-bonded samples at the molecular level. Among the many H-bonded systems, liquid methanol holds a key place in that it exhibits a very simple H-bond network where, on average, each molecule acts as a single H-bond donor and, at the same time, as a single H-bond acceptor. Here we report on the IR spectra emerging from a series of state-of-the-art ab initio molecular dynamics simulations of bulk liquid methanol under the action of static and homogeneous electric fields. In addition, the same analysis is here conducted in the absence of the external field and for different temperatures. Although some electric-field-induced effects resemble the response of other polar liquids (such as the global contraction of the IR spectrum upon field exposure), it turns out that, distinctly from water, the “electrofreezing” phenomenon is unlikely to happen in liquid methanol. Finally, we provide atomistic analyses magnifying the completely different nature of electric-field- and temperature-induced effects on bulk liquid methanol and on its vibrational response.

1,3,5-Triaza-7-Phosphaadamantane (PTA) as a 31P NMR Probe for Organometallic Transition Metal Complexes in Solution

Due to the rigid structure of 1,3,5-triaza-7-phosphaadamantane (PTA), its ³¹P chemical shift solely depends on non-covalent interactions in which the molecule is involved. The maximum range of change caused by the most common of these, hydrogen bonding, is only 6 ppm, because the active site is one of the PTA nitrogen atoms. In contrast, when the PTA phosphorus atom is coordinated to a metal, the range of change exceeds 100 ppm. This feature can be used to support or reject specific structural models of organometallic transition metal complexes in solution by comparing the experimental and Density Functional Theory (DFT) calculated values of this ³¹P chemical shift. This approach has been tested on a variety of the metals of groups 8-12 and molecular structures. General recommendations for appropriate basis sets are reported.

Modeling of Solute-Solvent Interactions Using an External Electric Field-From Tautomeric Equilibrium in Nonpolar Solvents to the Dissociation of Alkali Metal Halides

An implicit account of the solvent effect can be carried out using traditional static quantum chemistry calculations by applying an external electric field to the studied molecular system. This approach allows one to distinguish between the effects of the macroscopic reaction field of the solvent and specific solute-solvent interactions. In this study, we report on the dependence of the simulation results on the use of the polarizable continuum approximation and on the importance of the solvent effect in nonpolar solvents. The latter was demonstrated using experimental data on tautomeric equilibria between the pyridone and hydroxypyridine forms of 2,6-di-tert-butyl-4-hydroxy-pyridine in cyclohexane and chloroform.

Nuclear Quantum Effects Largely Influence Molecular Dissociation and Proton Transfer in Liquid Water under an Electric Field

Proton transfer in liquid water controls acid-base chemistry, crucial enzyme reactions, and the functioning of fuel cells. Externally applied static electric fields in water are capable of dissociating molecules and transferring protons across the H-bond network. However, the impact of nuclear quantum effects (NQEs) on these fundamental field-induced phenomena has not yet been reported. By comparing state-of-the-art ab initio molecular dynamics (AIMD) and path integral AIMD simulations of water under electric fields, I show that quantum delocalization of the proton lowers the molecular ionization threshold to approximately one-third. Moreover, also the water behavior as a protonic semiconductor is considerably modified by the inclusion of NQEs. In fact, when the quantum nature of the nuclei is taken into account, the proton conductivity is ∼50% larger. This work proves that NQEs sizably affect the protolysis phenomenon and proton transfer in room-temperature liquid water.

Solvent Organization and Rate Regulation of a Menshutkin Reaction by Oriented External Electric Fields are Revealed by Combined MD and QM/MM Calculations

When and how do external electric fields (EEFs) lead to catalysis in the presence of a (polar or nonpolar) solvent? This is the question that is addressed here using a combination of molecular dynamics (MD) simulations, quantum mechanical/molecular mechanical calculations with EEF, and quantum mechanical/(local) electric field calculations. The paper focuses on a model reaction, the Menshutkin reaction between CH3I and pyridine in three solvents of varying polarity. Using MD simulations, we find that the EEF causes the solvent to undergo organization; the solvent molecules gradually align with the applied field as the field strength increases. The collective orientation of the solvent molecules modifies the electrostatic environment around the Menshutkin species and induces a global electric field pointing in the opposite direction of the applied EEF. The combination of these two entangled effects leads to partial or complete screening of the EEF, with the extent of screening being proportional to the polarity/polarizability of the solvent. Nevertheless, we find that catalysis of the Menshutkin reaction inevitably emerges once the EEF exceeds the opposing field of the organizing solvent, i.e., once polarization of the Menshutkin complex is observed to set in. Overall, our analysis provides a lucid and pictorial interpretation of the behavior of solutions in the presence of EEFs and indicates that EEF-mediated catalysis should, in principle, be feasible in bulk setups, especially for nonpolar and mildly polar solvents. By application of the charge-transfer paradigm, it is shown that the emergence of OEEF catalysis in solution can be generalized to other reactions as well.

Possibility of realizing superionic ice VII in external electric fields of planetary bodies

In a superionic (SI) ice phase, oxygen atoms remain crystallographically ordered while protons become fully diffusive as a result of intramolecular dissociation. Ice VII's importance as a potential candidate for a SI ice phase has been conjectured from anomalous proton diffusivity data. Theoretical studies indicate possible SI prevalence in large-planet mantles (e.g., Uranus and Neptune) and exoplanets. Here, we realize sustainable SI behavior in ice VII by means of externally applied electric fields, using state-of-the-art nonequilibrium ab initio molecular dynamics to witness at first hand the protons' fluid dance through a dipole-ordered ice VII lattice. We point out the possibility of SI ice VII on Venus, in its strong permanent electric field.

Massive generation of metastable bulk nanobubbles in water by external electric fields

Nanobubbles (NBs) are nanoscopic gaseous domains than can exist on solid surfaces or in bulk liquids. They have attracted substantial attention due to their long-time (meta)stability and a high potential for real-world applications. Using an approach not previously investigated, we exploit surface-electrostatic NB formation and stabilization via application of external electric fields in gas-liquid systems, with the marked result of massively increased gas uptake into the liquid in NB form. The de facto gas solubility enhancement (over many months) ranges from 2.5-fold for oxygen to 30-fold for methane vis-à-vis respective Henry's law values for gas solubility; the more hydrophobic the gas, the more spectacular the increase. Molecular dynamics simulations reveal that the origin of NBs' movement lies in dielectrophoresis, while substantial NB stabilization arises from a surface-polarization interaction.

Polarizable MD simulations of ionic liquids: How does additional charge transfer change the dynamics?

A new model for treating charge transfer in ionic liquids is developed and applied to 1-ethyl-3-methylimidazolium tetrafluoroborate. The model allows for us to examine the roles of charge transfer, polarizability, and charge scaling effects on the dynamics of ionic liquids.

Ab initio spectroscopy of water under electric fields

Whereas a broad range of literature exists on the spectroscopy of water in disparate conditions, infrared (IR) and Raman spectra of water subjected to electric fields have never extensively been investigated so far. Based on ab initio molecular dynamics simulations, here we present IR and Raman spectra of bulk liquid water under the effect of static electric fields. A contraction of the entire frequency range is recorded upon increasing the field intensity both in the IR and in the Raman spectra. Whilst the OH stretching band is progressively shifted toward lower frequencies - indicating a field-induced strengthening of the H-bond network - all the other bands are up-shifted by the field. Furthermore, an evident modification of the librational mode band appears in all the spectra. Finally, the order-maker action of the field emerges also from the increase of the water orientational tetrahedral order. Upon field exposure, the water structure becomes more "ice like".

The influence of polarizability and charge transfer on specific ion effects in the dynamics of aqueous salt solutions

The diffusion rates for water molecules in salt solutions depend on the identity of the ions, as well as their concentration. Among the alkali metal ions, cesium and potassium increase and sodium strongly decreases the diffusion constant of water. The origin of the difference can be understood by examining the simulation results using different potential models. In this work, aqueous solutions of salts are simulated with a variety of models. Commonly used non-polarizable models, which otherwise reproduce many experimental properties, do not capture the trend in the diffusion constant, while models which include polarization and/or charge transfer interactions do. For the non-polarizable models, the diffusion constant decreases too strongly with salt concentration. The changes in the water diffusion constant with increasing salt concentration match the diffusion constant of the ion. The ion diffusion constant is dependent on the residence time for water in the ion solvation shell. The non-polarizable models over-estimate the residence time, relative to the translational diffusion constant and so tend to under-estimate the ion and water diffusion constants.

Synthesis of (d)-erythrose from glycolaldehyde aqueous solutions under electric field

A Miller-like numerical experiment demonstrates that ubiquitarious molecules such as water and glycolaldehyde can synthesize (d)-erythrose, one of the direct precursors of ribose.

Novel electrochemical route to cleaner fuel dimethyl ether

Methanol, the simplest alcohol, and dimethyl ether, the simplest ether, are central compounds in the search for alternative “green” combustion fuels. In fact, they are generally considered as the cornerstones of the envisaged “Methanol Economy” scenario, as they are able to efficiently produce energy in an environmentally friendly manner. However, despite a massive amount of research in this field, the synthesis of dimethyl ether from liquid methanol has never so far been reported. Here we present a computational study, based on ab initio Molecular Dynamics, which suggests a novel synthesis route to methanol dehydration – leading thus to the dimethyl ether synthesis – through the application of strong electric fields. Besides proving the impressive catalytic effects afforded by the field, our calculations indicate that the obtained dimethyl ether is stable and that it can be progressively accumulated thanks to the peculiar chemical pathways characterising the methanol reaction network under electric field. These results suggest that the experimental synthesis of dimethyl ether from liquid methanol could be achieved, possibly in the proximity of field emitter tips.

Ionic diffusion and proton transfer in aqueous solutions of alkali metal salts

We report on a series of ab initio molecular dynamics investigations on LiCl, NaCl, and KCl aqueous solutions under the effect of static electric fields.

One-step electric-field driven methane and formaldehyde synthesis from liquid methanol

By means of state-of-the-art computational approaches, a new fundamental chemical reaction, involving formaldehyde and methane, has been observed when an electric field is applied to liquid methanol.

The reaction pathways connecting methanol to methane and formaldehyde are among the most emblematic in chemistry because of their outstanding interest in the fields of energy, synthesis, and bio- and geo-chemistry. Despite of its fundamental nature, the one-pot synthesis of formaldehyde and methane stemming from methanol has never been reported before. Here we present a study, based on ab initio molecular dynamics and free-energy methods, in which the simultaneous oxidation and reduction (i.e., the disproportionation) of liquid methanol into methane and formaldehyde has been achieved at ambient temperature through the application of a static electric field. Because strong electric fields can be generated in the proximity of field emitter tips, this finding shows that the challenge of experimentally disproportionating methanol into formaldehyde and methane could be attempted. We show that the methanol “solvent” molecules play a major role in this process and that the chemical pathway connecting methanol to the detected products in the bulk liquid phase is very different from its reproduced gas-phase counterpart. Finally, we demonstrate that switching on an external electric field drastically modifies the reaction network of methanol, lowering some activation barriers, stabilizing the methane and formaldehyde products, and opening otherwise difficult-to-achieve chemical routes.

Ab initio molecular dynamics study of an aqueous NaCl solution under an electric field

Ab initio molecular dynamics simulations of salty water under an electric field reveal two regimes of the relative mobilities of chlorine and sodium ions. When water dissociation and proton transfer are actived at strong field intensities, the presence of the ions hinders the efficiency of the proton transfer mechanism.

Fluctuations of local electric field and dipole moments in water between metal walls

We examine the thermal fluctuations of the local electric field Ek (loc) and the dipole moment μk in liquid water at T = 298 K between metal walls in electric field applied in the perpendicular direction. We use analytic theory and molecular dynamics simulation. In this situation, there is a global electrostatic coupling between the surface charges on the walls and the polarization in the bulk. Then, the correlation function of the polarization density pz(r) along the applied field contains a homogeneous part inversely proportional to the cell volume V. Accounting for the long-range dipolar interaction, we derive the Kirkwood-Fröhlich formula for the polarization fluctuations when the specimen volume v is much smaller than V. However, for not small v/V, the homogeneous part comes into play in dielectric relations. We also calculate the distribution of Ek (loc) in applied field. As a unique feature of water, its magnitude |Ek (loc)| obeys a Gaussian distribution with a large mean value E0 ≅ 17 V/nm, which arises mainly from the surrounding hydrogen-bonded molecules. Since |μk|E0 ∼ 30kBT, μk becomes mostly parallel to Ek (loc). As a result, the orientation distributions of these two vectors nearly coincide, assuming the classical exponential form. In dynamics, the component of μk(t) parallel to Ek (loc)(t) changes on the time scale of the hydrogen bonds ∼5 ps, while its smaller perpendicular component undergoes librational motions on time scales of 0.01 ps.

A matter of quantum voltages

Voltages inside matter are relevant to crystallization, materials science, biology, catalysis, and aqueous chemistry. The variation of voltages in matter can be measured by experiment, however, modern supercomputers allow the calculation of accurate quantum voltages with spatial resolutions of bulk systems well beyond what can currently be measured provided a sufficient level of theory is employed. Of particular interest is the Mean Inner Potential (V(o))--the spatial average of these quantum voltages referenced to the vacuum. Here we establish a protocol to reliably evaluate V(o) from quantum calculations. Voltages are very sensitive to the distribution of electrons and provide metrics to understand interactions in condensed phases. In the present study, we find excellent agreement with measurements of V(o) for vitrified water and salt crystals and demonstrate the impact of covalent and ionic bonding as well as intermolecular/atomic interactions. Certain aspects in this regard are highlighted making use of simple model systems/approximations. Furthermore, we predict V(o) as well as the fluctuations of these voltages in aqueous NaCl electrolytes and characterize the changes in their behavior as the resolution increases below the size of atoms.