Search

Thermophoresis has emerged as a powerful tool for characterizing and manipulating colloids at the nano- and micro-scales due to its sensitivity to colloid–solvent interactions. The use of surfactants enables the tailoring of surface chemistry on colloidal particles and the tuning of interfacial interactions. However, the microscopic mechanisms underlying thermophoresis in surfactant solutions remain poorly understood due to the complexity of multiscale interaction coupling. To achieve a more fundamental understanding of the roles of surfactants, we investigated the thermophoretic behavior of silica beads in both ionic and nonionic surfactant solutions at various background temperatures. We provide a complete mechanistic picture of the effects of surfactants on interfacial interactions through mode-coupling analysis of both electrophoretic and thermophoretic experiments. Our results demonstrate that silica thermophoresis is predominantly governed by the dissociation of silanol functional groups at silica–water interfaces in nonionic surfactant solutions, while in ionic surfactant solutions, the primary mechanism driving silica thermophoresis is the adsorption of ionic surfactants onto the silica surface. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

In recent years, ionic microgels have garnered much attention due to their unique properties, especially their stimulus-sensitive swelling behavior. The tunable response of these soft, permeable, compressible, charged colloidal particles is increasingly attractive for applications in medicine and biotechnologies, such as controlled drug delivery, tissue engineering, and biosensing. The ability to model and predict variation of the osmotic pressure of a single microgel with respect to changes in particle properties and environmental conditions proves vital to such applications. In this work, we apply both nonlinear Poisson–Boltzmann theory and molecular dynamics simulation to ionic microgels (macroions) in the cell model to compute density profiles of microions (counterions, coions), single-microgel osmotic pressure, and equilibrium swelling ratios of spherical microgels whose fixed charge is confined to the macroion surface. The basis of our approach is an exact theorem that relates the electrostatic component of the osmotic pressure to the microion density profiles. Close agreement between theory and simulation serves as a consistency check to validate our approach. We predict that surface-charged microgels progressively deswell with increasing microgel concentration, starting well below close packing, and with increasing salt concentration, in qualitative agreement with experiments. Comparison with previous results for microgels with fixed charge uniformly distributed over their volume demonstrates that surface-charged microgels deswell more rapidly than volume-charged microgels. We conclude that swelling behavior of ionic microgels in solution is sensitive to the distribution of fixed charge within the polymer-network gel and strongly depends on bulk concentrations of both microgels and salt ions. [ABSTRACT FROM AUTHOR]

Molecular simulations in an open environment, involving ion exchange, are necessary to study various systems, from biosystems to confined electrolytes. However, grand-canonical simulations are often computationally demanding in condensed phases. A promising method [L. Belloni, J. Chem. Phys. 151, 021101 (2019)], one of the hybrid nonequilibrium molecular dynamics/Monte Carlo algorithms, was recently developed, which enables efficient computation of fluctuating number or charge density in dense fluids or ionic solutions. This method facilitates the exchange through an auxiliary dimension, orthogonal to all physical dimensions, by reducing initial steric and electrostatic clashes in three-dimensional systems. Here, we report the implementation of the method in LAMMPS with a Python interface, allowing facile access to grand-canonical molecular dynamics simulations with massively parallelized computation. We validate our implementation with two electrolytes, including a model Lennard-Jones electrolyte similar to a restricted primitive model and aqueous solutions. We find that electrostatic interactions play a crucial role in the overall efficiency due to their long-range nature, particularly for water or ion-pair exchange in aqueous solutions. With properly screened electrostatic interactions and bias-based methods, our approach enhances the efficiency of salt-pair exchange in Lennard-Jones electrolytes by approximately four orders of magnitude, compared to conventional grand-canonical Monte Carlo. Furthermore, the acceptance rate of NaCl-pair exchange in aqueous solutions at moderate concentrations reaches about 3% at the maximum efficiency. [ABSTRACT FROM AUTHOR]

For industrial applications of self-assembled wormlike micelles, measurement and characterization of a micellar material's microstructure and rheology are paramount for the development and deployment of new high-performing and cost-effective formulations. Within this workflow, there are significant bottlenecks associated with experimental delays and a lack of transferability of results from one chemistry to another. In this work, we outline a process to predict microscopic and thermodynamic characteristics of wormlike micelles directly from rheological data by combining a more robust and efficient fitting algorithm with a recently published constitutive model called the Toy Shuffling model [J. D. Peterson and M. E. Cates, J. Rheol. 64, 1465–1496 (2020) and J. D. Peterson and M. E. Cates, J. Rheol. 65, 633–662 (2021)]. To support this work, linear rheology measurements were taken for 143 samples comprising a common base formulation of commercial sodium lauryl ether sulfate, cocamidopropyl betaine, and salt (NaCl). The steady state zero shear viscosity evident in linear rheology was measured in duplicate via direct steady and oscillatory shear experiments. Fitting the collected data to the model, we found trends in the microstructural and thermodynamic characteristics that agree with molecular dynamics simulations. These trends validate our new perspective on the parameters that inform the study of the relationship between chemical formulation and rheology. This work, when implemented at scale, can potentially be used to inform and test strategies for predicting self-assembled micellar structures based on chemical formulation. [ABSTRACT FROM AUTHOR]

We present the coupling of two frameworks—the pseudo-open boundary simulation method known as constant potential molecular dynamics simulations (CμMD), combined with quantum mechanics/molecular dynamics (QMMD) calculations—to describe the properties of graphene electrodes in contact with electrolytes. The resulting CμQMMD model was then applied to three ionic solutions (LiCl, NaCl, and KCl in water) at bulk solution concentrations ranging from 0.5 M to 6 M in contact with a charged graphene electrode. The new approach we are describing here provides a simulation protocol to control the concentration of electrolyte solutions while including the effects of a fully polarizable electrode surface. Thanks to this coupling, we are able to accurately model both the electrode and solution side of the double layer and provide a thorough analysis of the properties of electrolytes at charged interfaces, such as the screening ability of the electrolyte and the electrostatic potential profile. We also report the calculation of the integral electrochemical double layer capacitance in the whole range of concentrations analyzed for each ionic species, while the quantum mechanical simulations provide access to the differential and integral quantum capacitance. We highlight how subtle features, such as the adsorption of potassium graphene or the tendency of the ions to form clusters contribute to the ability of graphene to store charge, and suggest implications for desalination. [ABSTRACT FROM AUTHOR]

An ionic liquid (IL) is a salt in the liquid state that consists of a cation and an anion, one of which possesses an organic component. Because of their non-volatile property, these solvents have a high recovery rate, and, hence, they are considered as environment-friendly green solvents. It is necessary to study the detailed physicochemical properties of these liquids for designing and processing techniques and find suitable operating conditions for IL-based systems. In the present work, the flow behavior of aqueous solutions of an imidazolium-based IL, 1-methyl-3-octylimidazolium chloride, is investigated, where the dynamic viscosity measurements indicate non-Newtonian shear thickening behavior in the solutions. Polarizing optical microscopy shows that the pristine samples are isotropic and transform into anisotropic after shear. These shear thickened liquid crystalline samples change into an isotropic phase upon heating, which is quantified by the differential scanning calorimetry. The small angle x-ray scattering study revealed that the pristine isotropic cubic phase of spherical micelles distort into non-spherical micelles. This has provided the detailed structural evolution of mesoscopic aggregates of the IL in an aqueous solution and the corresponding viscoelastic property of the solution. [ABSTRACT FROM AUTHOR]

Polarimetric angle-resolved second-harmonic scattering (AR-SHS) is an all-optical tool enabling the study of unlabeled interfaces of nano-sized particles in an aqueous solution. As the second harmonic signal is modulated by interference between nonlinear contributions originating at the particle's surface and those originating in the bulk electrolyte solution due to the presence of a surface electrostatic field, the AR-SHS patterns give insight into the structure of the electrical double layer. The mathematical framework of AR-SHS has been previously established, in particular regarding changes in probing depth with ionic strength. However, other experimental factors may influence the AR-SHS patterns. Here, we calculate the size dependence of the surface and electrostatic geometric form factors for nonlinear scattering, together with their relative contribution to the AR-SHS patterns. We show that the electrostatic term is stronger in the forward scattering direction for smaller particle sizes, while the ratio of the electrostatic to surface terms decreases with increasing size. Besides this competing effect, the total AR-SHS signal intensity is also weighted by the particle's surface characteristics, given by the surface potential Φ0 and the second-order surface susceptibility χ s , 2 2 . The weighting effect is experimentally demonstrated by comparing SiO2 particles of different sizes in NaCl and NaOH solutions of varying ionic strengths. For NaOH, the larger χ s , 2 2 values generated by deprotonation of surface silanol groups prevail over the electrostatic screening occurring at high ionic strengths; however, only for larger particle sizes. This study establishes a better connection between the AR-SHS patterns and surface properties and predicts trends for arbitrarily-sized particles. [ABSTRACT FROM AUTHOR]

The electric conductivity of ionic solutions is well understood at low ionic concentrations of up to a few millimolar but becomes difficult to unravel at higher concentrations that are still common in nature and technological applications. A model for the conductivity at high concentrations was recently put forth for monovalent electrolytes at low electric fields. The model relies on applying a stochastic density-functional theory and using a modified electrostatic pair-potential that suppresses unphysical, short-range electrostatic interactions. Here, we extend the theory to multivalent ions as well as to high electric fields where a deviation from Ohm's law known as the Wien effect occurs. Our results are in good agreement with experiments and recent simulations. [ABSTRACT FROM AUTHOR]

We have studied experimentally the electrical conductance–voltage curves of negatively and positively charged conical nanopores bathed in ionic solutions with monovalent, divalent, and trivalent cations at electrochemically and biologically relevant ionic concentrations. To better understand the interaction between the pore surface charge and the mobile ions, both single salts and salt mixtures have been considered. We have paid attention to the effects on the conductance of the cation valency, the pore charge asymmetry, and the pore charge inversion phenomena due to trivalent ions, both in single salts and salt mixtures. In addition, we have described how small concentrations of multivalent ions can tune the nanopore conductance due to monovalent majority ions, together with the effect of these charges on the additivity of ionic conductance and fluoride-induced negative differential conductance phenomena. This compilation and discussion of previously presented experimental data offers significant insights on the interaction between fixed and mobile charges confined in nanoscale volumes and should be useful in establishing and checking new models for describing ionic transport in the vicinity of charged surfaces. [ABSTRACT FROM AUTHOR]

Dynamic fluctuations in the hydrogen-bond network of water occur from femto- to nanosecond timescales and provide insight into the structural/dynamical aspects of water at ion–water interfaces. Employing terahertz spectroscopy assisted with molecular dynamics simulations, we study aqueous chloride solutions of five monovalent cations, namely, Li, Na, K, Rb, and Cs. We show that ions modify the behavior of the surrounding water molecules and form interfacial layers of water around them with physical properties distinct from those of bulk water. Small cations with high charge densities influence the kinetics of water well beyond the first solvation shell. At terahertz frequencies, we observe an emergence of fast relaxation processes of water with their magnitude following the ionic order Cs > Rb > K > Na > Li, revealing an enhanced population density of weakly coordinated water at the ion–water interface. The results shed light on the structure breaking tendency of monovalent cations and provide insight into the properties of ionic solutions at the molecular level. [ABSTRACT FROM AUTHOR]

We characterize a concentrated 7.3 m CaCl2 solution, combining neutron diffraction with chloride isotopic substitution (Cl-NDIS) in null water and molecular dynamics (MD) simulations. We elucidate the solution structure, thermodynamic properties, and extent of ion pairing previously suggested as concentration-dependent and often not observed at lower concentrations. Our Cl-NDIS measurements designate the solvent-shared ion pairing as dominant and the contact ion pairing (CIP) as insignificant even under conditions close to the solubility limit. The MD models parameterized against neutron diffraction with calcium isotopic substitution (Ca-NDIS) overestimate CIP despite successfully reproducing most of the Cl-NDIS signal. This drawback originates from the fact that Ca2+–Cl− interactions were primarily "hidden" in the Ca-NDIS signal due to overlapping with Ca2+–Ow and Ca2+–Hw contributions to the total scattering. Contrary, MD models with moderate CIP and possessing generally good performance at high concentrations fail to reproduce the NDIS measurements accurately. Therefore, the electronic polarization, introduced in most of the recent MD models via scaling ionic charges, resolves some but not all parameterization drawbacks. We conclude that despite improving the quality of MD models "on average," the question "which model is the best" has not been answered but replaced by the question "which model is better for a given research." An overall "good" model can still be inappropriate or, in some instances, "bad" and, unfortunately, produce erroneous results. The accurate interpretation of several NDIS datasets, complemented by MD simulations, can prevent such mistakes and help identify the strengths, weaknesses, and convenient applications for corresponding computational models. [ABSTRACT FROM AUTHOR]

Thermodiffusion is the migration of a species due to a temperature gradient and is the driving phenomenon in many applications ranging from early cancer detection to uranium enrichment. Molecular dynamics (MD) simulations can be a useful tool for exploring the rather complex thermodiffusive behavior of species, such as proteins and ions. However, current MD models of thermodiffusion in aqueous ionic solutions struggle to quantitatively predict the Soret coefficient, which indicates the magnitude and direction of species migration under a temperature gradient. In this work, we aim to improve the accuracy of MD thermodiffusion models by assessing how well different water models can recreate thermodiffusion in a benchmark aqueous NaCl solution. We tested four of the best available rigid non-polarizable water models (TIP3P-FB, TIP4P-FB, OPC3, and OPC) and the commonly used TIP3P and SPC/E water models for their ability to predict the inversion temperature and Soret coefficient in 0.5, 2, and 4M aqueous NaCl solutions. Each water model predicted a noticeably different ion distribution yielding different inversion temperatures and magnitudes of the Soret coefficient. By comparing the modeled Soret coefficients to published experimental values, we determine TIP3P-FB to be the water model that best recreates thermodiffusion in aqueous NaCl solutions. Our findings can aid future works in selecting the most accurate rigid non-polarizable water model, including water and ion parameters for investigating thermodiffusion through MD simulations. [ABSTRACT FROM AUTHOR]

When selecting a solvent for a given solute, the strongly held idiom "like dissolves like", meaning that polar solvents are used for polar solutes, is often used. This idea has resulted from the concept that most molecular solvents are homogeneous. In a deep eutectic solvent (DES), however, both components can be ionic or non-ionic, polar or non-polar. By tuning the components, DESs can solubilize a wide variety of solutes, often mixing hydrophobic and hydrophilic components, and the mixture can be designed to control phase behavior. The liquids often contain significant short-length order, and preferential solvation of one component often occurs. The addition of small polar molecules such as water or alcohols results in non-homogeneous liquids, which have significantly decreased viscosity and increased ionic conductivity. Accordingly, the areas covered in this special issue focus on structure and dynamics, solvation, the mobility of charged species, and the ability to obtain controllable phase behavior by adding polar diluents or using hydrophobic DESs. [ABSTRACT FROM AUTHOR]

The most intriguing hypothesis explaining many water anomalies is a metastable liquid–liquid phase transition (LLPT) at high pressure and low temperatures, experimentally hidden by homogeneous nucleation. Recent infrared spectroscopic experiments showed that upon addition of hydrazinium trifluoroacetate to water, the supercooled ionic solution undergoes a sharp, reversible LLPT at ambient pressure, possible offspring of that in pure water. Here, we calculate the temperature-dependent signature of the OH-stretching band, reporting on the low/high density phase of water, in neat water and in the same experimentally investigated ionic solution. The comparison between the infrared signature of the pure liquid and that of the ionic solution can be achieved only computationally, providing insight into the nature of the experimentally observed phase transition and allowing us to investigate the effects of ionic compounds on the high to low density supercooled liquid water transition. We show that the experimentally observed crossover behavior in the ionic solution can be reproduced only if the phase transition between the low- and high-density liquid states of water is coupled to a mixing–unmixing transition between the water component and the ions: at low temperatures, water and ions are separated and the water component is a low density liquid. At high temperatures, water and ions get mixed and the water component is a high-density liquid. The separation at low temperatures into ion-rich and ion-poor regions allows unveiling the polyamorphic nature of liquid water, leading to a crossover behavior resembling that observed in supercooled neat water under high pressure. [ABSTRACT FROM AUTHOR]

Living organisms can sense extracellular forces via mechanosensitive ion channels, which change their channel conformations in response to external pressure and regulate ion transport through the cell membrane. Such pressure-regulated ion transport is critical for various biological processes, such as cellular turgor control and hearing in mammals, but has yet to be achieved in artificial systems using similar mechanisms. In this work, we construct a nanoconfinement by reversibly blocking a single nanopore with a nanoparticle and report anomalous and ultra-mechanosensitive ionic transport across the resulting nanoconfinement upon assorted mechanical and electrical stimuli. Our observation reveals a suppressed ion conduction through the system as the applied pressure increases, which imitates certain behaviors of stretch-inactivated ion channels in biological systems. Moreover, pressure-induced ionic current rectification is also observed despite the high ionic concentration of the solution. Using a combined experimental and simulation study, we correlate both phenomena to pressure-induced nanoparticle rotation and the resulting physical structure change in the blocked nanopore. This work presents a mechanosensitive nano-confinement requiring minimal fabrication techniques and provides new opportunities for bio-inspired nanofluidic applications. [ABSTRACT FROM AUTHOR]

Interactions among ions and their specific interactions with macromolecular solutes are known to play a central role in biomolecular stability. However, similar effects in the conformational stability of protein loops that play functional roles, such as binding ligands, proteins, and DNA/RNA molecules, remain relatively unexplored. A well-characterized enzyme that has such a functional loop is Escherichia coli dihydrofolate reductase (ecDHFR), whose so-called M20 loop has been observed in three ordered conformations in crystal structures. To explore how solution ionic strengths may affect the M20 loop conformation, we proposed a reaction coordinate that could quantitatively describe the loop conformation and used it to classify the loop conformations in representative ecDHFR x-ray structures crystallized in varying ionic strengths. The Protein Data Bank survey indicates that at ionic strengths (I) below the intracellular ion concentration-derived ionic strength in E. coli (I ≤ 0.237M), the ecDHFR M20 loop tends to adopt open/closed conformations, and rarely an occluded loop state, but when I is >0.237M, the loop tends to adopt closed/occluded conformations. Distance-dependent electrostatic potentials around the most mobile M20 loop region from molecular dynamics simulations of ecDHFR in equilibrated CaCl2 solutions of varying ionic strengths show that high ionic strengths (I = 0.75/1.5M) can preferentially stabilize the loop in closed/occluded conformations. These results nicely correlate with conformations derived from ecDHFR structures crystallized in varying ionic strengths. Altogether, our results suggest caution in linking M20 loop conformations derived from crystal structures solved at ionic strengths beyond that tolerated by E. coli to the ecDHFR function. [ABSTRACT FROM AUTHOR]

Thermodiffusion in liquids (the Soret effect) has several unusual properties. In particular, transport can occur with or against a temperature gradient depending on the case. Numerous empirical correlations have been proposed with mixed success or range of applicability. Here, we show that physicochemical mechanics, derived from the Smoluchowski equation as a description of diffusive transport phenomena, is in accord with the experimental and simulated thermodiffusion data from colloidal beads and biomacromolecules to ionic solutions and ultracold fluid mixtures. It yields a simple formula for the Soret coefficient ST based on the reference molar entropy including non-ideality. Hydrodynamic and local non-equilibrium effects are discussed but not included as these are apparently not a major contribution for the wide range of solutes under the near-equilibrium experimental conditions considered here. [ABSTRACT FROM AUTHOR]

This work applies a molecular theory to study the formation of lateral self-assembled aggregates in mixed brushes composed of polyanion and polycation chains. In order to overcome the well-known limitations of mean-field electrostatics to capture polyelectrolyte complexation, the formation of ion pairs between anionic and cationic groups in the polyelectrolytes is explicitly modeled in our theory as an association reaction. This feature is essential to capture the microphase separation of the mixed brush and the formation of lateral aggregates triggered by polyelectrolyte complexation. The effects of solution pH and ionic strength, surface coverage, and chain length on the morphology of the mixed brush are systematically explored. It is shown that increasing salt concentration leads to the rupture of polyelectrolyte complexes and the stabilization of the homogeneous, non-aggregated brush, providing that the formation of ion pairs between the polyelectrolytes and the salt ions in solution is explicitly accounted for by the theory. The inclusion of ion-pairing association reactions between oppositely charged polyelectrolytes within a mean-field description of electrostatics emerges from this work as a useful and simple theoretical approach to capture the formation of polyelectrolyte complexes and their responsiveness to solution ionic strength and pH. [ABSTRACT FROM AUTHOR]

In this work, a force field for several ions in water is proposed. In particular, we consider the cations Li+, Na+, K+, Mg2+, and Ca2+ and the anions Cl− and S O 4 2 − . These ions were selected as they appear in the composition of seawater, and they are also found in biological systems. The force field proposed (denoted as Madrid-2019) is nonpolarizable, and both water molecules and sulfate anions are rigid. For water, we use the TIP4P/2005 model. The main idea behind this work is to further explore the possibility of using scaled charges for describing ionic solutions. Monovalent and divalent ions are modeled using charges of 0.85 and 1.7, respectively (in electron units). The model allows a very accurate description of the densities of the solutions up to high concentrations. It also gives good predictions of viscosities up to 3 m concentrations. Calculated structural properties are also in reasonable agreement with the experiment. We have checked that no crystallization occurred in the simulations at concentrations similar to the solubility limit. A test for ternary mixtures shows that the force field provides excellent performance at an affordable computer cost. In summary, the use of scaled charges, which could be regarded as an effective and simple way of accounting for polarization (at least to a certain extend), improves the overall description of ionic systems in water. However, for purely ionic systems, scaled charges will not adequately describe neither the solid nor the melt. [ABSTRACT FROM AUTHOR]

Electrochemistry is a highly versatile part of chemical research which is involved in many of the processes in the field of micromotion. Its input has been crucial from the synthesis of microstructures to the explanation of phoretic mechanisms. However, using electrochemical effects to propel artificial micromotors is still to be achieved. Here, we show that the forces generated by electrochemical reactions can not only create active motion, but they are also strong enough to overcome the adhesion to the substrate, caused by the increased ionic strength of the solutions containing the ions of more noble metals themselves. The galvanic replacement of copper by platinum ions is a spontaneous process, which not only provides a sufficiently strong electromotive force to propel the Janus structures but also results in asymmetric Pt-hatted structures, which can be further used as catalytic micromotors. [ABSTRACT FROM AUTHOR]

We investigate using molecular dynamics simulations the dynamical and structural properties of LiCl:6H2O aqueous solution upon supercooling. This ionic solution is a glass forming liquid of relevant interest in connection with the study of the anomalies of supercooled water. The LiCl:6H2O solution is easily supercooled and the liquid state can be maintained over a large decreasing temperature range. We performed simulations from ambient to 200 K in order to investigate how the presence of the salt modifies the behavior of supercooled water. The study of the relaxation time of the self-density correlation function shows that the system follows the prediction of the mode coupling theory and behaves like a fragile liquid in all the range explored. The analysis of the changes in the water structure induced by the salt shows that while the salt preserves the water hydrogen bonds in the system, it strongly affects the tetrahedral hydrogen bond network. Following the interpretation of the anomalies of water in terms of a two-state model, the modifications of the oxygen radial distribution function and the angular distribution function of the hydrogen bonds in water indicate that LiCl has the role of enhancing the high density liquid component of water with respect to the low density component. This is in agreement with recent experiments on aqueous ionic solutions. [ABSTRACT FROM AUTHOR]

We develop a nonlocal Nernst-Planck model for reaction and diffusion in multicomponent ionic systems. We apply the model to the one-dimensional liquid junction problem, in which two electrolytic solutions of different ionic concentrations are brought into contact via a permeable membrane. Transport of ions through the membrane induces an electric field which is modeled using two separate nonlocal conditions: charge conservation and Gauss' law. We inve1064s that the net charge at each point in the domain is zero, by considering four different initial scenarios. Charge conservation and Gauss' law yield similar results for most practical scenarios in which the initial condition satisfies strict electroneutrality. However, Gauss' law has two important advantages over charge conservation: (i) it is numerically more stable and can be applied even when the concentration of all the charged species drops to zero and (ii) computationally, it is significantly cheaper. Further, this study provides insights on the prescription of electroneutrality conditions necessary to handle the physics of evolving charges in nonlocal peridynamic models that are aimed at modeling nonlocal reaction-diffusion or corrosion-type processes. [ABSTRACT FROM AUTHOR]

Despite considerable efforts over more than two decades, our knowledge of the interactions in electrolyte solutions is not yet satisfactory. Not even one of the most simple and important aqueous solutions, NaCl(aq), escapes this assertion. A requisite for the development of a force field for any water solution is the availability of a good model for water. Despite the fact that TIP4P/2005 seems to fulfill the requirement, little work has been devoted to build a force field based on TIP4P/2005. In this work, we try to fill this gap for NaCl(aq). After unsuccessful attempts to produce accurate predictions for a wide range of properties using unity ionic charges, we decided to follow recent suggestions indicating that the charges should be scaled in the ionic solution. In this way, we have been able to develop a satisfactory non-polarizable force field for NaCl(aq). We evaluate a number of thermodynamic properties of the solution (equation of state, maximum in density, enthalpies of solution, activity coefficients, radial distribution functions, solubility, surface tension, diffusion coefficients, and viscosity). Overall the results for the solution are very good. An important achievement of our model is that it also accounts for the dynamical properties of the solution, a test for which the force fields so far proposed failed. The same is true for the solubility and for the maximum in density where the model describes the experimental results almost quantitatively. The price to pay is that the model is not so good at describing NaCl in the solid phase, although the results for several properties (density and melting temperature) are still acceptable. We conclude that the scaling of the charges improves the overall description of NaCl aqueous solutions when the polarization is not included. [ABSTRACT FROM AUTHOR]

Ionic solutions are often regarded as fully dissociated ions dispersed in a polar solvent. While this picture holds for dilute solutions, at higher ionic concentrations, oppositely charged ions can associate into dimers, referred to as Bjerrum pairs. We consider the formation of such pairs within the nonlinear Poisson-Boltzmann framework and investigate their effects on bulk and interfacial properties of electrolytes. Our findings show that pairs can reduce the magnitude of the dielectric decrement of ionic solutions as the ionic concentration increases. We describe the effect of pairs on the Debye screening length and relate our results to recent surface-force experiments. Furthermore, we show that Bjerrum pairs reduce the ionic concentration in bulk electrolyte and at the proximity of charged surfaces, while they enhance the attraction between oppositely charged surfaces. [ABSTRACT FROM AUTHOR]

The motion of ionic solutes and charged particles under the influence of an electric field and the ensuing hydrodynamic flow of the underlying solvent is ubiquitous in aqueous colloidal suspensions. The physics of such systemsis described by a coupled set of differential equations, along with boundary conditions, collectively referred to as the electrokinetic equations. Capuani et al. [J. Chem. Phys. 121, 973 (2004)] introduced a lattice-based method for solving this system of equations, which builds upon the lattice Boltzmann algorithm for the simulation of hydrodynamic flow and exploits computational locality. However, thus far, a description of how to incorporate moving boundary conditions into the Capuani scheme has been lacking. Moving boundary conditions are needed to simulate multiple arbitrarily moving colloids. In this paper, we detail how to introduce such a particle coupling scheme, based on an analogue to the moving boundary method for the pure lattice Boltzmann solver. The key ingredients in our method are mass and charge conservation for the solute species and a partial-volume smoothing of the solute fluxes to minimize discretization artifacts. We demonstrate our algorithm's effectiveness by simulating the electrophoresis of charged spheres in an external field; for a single sphere we compare to the equivalent electro-osmotic (co-moving) problem. Our method's efficiency and ease of implementation should prove beneficial to future simulations of the dynamics in a wide range of complex nanoscopic and colloidal systems that were previously inaccessible to lattice-based continuum algorithms. [ABSTRACT FROM AUTHOR]

We employ a field-theoretical variational approach to study the behavior of ionic solutions in the grand canonical ensemble. To describe properly the hardcore interactions between ions, we use a cutoff in Fourier space for the electrostatic contribution of the grand potential and the Carnahan- Starling equation of state with a modified chemical potential for the pressure one. We first calibrate our method by comparing its predictions at room temperature with Monte Carlo results for excess chemical potential and energy. We then validate our approach in the bulk phase by describing the classical "ionic liquid-vapor" phase transition induced by ionic correlations at low temperature, before applying it to electrolytes at room temperature confined to nanopores embedded in a low dielectric medium and coupled to an external reservoir of ions. The ionic concentration in the nanopore is then correctly described from very low bulk concentrations, where dielectric exclusion shifts the transition up to room temperature for sufficiently tight nanopores, to high concentrations where hardcore interactions dominate which, as expected, modify only slightly this ionic "capillary evaporation". [ABSTRACT FROM AUTHOR]

The translational diffusivity of water in solutions of alkali halide salts depends on the identity of ions, exhibiting dramatically different behavior even in solutions of similar salts of NaCl and KCl. The water diffusion coefficient decreases as the salt concentration increases in NaCl. Yet, in KCl solution, it slightly increases and remains above bulk value as salt concentration increases. Previous classical molecular dynamics simulations have failed to describe this important behavior even when polarizable models were used. Here, we show that inclusion of dynamical charge transfer among water molecules produces results in a quantitative agreement with experiments. Our results indicate that the concentration-dependent diffusivity reflects the importance of many-body effects among the water molecules in aqueous ionic solutions. Comparison with quantum mechanical calculations shows that a heterogeneous and extended distribution of charges on water molecules around the ions due to ion-water and also water-water charge transfer plays a very important role in controlling water diffusivity. Explicit inclusion of the charge transfer allows us to model accurately the difference in the concentration-dependent water diffusivity between Na+ and K+ ions in simulations, and it is likely to impact modeling of a wide range of systems for medical and technological applications. [ABSTRACT FROM AUTHOR]

The 2D Raman-THz spectrum of the MgCl2 solution was simulated using the molecular dynamics simulation and the stability matrix method and compared with that of the pure water. The 2D Raman-THz signal provides more information on the ion effects on the collective water motion than the conventional 1D signal. The presence of MgCl2 suppresses the cross peak of water between the hydrogen bond bending and the other intermolecular vibrational mode, which clearly illustrates that the water hydrogen bending motion is affected by the confining effect of the ions. Our theoretical work thus demonstrates that the 2D Raman-THz technique can become a valuable nonlinear vibrational probe for the molecular dynamics in the ionic solutions. [ABSTRACT FROM AUTHOR]

Numerical properties of the smooth particle mesh Ewald (SPME) sum [U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee, and L. G. Pedersen, J. Chem. Phys.103, 8577 (1995)] have been investigated by molecular dynamics simulation of ionic solutions and dipolar fluids. Scaling dependence of execution time on the number of particles at optimal performance have been determined and compared with the corresponding data of the standard Ewald (SE) sum. For both types of systems and over the range from N = 10³ to 105 particles, the SPME sum displays a sub φ(N ln N) complexity, whereas the SE sum possesses an φ(N 3/2) complexity. The breakeven of the simulation times appears at (10³) particles, and the SPME sum is ≈20 times faster than the SE sum at 105 particles. Furthermore, energy truncation error and the energy and force execution time of the reciprocal space evaluation as function of the number of particles and the convergence parameters of the SPME sum have been determined for both types of systems containing up to 106 particles. [ABSTRACT FROM AUTHOR]

In this work the system formed by charged (ionic) microgels in the presence of monovalent salt is investigated by solving numerically the Ornstein-Zernike integral equations within the Hypernetted-Chain approximation. The ionic density profiles, effective interaction between microgel particles, and the effective charge of the particles are calculated. In addition to the electrostatic interaction, the excluded-volume repulsion between the microgel particle and the ions is also explicitly taken into account. Although this steric interaction is irrelevant in the swollen state (when the packing fraction of the polymer network is low), it becomes a very important contribution close to the de-swollen state, hindering the counterion penetration inside the microgel mesh. The theoretical predictions indicate that the ionic density profiles are strongly affected by the degree of swelling, going from a volumetric distribution of counterions in the swollen state to a surface accumulation outside the particle that becomes more important as the particle shrinks. The electrostatic effective interaction between pairs of microgel particles is shown to be the result of a complex interplay between electrostatic and depletion effects that strongly depend on the bare charge density of the particle. For sufficiently charged microgel particles, the steric exclusion leads to a less efficient screening of the microgel charge near the de-swollen configuration, and so to a significant increase of the effective charge of the microgel particle. [ABSTRACT FROM AUTHOR]

In this, the second part of our series on the dielectric spectrum symmetrical broadening of water, we consider ionic aqueous solutions. If in Part I, dipole-dipole interaction was the dominant feature, now ion-dipole interplay is shown to be the critical element in the dipole-matrix interaction. We present the results of high-frequency dielectric measurements of different concentrations of NaCl/KCl aqueous solutions. We observed Cole-Cole broadening of the main relaxation peak of the solvent in the both electrolytes. The 3D trajectory approach (described in detail in Part I) is applied in order to highlight the differences between the dynamics and structure of solutions of salts on one hand and dipolar solutes on the other hand. [ABSTRACT FROM AUTHOR]

The present investigations continue previous research on transport in aqueous ionic solutions through carbon nanotubes. Specifically, the effects of the nanotube radius, solute concentration, and applied external electric fields on the transport properties are investigated in terms of mobilities, currents, and pairing times of the solute ions. The simulated transport features are corroborated with general theoretical results of nanofluidics (such as the linear log-log regime of the nanochannel conductance as function of the solute concentration and the current-voltage curve of the channel). Discontinuities in the partial ionic currents are explained on the basis of a recent theoretical model of quantized ionic conductance in nanopores, developed by Zwolak et al. Correlations between the structural and dynamic properties are established, linking causally the highly structured spatial density profiles, the ion pairing phenomenon and the ionic currents. [ABSTRACT FROM AUTHOR]

The reported work extends previously published research on transport in aqueous ionic solutions through carbon nanotubes. Specifically, the effects of the nanotube radius, solute concentration, and applied external electric fields on the solution structuring are investigated in terms of spatial density distributions, pair distribution functions, and electrostatic potential profiles. Several simulated structural features are consistent with general theoretical results of nanofluidics and can be interpreted fairly well with respect to these (such as the Donnan-type voltages established at the channel apertures depending on the logarithm of the maximum ion concentration). The simulated properties are based on averages over the largest data collection times reported in the literature (0.8 μs), providing accurate estimates of the measured quantities. [ABSTRACT FROM AUTHOR]

We use classical density functional theory (DFT) and an explicit solvent description to investigate the wetting and drying behavior of ionic solutions in contact with a charged solid substrate. The solvent is modeled by dipolar hard spheres, while the monovalent ions are oppositely charged hard spheres; cohesion is ensured by a Yukawa attraction between all three species. The free energy functional describing the inhomogeneous solution includes the best available fundamental measure description of excluded volume correlations in a ternary mixture of hard spheres, whereas all electrostatic and cohesive interactions are treated within the mean-field approximation. We find both first and second order wetting transitions which are rather little affected by ions at low and moderate concentrations, compared to the wetting behavior of the pure solvent. A novel drying scenario is predicted, where complete drying is prevented by the electrostatic attraction between a positively charged substrate and the anions, while near a negatively charged substrate, a first order transition with a predrying line is observed. The various scenarios are surprisingly similar to our previous predictions based on a “semiprimitive” model where the solvent particles carry no dipole, but the ion-ion interactions are reduced by a local dielectric permittivity. [ABSTRACT FROM AUTHOR]

Bulk properties of ionic liquid crystals are investigated using density functional theory. The liquid crystal molecules are represented by ellipsoidal particles with charges located in their center or at their tails. Attractive interactions are taken into account in terms of the Gay–Berne pair potential. Rich phase diagrams involving vapor, isotropic and nematic liquid, as well as smectic phases are found. The dependence of the phase behavior on various parameters such as the length of the particles and the location of charges on the particles is studied. [ABSTRACT FROM AUTHOR]

In this paper, we perform molecular dynamics simulations of an interface containing charged functional groups of different valences in contact with 2:1 ionic solution. We take into account both the finite sizes of the ions in solution and the functional groups but we neglect the structural details of the solvent (primitive model). We show that the distribution of ions and the electrostatic properties of the system depend strongly on the valence of the interfacial charged groups. In the case of surfaces containing well-separated charged interfacial groups, we observe counterion binding at these groups induced by electrostatic interactions. A detailed analysis of the potential of mean force between interfacial charged groups and ions reveals significant features not anticipated by present theories of electrolytes near interfaces. Overall, our results show that, in primitive models of the ion-interface interaction, not only the ionic size and valence are important but the size and valence of the interfacial charged groups also have a significant impact. [ABSTRACT FROM AUTHOR]

Molecular dynamics simulations of ionic solutions depend sensitively on the force fields employed for the ions. To resolve the fine differences between ions of the same valence and roughly similar size and in particular to correctly describe ion-specific effects, it is clear that accurate force fields are necessary. In the past, optimization strategies for ionic force fields either considered single-ion properties (such as the solvation free energy at infinite dilution or the ion-water structure) or ion-pair properties (in the form of ion-ion distribution functions). In this paper we investigate strategies to optimize ionic force fields based on single-ion and ion-pair properties simultaneously. To that end, we simulate five different salt solutions, namely, CsCl, KCl, NaI, KF, and CsI, at finite ion concentration. The force fields of these ions are systematically varied under the constraint that the single-ion solvation free energy matches the experimental value, which reduces the two-dimensional {σ,[variant_greek_epsilon]} parameter space of the Lennard-Jones interaction to a one dimensional line for each ion. From the finite-concentration simulations, the pair potential is extracted and the osmotic coefficient is calculated, which is compared to experimental data. We find a strong dependence of the osmotic coefficient on the force field, which is remarkable as the single-ion solvation free energy and the ion-water structure remain invariant under the parameter variation. Optimization of the force field is achieved for the cations Cs+ and K+, while for the anions I- and F- the experimental osmotic coefficient cannot be reached. This suggests that in the long run, additional parameters might have to be introduced into the modeling, for example by modified mixing rules. [ABSTRACT FROM AUTHOR]

Most aqueous biological and technological systems contain solvated ions. Atomistic explicit-water simulations of ionic solutions rely crucially on accurate ionic force fields, which contain most commonly two adjustable parameters: the Lennard-Jones diameter and the interaction strength. Assuming these parameters to be properly optimized, the plethora of parameters one finds in the literature for one and the same ion is surprising. In principle, the two parameters should be uniquely determined by matching two ionic properties obtained for a particular water model and within a given simulation protocol with the corresponding experimental observables. Traditionally, ion parameters were chosen in a somewhat unsystematic way to reproduce the solvation free energy and to give the correct ion size when compared with scattering results. Which experimental observable one chooses to reproduce should in principle depend on the context within which the ionic force field is going to be used. In the present work we suggest to use the solvation free energy in conjunction with the solvation entropy to construct thermodynamically sound force fields for the alkali and halide ions for the simulation of ion-specific effects in aqueous environment. To that end we determine the solvation free energy and entropy of both cations and anions in the entire relevant parameter space. As an independent check on the quality of the resulting force fields we also determine the effective ionic radius from the first peak of the radial ion-water distribution function. Several difficulties during parameter optimization are discussed in detail. (i) Single-ion solvation depends decisively on water-air surface properties, which experimentally becomes relevant when introducing extrathermodynamic assumptions on the hydronium (H3O+) solvation energy. Fitting ion pairs circumvents this problem but leaves the parameters of one reference ion (here we choose chloride) undetermined. (ii) For the halides the problem is almost underdetermined, i.e., there is a whole set of degenerate parameters that equally well describe, e.g., chloride and bromide ions. (iii) For the heavy cations the problem is overdetermined, i.e., no combination of Lennard-Jones parameters is able to reproduce simultaneously energy and entropy of solvation. We discuss various possibilities to deal with these problems and finally present an optimized force field for the halide anions that reproduces the free energy and the entropy of solvation. For the alkali metal cations there is no unambiguous choice of parameters. Therefore, we give three different parameter sets for every ion with a small, intermediate, or large Lennard-Jones interaction strength, where the Lennard-Jones diameters are optimized to reproduce the solvation free energy. The ionic radius is reproduced with acceptable accuracy by this optimization strategy, meaning that the proposed force fields are reliable beyond the target observables (i.e., free energy and entropy of solvation). [ABSTRACT FROM AUTHOR]

We consider a new qualitative approach for treating theoretically the solvation of single-atomic ionic solutes in binary mixtures of polar and nonpolar aprotic solvents. It is based on the implicit continuum electrostatic model of the solvent mixture involving distance-dependent dielectric permittivity [variant_greek_epsilon](R) (where R is the distance from the ion) and local concentrations C1(R) and C2(R) of the solvent ingredients. For a given R, the condition for local thermodynamic equilibrium provides the transcendental equation for explicitly establishing the permittivity and concentration profiles. Computations performed with real Cl- and model Cl+ ions as solutes in benzene/DMSO mixtures are compared with the molecular dynamics simulations of the same systems. A significant discrepancy of molecular and continuum results is revealed for the concentration profiles in the close vicinity of the ion boundary, although the general trends are similar. The continuum methodology cannot account for the formation of rigid solvent structures around ions, which is most significant for the case of Cl+. Such defect, however, proves to become of less importance in calculations of the solvation free energy, which are quite satisfactory for Cl- ion. Free energy calculations for Cl+ are less successful in the range of low DMSO concentration. [ABSTRACT FROM AUTHOR]

In the present work, we have performed quantum chemical calculations to determine preferable species among the ionic complexes that are present in ambient water due to the autodissociation of water molecule. First, we have formulated the relative population of the hydrated complexes with respect to the bare ion (H3O+ or OH-) in terms of the solvation free energies of the relevant molecules. The solvation free energies for various ionic species (H3O+, H5O2+, H7O3+, H9O4+ or OH-, H3O2-, H5O3-, H7O4-, H9O5-), categorized as proton or hydroxide ion in solution, have been computed by employing the QM/MM-ER method recently developed by combining the quantum mechanical/molecular mechanical (QM/MM) approach with the theory of energy representation (ER). Then, the computed solvation free energies have been used to evaluate the ratio of the populations of the ionic complexes to that of the bare ion (H3O+ or OH-). Our results suggest that the Zundel form, i.e., H5O2+, is the most preferable in the solution among the cationic species listed above though the Eigen form (H9O4+) is very close to the Zundel complex in the free energy, while the anionic fragment from water molecules mostly takes the form of OH-. It has also been found that the loss of the translational entropy of water molecules associated with the formation of the complex plays a role in determining the preferable size of the cluster. [ABSTRACT FROM AUTHOR]

The interaction of water with room temperature ionic liquid (RTIL) [bmim][PF6] has been studied in [bmim][PF6]/TX-100/water ternary microemulsions by solvent and rotational relaxation of coumarin 153 (C-153) and coumarin 490 (C-490). The rotational relaxation and average solvation time of C-153 and C-490 gradually decrease with increase in water content of the microemulsions. The gradual increase in the size of the microemulsion with increase in w0 (w0=[water]/[surfactant]) is evident from dynamic light scattering measurements. Consequently the mobility of the water molecules also increases. In comparison to pure water the retardation of solvation time in the RTIL containing ternary microemulsions is very less. The authors have also reported the solvation time of C-490 in neat [bmim][PF6]. The solvation time of C-490 in neat [bmim][PF6] is bimodal with time constants of 400 ps and 1.10 ns. [ABSTRACT FROM AUTHOR]

A theoretical expression of the complex dielectric constant attributed to space-charge polarization has been derived under an electric field calculated using Poisson’s equation considering the effects of bound charges on ions. The frequency dependence of the complex dielectric constant of chlorobenzene solutions doped with tetrabutylammonium tetraphenylborate (TBATPB) has been analyzed using the theoretical expression, and the impact of the bound charges on the complex dielectric constant has been clarified quantitatively in comparison with a theory that does not consider the effect of the bound charges. The Stokes radius of TBA+(=TPB-) determined by the present theory shows a good agreement with that determined by conductometry in the past; hence, the present theory should be applicable to the direct determination of the mobility of ion species in an electrolytic solution without the need to measure ionic limiting equivalent conductance and transport number. [ABSTRACT FROM AUTHOR]

The relaxation of the electronic spin S of a paramagnetic metal ion with fully quenched orbital angular momentum in its ground state is investigated in an external magnetic field through a systematic study of the time correlation functions governing the evolution of the statistical operator (density matrix). Let ω0 be the Larmor angular frequency of S. When the relaxation is induced by a time-fluctuating perturbing Hamiltonian hH1(t) of time correlation τc, it is demonstrated that after a transient period the standard Redfield approximation is relevant to calculate the evolution of the populations of the spin states if |H1|2τc2/(1+ω02τc2)<1 and that this transient period becomes shorter than τc at sufficiently high field for a zero-field splitting perturbing Hamiltonian. This property, proven analytically and confirmed by numerical simulation, explains the surprising success of several simple expressions of the longitudinal electronic relaxation rate 1/T1e derived from the Redfield approximation well beyond its expected validity range |H1|τc<1. It has favorable practical consequences on the interpretation of the paramagnetic relaxation enhancement of nuclei used for structural and dynamic studies. [ABSTRACT FROM AUTHOR]

A molecular dynamics simulation model for an electroactive interface in which a metallic electrode is maintained at a preset electrical potential is described. The model, based on earlier work of Siepmann and Sprik [J. Chem. Phys. 102, 511 (1995)], uses variable charges whose magnitudes are adjusted on-the-fly according to a variational procedure to maintain the constant potential condition. As such, the model also allows for the polarization of the electrode by the electrolyte, sometimes described by the introduction of image charges. The model has been implemented in a description of an electrochemical cell as a pair of parallel planar electrodes separated by the electrolyte using a two-dimensional Ewald summation method. The method has been applied to examine the interfacial structure in two ionic liquids, consisting of binary mixtures of molten salts, chosen to exemplify the influences of dissimilar cation size and charge. The stronger coordination of the smaller and more highly charged cations by the anions prevents them from approaching even the negatively charged electrode closely. This has consequences for the capacitance of the electrode and will also have an impact on the rates of electron transfer processes. The calculated capacitances exhibit qualitatively the same dependence on the applied potential as has been observed in experimental studies. [ABSTRACT FROM AUTHOR]

We have carried out a series of molecular dynamics simulations to investigate the hydrogen bond and residence dynamics of X-–water (X=F, Cl, and I) and pairs in aqueous solutions at a temperature of 673 K. The calculations are done at six different water densities ranging from 1.0 to 0.15 g cm-3. The hydrogen bonds are defined by using a set of configurational criteria with respect to the anion(oxygen)–oxygen and anion(oxygen)–hydrogen distances and the anion(oxygen)–oxygen–hydrogen angle for an anion(water)–water pair. The F-–water hydrogen bonds are found to have a longer lifetime than all other hydrogen bonds considered in the present study. The lifetime of Cl-–water hydrogen bonds is shorter than that of F-–water hydrogen bonds but longer than the lifetime of water–water hydrogen bonds. The lifetimes of I-–water and water–water hydrogen bonds are found to be very similar. Generally, the lifetimes of both anion–water and water–water hydrogen bonds are found to be significantly shorter than those found under ambient conditions. In addition to hydrogen bond lifetimes, we have also calculated the residence times and the orientational relaxation times of water molecules in ion(water) hydration shells and have discussed the correlations of these dynamical quantities with the observed dynamics of anion(water)-water hydrogen bonds as functions of the ion size and density of the supercritical solutions. [ABSTRACT FROM AUTHOR]

Despite extensive experimental, theoretical, and simulation efforts, a unified description of ionic phase transitions and criticality has not yet emerged. In this work, we investigate the phase behavior of the restricted primitive model of elec0trolyte solutions on the simple cubic lattice using grand canonical Monte Carlo simulations and finite-size scaling techniques. The phase diagram of the system is distinctly different from its continuum-space analog. We find order-disorder transitions for reduced temperatures T*≤=0.51, where the ordered structures resemble those of the NaCl crystal. The order-disorder transition is continuous for 0.15≤=T*≤=0.51 and becomes first order at lower temperatures. The line of first-order transitions is a line of three-phase coexistence between a disordered and two ordered phases. The line of continuous, second-order transitions meets this line of triple points at a tricritical point at T*≃0.1475. We locate the line of continuous transitions, and the line of triple points using finite-size scaling techniques. The tricritical temperature is estimated by extrapolation of the size-dependent tricritical temperatures obtained from a sixth-order Landau expansion of the free energy. Our calculated phase diagram is in qualitative agreement with mean-field theories. [ABSTRACT FROM AUTHOR]

We perform molecular dynamics simulations to study the ionic solvation and association behavior in concentrated aqueous LiCl solutions at ambient conditions, including consideration of expected signatures of ion pairing that might be found in neutron diffraction experiments with isotopic substitution. The ten possible pair radial distribution functions that define the microstructure of the systems are determined and used to assess the first-order difference of the neutron-weighted correlation functions for these solutions in heavy and null water. Then, both sets of correlation functions are applied to the interpretation of the ion’s local environment in terms of the location of the relevant peaks and the penetration of ions into the counterion solvation shells as a signature of ion-pair formation. Finally, we illustrate how first-order difference experiments involving null and heavy water might be used to assess the magnitude of the Mv+-Xv- ion-pair formation for a salt Mv+Xnv- in an aqueous solution, provided the significant experimental challenges in these studies could be overcome. [ABSTRACT FROM AUTHOR]

1H, 23Na, 35Cl, 79Br, and 81Br NMR chemical shifts (δ) and signal half widths (Δ1/2) have been measured in aqueous electrolyte mixtures [tetrahydrofuran/H2O/NaCl and 3-methylpyridine (3MP)/H2O/NaBr] at different mass fractions of salt (X) in the one-phase region, close to their lower critical solution points (TCL). Discontinuous changes in slope of δ=f(X) and Δ1/2=f(X) have been found in 23Na and 81Br NMR spectra of 3MP/water/NaBr solution at X≈0.1 and T=301 K. The dependency of 1H NMR signals of 3MP is continuous over the whole investigated range of X=0.002–0.2, whereas changes in the slope of H2O chemical shifts are hardly noticeable. In the two-phase region, i.e., at T>TCL, a doubling of all NMR signals has been observed. The sensitivity of NMR parameters depends more on composition of solution for anions (Cl- and Br-) than for cations (Na+). A very strong relaxation effect for 81Br nuclei with relaxation rates reaching 14 000 s-1 was observed. The results are interpreted in terms of ion-molecular clustering and changes in coherency of dipole configurations of water molecules during supramolecular restructuring of solutions. [ABSTRACT FROM AUTHOR]

Interactions between salt ions and lipid components of biological membranes are essential for the structure, stability, and functions of the membranes. The specific ionic composition of aqueous buffers inside and outside of the cell is known to differ considerably. To model such a situation we perform atomistic molecular-dynamics (MD) simulations of a single-component phosphatidylcholine lipid bilayer which separates two aqueous reservoirs with and without NaCl salt. To implement the difference in electrolyte composition near two membrane sides, a double bilayer setup (i.e., two bilayers in a simulation box) is employed. It turns out that monovalent salt, being in contact with one leaflet only, induces a pronounced asymmetry in the structural, electrostatic, and dynamical properties of bilayer leaflets after 50 ns of MD simulations. Binding of sodium ions to the carbonyl region of the leaflet which is in contact with salt results in the formation of “Na-lipids” complexes and, correspondingly, reduces mobility of lipids of this leaflet. In turn, attractive interactions of chloride ions (mainly located in the aqueous phase close to the water-lipid interface) with choline lipid groups lead to a substantial (more vertical) reorientation of postphatidylcholine headgroups of the leaflet adjoined to salt. The difference in headgroup orientation on two sides of a bilayer, being coupled with salt-induced reorientation of water dipoles, leads to a notable asymmetry in the charge-density profiles and electrostatic potentials of bilayer constitutes of the two leaflets. Although the overall charge density of the bilayer is found to be almost insensitive to the presence of salt, a slight asymmetry in the charge distribution between the two bilayer leaflets results in a nonzero potential difference of about 85 mV between the two water phases. Thus, a transmembrane potential of the order of the membrane potential in a cell can arise without ionic charge imbalance between two aqueous compartments. [ABSTRACT FROM AUTHOR]