Formation of mesoscopic water networks in aqueous systems

Molecular-scale understanding of diluent effects on ligand assembly for metal ion separations

Use of metal-selective ligands in solvent extraction is instrumental in extraction of critical materials and recycling, yet, diluent effects on extraction performance are not well understood. Experimental and empirical solvent parameters have been proposed to correlate with extraction performance, but are often inadequate predictors. We follow the hypothesis that the diluents' primary influence on extraction efficiency is whether or not it hinders assembly of the bulky extracting ligands into a geometry necessary for metal complexation. This behavior is readily accessible with molecular dynamics (MD), where the atomistic description of molecules can be applied to arbitrary extractant-solvent molecules and their mixtures. Several simulated quantities are considered, from both pairwise and graph theoretical analyses, and compared to experimental distribution ratio data for americium extraction by TODGA in a series of inert, non-interacting diluents. These simple properties, especially the formation of closed triplets corresponding to the 3 : 1 ligand : metal stoichiometric solvate, suggest a potential predictive power of this approach. This methodology provides a path forward to comprehensively understand and predict diluent effects in more complex systems involving different extracting ligands and multi-component diluent mixtures.

Modeling Henry's law and phase separations of water-NaCl-organic mixtures with solvation and ion-pairing

Empirical measurements of solution vapor pressure of ternary acetonitrile (MeCN) H2O-NaCl-MeCN mixtures were recorded, with NaCl concentrations ranging from zero to the saturation limit, and MeCN concentrations ranging from zero to an absolute mole fraction of 0.64. After accounting for speciation, the variability of the Henry's law coefficient at vapor-liquid equilibrium (VLE) of MeCN ternary mixtures decreased from 107% to 5.1%. Solute speciation was modeled using a mass action solution model that incorporates solute solvation and ion-pairing phenomena. Two empirically determined equilibrium constants corresponding to solute dissociation and ion pairing were utilized for each solute. When speciation effects were considered, the solid-liquid equilibrium of H2O-NaCl-MeCN mixtures appear to be governed by a simple saturation equilibrium constant that is consistent with the binary H2O-NaCl saturation coefficient. Further, our results indicate that the precipitation of NaCl in the MeCN ternary mixtures was not governed by changes in the dielectric constant. Our model indicates that the compositions of the salt-induced liquid-liquid equilibrium (LLE) boundary of the H2O-NaCl-MeCN mixture correspond to the binary plateau activity of MeCN, a range of concentrations over which the activity remains largely invariant in the binary water-MeCN system. Broader comparisons with other ternary miscible organic solvent (MOS) mixtures suggest that salt-induced liquid-liquid equilibrium exists if: (1) the solution displays a positive deviation from the ideal limits governed by Raoult's law; and (2) the minimum of the mixing free energy profile for the binary water-MOS system is organic-rich. This work is one of the first applications of speciation-based solution models to a ternary system, and the first that includes an organic solute.

Dynamic Community Detection Decouples Multiple Time Scale Behavior of Complex Chemical Systems

Although community or cluster identification is becoming a standard tool within the simulation community, traditional algorithms are challenging to adapt to time-dependent data. Here, we introduce temporal community identification using the Δ-screening algorithm, which has the flexibility to account for varying community compositions, merging and splitting behaviors within dynamically evolving chemical networks. When applied to a complex chemical system whose varying chemical environments cause multiple time scale behavior, Δ-screening is able to resolve the multiple time scales of temporal communities. This computationally efficient algorithm is easily adapted to a wide range of dynamic chemical systems; flexibility in implementation allows the user to increase or decrease the resolution of temporal features by controlling parameters associated with community composition and fluctuations therein.

Hierarchical phenomena in multicomponent liquids: simulation methods, analysis, chemistry

Complex, multicomponent, solutions have often been studied solely through the lens of specific applications of interest. Yet advances to both simulation methodologies (enhanced sampling, etc.) and analysis techniques (network analysis algorithms and others), are creating a trove of data that reveal transcending characteristics across vast compositional phase space. This perspective discusses technical considerations of the reliable and accurate simulations of complex solutions, followed by the advances to analysis algorithms that elucidate coupling of different length and timescale behavior (hierarchical phenomena). The different manifestations of hierarchical phenomena are presented across an array of solution environments, emphasizing fundamental and ongoing science questions. With a more advanced molecular understanding in hand, a quintessential application (solvent extraction) is discussed, where significant opportunities exist to re-imagine the technical scope of an established technology.

Tuning Water Networks via Ionic Liquid/Water Mixtures

Water in nanoconfinement is ubiquitous in biological systems and membrane materials, with altered properties that significantly influence the surrounding system. In this work, we show how ionic liquid (IL)/water mixtures can be tuned to create water environments that resemble nanoconfined systems. We utilize molecular dynamics simulations employing ab initio force fields to extensively characterize the water structure within five different IL/water mixtures: [BMIM + ][BF 4 - ], [BMIM + ][PF 6 - ], [BMIM + ][OTf - ], [BMIM + ][NO 3 - ]and [BMIM + ][TFSI - ] ILs at varying water fraction. We characterize water clustering, hydrogen bonding, water orientation, pairwise correlation functions and percolation networks as a function of water content and IL type. The nature of the water nanostructure is significantly tuned by changing the hydrophobicity of the IL and sensitively depends on water content. In hydrophobic ILs such as [BMIM + ][PF 6 - ], significant water clustering leads to dynamic formation of water pockets that can appear similar to those formed within reverse micelles. Furthermore, rotational relaxation times of water molecules in supersaturated hydrophobic IL/water mixtures indicate the close-connection with nanoconfined systems, as they are quantitatively similar to water relaxation in previously characterized lyotropic liquid crystals. We expect that this physical insight will lead to better design principles for incorporation of ILs into membrane materials to tune water nanostructure.

Continuous noninvasive glucose monitoring; water as a relevant marker of glucose uptake in vivo

With diabetes set to become the number 3 killer in the Western hemisphere and proportionally growing in other parts of the world, the subject of noninvasive monitoring of glucose dynamics in blood remains a "hot" topic, with the involvement of many groups worldwide. There is a plethora of techniques involved in this academic push, but the so-called multisensor system with an impedance-based core seems to feature increasingly strongly. However, the symmetrical structure of the glucose molecule and its shielding by the smaller dipoles of water would suggest that this option should be less enticing. Yet there is enough phenomenological evidence to suggest that impedance-based methods are truly sensitive to the biophysical effects of glucose variations in the blood. We have been trying to answer this very fundamental conundrum: "Why is impedance or dielectric spectroscopy sensitive to glucose concentration changes in the blood and why can this be done over a very broad frequency band, including microwaves?" The vistas for medical diagnostics are very enticing. There have been a significant number of papers published that look seriously at this problem. In this review, we want to summarize this body of research and the underlying mechanisms and propose a perspective toward utilizing the phenomena. It is our impression that the current world view on the dielectric response of glucose in solution, as outlined below, will support the further evolution and implementation toward practical noninvasive glucose monitoring solutions.

Widom line, dynamical crossover, and percolation transition of supercritical oxygen via molecular dynamics simulations

Supercritical oxygen, a cryogenic fluid, is widely used as an oxidizer in jet propulsion systems and is therefore of paramount importance in gaining physical insights into processes such as transcritical and supercritical vaporization. It is well established in the scientific literature that the supercritical state is not homogeneous but, in fact, can be demarcated into regions with liquid-like and vapor-like properties, separated by the "Widom line." In this study, we identified the Widom line for oxygen, constituted by the loci of the extrema of thermodynamic response functions (heat capacity, volumetric thermal expansion coefficient, and isothermal compressibility) in the supercritical region, via atomistic molecular dynamics simulations. We found that the Widom lines derived from these response functions all coincide near the critical point until about 25 bars and 15-20 K, beyond which the isothermal compressibility line begins to deviate. We also obtained the crossover from liquid-like to vapor-like behavior of the translational diffusion coefficient, shear viscosity, and rotational relaxation time of supercritical oxygen. While the crossover of the translational diffusion coefficient and shear viscosity coincided with the Widom lines, the rotational relaxation time showed a crossover that was largely independent of the Widom line. Further, we characterized the clustering behavior and percolation transition of supercritical oxygen molecules, identified the percolation threshold based on the fractal dimension of the largest cluster and the probability of finding a cluster that spans the system in all three dimensions, and found that the locus of the percolation threshold also coincided with the isothermal compressibility Widom line. It is therefore clear that supercritical oxygen is far more complex than originally perceived and that the Widom line, dynamical crossovers, and percolation transitions serve as useful routes to better our understanding of the supercritical state.

What is the primary mover of water dynamics?

We present evidence that the microscopic origin of both the excess wing and the main relaxation process of pure water is the same.

Two-dimensional percolation at the free water surface and its relation with the surface tension anomaly of water

The percolation temperature of the lateral hydrogen bonding network of the molecules at the free water surface is determined by means of molecular dynamics computer simulation and identification of the truly interfacial molecules analysis for six different water models, including three, four, and five site ones. The results reveal that the lateral percolation temperature coincides with the point where the temperature derivative of the surface tension has a minimum. Hence, the anomalous temperature dependence of the water surface tension is explained by this percolation transition. It is also found that the hydrogen bonding structure of the water surface is largely model-independent at the percolation threshold; the molecules have, on average, 1.90 ± 0.07 hydrogen bonded surface neighbors. The distribution of the molecules according to the number of their hydrogen bonded neighbors at the percolation threshold also agrees very well for all the water models considered. Hydrogen bonding at the water surface can be well described in terms of the random bond percolation model, namely, by the assumptions that (i) every surface water molecule can form up to 3 hydrogen bonds with its lateral neighbors and (ii) the formation of these hydrogen bonds occurs independently from each other.

Molecular structure and transport dynamics in perfluoro sulfonyl imide membranes

We report a detailed and comprehensive analysis from classical molecular dynamics simulations of the nanostructure of a model of hydrated perfluoro sulfonyl imide (PFSI) membrane, a polymeric system of interest as a proton conductor in polymer electrolyte membrane fuel cells. We also report on the transport dynamics of water and hydronium ions, and water network percolation in this system. We find that the water network percolation threshold for PFSI, i.e. the threshold at which a consistent spanning water network starts to develop in the membrane, is found to occur between hydration levels (λ) 6 and 7. The higher acidity of the sulfonyl imide acid group of PFSI compared to the sulfonic acid group in Nafion, as computationally characterized in our earlier ab initio study (Idupulapati et al 2010 J. Phys. Chem. A 114 6904-12), results in a larger fraction of 'free' hydronium ions at low hydration levels in PFSI compared to Nafion. However, the calculated diffusion coefficients of the H(3)O(+) ions and H(2)O molecules as a function the hydration level are observed to be almost the same as that of Nafion, indicating similar conductivity and consistent with experimental data.

Atomistic simulation of water percolation and proton hopping in Nafion fuel cell membrane

We have performed a detailed analysis of water clustering and percolation in hydrated Nafion configurations generated by classical molecular dynamics simulations. Our results show that at low hydration levels H(2)O molecules are isolated and a continuous hydrogen-bonded network forms as the hydration level is increased. Our quantitative analysis has established a hydration level (λ) between 5 and 6 H(2)O/SO(3)(-) as the percolation threshold of Nafion. We have also examined the effect of such a network on proton transport by studying the structural diffusion of protons using the quantum hopping molecular dynamics method. The mean residence time of the proton on a water molecule decreases by 2 orders of magnitude when the λ value is increased from 5 to 15. The proton diffusion coefficient in Nafion at a λ value of 15 is about 1.1 × 10(-5) cm(2)/s in agreement with experiment. The results provide quantitative atomic-level evidence of water network percolation in Nafion and its effect on proton conductivity.

Percolation Threshold of Water in Ideal Binary Mixture

Percolation transitions of water and solute molecules in an ideal completely miscible aqueous solution are studied by computer simulations. Three concentration ranges with different kinds of water and solute clustering can be distinguished. An infinite water network, present in pure liquid water, breaks at some solute concentration. Accordingly, the infinite network of solute molecules, present in pure liquid solute, breaks upon addition of water. At ambient and high temperatures, there is a concentration range were both components are below their respective percolation thresholds and there are no water or solute infinite networks. The existence of such concentration range is a characteristic of a completely miscible binary mixture. Upon supercooling, the threshold concentrations decrease and simultaneous existence of the infinite networks of water and of solute becomes possible. The presence of two inter-penetrating infinite networks of molecules may be one of the conditions required for the liquid-liquid transition of one-component isotropic fluid.

Demixing transition of the aqueous solution of amyloidogenic peptides: a REMD simulation study

The aggregation of amyloidogenic peptides in liquid water is studied at various temperatures by replica exchange molecular dynamics (REMD) simulations. The formation of a peptide aggregate upon decreasing the temperature reveals features typical for a first-order demixing phase transition, which is smeared out due to the finite size of the simulation box. Various properties of the ensemble of peptides were used to describe the temperature-induced demixing phase transition, which was found to occur at about 375 K. The hydrational and volumetric properties of the peptides and their aggregates are analyzed.

Percolation and three-dimensional structure of supercritical water

It is well established that at ambient and supercooled conditions water can be described as a percolating network of H bonds. This work is aimed at identifying, by neutron diffraction experiments combined with computer simulations, a percolation line in supercritical water, where the extension of the H-bond network is in question. It is found that in real supercritical water liquidlike states are observed at or above the percolation threshold, while below this threshold gaslike water forms small, sheetlike configurations. Inspection of the three-dimensional arrangement of water molecules suggests that crossing of this percolation line is accompanied by a change of symmetry in the first neighboring shell of molecules from trigonal below the line to tetrahedral above.

Percolation transition in fluids: scaling behavior of the spanning probability functions

We examine different spanning probability functions (wrapping and crossing) near the percolation threshold of a supercritical square-well fluid and determine the threshold values of these probabilities, which may be universal for all fluids. It is shown that for a continuous system, over a wide range of system size, the wrapping probabilities can be described by universal scaling functions, whereas the crossing probabilities do not show such universal behavior over the same range of system size. The obtained universal functions for the wrapping probabilities can be used for an estimation of the percolation threshold in fluids in general. The results for the crossing probabilities allow us then to characterize large clusters in real fluids.

Percolation Transition in Supercritical Water: A Monte Carlo Simulation Study

Computer simulations of water have been performed on the canonical ensemble at 15 different molecular number densities, ranging from 0.006 to 0.018 A-3, along the supercritical isotherm of 700 K, in order to characterize the percolation transition in the system. It is found that the percolation transition occurs at a somewhat higher density than what is corresponding to the supercritical extension of the boiling line. We have shown that the fractal dimension of the largest cluster and the probability of finding a spanning cluster are the most appropriate properties for the location of the true percolation threshold. Thus, percolation transition occurs when the fractal dimension of the largest cluster reaches 2.53, and the probability of finding a cluster that spans the system in at least one dimension and in all the three dimensions reaches 0.97 and 0.65, respectively. On the other hand, the percolation threshold cannot be accurately located through the cluster size distribution, as it is distorted by appearance of clusters crossing the finite simulated system even far below the percolation threshold. The structure of the largest water cluster is dominated by a linear, chainlike arrangement, which does not change noticeably until the largest cluster becomes infinite.

Effect of confinement on the liquid-liquid phase transition of supercooled water

We report on an observation of the phase transition between two liquid phases of supercooled confined water in simulations. The temperature of the liquid-liquid transition of water at zero pressure slightly decreases due to confinement in the hydrophobic pore. The hydrophilic confinement affects this temperature in the opposite direction and shifts the critical point of the liquid-liquid transition to a higher pressure. As a result, in a strongly hydrophilic pore the liquid-liquid phase transition becomes continuous at zero pressure, indicating the shift of its critical point from negative to a positive pressure. These findings indicate that experimental studies of water confined in the pores of various hydrophobicity/hydrophilicity may clarify the location of the liquid-liquid critical point of bulk water.