Hydration dynamics model systems

For the model system considered in Eikerling et al., i the chemical composition and water content are fixed. Only minimal hydration could be considered. A more recently begun work aims explicitly at the understanding of structural correlations and dynamics at acid-functionalized interfaces between polymer and water in PEMs. It directly addresses the question of... 

The structural and dynamic properties of water may be affected by both purely geometrical confinement and/or interaction forces at the interface. Therefore, a detailed description of these properties must take into account the nature of the substrate and its affinity to form bonds with water molecules, as well as the hydration level or number of water layers. In order to discriminate between these effects, reliable model systems exhibiting hydrophilic or hydrophobic interactions with water are required. This appears to be the appropriate strategy to permit some understanding of the behavior of water close to a biological macromolecule, as presented in the following. 

When partially hydrated samples are cooled down to 77 K, no crystallization peak is detected by differential thermal analysis. The x-ray and neutrons show that an amorphous form is obtained and its structure is different from those of low-and high-density amorphous ices already known [5]. Samples with lower levels of hydration corresponding to one monolayer coverage of water molecules are under investigation. This phenomenon looks similar in both hydrophilic and hydrophobic model systems. However, in order to characterize more precisely the nature of the amorphous phase, the site-site partial correlation functions need to be experimentally obtained and compared with those deduced from molecular dynamic simulations. 

In order to reveal the origin of the apparent transition from exclusion mechanism to some kind of adsorption of smaller ions at hi concentrations of the mixture, a model system, composed of calcium chloride and hydrochloric acid at varying concentrations, has been examined in detail [170-172]. This mixture is environmentally benign, allows preparation of highly concentrated solutions, and has the advantage that each of the three constituent ions, Ca, H, and Cl, can be easily determined by titration. The test mixture has been studied in both static and dynamic experiments in combination with the hydrated nanoporous polystyrene NN-381. 

An interesting combined use of discrete molecular and continuum techniques was demonstrated by Floris et al.181,182 They used the PCM to develop effective pair potentials and then applied these to molecular dynamics simulations of metal ion hydration. Another approach to such systems is to do an ab initio cluster calculation for the first hydration shell, which would typically involve four to eight water molecules, and then to depict the remainder of the solvent as a continuum. This was done by Sanchez Marcos et al. for a group of five cations 183 the continuum model was that developed by Rivail, Rinaldi et al.14,108-112 (Section III.2.ii). Their results are compared in Table 14 with those of Floris et al.,139 who used a similar procedure but PCM-based. In... 

It is important to propose molecular and theoretical models to describe the forces, energy, structure and dynamics of water near mineral surfaces. Our understanding of experimental results concerning hydration forces, the hydrophobic effect, swelling, reaction kinetics and adsorption mechanisms in aqueous colloidal systems is rapidly advancing as a result of recent Monte Carlo (MC) and molecular dynamics (MO) models for water properties near model surfaces. This paper reviews the basic MC and MD simulation techniques, compares and contrasts the merits and limitations of various models for water-water interactions and surface-water interactions, and proposes an interaction potential model which would be useful in simulating water near hydrophilic surfaces. In addition, results from selected MC and MD simulations of water near hydrophobic surfaces are discussed in relation to experimental results, to theories of the double layer, and to structural forces in interfacial systems. 

The electrolyte concentration is very important when it comes to discussing mechanisms of ion transport. Molar conductivity-concentration data show conductivity behaviour characteristic of ion association, even at very low salt concentrations (0.01 mol dm ). Vibrational spectra show that by increasing the salt concentration, there is a change in the environment of the ions due to coulomb interactions. In fact, many polymer electrolyte systems are studied at concentrations greatly in excess of 1.0 mol dm (corresponding to ether oxygen to cation ratios of less than 20 1) and charge transport in such systems may have more in common with that of molten salt hydrates or coulomb fluids. However, it is unlikely that any of the models discussed here will offer a unique description of ion transport in a dynamic polymer electrolyte host. Models which have been used or developed to describe ion transport in polymer electrolytes are outlined below. 

Because of the complexity of hydrated PEMs, a full atomistic modeling of proton transport is impractical. The generic problem is a disparity of time and space scales. While elementary molecular dynamics events occur on a femtosecond time scale, the time interval between consecutive transfer events is usually 3 orders of magnitude greater. The smallest pore may be a few tenth of nanometer while the largest may be a few tens of nanometers. The molecular dynamics events that protons transfer between the water filled pores may have a timescale of 100-1000 ns. This combination of time and spatial scales are far out of the domain for AIMD but in the domain of MD and KMC as shown in Fig. 2. Because of this difficulty, in the models the complexity of the systems is restricted. In fact in many models the dynamics of excess protons in liquid water is considered as an approximation for proton conduction in a hydrated Nation membrane. The conformations and energetics of proton dissociation in acid/water clusters were also evaluated as approximations for those in a Nation membrane.16,19 20 22 24 25... 

Recently attempts have been made to set up computer models of salt solutions (Kistenmacher et al., 1974 Fromm et al., 1975 Watts et al., 1974). A molecular dynamic study has been reported for a system comprising 198 H20 + 9Li+Cl- (Heinzinger and Vogel, 1974). The results show good agreement with experiment. For example, a water molecule hydrating Li+ rotates faster than water molecules in bulk (Hertz et al., 1971). 

Molecular modeling using either Monte-Carlo simulations or molecular dynamics is used to apply molecular mechanics energy minimizations to very complex systems [348]. In complex flexible molecules such as proteins or nucleic adds, the number of variable parameters, i.e., bond torsion angles, is such that the global search for energy minima becomes impossible The same problem occurs with theoretical calculations of water structure in aqueous solutions or in heavily hydrated crystals. 

In this work, we have approaehed the understanding of proton transport with two tasks. In the first task, deseribed above, we have sought to identify the moleeular-level stmeture of PFSA membranes and their relevant interfaees as a funetion of water content and polymer architecture. In the second task, described in this Section, we explain our efforts to model and quantify proton transport in these membranes and interfaces and their dependence on water content and polymer architecture. As in the task I, the tool employed is molecular dynamics (MD) simulation. A non-reactive algorithm is sufficient to generate the morphology of the membrane and its interfaces. It is also capable of providing some information about transport in the system such as diffusivities of water and the vehicular component of the proton diffusivity. Moreover, analysis of the hydration of hydronium ion provides indirect information about the structural component of proton diffusion, but a direct measure of the total proton diffusivity is beyond the capabilities of a non-reactive MD simulation. Therefore, in the task II, we develop and implement a reactive molecular dynamics algorithm that will lead to direct measurement of the total proton diffusivity. As the work is an active field, we report the work to date. 

---