Dynamics and Spectral Response of Water Molecules around Tetramethylammonium Cation

A Water-in-Salt Electrolyte for Room-Temperature Fluoride-Ion Batteries Based on a Hydrophobic-Hydrophilic Salt

Realizing room-temperature, efficient, and reversible fluoride-ion redox is critical to commercializing the fluoride-ion battery, a promising post-lithium-ion battery technology. However, this is challenging due to the absence of usable electrolytes, which usually suffer from insufficient ionic conductivity and poor (electro)chemical stability. Herein we report a water-in-salt (WIS) electrolyte based on the tetramethylammonium fluoride salt, an organic salt consisting of hydrophobic cations and hydrophilic anions. The new WIS electrolyte exhibits an electrochemical stability window of 2.47 V (2.08-4.55 V vs Li+/Li) with a room-temperature ionic conductivity of 30.6 mS/cm and a fluoride-ion transference number of 0.479, enabling reversible (de)fluoridation redox of lead and copper fluoride electrodes. The relationship between the salt property, the solvation structure, and the ionic transport behavior is jointly revealed by computational simulations and spectroscopic analysis.

Hydrophilic and Apolar Hydration in Densely Grafted Cationic Brushes and Counterions with Large Mobilities

We employ an all-atom molecular dynamics (MD) simulation framework to unravel water microstructure and ion properties for cationic [poly(2-(methacryloyloxy)ethyl) trimethylammonium chloride] (PMETAC) brushes with chloride ions as counterions. First, we identify locally separate water domains (or first hydration shells) each around {N(CH3)3}+ and the C═O functional groups of the PMETAC chain and one around the Cl- ion. These first hydration shells around the respective moieties overlap, and the extent of the overlap depends on the nature of the species triggering it. Second, despite the overlap, the water molecules in these domains demonstrate disparate properties dictated by the properties of the atoms and groups around which they are located. For example, the presence of the methyl groups makes the {N(CH3)3}+ group trigger apolar hydration as evidenced by the corresponding orientation of the dipole of the water molecules around the {N(CH3)3}+ moiety. These water molecules around the {N(CH3)3}+ group also have enhanced tetrahedrality compared to the water molecules constituting the hydration layer around the C═O group and the Cl- counterion. Our simulations also identify that there is an intervening water layer between the Cl- ion and {N(CH3)3}+ group: this layer prevents the Cl- ion from coming very close to the {N(CH3)3}+ group. As a consequence, there is a significantly large mobility of the Cl- ions inside the PMETAC brush layer. Furthermore, the C═O group of the polyelectrolyte (PE) chain, due to the partial negative charge on the oxygen atom and the specific structure of the PMETAC brush system, demonstrates strongly hydrophilic behavior and enforces a specific dipole response of water molecules analogous to that experienced by water around anionic species of high charge density. In summary, our findings confirm that PMETAC brushes undergo hydrophilic hydration at one site and apolar hydration at another site and ensure large mobility of the supported Cl- counterions.

Structural and dynamical properties of concentrated alkali- and alkaline-earth metal chloride aqueous solutions

Concentrated ionic aqueous electrolytes possess a diverse array of applications across various fields, particularly in the field of energy storage. Despite extensive examination, the intricate relationships and numerous physical mechanisms underpinning diverse phenomena remain incompletely understood. Molecular dynamics simulations are employed to probe the attributes of aqueous solutions containing LiCl, NaCl, KCl, MgCl2, and CaCl2, spanning various solute fractions. The primary emphasis of the simulations is on unraveling the intricate interplay between these attributes and the underlying physical mechanisms. The configurations of cation-Cl- and Cl--Cl- pairs within these solutions are disclosed. As the solute fraction increases, consistent trends manifest regardless of solute type: (i) the number of hydrogen bonds formed by the hydration water surrounding ions decreases, primarily attributed to the growing presence of counter ions in proximity to the hydration water; (ii) the hydration number of ions exhibits varying trends influenced by multiple factor; and (iii) the diffusion of ions slows down, attributed to the enhanced confinement and rebound of cations and Cl- ions from the surrounding atoms, concurrently coupled with the changes in ion vibration modes. In our analysis, we have, for the first time, clarified the reasons behind the slowing down of the diffusion of the ions with increasing solute fraction. Our research contributes to a better understanding and manipulation of the attributes of ionic aqueous solutions and may help designing high-performance electrolytes.

Direct Correlation between Short-Range Vibrational Spectral Diffusion and Localized Ion-Cage Dynamics of Water-in-Salt Electrolytes

The molecular dynamics simulations of a "water-in-salt" electrolyte, lithium bis(trifluoromethyl sulfonyl) imide (LiNTf₂), with a varying concentration range of 3 to 20 m were performed to establish a direct connection between a dynamic property like the ion-cage lifetime with the short-range vibrational stretching frequency shift of the used probe, HOD. The properties reported here are compared to that obtained from experiments performed at the same concentrations. The time-series wavelet transform was adopted as a preferable mathematical tool for calculating the instantaneous fluctuating frequencies of the probe O-D stretch mode and the concentration-dependent vibrational stretch spectral signature based on the variable functions associated with a particular chemical bond derived from classical molecular dynamics trajectories. The decay time constants of frequency fluctuations and the lifetime of the ion cage (τ_IC) were estimated as a function of salt concentration. Herein, we emphasize the correlation between the slowest time constant (τ₃) of the decay of O-D stretch frequency fluctuations and the timescales associated with the lifetime of ion cages (τ_IC). The results exhibit that the existing relationships were also concentration-dependent. Therefore, this study highlights the connection between the ionic motions that regulate the overall system dynamics with the short-range vibrational frequency shift of the used probe, which was used similar to experiments. It also provides an understanding of the interionic interactions and the dynamical and spectral properties of the electrolytic mixtures. We establish a direct correlation between short-range frequency profile and localized ion-cage lifetime, which can fill the gap of understanding between viscosity, vibrational frequency, and ion-cage dynamics of electrolytes.

Molecular Simulation-Guided Spectroscopy of Imidazolium-Based Ionic Liquids and Effects of Methylation on Ion-Cage and -Pair Dynamics

Classical molecular dynamics simulations were performed to assess an atomistic interpretation of the ion-probe structural interactions in two typical ionic liquids (ILs), 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [BMIm][NTf₂] and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide [BDimIm][NTf₂] through computational ultrafast spectroscopy. The nitrile stretching vibrations of the thiocyanate anion, [SCN]^-, serve as the local mode of the ultrafast system dynamics within the imidazolium-based ionic liquid environment. The wavelet transform of classical trajectories determines the time-varying fluctuating frequencies and the stretch spectral signatures of SCN^- in the normalized distribution. However, computational modeling of the two-dimensional (2D) spectra from the wavelet-derived vibrational frequencies yields time evolution of the local molecular structure along with the varied time-dependent dynamics of the spectral diffusion process. We calculated the frequency-frequency correlation functions (FFCFs), time correlations associated with the ion-pair and -cage dynamics, and mean square displacements as a function of time, depicting diffusive dynamics. The calculated results based on the pair correlation functions and the distribution of atomic density suggest that the hydrogen and methylated carbon at the two-position of the imidazolium ring of [BMIm] and [BDimIm] cations, respectively, strongly interact with the probe through the N of the thiocyanate anion rather than the S atom. The center-of-mass center-of-mass (COM-COM) cation-probe radial distribution functions (RDFs) in conjunction with the site-specific structural analysis further reveal well-structured interactions of the thiocyanate ion and [BMIm]⁺ cation rather than the [BDimIm] cation. In contrast, the anion-probe COM-COM RDFs depict weak interactive associations within the vibrational probe [SCN]^- and [NTf₂]^- ions. Methylation at the two-position of the imidazolium ring predicts slower structural reorganization and breaking and reformation dynamics of the ion pairs and cages within the ionic liquid framework.

Multiple Ensembles of the Hydrogen-bonded Network in Ethylammonium Nitrate versus Water from Vibrational Spectral Dynamics of SCN- Probe

We performed classical molecular dynamics simulations to monitor the structural interactions and ultrafast dynamical and spectral response in the protic ionic liquid, ethylammonium nitrate (EAN) and water using the nitrile stretching mode of thiocyanate ion (SCN-) as the vibrational probe. The normalized stretch frequency distribution of nitrile stretch of SCN- attains an asymmetric shape in EAN, indicating the existence of more than one hydrogen-bonding environment in EAN. We computed the 2D IR spectrum from classical trajectories, applying the response function formalism. Spectral diffusion dynamics in EAN undergo an initial rattling of the SCN - inside the local ion-cage occurring at a timescale of 0.10 ps, followed by the breakup of the ion-cage activating molecular diffusion at 7.86 ps timescale. In contrast, the dynamics of structural reorganization occur at a timescale of 0.58 ps in H 2 O. Hence, the time dependence of the frequency-frequency correlation function decay hints at the local molecular structure and ultrafast ion dynamics of the SCN - probe. The loss of frequency correlation read from the peak shape changes in the 2D correlation spectrum as a function of waiting time is faster in H 2 O than in EAN due to the enhanced structural ordering and higher viscosity of the latter. We provide an atomic-level interpretation of the solvation environment around SCN - in EAN and water, which indicates the multiple ensembles of the hydrogen bond network in EAN.

Ionic Dynamics and Vibrational Spectral Diffusion of a Protic Alkylammonium Ionic Salt through Intrinsic Cationic N-H Vibrational Probe from FPMD Simulations

We employed density functional theory (DFT)-based molecular dynamics simulations to explore the structure, dynamics, and spectral properties of the protic ionic entity trimethylammonium chloride (TMACl). Structural investigations include calculating the site-site radial distribution functions (RDFs), the distribution of constituent cations and anions in three-dimensional space, and combined distribution functions of the hydrogen-bonded pair RDF versus angle, revealing the structural characteristics of the ionic solvation and the intermolecular interactions within ions. Further, we determined the instantaneous vibrational stretching frequencies of the intrinsic N-H stretch probe modes by applying the time-series wavelet method. The associated ionic dynamics within the protic ionic compound were investigated by calculating the time-evolution of the fluctuating frequencies and the frequency-time correlation functions (FFCFs). The time scale related to the local structural relaxation process and the average hydrogen bond lifetime, ion cage dynamics, and mean squared displacement were investigated. The faster decay component of the FFCFs, depicting the intermolecular motion of intact hydrogen bonds in TMACl, is 0.07 ps for the Perdew-Burke-Ernzerhof (PBE)-based simulation and 0.06 ps for the PBE-D2 representation. The slower time scale of the longer picosecond decay time component of PBE and PBE-D2 representations are 3.13 and 2.87 ps, respectively. These picosecond time scales represent more significant fluctuations of the hydrogen-bonding partners in the ionic entity and hydrogen-bond jump events accompanied by large angular jumps. The longest picosecond time scales represent structural relaxation, including large angular jumps and ion-pair dynamics. Also, ion cage lifetimes correlate with the slowest time scale of the associated dynamics of vibrational spectral diffusion despite the type of DFT functional. This study benchmarks DFT treatments of the exchange-correlation functional with and without the van der Waals (vdW) dispersion correction scheme. The inclusion of vdW interactions to the PBE functional represents a less structured state of the ionic entity and faster dynamics of the molecular motions relative to the one predicted by the PBE system. All the results illustrate the necessity of accurately describing the Coulomb interactions, vdW dispersive interactive forces, and localized hydrogen bonds required to sustain the energetic balance in this ionic salt.

Microheterogeneity-Induced Vibrational Spectral Dynamics of Aqueous 1-Alkyl-3-methylimidazolium Tetrafluoroborate Ionic Liquids of Different Cationic Chain Lengths

We have monitored the impacts of an increment in the alkyl chain length of the imidazolium-based tetrafluoroborate ionic liquids on the local deuteroxyl probe modes of interest. For this study, we have taken 1-ethyl-3-methylimidazolium tetrafluoroborate [EMIm][BF₄], 1-butyl-3-methylimidazolium tetrafluoroborate [BMIm][BF₄], 1-octyl-3-methylimidazolium tetrafluoroborate [OMIm][BF₄], and 1-decyl-3-methylimidazolium tetrafluoroborate [DMIm][BF₄] ionic liquid solutions with 5% HOD in H₂O as the vibrational reporter of the associated ultrafast system dynamics. Classical molecular dynamics (MD) simulations were employed to determine molecular structure and dynamic properties, while the spectral profiles were derived by applying the wavelet analysis of classical trajectories. Spatial distribution functions reveal the heterogeneity within the molecular structures of the ionic liquids (ILs) with varying alkyl chain lengths. The intense position of the spectral peak, the frequency corresponding to the shoulder peak, and the spectral linewidth of the O-D stretch distribution are not influenced by the increment in the cationic chain length. In addition, the ionic liquid (IL) [BMIm][BF₄] exhibits a notable trend; the dynamic timescales are longer than the other studied systems. Therefore, we have performed the Voronoi decomposition analysis of the ionic and the polar-apolar domains, symmetrically increasing the length of alkyl chains on the IL cations. Domain analysis reveals structural microheterogeneity; the anions form discrete domains, and the ionic liquid constituting cations form continuous domains irrespective of the alkyl chain length on the imidazolium cations. Therefore, this computational ultrafast spectroscopy study aids in forming a molecular-level picture of the ionic liquid cations and anions in the liquid phase, providing a detailed interpretation of the spectral properties of the probe stretching vibrations.

How Far Is "Bulk Water" from Interfaces? Depends on the Nature of the Surface and What We Measure

Using systematic molecular dynamics (MD) simulations, we revisit the question: At what distance from an interface do the properties of "bulk water" get recovered? We have considered three different kinds of interfaces: nonpolar (hydrophobic; isooctane-water interface), charged (negative; AOT bilayer), and polar (zwitterionic; POPC bilayer). In order to interrogate the extent of perturbation of the interfacial water molecules as a function of the distance from the interface, we utilize a diverse range of structural and dynamical parameters. To capture the structural perturbations, we look into local density (translational order), local tetrahedral order parameter, and dipolar orientation of the water molecules. We also explore the anisotropic diffusion of the water molecules in the direction perpendicular to the interface as well as the planar diffusion parallel to the interface in a distance dependent manner. In addition, the orientational time correlation functions have been computed to understand the extent of slowdown in the rotational dynamics. As expected, the electrostatic field emanating from the charged AOT interface seems to have the highest long-range effect on the orientational order and dynamics of the water molecules, whereas specific interactions like hydrogen bonding and electrostatic interaction lead to significant trapping and kinetic slowdown for both AOT and POPC (zwitterionic) very close to the interface. Our analysis highlights that not only the length-scale of perturbation depends on the nature of the interfaces and specific interactions but also the type of water property that we measure/calculate. Different water properties seem to have widely different length-scale of perturbation. Orientational order parameters seem to be perturbed to a much longer length-scale as compared to translational order parameters. The global orientational order of water can be perturbed even up to ∼4-5 nm near the negatively charged AOT surface in the absence of any extra electrolyte. This observation has significant implication toward the interpretation of experimental measurements as well since different spectroscopic techniques would probe different parameters or water properties with possible mutual disagreement and inconsistency between different types of measurements. Thus, our study provides a broader and unifying perspective toward the aspect of "context dependent" structural and dynamical perturbation of "interfacial water".

Size-Dependent Order-Disorder Crossover in Hydrophobic Hydration: Comparison between Spherical Solutes and Linear Alcohols

Theory and computer simulation studies have predicted that water molecules around hydrophobic molecules should undergo an order-disorder transition with increasing solute size around a 1 nm length scale. Some theories predict the formation of a clathrate-like ordered structure around smaller hydrophobic solutes (<1 nm) and the formation of disordered vapor-liquid interfaces around larger solutes (>1 nm) and surfaces. Experimental validation of these predictions has often been elusive and contradictory. High-resolution Raman spectroscopy has detected that water around small hydrophobic solutes shows a signature similar to that of bulk water at lower temperature (increased ordering and a stronger hydrogen-bonded network). Similarly, water around larger solutes shows an increasing population of dangling OH bonds very similar to higher temperature bulk water. Thus, the solute size dependence of the structure and dynamics of water around hydrophobic molecules seems to have an analogy with the temperature dependence in bulk water. In this work, using atomistic classical molecular dynamics (MD) simulations, we have systematically investigated this aspect and characterized this interesting analogy. Structural order parameters including the tetrahedral order parameter (Q), hydrogen bond distribution, and vibrational power spectrum highlight this similarity. However, in contrast to the experimental observations, we do not observe any length-dependent crossover for linear hydrophobic alcohols (n-alkanols) using classical MD simulations. This is in agreement with earlier findings that linear alkane chains do not demonstrate the length-dependent order-disorder transition due to the presence of a sub-nanometer length scale along the cross section of the chain. Moreover, the collapsed state of linear hydrocarbon chains is not significantly populated for smaller chains (number of carbons below 20). In the context of our computational results, we raise several pertinent questions related to the sensitivity of various structural and dynamical parameters toward capturing these complex phenomena of hydrophobic hydration.

Amphiphilicity of Intricate Layered Graphene/g-C3N4 Nanosheets

The hybrid heterostructure of the tri-s-triazine form of graphitic carbon nitride (g-C₃N₄), a stable two-dimensional material, results from intricate layer formation with graphene. In this material, g-C₃N₄, an amphiphilic material, stabilizes Pickering emulsions as an emulsifier and can effectively disperse graphene. Due to the various technological applications of the hybrid nanosheets in an aqueous environment, it is essential to study the interaction of water molecules with graphene and g-C₃N₄ (Gr/g-C₃N₄)-combined heterostructure. Although few studies have been performed signifying the water orientation in the interfacial layer, we find that there is a lack of detailed studies using various dynamical and structural properties of the interfacial water molecules. The interface of the Gr/g-C₃N₄ hybrid structure, one of the rarely found amphiphilic interfaces (on the g-C₃N₃ side), is appropriate for exploring the water affinity due to the availability of heterogeneous interfacial aqueous interactions. We adopted classical molecular dynamics simulations using two models for water molecules to study the structure and dynamics of an aqueous interface. We have correlated the structural properties to dynamics and spectral properties to understand the overall behavior of the amphiphilic interface. Our results branch into two significant hydrogen bond (HB) properties in HB count and HB strength among the water molecules in the different layers. The HB counts in the different layers of water are correlated using the average distance distribution (P_rO4), tetrahedral order parameters, HB donor/acceptor count, and total HBs per water molecule. A conspicuous difference is found in the HB count and related dynamics of the system. The HB lifetime and diffusion coefficient hint at the equivalent strength of HBs in the different layers. All the findings conclude that the amphiphilicity of the Gr/g-C₃N₄ interface can help in understanding various interfacial physical and chemical processes.

Aqueous Affinity and Interfacial Dynamics of Anisotropic Buckled Black Phosphorous

The structure of black phosphorous (BP) is similar to the honeycomb arrangement of graphene, but the layered BP is found to be buckled and highly anisotropic. The buckled surface structure affects interfacial molecule mobility and plays a vital role in various nanomaterial applications. The BP is also known for wettability, droplet formation, stability, and hydrophobicity in the aqueous environment. However, there is a gap concerning the structural and dynamical behavior of water molecules, which is available in abundance for other monoatomic and polyatomic two-dimensional (2D) materials. Motivated by the technological importance, we try to bridge the gap by explaining the surface anisotropy-facilitated behavior of water molecules on bilayer BP using classical and first principles molecular dynamics (MD) simulations. From our classical MD study, we find three distinct layers of water molecules. The water layer closest to the interface is L1, followed by L2 and L3/bulk perpendicular to the BP surface. Water molecules in the L1 layer experience some structural disintegration in hydrogen bond (HB) phenomena compared to the bulk. There is a loss of HB donor-acceptor count per water molecule. The average HB count decreases because of an elevated rate of HB formation and deformation; this would affect the dynamic properties in terms of HB lifetime. Therefore, we observe the reduced lifetime of HB in the layer in close contact with BP, which again complements our finding on the diffusion coefficient of water molecules in distinct layers. Water diffuses relatively faster with diffusion coefficient 3.25 × 10^-9 m² s^-1 in L1, followed by L2 and L3. The BP layer shows moderate hydrophobic nature. Our results also indicate the anisotropic behavior as the diffusion along the x-direction is faster than that along the y-direction. The gap in the slope of the x and y components of mean-squared displacement (MSD) complements the pinning effect in an aqueous environment. We observe blue-shifted and red-shifted libration and O-H stretching modes from the calculated power spectra for the L1 water molecules compared to the L2 and L3 molecules from first principles MD simulations. Our analysis may help understand the physical phenomena that occur during the surface wetting of the predroplet formation process observed experimentally.

Dynamics of Ionic Liquid through Intrinsic Vibrational Probes Using the Dispersion-Corrected DFT Functionals

First principles molecular dynamics simulations have been utilized to study the spectral properties of the protic ionic liquid, methylammonium formate (MAF). All simulations were performed using density functional theory (DFT) and various van der Waals-corrected exchange-correlation functionals. We calculated the vibrational stretch frequency distributions, determined the time-frequency correlations of the intrinsic vibrational probes, the N-H and C-O modes in MAF, and the frequency-structure correlations. We also estimated the average hydrogen-bond lifetimes and orientation dynamics to capture the ultrafast spectral response. The spectroscopic signature of the N-H stretching vibrations using the Becke-Lee-Yang-Parr (BLYP) and Perdew-Burke-Ernzerhof (PBE) functionals displays a spectral shift in the lower frequency side, suggesting stronger hydrogen-bonding interactions represented by the gradient approximation functionals than the van der Waals (vdW)-corrected simulations. The carboxylate frequency profiles with the dispersion-corrected representations are almost similar without a significant difference in the normalized distributions. Besides, the COO stretching frequencies at the peak maxima positions of the PBE functionals exhibit a lesser deviation from the experimental data. Spectral diffusion dynamics of the intrinsic vibrational probes on the cationic and anionic sites of the ionic liquid proceed through a short time relaxation of the intact hydrogen bonds followed by an intermediate time constant and a longer time decay indicating the switchover of hydrogen bonds. Dispersion-corrected atom-centered one-electron potential (DCACP) correction added to the BLYP system slows down the picosecond time scales of frequency correlation and the time constants of rotational motion, lengthening the overall system dynamics. The observed trends in the time-dependent decays of frequency fluctuations and the orientation autocorrelation functions correlate with the structural interactions in liquid MAF and hydrogen-bond dynamics. In this study, we examine the predictions made by different density functional treatments comparing the results of the uncorrected BLYP and PBE representations with the semiempirical vdW methods of Grimme and matching our calculated data with the experimental observations.

Insignificant Effect of Temperature on the Structure and Angular Jumps of Water near a Hydrophobic Cation

The ambiguity in the behavior of water molecules around hydrophobic solutes is a matter of interest for many studies. Motivated by the earlier results on the dynamics of water molecules around tetramethylammonium (TMA) cation, we present the effect of temperature on the structure and angular jumps of water due to hydrophobicity using first principles molecular dynamics simulations. The average intermolecular distance between the central oxygen and four nearest neighbors is found to be the highest for water molecules in the solvation shell of TMA at 400 K, followed by the same at 330 K. The hydrogen bond (HB) donor-acceptor count, HB per water molecule, and tetrahedral order parameter suggests the loss of tetrahedrality in the solvation shell. Elevated temperature affects the tetrahedral parameter in local regions. The HB jump mechanism is studied for methyl hydrogen and water molecules in the solvation shell. Observations hint at the presence of dangling water molecules in the vicinity of the hydrophobic cation, and no evidence is found for the enhanced structural ordering of nearby water molecules.

Conformation-induced vibrational spectral dynamics of hydrogen peroxide and vicinal water molecules

We studied the conformation-induced spectral response of water molecules due to site-specific structural alterations of solvated hydrogen peroxide (H₂O₂) employing DFT-based first principles molecular dynamics (FPMD) simulations. Wavelet transform was used to determine the time-dependent frequencies of the hydroxyls of water molecules and the O-H stretch modes of H₂O₂. Shifts in the vibrational frequency of the hydrogen-bonded hydroxyls inside the solvation shell of H₂O₂ support multiple distinctive hydrogen bonding environments. This paper classifies two distinct hydrogen bond types inside the O-O_W solvation shell of H₂O₂, and the dynamical calculations provide a quantitative estimation of the relative hydrogen bond strength. We ascertain the reason for not observing the escape of water molecules from the hydrogen peroxide hydration shell, unlike the solvation shell of ionic solutions and neutral solutes. Besides, we provide a comprehensive analysis of the spectral shifts in the normalized frequency distribution, the time-dependent decay of frequency-frequency correlation functions, and the hydrogen bond length scale fluctuations. We also quantify the relative contribution of the cisoid and transoid conformers affecting the vibrational spectral signature of the vicinal water molecules. While the transoid conformers promote the hydrogen bonding interactions through the oxygen site (OH_W), the cisoid conformers facilitate hydrogen peroxide-water hydrogen bond formation through the hydrogen site (HO_W). These non-identical hydrogen bond associations stabilize hydrogen peroxide in water.

Conformational dynamics of aqueous hydrogen peroxide from first principles molecular dynamics simulations

We performed first principles molecular dynamics simulations of a relatively dilute aqueous hydrogen peroxide (H2O2) solution to examine its structural alterations and relevant dynamics upon solvation. The internal rotation of the OH groups about the O-O bond facilitates the flexible structure of H2O2. Structural calculations reveal dihedral angle fluctuations in the aqueous solution. Water molecules make stronger hydrogen bonds through the hydrogen atom of the solute than the oxygen atom leading to distinct hydrogen bonding configurations inside the first solvation shell. Time-dependent dihedral angle alterations result in conformational changes and the normalized dihedral angle distribution plot displays characteristic peaks at ∼100-120° and ∼230°, illustrating various conformational states. Within the simulation time, flexibility-induced interconversion of hydrogen peroxide gives rise to several cisoid and transoid conformers. In this study, we examine the relative population of the associated conformational states and the lifetime of the cisoid and transoid conformers from the torsion angle variations. We also determine the free energy landscape of the rotational isomerization process in H2O2 and explore two distinct energy barriers during such interconversion.

Solvation Structures of Tetraethylammonium Bromide and Tetrafluoroborate in Aqueous Binary Solvents with Ethanol, Trifluoroethanol, and Acetonitrile

The solvation structures of tetraethylammonium bromide and tetrafluoroborate (TEABr and TEABF₄) in aqueous binary solvents with ethanol (EtOH), 2,2,2-trifluoroethanol (TFE), and acetonitrile (AN) have been clarified by molecular dynamics (MD) simulations. In addition, ¹H and ¹³C NMR chemical shifts of the H and C atoms within TEA⁺ in the binary solvents have been measured as a function of the mole fraction of the organic solvent, x_OS. The variations of the chemical shifts with an increase in x_OS were interpreted according to the solvation structures of TEA⁺, Br^-, and BF₄^- obtained from the MD simulations. It has been found that TEABF₄ at 130 mmol dm^-3 cannot be dissolved into the EtOH and TFE solvents above x_OS ≈ 0.7 and 0.6, respectively, while TEABr can be done in both solvents. Interestingly, TEABr and TEABF₄ at the concentration can be dissolved in the AN solvents over the entire x_OS range. The solvation of TEA⁺, Br^-, and BF₄^- in each solvent has been discussed in terms of the electrostatic force, the weak hydrogen bond of C-H···F-C, and the dipole-dipole interaction.