Next Article in Journal
Multiphysics Modeling and Analysis of a Solar Desalination Process Based on Vacuum Membrane Distillation
Next Article in Special Issue
Gradient NMR Method for Studies of Water Translational Diffusion in Plants
Previous Article in Journal
Development of a Glucose Sensor Based on Glucose Dehydrogenase Using Polydopamine-Functionalized Nanotubes
Previous Article in Special Issue
Ionic Mobility in Ion-Exchange Membranes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Molecular and Ionic Diffusion in Ion Exchange Membranes and Biological Systems (Cells and Proteins) Studied by NMR

by
Vitaliy I. Volkov
1,2,*,
Alexander V. Chernyak
1,2,
Irina A. Avilova
1,
Nikita A. Slesarenko
1,
Daria L. Melnikova
3 and
Vladimir D. Skirda
3
1
Institute of Problems of Chemical Physics RAS, 142432 Chernogolovka, Russia
2
Scientific Center in Chernogolovka RAS, 142432 Chernogolovka, Russia
3
Institute of Physics, KazanFederal University, 420008 Kazan, Russia
*
Author to whom correspondence should be addressed.
Membranes 2021, 11(6), 385; https://doi.org/10.3390/membranes11060385
Submission received: 31 March 2021 / Revised: 30 April 2021 / Accepted: 12 May 2021 / Published: 24 May 2021

Abstract

:
The results of NMR, and especially pulsed field gradient NMR (PFG NMR) investigations, are summarized. Pulsed field gradient NMR technique makes it possible to investigate directly the partial self-diffusion processes in spatial scales from tenth micron to millimeters. Modern NMR spectrometer diffusive units enable to measure self-diffusion coefficients from 10−13 m2/s to 10−8 m2/s in different materials on 1 H, 2 H, 7 Li, 13 C, 19 F, 23 Na, 31 P, 133 Cs nuclei. PFG NMR became the method of choice for reveals of transport mechanism in polymeric electrolytes for lithium batteries and fuel cells. Second wide field of application this technique is the exchange processes and lateral diffusion in biological cells as well as molecular association of proteins. In this case a permeability, cell size, and associate lifetime could be estimated. The authors have presented the review of their research carried out in Karpov Institute of Physical Chemistry, Moscow, Russia; Institute of Problems of Chemical Physics RAS, Chernogolovka, Russia; Kazan Federal University, Kazan, Russia; Korea University, Seoul, South Korea; Yokohama National University, Yokohama, Japan. The results of water molecule and Li+, Na+, Cs+ cation self-diffusion in Nafion membranes and membranes based on sulfonated polystyrene, water (and water soluble) fullerene derivative permeability in RBC, casein molecule association have being discussed.

1. Cation-Exchange Membranes. Structure, Hydration, Ionic, and Molecular Mobility

Ion-exchange membranes are widely applied for modern electrochemical technologies and separation processes. New materials design requires an electro mass transfer investigation. This research is mainly concerned about macroscopic transport processes [1,2,3,4,5,6,7]. However, ion and molecular translation microscopic mobilities have to be investigated for membrane selectivity mechanism understanding.
Of most interest is the relationship between the following fundamentally important characteristics that determine the ion and molecular transport:
  • The nanoscale structure of ion transport channels. The structure and dynamics of polymer matrix at the submicro level from several tenths of nanometer (sizes of solvated ions and molecules) to several nanometers or several tens of nanometers (characteristic lateral dimensions and lengths of ionic channels), determine the selective ion transport because these structural units form transport path for ion transfer by macroscopic distances. Studying the nanostructure opens up the prospects for targeted synthesis of ion exchange polymer, insofar as their preparation is accompanied by the formation of the nanostructure.
  • The type of interaction of mobile ions and hydration water molecules with functional groups. Data on the structure of ionic complexes and on the mechanisms of interaction of ions and water molecules with the polymer matrix are necessary for understanding the mechanisms of selectivity of ion-exchange membranes and elementary steps of the diffusion transport of ions.
  • The elementary steps of diffusion of ions and molecules, which can be characterized by the lifetime of a species on functional group, the time of translational displacement, the partial diffusion coefficient on various spatial scales (if diffusion occurs in a heterogeneous medium).
The problem of elementary diffusion jumps logically follows from the aforesaid. Evidently, the time of elementary jump and the height of the potential barrier overcome by a moving species are largely determined by the geometry of diffusion channels and the structure of hydrate ionic complexes. This information is necessary for both the elaboration of adequate transport models and the targeted synthesis of high-performance ion exchange polymers.
The knowledge of structure and dynamics in a different spatial scale and in a broad band of molecular motion frequencies may be obtained by NMR directly.
NMR spectroscopy The most popular method is 1H NMR, which was used to study Dowex 50 W, CU-2 sulfonate cation exchanger resin and the corresponding membranes MC-40, cation exchange membranes based on polyethylene and sulfonated grafted polystyrene MSC [8,9,10,11,12,13,14,15,16,17] and perfluorinated cation-exchange membranes [18,19,20,21,22,23,24,25,26]. To date, techniques have been developed for recording high-resolution NMR spectra and the main factors that determine the chemical shift of water protons in granulated sulfonate cation-exchangers and ion-exchange membranes have been elucidated. The required information can also be obtained from the solid-state high-resolution NMR spectroscopy data [23,24]. Information on hydration of ionic channels in membranes is of fundamental importance for understanding the mechanism of migration of cations and water molecules.
Alkaline metal cations of lithium, sodium, and cesium were studied by NMR on 7 Li, 23 Na and 133 Cs nuclei in cation-exchange membranes [16,17,20,21,27,28,29,30,31,32] and in sulfonated polystyrene salts [33,34]. Some qualitative data about ionogenic group-cation interaction and cation motion were obtained.
Pulse NMR methods NMR relaxation techniques were for the first time applied for local cationic and water molecules mobility characterization in polymeric electrolytes more than 50 years ago [10]. Spin–lattice and spin–spin relaxation times measurements on 1 H, 7 Li, 19 F NMR nuclei were performed in Nafion and MF-4SC (Russian Nafion type membrane) membranes [35,36,37,38,39,40,41]. Unfortunately, the numerical calculation of correlation times is hard work because of wide molecular motion frequency distribution.
The study of the metal ion mobility is associated with even more serious complications. The most serious obstacle is the absence of theoretical works to serve as the base for studying the region of diffusion motion with the characteristic correlation times longer than ωo−1 (where ω0 is the NMR frequency, usually, ~109 Hz). To date, 7 Li and 23 Na NMR relaxation works have been performed dealing with the mobilities of lithium and sodium cations in the CU-2 type sulfonate cation exchangers and the corresponding membranes, perfluorinated sulfonate cation-exchange membranes [31,32,33,34,39,40,41].
The pulsed field gradient NMR method, [42,43] which makes it possible to directly measure the diffusion coefficients of protons and other ions in heterogeneous media, is free from these drawbacks. The number of studies of this type substantially increased in recent decades, [5,6,19,38,39,40,41,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83] which was associated with the increased interest in the problems of ionic mobility in polyelectrolytes.
To summarize the foregoing, the following points should be underlined. The magnetic resonance techniques and, especially, NMR methods provide the unique possibility of acquiring detailed information on the state of molecules and ions, the local molecular and ionic mobility, and the diffusion on the spatial scale from several tenths of nanometer to several millimeters. The advantages of NMR spectroscopy also involve the possibility of studying one and the same sample under conditions resembling the service conditions by several methods simultaneously, which makes it possible to compare the results of different measurements and unambiguously interpret them. For these studies to be performed the experimental procedures should be worked out and the problems associated with theoretical quantitative description of data should be solved.
The NMR methods are especially attractive for acquiring detailed information on the ion and molecular transport in polymer electrolytes. The modern level of experimental research instruments allows one to study both elementary processes and macroscopic transfer under the service conditions of electrochemical systems. The successful introduction of NMR methods into the research and technological practice is limited by the lack of publications devoted, first of all, to demonstration of the potential of experimental NMR techniques in this research field.
Despite the considerable number of NMR studies of polymer electrolytes, the reviews on this subject are scarce. In the present review, the experimental results obtained by NMR methods on the ion and water molecular transport in polymer ion exchangers carried out in Russia and abroad are analyzed and generalized. From our point of view, such an analysis will demonstrate the potential of modern NMR methods and help to reveal some fundamental features of ion and molecular transport in ion exchange membranes at the molecular level.
The main results of NMR studies in ion-exchange membranes are discussed. Attention is focused on the potential of NMR techniques in solving particular problems in relation to the most widely known ion exchangers. The most thorough studies were carried out for perfluorinated membranes. Using these membranes as examples, an attempt is made to find the relationship between the polymer matrix structure, the ion hydration details, and the diffusion mobility of ions and molecules on different spatial scales and then to apply this information for revealing details of the ion transport mechanism in ion-exchange membranes.

1.1. Ion-Exchange Membranes, Nanochannel Structure

In this review we consider mainly two types of well-known cation exchange membranes.
The first type is perfluorinated sulfonic cation exchange membrane Nafion and its Russian analog MF-4SC, carboxyl cation exchange membrane F-4CF. The second one is membranes on the basis of sulfonated polystyrenes.
The main studies with the use of NMR methods were carried out on perfluorinated sulfonate cation-exchange membranes Nafion and MF-4SC (Russian analogue of Nafion). At present, the nanostructure of perfluorinated sulfonate cation-exchange membranes has been studied in sufficient detail. Several models were proposed for describing the transport channels [5,84,85,86,87,88,89,90,91,92,93,94,95] of which the Gierke model [5,6,89,90,91] was used most widely (Figure 1). This model is also applied for describing of inorganic channels in polystyrene sulfonic cation exchange resins and membranes.
This model is based on the concept that water molecules and counterions form spherical clusters with walls built of sulfo groups. It was proposed that the clusters are connected by channels. The limiting step of diffusion is ion and molecular transport in the channels, which were not observed experimentally. On the basis of more detailed small angle X-ray scattering investigation, a group of authors [5,84,85,95] proposed the channel model of perfluorinated membranes (Figure 1 and Figure 2).
The detailed Nafion structure model is based on the data of small-angle X-Ray scattering [84,85] Mössbauer spectroscopy [93] and standard porosimetry [94] and shown in Figure 3. It was found that associates including sulfo groups, counterions, and water molecules are formed in the amorphous part of the membrane. These associates form the ion transport channels with the volume of approximately a fourth of the total membrane volume [96,97].
The key point of this structure model is the comparison of macroscopic lithium cations and water molecules self-diffusion coefficients with calculated microscopic self-diffusion coefficients (Table 1). From NMR relaxation data the correlation times τi of water molecules and Li+ cations jumping were calculated. Self-diffusion coefficients of water molecules and lithium cations were estimated from Einstein equation as l2/6τi where l is the jump length of water molecule (0.3 nm) or Li+ cation (0.7 nm, the average distance between SO3 groups).
The good agreement between experimental macroscopic and calculated self-diffusion coefficients is observed. It is confirmed the conclusion from channel structure model that macroscopic water and ion transfer is controlled by microscopic particle jumping.
Therefore, this interpretation of 1 H and 7 Li relaxation data is conformed to Nafion structure model (Figure 3).
Another extensive class of membranes is sulfonate cation-exchange membranes based on styrene-divinylbenzene. Both heterogeneous and homogeneous membranes were studied as well as gel and macroporous sulfonate cation exchangers CU-2 and CU-23 and their foreign analogues. The ESR studies of Cu2+ cations used as paramagnetic probes revealed two domains containing hydrated copper (II) complexes with different mobility [100]. The dependences of the correlation times and diffusion coefficients of water on the moisture content [44,45] were consistent with the assumption that the structure of transport channels in such ion exchangers is described by the model in Figure 1, which is widely used for ionic transport in these ion-exchangers. Cluster-channel structure model is also suitable for explanation of electro mass transfer in membranes (MSC) based on polyethylene and sulfonated grafted polystyrene [16,17].
It should be concluded that ion transport is related in many respects to the mobility of water molecules and is determined by the character of membrane hydration. To reveal the mechanisms of ionic conduction, it is important to correlate the character of hydration of ionogenic groups with the translational mobility of water molecules and ions and the ionic conductivity of membranes. Such studies should be carried out for membranes with a well-known structure of transport channels.

1.2. Specific Features of Cation Hydration. Mechanism of Cation—Functional Group Interaction

For a cation hydration characterization in cation-exchangers, 1H NMR spectra of water molecules are analyzed. The first calculation of one charge cation hydration numbers were carried out on sulfonate cation exchange resin Dowex 50 W in the beginning of 1970-th just after spectrometers of NMR started to be on sale.

1.2.1. Cation Hydration

The Dependence of 1H Chemical Shift on the Humidity

Acid ionic form
The most thoroughly studied membranes are perfluorinated sulfonate cation-exchange membranes Nafion, MF-4SC [18,19,20,21,22,23,25], MF-4SC membranes modified with inorganic dopants [21] and perfluorinated carboxylic membranes F-4CF in various ionic forms [20,21]. Membranes 11 00385 i001
The typical NMR spectra of protons in a perfluorinated sulfonate cation exchange membrane in its Li+, Na+, Cs+ forms and of hydrated H+ cations measured at various relative humidity of the environment are rather narrow singlet lines, indicating high mobility of water molecules and H+ counter-ions in membranes. Some examples of the evolution of 1H NMR spectra at different humidity and temperature are shown in Figure 4 and Figure 5.
The H+ ionic form 1H NMR spectrum is a singlet, which position is shifted to lower magnetic fields relatively to bulk water signal. The NMR line is rather narrow even at low humidity and below 0 °C which indicates high proton mobility at these conditions. The chemical shift value depends on water content and temperature (Figure 4 and Figure 5).
Counter ion H+ hydration numbers calculation in sulfonate cation exchange resins and membranes [8,9,13,14,15,16,25], as well as in aqueous solutions of sulfonic acids [101,102,103,104] on the basis of 1H chemical shift temperature and moisture content dependences analysis, was carried out. The absolute value of water molecules per sulfonate group (λ) have to be obtained in order to calculate hydration number correctly. As a rule membrane moisture content (λg) is determined by gravimetric technique. In our case during this procedure the sample of membrane dries on phosphoric anhydride or at temperature 110 °C until the constant mass. After this drying 1H NMR line remains narrow (spectrum 7 in Figure 4) which is evidence of high proton mobility. Therefore it was assumed that in perfluorinated sulfonate cation membranes two residual water molecules per H+ cation are persisted and form hydroxonium ion [H5O2]+ [13,14,15,18,19,25,105]. The absolute water content was determined by 1 H NMR spectroscopy technique directly and residual water content is 1.50 ± 0.50 in Nafion 115 membrane after drying to constant weight at 100 °C [106]. Following this paper, we have measured the absolute water content in Nafion 117 membrane dried to the constant weight at 110 °C or equilibrated with P2O5 at room temperature. 1 H NMR spectra of the membrane at 10% and 75% humidity are shown in Figure 6. Cyclohexane solution in carbon tetrachloride inserted in the NMR sample tube was used as standard (line 1 in Figure 6). Line 2 in the Figure 6 is the averaged signal of water molecule protons and hydrated H+ cation.
The amount of residual water molecules, calculated as the average value from 20 measurements for different drying procedures (drying at 110 °C and with P2O5) was λ0 = 1.9 ± 0.4.
The observed 1H NMR line is the average signal from H+, hydrated H+(H2O)h water molecules, and water molecules of the next hydration spheres. Due to the fast (compared to the chemical shift difference between these three proton positions) molecular exchange the only NMR line with an average chemical shift is observed as it is shown in Figure 4 and Figure 5. The resulting chemical shift δ is a superposition of chemical shifts of H+ and hydrated water molecules δc in H+(H2O)h and water protons δH2O—Equation (1), where pc and pH2O are the respective relative protons fractions.
δ = p c δ c + p H 2 O δ H 2 O
It is assumed that H+ ion forms stable aqua complex with humidity independent hydration number h. The chemical shift (δ-δH2O) is approximated by Equation (2).
δ δ H 2 O = 2 h + 1 δ c δ H 2 O 2 λ g + h + 1
The dependences of chemical shift on the humidity are shown in Figure 7. This dependence is described well by Equation (2) when h is 2 (curve 1, Figure 7). This h value is equal to the residual water amount λ0 = 1.9 ± 0.4 which was obtained above by the direct NMR measurement. Therefore, it may be concluded that after membrane drying until the constant weight two hydrated water molecules form hydroxonium cation H5O2+. This conclusion is agreed with sulfonic cation exchanger resins and MF-4SC membrane H+ ionic form 1H NMR investigation [8,9,13,14,15,16].
Some authors [107,108] have proposed that in dry Nafion H+ ionic forms the hydronium ion H3O+ (h = 1) is formed. According to the direct NMR measurements carried by us and the authors of [106] the amount of residual water λ0 is about two and the dependence of chemical shift on humidity is approximated by Equation 1.2 at h = 2 (curve 1 in Figure 7) that corresponds to hydroxonium ion. Assumption of one residual water molecule (hydronium ion, h = 1) contradicts to the direct water content NMR measurements and Equation (2) (curve 2, Figure 7). Therefore, the assumption of hydroxonium ion formation at low water content looks more probable.
Two residual water molecules and (H5O2)+ hydrated ion formation at low water content are typical for all sulfonic cation exchangers: sulfonic cation exchange resins, MC, and MSC membranes.
Salt ionic forms
As the moisture content decreases, the NMR spectral lines shift downfield for membranes in the acidic form and upfield for membranes in the salt form relative to those of free water. Such changes in the chemical shifts of water protons cause by a destruction of a hydrogen-bond network in the salt-form membranes and strengthening of hydrogen bonds as a result of hydration of H+ counter-ions. Dependences of 1H chemical shift on water content in different ionic forms of sulfonated and carboxylic membranes are shown in Figure 8 and Figure 9.
The moisture content measured by gravimetry is expressed as the number of water molecules (λ) per ionogenic group. The 1 H chemical shifts sharply change starting from a certain critical value λ because at low water content, all water molecules occupy the first hydration shell of the cation and are strongly affected by the cation. The hydration energy of H+, Li+, and Na+ cations is higher and the hydration energy of Cs+ is lower than the hydrogen bond energy; hence, the changes in the chemical shifts of water protons are not so pronounced in Cs films (curve 4, Figure 8 and curve 8, Figure 9). In the acidic form of sulfonate cation exchange membranes two remaining water molecules forms the stable cation H5O2+ with the counter-ion H+. Therefore the observed dependence of the chemical shift on the real fraction of water protons in the membrane 1/(2λ + 5) is a straight line (line 9 in Figure 9) [19].
The calculated hydration numbers (the average number of water molecules in the hydration shell of an ion (h)) for H+, Li+, Na+ and Cs+ cations are as follows [19,20], Table 2:
The hydration numbers for Li+, Na+ and Cs+ were also assessed from the dependences of 7 Li, 23 Na, and 133 Cs chemical shifts on the moisture content. These values agree well with the data of proton magnetic resonance.
The hydration numbers of cations in carboxylic membranes are smaller than in sulfonate cation exchange membranes. This is probably caused by the fact that in F-4SC membranes, the hydration shell of cations includes, first of all, carboxyl groups, whereas in MF-4SC membranes even at λ ~ h the local environment of a cation contains only water molecules.
It should be noted that the hydration numbers of doubly charged ions substantially exceed those of singly charged ions [19]. In the membranes with a different polymer matrix and, correspondingly, a different structure of transport channels, for instance, styrene–divinylbenzene systems, the hydration is of the same nature.
This fact and the weak dependence of the chemical shift in electrolyte solutions on the type of anion suggest that in cation-exchange membranes, water molecules interact mainly with cations. The fraction of broken hydrogen bonds increases with increase in the hydration energy of cations [29]. For λ > h, the cation and the functional group form a solvent-separated ion pair and the hydration shells of ionogenic groups overlap (Figure 10a). For a low moisture content λ < h, a network of hydrogen bonds between water molecules is broken, thus hampering the ion jumping between the neighboring ionogenic groups (Figure 10b).

Temperature Dependence of 1H Chemical Shift

Acid ionic form
Counter ion H+ hydration number may be also calculated from 1 H chemical shift temperature variation. The 1 H NMR spectra in Nafion 117 membrane at different temperature is shown in Figure 5. In acid ionic form of membrane lines are shifted to lower magnetic fields, but line width increases with the lowering temperature. The narrow line width at temperature below 0 °C indicates high proton mobility at low temperature. The temperature dependences of chemical shift at different water contents are shown in Figure 11. These dependences are straight lines until λ = 7–8 in the whole temperature range from −60 °C to + 50 °C (curves 4-9 in Figure 11) and at λ = 12 in temperature range above −25 °C (curve 3 in Figure 11) and at λ = 17.5 above −10 °C (curve 2 in Figure 11). The slope of these lines increases with the increasing of λ; the slope of the lines is highest for bulk water (curve 1 in Figure 11).
Follow to 1 H NMR temperature investigation of acidic aqueous solutions [101,102,103] and sulfonic cation exchange resins [8,9,14,15,16] in H+ ionic form the hydration numbers were calculated from Equation (3):
h = λ 0.5 + λ d δ dT d δ H 2 O dT ,
where λ = λg + λ0. Hydration numbers h of H+ cation in H+ ionic form Nafion membranes at different λ are listed in Table 3.
Hydration numbers are differing with humidity variation. At low water content (λ is about 2–4) h is closed to 2, which agrees with previous result, obtained from proton chemical shift humidity dependences. With following λ increasing h also increases until h ≈ 4.
Hereby in spite of rather crude approaches of proton chemical shift humidity and temperature dependences interpretation some important conclusions may be done.
At low water content H+ cation strongly binds two water molecules which are not able to be desorbed at high temperature and vacuum drying forming hydroxonium H5O2+ ion. With humidity increasing the more water rich hydrated complex H9O4+ is formed. The same hydration peculiarities are observed for all sulfonic cation exchangers.
Salt ionic forms
From 1 H chemical shift temperature dependences hydration numbers of Li+, Na+, and Cs+ ions may be calculated. Relevant techniques were applied to sulfonic cation exchangers [8,9,14,15,16] and salt aqueous solutions [101,102,103].
Temperature chemical shift dependences at different water contents are shown in Figure 12. These dependences are straight lines, which slope increased following by λ increasing; the line slope is highest for bulk water (curve 10 in Figure 12).
Following 1 H NMR temperature investigation of sulfonic cation exchange resins and membranes [8,9,14,15,16] as well as salt aqueous solutions [101,102,103] in salt ionic form the hydration numbers h were calculated from Equation (4):
h = λ 1 d δ dT d δ H 2 O dT ,
where λ is the number of water molecules per sulfonate group; δ is the measured 1H chemical shift; δH2O is the bulk water 1 H chemical shift. Hydration numbers h of Li+ cation in appropriate ionic form of Nafion membranes at different humidity are listed in Table 4.
Hydration numbers are differed with humidity variation. At high water content λ > 10.7, hydration number h about 4–6, this value is closed to lithium cation hydration number in dilute lithium salt aqueous solutions (h = 5–6), where Li+ forms a separate ionic pair with anion. With reducing water content Li+ hydration number is decreased and oxygen atom of sulfonate group replaces the oxygen of water molecule; thereby a contact ionic pair is created. Sodium cation hydration number h is 6 ± 1 (λ = 10 at 98% RH) which is the same as h in dilute sodium salt aqueous solutions [103,104] so the separate ionic pair Na+—SO3 group is formed at these conditions. For Cs+ cation hydration number h is 1 ± 0.2 (λ = 4 at 98% RH) this value is less compare with h in aqueous solution [103] and in Dowex 50W (cation resin) [8] (h = 3–4). It means that even at maximum water content Cs+ interact with Nafion SO3 group directly and contact ionic pair is forming. This cesium cation hydration particularity is due to the fact that compare to lithium and sodium cations Cs+ hydration energy is less than hydration bond energy between water molecules [110].
Hereby, in spite of rather crude approaches of proton chemical shift temperature dependences interpretation, some important conclusions may be done.
At low water content Li+ and Na+ cations are in direct contact with sulfonate groups, while at high water content water molecules are built in between cation and SO3 groups. Cesium cation interacts with sulfonate group directly at any water content.
Hydration numbers h of Li+, Na+ Cs+ in MSC membrane are shown in Table 5, hydration numbers for same cations in equimolar aqueous salt chloride solutions and membrane water uptakes are also shown for comparison. The h values of Li+, Na+, and Cs+ in MSC membranes at high humidity is practically equal to salt solution ones. The crystallography radii and Stokes–Einstein hydrodynamic ion radii are also presented in Table 5. We have calculated hydrodynamic ion radii from Stokes–Einstein equation on the basis of ionic diffusion coefficient in chloride aqueous solution. Ion diffusion coefficient concentration dependences were approximated to infinite dilute concentration.
Li, 23 Na, and 133 Cs NMR Spectroscopy. Mechanisms of Cation—Ionogenic Group Interaction7 Li, 23 Na, and 133 Cs NMR spectra of in MSC membranes are shown in Figure 13, for example. These spectra are represented by narrow lines that belong to Li+, Na+, and Cs+ cations in appropriate ionic forms.
Figure 14 shows the dependences of chemical shifts in 7 Li NMR spectra on the moisture content in the Li+ form of MF-4SC, Nafion 117, and carboxylic F-4CF membranes.
For water content λ > 4, the chemical shifts depend slightly on the membrane type and the moisture content; for λ < 4, the 7Li chemical shift for sulfonate cation exchange membranes to decrease. The half-widths of 7Li NMR lines at λ > 3–4 are almost the same for all membranes, while their line width increase is observed at λ < 2–3 which is due to the longer lifetimes of Li+ ions located on ionogenic groups.
The dependences of chemical shifts on the moisture content in the Na+ form of MF-4SC and F-4CF membranes (Figure 15a) are similar to the corresponding dependences for their Li+ form. For λ > 4–5, the 23 Na chemical shifts are close for sulfonate cation-exchange and carboxylic membranes and start to dramatically decrease at λ < 4–5. As in the case of Li+ ions, this can be attributed to the formation of ionogenic group-cation contact pairs, which sharply reduces the sodium ion mobility.
Cs+ cation 113 Cs chemical shift humidity dependence is not so sharp in opposite to 7 Li and 23 Na ones for Li+ and Na+ cations (Figure 15b). These signals in carboxylic membranes shift in high field compare to sulfonate membranes (Figure 15).
The small influence of water content on 133 Cs NMR line shift and width is due to low ion hydration energy Cs+ ion in comparison with Li+ and Na+ hydration energies. Cesium ion hydration energy is less than water molecule hydrogen bond energy. Therefore, water molecules make own hydrogen bond network and Cs+ cations forming a contact ionic pair with charge groups even at high humidity. For this reason, cesium ions’ immediate surrounding does not depend on humidity which is followed by slight humidity dependence of 133 Cs chemical shift and line width.
Line width and chemical shift of Li+ and Na+ nuclei humidity behavior may be explained as follows. At high water content (λ >> h) a hydrate cation shell similar to a hydrate cation shell in salt aqueous solutions and cation-ionogenic group ionic pair is separated. At low water content when λ < h sulfonate or carboxylic group oxygen atoms may enter to cation sphere forming contact ionic pairs. In this case cation surroundings symmetry is disturbed which accompanied with the chemical shift and line width increasing owing to slowing-down cation mobility [20,21].
The chemical shift of alkali metal nuclei can be calculated by the formula:
δ i = p c δ i c + p b δ i b ,
where pc and pb are the relative fractions of contact and solvent-separated ion pairs, respectively; δic, δib are the chemical shifts of ions forming either the contact or solvent-separated ion pairs [21] Using this equation, the dependences of the relative proportion of contact ion pairs on the moisture content were calculated (Figure 16).
The fraction of contact ion pairs was found to vary in the following series for both sulfonate cation-exchange and carboxylic membranes:
Cs+ > Na+ > Li+.

1.2.2. Water Behavior at Temperature Lower 0 °C

Hydration particularities control the structural and dynamic water behavior. On the basis of DSC measurements water in membrane nanochannels usually is divided on “unbounded” and “bounded” water. Unbounded water forms ice phase at temperature below 0 °C, while bounded water is mobile at freezing temperature [25,111,112]. On the 1H NMR data water hydrogen bond network in Nafion membrane is destroyed. Therefore, water molecules are not able to organize ice phase and their mobility is high at temperatures below 0 °C. This proposal explains the narrow proton NMR line width at freezing temperatures. Let us discuss this phenomenon in more details. The observed 1H NMR lines are caused by protons of H+ and mobile water molecules. The mobile water molecules amount proportional to NMR line area. The dependences of mobile water molecule amount on temperature at different water content are shown in Figure 17. The number of mobile water molecules is temperature independent until λ about 8 (curves 1–5 in Figure 17), the similar result was obtained by A. Guillermo and co-workers [74]. With λ increasing (λ > 9), the mobility of some water molecules at temperature below 0 °C is decreased, anyway about 8–9 water molecule per sulfonate group possess high mobility independently on temperature (curve 6 in Figure 17).
This effect may be explained, keeping in mind the existence of macroscopic dimension pores with low sulfonated group concentration in perfluorinated membranes [113]. The first water molecules are sorbed in nanochannels where ionogenic group concentration is high, after occupying these sites (about 9 molecules per sulfonated group) the other water molecules fill macropores. In macropores the water structure is similar to water structure in bulk, so this water is able to form ice phase at low temperature. Because of low ice molecule mobility their NMR line is very wide (about 22 kHz) it is impossible to observe ice molecules by high resolution NMR technique. It is also very important to mention that in spite of the fact that water does not freeze at temperatures below 0 °C the DSC thermogram peak is observed (as it is shown for example in Figure 18 [38]) which is usually interpreted as the freezing of free water [111,112,114,115].
The same water behavior is also observed for other ionic forms and type of membranes.
As it is shown in Figure 19, water in MSC membrane doesn’t frees until λ < 18 water molecules per sulfonate group [16].
The same phenomenon is observed for salt membrane forms.
In Figure 20 the independent of amount of mobile water molecules and Li+ cations on temperature is illustrated.
This disagreement between DSC and NMR results may be explained by the 1 H spin–relaxation data. Let us investigate water mobility at freezing condition in more details using NMR pulsed field gradient and NMR relaxation techniques.

1.3. NMR Relaxation. Local Mobility of Molecules and Ions

H NMR Relaxation, Local Proton Motion

Despite considerable progress in the research of ion and molecular transport in polymer electrolytes, the micro-scale transfer mechanisms are still unclear. As was stated above, the NMR methods furnish a lot of information on the local motion of ions, solvent molecules, and polymer matrices. The times of spin–lattice (T1) and spin–spin (T2) relaxation and the width of NMR line are very sensitive to local mobility. Local proton motion in Nafion membrane may be characterized by 1H spin relaxation data. The dependences of spin–lattice (T1) and spin–spin relaxation (T2) times on temperature and humidity are usually analyzed [35,36,37,38,41,105]. The detailed quantitative analysis based on Bloembergen, Purcell, Pound (BPP) magnetic dipole-dipole relaxation mechanism for distribution of correlation time in H+ and Li+ ionic form sulfo cationic perfluorinated membrane was carried out [31,37,38,105].
The dependences of spin–lattice relaxation rate (T1−1) on temperature in H+ ionic form of Nafion 117 for different water content are shown in Figure 21 and Figure 22. These dependences show wide maximum in low temperature region and maximum or shoulder in high temperature region. Spin–lattice relaxation rate low temperature maximum is shifted to high temperature region with decrease of membrane humidity.
The authors of [26] interpreted the T1−1(T) shape by two phases proton relaxation as a result of two different types of water motion in Nafion membrane. This interpretation is not unequivocal. For more defined interpretation the mutual analysis of proton spin–lattice (T1) and spin–spin (T2) relaxation processes is necessary. This investigation has been carried out in Li+ ionic form of perfluorinated sulfo cation exchange membrane MF-4SC [31,37,38,105]. The temperature dependences of T1−1(T) and T2−1(T) in samples with different water content were measured.
To interpret the measured relaxation times, it is necessary to determine the types of interactions and molecular motions.
For nuclei with spin 1/2, the principal relaxation interactions are the magnetic homo- and heteronuclear dipole–dipole couplings. The NMR relaxation for nuclei with spin > 1/2, e.g., 7 Li, is mainly determined by two mechanisms:
  • quadrupole relaxation due to interactions between the quadrupole moment of a nucleus with the fluctuating electric field gradient induced by the charge distribution around this nucleus;
  • dipole–dipole relaxation caused by random fluctuations of magnetic moments of nuclei.
It was shown that the magnetic dipole–dipole coupling is the main interaction inducing relaxation of water protons. The magnetic relaxation of 7 Li nuclei is modulated by the 7 Li−1 H magnetic dipole–dipole coupling and the quadrupole coupling between lithium nuclei and the electric field gradient. The 1H relaxation was described in terms of the Bloembergen–Purcell–Pound model, which takes into account the distribution of correlation times in the low-temperature region. The distribution of correlation times are either Gaussian or rectangular [31,37,38]. The processes of spin–lattice relaxation of 7 Li nuclei were determined by rapid librations of water molecules within the cation hydration shell, which allowed their lifetime to be calculated [38]. In contrast to the kinetics of longitudinal magnetization decay, the kinetics of transverse magnetization decay was biexponential.
Detailed analysis of spin–spin relaxation processes was carried out. The spin–spin relaxation was assumed to be modulated by slower motions as compared with librations of water molecules within the hydration shell of a lithium cation. These motions are associated with variations of the electric field gradient on a 7 Li nucleus as the latter approaches (or moves away from) a sulfo group. Based on these data, the lifetime of lithium cations on sulfo groups was calculated. The algorithms of calculation of correlation times of water molecules, the lifetimes of lithium hydration shells, and lithium cations located on sulfonic groups for MF-4SC membranes can be found in several publications [31,37,38]. Below, we present the main conclusions made in these publications.
Figure 23 and Figure 24 show temperature dependences of the longitudinal and transverse relaxation rates of 1 H and 7 Li nuclei.
Water molecule 1 H relaxation rates show two maximums of R1 and minimum of R2 (Figure 23). The minimum R2 is an unusual shape because for two phase relaxation process (the phase populations are temperature independent) R2(T) increasing with decreasing temperature smoothly without extremum [105]. The only way to explain a spin–lattice and spin–spin relaxation together, it was proposed, that the patterns of motion of water molecules are different in the low- (T < 280 K) and high-temperature (T > 320 K) ranges. In this case the two phases relaxation is also taken place, but the phase populations pl and ph are temperature dependent.
1 T 1 T = p l T T 1 l T + p h T T 1 h T
This was connected with the fact that at low temperatures, water molecules form associates that are highly mobile at T < 0 °C because in membranes, water molecules interact with charged ionic groups and cannot form the solid phase (ice). In the high-temperature region, water molecules can move as in electrolyte solutions. In the intermediate temperature range, both phases coexist. Figure 25 shows the temperature dependence of the high-temperature ph (solid line) and low temperature pl (dashed line) phase populations. The transition from one phase to another is accompanied by absorption or released energy due to the formation (rupture) of hydrogen bonds. Note that the observed thermal DSC effect was not connected with a freezing (melting) of free water, because according to NMR data, the number of mobile water molecules did not change upon the phase transition [25,38]. Associate formation, but not water crystallization shows DSC peak.
As it will be shown, the assumption of water molecule association below 0 °C, have explained the shape of self-diffusion coefficients temperature dependences.

1.4. Pulsed Field Gradient NMR. Diffusion and Ionic Conductivity

Pulsed field gradient NMR is very attractive for ion and molecular diffusion investigation in ion-exchangers because it is a direct technique of partial self-diffusion coefficient and diffusant relative content characterization [14,15,16,17,20,21,24,25,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,105,109,115]. Stimulated echo sequences are generally used for measurements [42,43]. We have reviewed the water molecule and alkali metal ion diffusion data, which are compared with ionic conductivity results [16,17,21,25,55,56,57,58,59,60,61,105]. Some details of pulsed field gradient NMR (PFG NMR) technique and self-diffusion coefficients calculations are given in Appendix A.

1.4.1. Pulsed Field Gradient 1 H NMR, Water Molecule Self-Diffusion and Ionic Conductivity

As it is shown in Figure 26, water molecule translational mobility controls an ionic conductivity.
Figure 27 shows the dependence of diffusion coefficients of water molecules on the moisture content in perfluorinated MF-4SC and F-4CF membranes of various ionic forms. The diffusion coefficients of water molecules (and hydrated H+ counter-ions for a membrane in its acidic form) vary in the following sequence of ionic forms for sulfonate cation exchange membranes at the same water content:
H+ > Cs+ > Na+ > Li+
The identical series of water diffusion coefficients are typical for carboxylic membranes in the salt form. However, in contrast to MF-4SC membranes, in the acidic form of carboxylic membranes, water and hydrated H+ cations behave much differently. In carboxylic membranes, the diffusion coefficients are two orders of magnitude smaller than in sulfonate cation-exchange membranes (see curves 5 and 1 in Figure 27a). In carboxylic membranes, restricted diffusion occurs, which is manifested as a decrease in the diffusion coefficients of water and hydrated H+ cations with increase in the diffusion time (Figure 27b). The size of restriction regions is about one µm and decreases with a decrease in the moisture content. Presumably, these restrictions are induced by the formation of internal hydrogen bonds between carboxylic groups of neighboring macromolecules.
The dependences of ionic conductivity of membranes on the moisture content have a similar form as the analogous dependences for diffusion coefficients (Figure 28) [20,21]. In the region λ < h, the conductivity of samples sharply drops down as the moisture content decreases. In sulfonate cation-exchange membranes, the conductivity varies in the following series at maximum water content H+ >> Na+ ≈ Li+ > Cs+, whereas in carboxylic membranes, the inverse dependence is observed. Sulfonate cation-exchange membranes in their Li+ and Na+ forms exhibit higher conductivity as compared with carboxylic membranes. For their Cs+ form, the inverse dependence was observed. This probably can be explained by the larger radius of cesium ions. The interaction of Cs+ with the membrane matrix and water molecules is fairly weak, even for low degrees of hydration. Moreover, for low degrees of hydration, hopping of cesium cations between the oxygen atoms of ionogenic groups is rather limited by the distance between the oxygen atoms, which is similar for membranes of both types, and not by the negative charge on oxygen ions. For carboxyl membranes, this distance may be smaller, due to a larger negative charge on oxygen atoms.
The conductivity of carboxylic membranes in their acidic form is four orders of magnitude lower than the conductivity of sulfonate membranes in the same form.
The ionic conductivity was calculated based on the Nernst–Einstein equation:
σ = ne 2 D kT
where n is the number of charge carriers, cm3; D is the self-diffusion coefficient, m2/s; e is the electron charge,1.9·10–19 C; k is the Boltzmann constant, 1.38·10–23 J/K; and T is the temperature.
The obtained results were compared with the results of direct measurements (Table 6).
For MF-4SC membranes in the Li+ form, the experimental value of conductivity and the value calculated from diffusion coefficients of lithium cations coincide to within the experimental error (6.2·10−3 and 6.5·10−3 S·cm−1, respectively). The ionic conductivity found from water diffusion coefficients is only several times higher than the conductivity measured experimentally. The similarity of these values suggests that the translational displacements of the cations and water molecules are correlated.
Temperature dependences of self-diffusion coefficients at different water content are shown in Figure 29a. The main features of the curves are different slopes at high and at low temperature regions [16,19,25,38,74,105,109,115]. High and low temperature parts were approximated by the Arrhenius equitation:
D = D 0 e E a R T
where Do is temperature independent, R is gas constant, T is absolute temperature, Ea is self-diffusion activation energy.
The corresponding activation energies for different humidities are shown in Figure 29b. At high humidity the activation energies do not depend on water content and are close to bulk water activation energy. At low water content, high and low temperature activation energies increase and become closer to each other. The changes of the slope in temperature dependences is not observed for λ less than four water molecules per sulfonate group (curve 6, Figure 29a, star in Figure 29b).
Proton conductivities temperature dependences [25,78,79,80,81,111,112] are very similar to the temperature dependences of self-diffusion coefficients. Temperature dependences of self-diffusion coefficients, protonic conductivities calculated from Nernst–Einstein Equation (7), and experimentally measured conductivities are shown in Figure 30 for λ ≈ 4–5 as an example.
The calculated conductivities are 3–4 times bigger compared to the experimental values. There are some reasons for this difference. On the one hand, the measured self-diffusion coefficient is characterized as average translation mobility of water molecules and hydrated H+ and it is impossible to distinguish H+ translational mobility, which controls a charge transfer. On the other hand, not every cation jump is accompanied by charge transfer. The changes in the conductivity and self-diffusion curve slopes are observed at the same temperatures. The increase in the slopes of the self-diffusion and proton conductivity curves at low temperatures is usually explained by the freezing of free (unbound) water at temperatures below 0 °C [81,111,112].
In our opinion this explanation contradicts to the temperature dependence of the mobile water amount. As it was found above water does not freeze until λ = 9–10 (Figure 17), nevertheless the change in the slope on temperature dependences is observed at λ less than this water content (Figure 30a). It is also very important to mention that in spite of the fact that water doesn’t freeze at temperatures below 0 °C the DSC thermogram peak is observed (as it is shown for example in Figure 18 [38]), which is usually interpreted as the freezing of free water [81,111,112,114].
This disagreement may be explained by the 1 H spin–relaxation data. As it was shown above at low temperature water molecules form associates, therefore the self-diffusion activation energy is higher compare to high temperature region, where the water structure is similar to bulk water one. The temperature dependence of the diffusion coefficients D in the MSC membrane held at different relative humidity is also linearized in the coordinates of the Arrhenius Equation (7) (Figure 31).
The main feature of the temperature dependences of the self-diffusion coefficients is an increase in the self-diffusion activation energy with a decrease in moisture content in the entire temperature range. This phenomenon is similar to Nafion self-diffusion coefficient temperature dependences. It should be also mentioned, that water in MSC membrane is mobile at freezing temperatures until λ = 18 [16].
Experimental diffusion coefficients and those calculated from Nernst–Eistein Equation (7) at different moisture contents are shown in Figure 32.
It is noteworthy that the dependences of the diffusion coefficients of H+ cations calculated from the conductivity and measured by NMR on the moisture content are similar. The values of diffusion coefficients decrease with decreasing relative humidity. The values of the diffusion coefficients determined by NMR are higher than those calculated based on the data on conductivity. This is quite typical for ion-exchange membranes and arises because the conductivity is limited by proton transfer in narrow channels (Figure 1) of the membranes, while NMR characterizes the averaged mobility of protons and water in larger pores and channels [16,17].

1.4.2. Pulsed Field Gradient 7 Li, 23 Na, 133 Cs NMR, Alkaline Metal Cation Self-Diffusion

Self-diffusion of water molecules and Li+, Na+, Cs+ cations were investigated in appropriate salt ionic form of Nafion 117 membranes.
Lithium cations and water molecules
Lithium cations and water molecules self-diffusion in Li+ Nafion ionic form were measured at different water contents. Spin–echo signal attenuations (diffusion decays) of 7 Li nuclei of Li+ cation and water molecule 1 H nuclei are shown in Figure 33.
Spin–echo attenuations were exponential approximated by Equation (A2) in three orders of magnitudes in the whole range of water contents. It indicates that translation mobility of water molecules and lithium cations was characterized as the only self-diffusion coefficient. Self-diffusion coefficients of lithium cation and water molecule dependences on water content λ are given in Figure 34.
These dependences shapes are similar. It may be concluded that translational motions of Li+ cation and hydrated water molecules are correlated as it was indicated in [105].
Temperature dependences of Li+ and water molecules self-diffusion coefficients at different moisture content are shown in Figure 35 and Figure 36.
The dependences Ds(T) were approximated by Arrhenius Equation (8).
Equation (8) describes well Li+ self-diffusion in temperature region from −10 °C to +30 °C (Figure 35) and water molecule self-diffusion from −40 °C to +30 °C (Figure 36). Activation energies at different water content λ summarize in Table 7.
Activation energies Ea increase with hydration degree decreasing. Lithium cation self-diffusion activation energies are essentially more compare to water molecule self-diffusion activation energies.
Sodium and cesium cation self-diffusion
The measurement of Na+ and Cs+ self-diffusion coefficients is more effortful compare to Li+ self-diffusion coefficient. Until now the sodium and cesium self-diffusion data in Nafion membrane are unknown. There are two reasons of these difficulties. First is low magnetic moment of 23 Na and 133 Cs nuclei. The second is the large quadrupole moment of these nuclei. The main spin relaxation mechanism of 23Na and 133Cs nuclei is quadrupole moment–gradient ligand electric field interaction. We have observed these nuclei diffusion decays in aqueous solution and in sulfocation-exchange membranes (MSC) based on polyethylene sulfonated grafted polystyrene [16,17], which maximum water content λ is more compared to Nafion membrane. At this condition cation water molecule shell is high-symmetric. Therefore, electric field gradient is low and nuclear spin–spin relaxation time is long enough to observe spin–echo signal. We also succeeded to get spin–echo attenuation in Nafion membrane at RH = 98%. The examples of 23 Na and 133 Cs decays of Na+ and Cs+ cations are given in Figure 37 and Figure 38.
Lithium, sodium, and cesium cation self-diffusion coefficient temperature dependences in Nafion membrane at maximum humidity (RH = 95%) in the temperature range from 20 °C to 80 °C are shown in Figure 39.
The dependences Ds(T) were approximated by Arrhenius Equation (8). Self-diffusion coefficients increase in the next sequence Li+ ≈ Na+ > Cs+. Let us consider cation self-diffusion behaviors in MSC membrane. Spin–echo attenuation (diffusion decay) of 7 Li+, 23 Na+, and 133 Cs+ is exponential in salt ionic form of MSC membrane. Diffusion decay is well approximated by Equation (A2)–Figure 40.
Diffusion coefficient temperature dependences are also linearized in the coordinates of the Arrhenius equation (Figure 41). Cation diffusion coefficients increase in a sequence Li+ ≈ Na+ < Cs+. This row is the same for cation diffusion coefficients of chloride aqueous solutions [17]. Cation diffusion activation energies are about 16–18 kJ/mol.
Self-diffusion coefficients at 20 °C, self-diffusion activation energies of Li+, Na+, and Cs+ cations in Nafion membrane are indicated in Table 8. The values of Li+, Na+, Cs+ self-diffusion coefficients and activation energies in MSC sulfonic cation-exchange membrane and dilute chloride aqueous solution are also given for comparison [17].
Self-diffusion activation energies of Li+ and Na+ cations in Nafion 117 membrane are about 20 kJ/mol which is close to water self-diffusion activation energy in these membranes [16,37,105], but for Cs+ activation energy is distinctly more compared with Li+ and Na+ cations. It should be mentioned that in MSC membrane cation self-diffusion activation energies are little less and Cs+ activation energy are even smaller than Li+ and Na+ compare with Nafion membrane. Cation self-diffusion coefficients are changed in the next rows Li+ ≤ Na+ > Cs+ in Nafion; Li+ < Na+ < Cs+ in MSC [17] and chloride aqueous solutions. Higher Cs+ activation energy and lower self-diffusion coefficient in Nafion compared to aqueous solution and MSC membrane may be explained in the following way. As it was mentioned above as opposed to Li+ and Na+, Cs+ cation interacts with Nafion SO3 group directly and a contact ionic pair is forming. This is the reason of low mobility and high self-diffusion activation energy of cesium cation in Nafion.

1.4.3. Lithium, Sodium, Cesium Cation Self-Diffusion and Ionic Conductivity

Ionic conductivities σc of Li+, Na+, and Cs+ were calculated in Nafion membranes on the basis of Nernst–Einstein Equation (7). Calculated conductivities σc was compared with experimental values σe measured by impedance spectroscopy. The data are listed in Table 9.
Ionic conductivities calculated from cation self-diffusion coefficients are appreciably more compared with experimental meanings. The following discrepancy reason may be the next: ionic conductivity is controlled only by cation transfer along applied electric field, but PFG NMR fixes all translational jumping of cation, for instance, cation coming and leaving of sulfonate group.
Temperature dependences of experimental and calculated from Nernst–Einstein equation ionic conductivities of different ionic forms of MSC membranes are shown in Figure 42.
The conductivity values at different humidity and activation energies calculated from Arrhenius equation are listed in Table 10.
Ionic conductivity of investigated membranes increases in the sequence Li+ < Na+ < Cs+ << H+. It should be noted that the diffusion coefficients of lithium, sodium, and cesium cations in MSC and in aqueous solutions change in the same sequence.
As it is shown in Figure 42, the calculated and experimental conductivity curve are similar. Conductivity activation energies of alkaline metal cations are close to each other. However, the calculated conductivity values are in one or two orders of magnitude more in comparison with experimental ones. This difference seems natural. Ionic transport in these membranes is realized through the system of channels and pores, which size is depended on polymeric matrix nature and hydration degree (Figure 1). Ionic conductivity is limited by the transport of ions in the narrow channels. They are usually called the “bottle neck” [116,117]. The diffusion coefficient measured by NMR in the first turn may be due to high mobility ions localized in wide pores [16].

1.4.4. Li+, Na+, and Cs+ Hydration and Diffusion in Chloride Aqueous Solutions

A dependence of 1H water molecule chemical shift on solution concentration is shown in Figure 43. Proton signal shifts to the high field with chloride concentration increasing. This fact is explained by destroying of hydrogen bonds between water molecules [58,103,104,105]. This phenomenon is stronger for CsCl solution since Cs+ possesses low polarizing properties (due to the large ion size), which causes a hydrogen bond system destroying [17].
As a result, water molecule translational mobility in CsCl aqueous solution should increase with increasing concentration, which is observed experimentally (curve 3 in Figure 44).
Water and alkaline metal cation self-diffusion coefficients increase in the row of cation atomic mass growing: Li+ < Na+ < Cs+. This increasing sequence credits with cation hydrated radius decreasing due to ion hydration energy reduction from Li+ to Cs+. In opposition to membrane, measured ionic self-diffusion coefficients in dilute salt aqueous solution are equal to diffusion coefficients calculated from ion conductivity data [118]. This indicates that conductivity and PFG NMR are measuring different kinds of ionic mobility in membranes. Self-diffusion coefficient obtained by PFG NMR is an average self-diffusion coefficient while ionic conductivity is restricted by ion motion in narrow channels.
It should be mentioned that Li+, Na+, and water diffusion coefficients are reduced greater compared to Cs+ ion for which water diffusion coefficient even increases with an increase in electrolyte concentration (curves 1′, 2′;1, 2 Figure 44 and curves 3′, 3 Figure 44). The water hydrogen bond network is destroying to a greater extent with increasing in Cs+ concentration. As contrasted to Cs+, hydration energy of Li+ and Na+ ions is more. Therefore, mobility of water molecules connected with these cations drops and water diffusion coefficient is reduced with salt concentration rising.
Now the question should be considered about the reason of the self-diffusion coefficient decreasing with solution concentration growing. Water molecules hydrate both cations and anions. Cation hydrate shell structure is rebuilt during ion motion. The neighboring ions obstruct to hydrated shell forming and consequently ionic mobility decreases. This obstruction factor rises with concentration growth or ion activity coefficient decreasing [17]. Therefore, the water and ion self-diffusion coefficients are reduced with solution concentration increasing.
In Figure 45 7 Li, 23 Na, 133 Cs nuclei chemical shift concentration dependences are shown. This dependence is stronger for 133Cs. Chemical shift increases sometimes (compared to 7 Li, 23 Na), while CsCl concentration is varied from 1 mol/L to 4 mol/L (curve 3, Figure 45).
The chemical shift of these nuclei is determined by nuclear quadrupole moment interaction with gradient of electric field induced by hydration molecules of water. Because high hydration energy lithium and sodium cation hydrated shells are steady and symmetric. The electric field symmetry and consequently a cation nuclear chemical shift only insignificantly varies with concentration changes (curves 1, 2 Figure 45). In contrast to lithium and sodium, cesium cation hydration shell is not steady therefore a symmetry of hydrated electric field on Cs+ decreases with cesium chloride solution concentration rising which accompanied by appreciable increasing of 133Cs chemical shift (curve 3, Figure 45).

1.4.5. Specific Features of Translational Mobility of Water Molecules and Ions at a Low Moisture Content

For water molecule, cation self-diffusion, and conductivity dependence on humidity at low water content in MF-4SC membranes explanation, the percolation theory was applied. The equation for diffusion (conductivity) coefficient dependence on the amount of water molecule per sulfonated group (λ) derived by professor S.F. Timashev [5]:
D = D 0 e x p p a λ 1 3 μ
The physical sense of parameter a, p, and µ explained in [5]. This equation describes well self-diffusion coefficients of water molecule and lithium cation humidity dependences at low water content (Figure 46).
For effective membrane application, high water and appropriately cation high translation mobility is very important. As mentioned above, the high proton (ionic) conductivity is explained by a continuous network of hydrogen bonds. Hydrogen bond destruction is attended by threshold water and cation mobility reduction. The high ion mobility at low humidity may be achieved by reduction of charge group spacing because in this case for hydrogen bond network forming the smaller value of water molecules per sulfonate group λ is necessary. As it is shown in Figure 47, water self-diffusion coefficients at low λ rise in an order of magnitude, while SO3 group concentration (ion-exchange capacity) increases in two times.
The next way to increase ionic mobility at low water content is to insert oxygen atoms in polymer matrix. This oxygen “bridge” forms additional hydrogen bonds for water molecule motions. For this goal achievement, Nafion membranes were modified by inorganic dopants. The membrane modification and membrane conductivity properties specify in [119].
The next way is to synthesize ion exchange polymers containing oxygen atoms in its structure. This idea was realized in aromatic sulfo-containing polyamides, which fragments are shown in Figure 48.
Figure 48 shows the diffusion coefficients of lithium cations in sulfonate membranes with different structures of transport channels and in isomeric disulfo-containing aromatic polyamides as a function of the moisture content [20,39,62]. In systems with amide groups, which serve as bridges, reinforcing the network of hydrogen bonds between the sulfo groups, the diffusion coefficients of Li+ ions are two orders of magnitude higher than in systems containing only sulfo groups: MF-4SC and macro pore sulfocation exchange resin on the basis of sulfonated polystyrene CU-23. The structure of ionogenic transport channels substantially affects the translational mobility of ions and molecules. Thus, although the exchange capacity of CU-23 is five times higher compared with MF-4SC, the diffusion coefficients of lithium cations in MF-4SC membranes with their uniform distribution of sulfo groups and regular structure of channels were found to be comparable with Li+ diffusion coefficients in CU-23 in which the transport routes are formed by ionogenic clusters connected by narrow channels (see curves 3, 4 in Figure 49). The translational mobility of Li+ in terephthalic salt (PA-2) having a regular rod-like structure of channels was shown to be higher than in isophthalic salt (PA-1) with a globular packing of polymer chains (see curves 1, 2 in Figure 49).

1.4.6. Diffusion of Saturated Monatomic Alcohols and Water–Alcohol Mixtures

Alcohol molecules’ transport through ion-exchange membranes is important for a water-alcohol separation processes development. Investigation of alcohol membrane transfer mechanism is also topical for methanol and ethanol fuel cell operation understanding. A lot of publications were devoted to diffusion of saturated monoatomic alcohols and their mixtures in Nafion [64,67,68,72] and MF-4SC membranes [47,51,73]. Diffusion of water-ethanol solutions through membranes based on mono- and disulfo-containing aromatic polyamides [51] and through composite polyacrylic acid-polysulfonate membranes [48] is studied.
Alcohols’ molecules sorbing take place mainly in ionogenic channels of sulfonate MF-4SC membranes. Water, methanol, ethanol, and propanol self-diffusion coefficient dependences on diffusant content are shown in Figure 50. At maximum water and alcohol contents, self-diffusion coefficient of water molecules is more compared to alcohol molecules, but at low membrane diffusant occupation, water self-diffusion coefficient became less than alcohol ones.
The diffusion coefficients and diffusion activation energies sharply depend on the ionic form of the membrane and the concentration dependences of diffusion coefficients of alcohols follow a threshold pattern (Figure 50 and Figure 51). Hence, it can be concluded that diffusion of alcohols in MF-4SC membranes is similar to water diffusion.
Translational transfer of alcohols and water molecules in MF-4SC membranes and monosulfo-containing polyamide films proceeds via the same ionogenic channels. For this reason, in these systems, the separation factors of water–alcohol mixtures are low [51]. Table 11 shows the partial diffusion coefficients of water and ethanol and the separation factors of their mixture in MF-4SC membranes, films of monosulfo-containing polyamide (PA-1) and disulfo-containing polyamide (PA-2), composite membranes based on polyacrylic acid and polysulfone (PAA-PSF). Apparently, for PA-2 and PAA-PSF, the separation factors are higher by several orders of magnitude, as compared with MF-4SC and PA-1. This is connected with the fact that in contrast to MF-4SC and PA-1, the water and ethanol molecules in PA-2 and PAA-PSF are located in different transport channels [48,51].

2. Conclusions

Nowadays, chemical power sources, based on sulfocation exchange membranes find wide application. The principal part of these systems is formed by a polymer ion exchange membrane that should have high ionic conductivity. The problem of revealing the mechanism of ion transport in these membranes becomes quite challenging. The last decade was characterized by active studies of the state and mobility of cations and molecules in ion-exchange membranes of various types by using the methods of NMR and pulsed field gradient NMR on various nuclei.
Therefore, the review of NMR technique applications is a problem of today. In this review paper the results of cation-exchange membrane investigation obtained by hetero nuclear high resolution NMR spectroscopy, NMR relaxation, pulsed field gradient NMR are discussed.
The main attention is given to interconnection of membrane diffusion channel nanostructure, cation hydration and water molecule and cation mobility in different spatial scales. NMR self-diffusion data are compared with ionic conductivity measurements.
The main parts of NMR investigations were carried out in Nafion (or Russian Nafion analog) MF-4SC sulfonate perfluorinated membranes. These membranes are most studied by different physical techniques and could be a model system for a wide set of ion-exchangers.
The comparison of local water molecule and Li+ cation mobility calculated from 1 H and 7 Li spin relaxation data with water and lithium cation self-diffusion coefficients measured by PFG NMR shows that macroscopic transfer is controlled by ion and molecular jumping near sulfonate groups. This result is conformed to Nafion channel structure model in Figure 3. Therefore, a cation hydration governs by ionic motion.
Hydration numbers of alkaline and alkaline–earth metal cations were calculated from water molecule 1H chemical shift temperature dependences. For Li+, Na+, Cs+ counter ions, the relative part of contact pairs cation-charge group dependently on humidity was measured by 7 Li, 23 Na,133 Cs NMR. Some conclusions about membrane selectivity mechanism to these ions were proposed.
It was definitely shown that in sulfonate cation–exchangers in acid ionic form the least hydration number is equal two and at low water content hydrated cation [H5O2]+ is formed.
In opposite to conception based on DSC data about water freezing in membranes below 0 °C it was shown that amount of mobile water molecules does not change at temperature variation in spite of DSC peak observing. On the basic of 1 H spin–relaxation data it was supposed that at freezing temperature water molecules form mobile associates, but not ice phase. This assumption explains water and cation self-diffusion and ion conductivity temperature dependences.
A comparison of ion conductivity calculated from cation self-diffusion coefficients with experimental values confirms the cluster-channel structural model for membranes based on sulfonated polystyrene (Figure 1) and channel structural model for Nafion membranes (Figure 3).
From our opinion these NMR results give opportunity to understand mechanism of ionic and molecular transport in ion-exchange membranes more deeply.

3. Biological Systems: Protein and Cell Membranes

Another subject of our paper is biological systems: protein and red blood cells. In Figure 52 the dependence of water self-diffusion coefficients on gelatin and bovine serum albumin concentration in aqueous solution is shown (black). This dependence is similar to the dependence of water self-diffusion coefficient on water content in ion exchangers (blue in Figure 52). It may be proposed that water behavior in ion exchange membrane ionogenic channels and in protein ionic channels has some analogy. The forming polymer structure from macromolecules also may be similar to synthetic and biological polymers.

3.1. Protein Association

At the present time, the method of nuclear magnetic resonance with a pulsed field gradient magnetic (PFG NMR) is applied to study the features of translational mobility of a wide variety of complex molecular systems. Particular attention is paid to study the structure, dynamics, and functional (protein–protein) interactions of biomolecules. Note that determination of the geometric properties of molecules from the measured values of the self-diffusion coefficients using the Stokes–Einstein relation is only a part of the capabilities of the NMR method.
The interaction of proteins with other macromolecules or small partner molecules plays an important role in most biological processes [120]. Often, such interactions are manifested in a conformational exchange that occurs on a much faster time scale than most experimental biophysical methods [121]. NMR spectroscopy techniques in combination with pulsed magnetic field gradient NMR have a unique ability to extract information about these interactions and are used more often to investigate protein systems of increasing complexity, including proteins with an internally disordered structure (IDP). A characteristic feature of such proteins is that sample an ensemble of rapidly interconverting alternative conformations ranging from random coils to more structured conformations with secondary structure and residual tertiary structure elements. [122]. At the same time, IDP form well-defined complexes when interacting with a partner [123].
Another feature of PGF NMR method is that the translational mobility of the molecule under interest is registered directly existing at its native conditions without a need for any labeling. In addition, there is also no need for continuous calibrating and comparing with reference samples, when one analyzes PFG NMR data.
Besides the structural and geometric parameters of the molecules, PFG NMR investigations can afford fundamental insights into molecular mechanisms underlying the protein functionality, including conductivity of ion channels of membrane proteins [124,125], dynamics, and binding of intact ligands [126,127]. The effect of the bicells association with molecules such as proteins, peptides, detergents, or drugs can not only be detected with PFG NMR, but also quantitatively characterized based on its results, since the method provides information on the spectrum of self-diffusion coefficients for all components of the molecular system.
The topology features of the macromolecule, which determine its ability to change its conformation to a greater or lesser extent in a particular solvent [128,129], significantly affect the nature of the translational mobility of macromolecules. This fact was confirmed by the results of studying the self-diffusion of globular proteins in aqueous suspensions [130].
When immersed in water, some protein chains assume a rigid globular conformation [131]. Self-diffusion of globular proteins in aqueous solutions in [130] was studied by PFG NMR in a fairly wide range of concentrations of macromolecules at different pH-values of the medium and at different temperatures. In particular, aqueous suspensions of myoglobin, bovine serum albumin, barstar, and lysozyme were studied. Using the method of constructing generalized concentration dependences of self-diffusion coefficient of macromolecules [132], authors of the paper [130] managed to obtain a generalized concentration dependence for proteins, which, in turn, gave certain grounds to speculate about the existence of some general regularities of self-diffusion of globular proteins in aqueous suspensions (Figure 53).
It was shown [130] that the generalized concentration dependence of the self-diffusion coefficients of proteins, like the universal concentration dependence of the self-diffusion coefficients of linear flexible-chain polymers, approaches the characteristic asymptotes with ϕ−0 in the limit of dilute solutions and with ϕ−3 in the region of concentrated solutions.
However, in the intermediate concentration range, the obtained generalized dependence of self-diffusion coefficients for proteins differs significantly from the analogous curve for polymers. In [130], this difference is explained by the peculiarities of the dynamic behavior of proteins, which are due to the specific globular conformation of polypeptide chains. At the same time, a comparison of the diffusion behavior of globular proteins and rigid Brownian particles (see, for example, [133,134]) showed that some differences are observed only in the region of concentrated solutions, at φ φ ^ [130].
Thus, in the case of aqueous solutions of globular proteins, the study of self-diffusion of macromolecules by PFG NMR, the construction of a generalized concentration dependence of the self-diffusion coefficients of proteins, and a comparison of this dependence with a similar curve for flexible-chain polymers made it possible to qualitatively describe the features of the translational mobility of polypeptide chains and to characterize the structure of globules in the studied solutions.
At this point in time, in addition to structured globular proteins, studies of proteins with an internally disordered structure, which are characterized by high flexibility and the presence of disordered regions of the polypeptide chain, are relevant, and at the same time retain their functionally active state.
Many studies have shown [135,136] that IDP is capable of forming ordered aggregates—amyloid fibrils. Such formations can be the cause of many serious diseases, such as neurodegenerative diseases of Alzheimer’s and Parkinson’s [137,138].
It should be noted that it is necessary to be able to adequately describe this “family” of proteins and to correctly correlate the physical properties of IDP with their function in order to better understand the structural and dynamic features of such proteins as a whole. For the study of such systems, the NMR method has a wide range of possibilities and is, in fact, the most informative physicochemical method for studying protein molecules.
Melnikova et al. [139] investigated the dependence of the translational diffusion of α-casein, a typical representative of intrinsically disordered proteins, on the protein concentration, diffusion time, and storage time of the sample using the PFG NMR method. In this work, the authors compared the concentration dependence of the self-diffusion coefficient of a representative of IDP with generalized curves as a first step towards understanding the features of IDP diffusion in a wide range of experimental conditions. The dependence of α-casein self-diffusion coefficient < D > on protein concentration, expressed as a volume fraction φ, is shown in Figure 54 (purple diamonds).
As can be seen from the Figure 54, the concentration dependence of the self-diffusion coefficients of α-casein from the side of low concentrations shows, at first, the same tendency as globular proteins. However, in concentrated solutions it crosses the curve for linear flexible-chain macromolecules and reaches the asymptotic < D s > φ 12 .
The authors of the work [139], found that it is the self-organization effect, which manifests itself in the α-casein/water system above a certain critical concentration, that is the reason for the anomalously strong concentration dependence of the average values of the self-diffusion coefficients of α-casein molecules in the region of high concentrations. Moreover, the obtained concentration dependence does not coincide with any of the known generalized concentration dependences of the average self-diffusion coefficients previously established for linear flexible-chain polymers, as well as for dendrimers and globular proteins.
To establish the effect of self-organization, the authors of the investigation analyzed the dependences of the forms of diffusion decays of the spin echo signal on the diffusion time. In regions of relatively high concentrations, anomalous behavior of the self-diffusion coefficient was found. Figure 55 shows diffusion attenuations for the 15% α-casein solution, acquired at td values of 200, 400, and 800 ms.
The graphs are presented in conventional coordinates (Figure 55a) and in coordinates corresponding to testing the fully limited diffusion mode [140,141] (Figure 55b). As can be seen from Figure 55a, with an increase in the diffusion time for a part of protein molecules characterized by small values of self-diffusion coefficient, the slope of diffusion attenuations, which corresponds to a decrease in the values of the self-diffusion coefficients for this part of molecules with an increase in the diffusion time td. Figure 55b shows the same diffusion attenuations, but in coordinates in which fully limited diffusion mode can be easily tested. In these coordinates, it is clearly seen that the slope of the diffusion attenuation component with the lowest self-diffusion coefficient remains unchanged for all td values, as shown by the dashed lines in Figure 55b. This proves that the self-diffusion coefficient Dmin is inversely proportional to the diffusion time, and the root-mean-square displacement of α-casein molecules remains constant < r2 > ~ td0 as follows from the Equation (A4).
The presence of such a dependence allows us to assert that in concentrated solutions of α-casein, processes are observed that are very similar to the processes of gelation. Similar effects were observed in [43], especially in gelatin gels [140], as well as in radiation cross-linked polybutadienes [142].
The size of constraints calculated by (A4) is equal to r 2 50   ±   5   nm . The obtained value significantly exceeds the size of the α-casein molecule itself, the hydrodynamic radius of which is in the range of 2.2–3.2 nm [143]. Thus, the effect of limited diffusion of α-casein can only be associated with the fact that several α-casein molecules interact with each other to form a supramolecular structure similar to a three-dimensional gel network. In this case, interactions through the hydrophobic regions of the protein, hydrogen, ionic and other non-covalent bonds can act as intermolecular bonds [144,145,146].
In addition, the authors of the article [139] demonstrated that in the studied concentrated solution there are two sub-assemblies of α-casein molecules with different characteristics of translational mobility: by the values of the self-diffusion coefficients and by their dependence on the diffusion time. Moreover, it was possible to show that a state of dynamic equilibrium is realized between the indicated sub-assemblies of α-casein molecules. In other words, molecular exchange occurs between α-casein molecules in the gel state and free molecules with a characteristic lifetime of α-casein molecules in the gel phase 3.5 ± 0.4 s. In this case, the stationary value of the fraction of α-casein molecules in the gel state is estimated at 0.93 ± 0.01.
It should be noted that as we know in the studies of systems with gelation by the PFG NMR (see, for example, [43]), no signs of exchange between free polymer molecules and in a gel state were found.
It should be recognized that the amount of experimental material concerning the transport properties of IDP in solutions and mixtures remains insufficient. As shown in a review of experimental results of studying the self-diffusion of polymers and globular proteins in solutions, NMR diffusion (PFG NMR) is one of the most effective methods for studying the structural and dynamic properties of molecules in condensed systems.

3.2. Water and Biological Active Substances Self-Diffusion in Red Blood Cells

The transport of substances into biological cells is carried out in two main ways: active and passive. Active transport is the transfer of a substance through a biomembrane, flowing against a concentration gradient (from a region of low concentration to a region of high concentration), i.e., with the expenditure of free energy.
Due to passive transport, substances are transferred along a concentration gradient from a high concentration region to a low concentration region without energy consumption (diffusion, osmosis). There is a wide variety of diffusion transfer methods: diffusion of fat-soluble substances through the lipid part of the membrane; transport of hydrophilic substances through the pores formed by membrane lipids and proteins; facilitated diffusion with the participation of special carrier molecules. In the lipid phase of the membrane, non-polar substances are readily soluble: organic and fatty acids, esters, since they have an increased affinity for the lipid bilayer.
At the same time, polar substances (inorganic salts, sugars) poorly pass through the lipid bilayer. However, for water, the value of permeability is about 10−6 m/s, which is quite large for a polar substance insoluble in lipids [147]. An explanation for this phenomenon was found with the discovery of aquaporins—integral membrane proteins that form pores in cell membranes, through which water molecules are selectively passed, allowing it to enter and leave the cell, while at the same time ions and other soluble substances are not allowed through (Figure 56) [148,149]. That is, there are two ways for water to enter the cell: through the lipid bilayer and through aquaporins. The schematic movement of water through the channel formed by aquaporin is shown in Figure 56 [150]. The aquaporin channel is shaped like an hourglass. Bulk water is located in the extracellular vestibule and the intracellular vestibule of the canal. Further, the vestibule turns into a narrow tunnel with a width of up to 2.8 Å, which is enough for the passage of one water molecule only, while big molecules or hydrated ions are unable to pass through it. It has now been ascertained that some aquaporins also allow other solutes to pass through (glycerin, urea, ammonia, etc.) [151].
Water molecules are very important in the functioning of biological cells. Water is a solvent for many substances; acts as a medium for chemical reactions and participates in them (hydrolysis); maintains the shape of the cell and participates in its thermoregulation. Therefore, it is necessary to study the characteristics of water exchange in biological cells. Red blood cells (RBCs, erythrocyte) are of interest for such studies, since they are the main component of blood and play a key role in the formation of rheological parameters of blood. The erythrocyte membrane also has general principles of the organization of biological membranes.

3.2.1. RBC Water Permeability Studied by Paramagnetic Doping Technique

A number of papers are devoted to the permeability study of humans, birds and animals RBCs membranes for water molecules using various techniques [152,153,154,155,156]. Studies of water exchange in RBCs have not lost their relevance for more than 50 years. So, at the beginning, the look of the researchers was aimed at establishing the pathways for transporting water into the cell, and the presence of water channels was only assumed [152,153]. The discovery of aquaporins has brought some clarity to this issue. The search for methods for assessing the permeability of cell walls is still important [154,155,156]. This is due to the peculiarities of obtaining RBCs, the short time of their “life” outside the organism, which is due to their rapid response to changes in the environment. The most methods provide only indirect information about metabolic processes (including the water permeability of membranes).
In the field of studying the permeability of RBCs membranes for water molecules, NMR techniques are widely used [156,157,158,159,160,161]. In practice, two methods are used to assess the permeability of biological membranes: paramagnetic doping NMR technique and pulsed field gradient NMR technique.
In the first case, paramagnetic salt solutions are introduced into the extracellular region of the cell suspension, which leads to a decrease in the relaxation times of the protons of water molecules. The exchange rate k = 1/τ is calculated as the difference between the rates of spin–spin relaxation of protons in intra- and extracellular water, where τ is intracellular water lifetime. Permeability is calculated by the formula: P = k V S , where V/S is the ratio of cell volume to its surface area [157].
This method is widely used in the works of Benga and Kuchel with colleagues [156]. The authors have evaluated the permeability of RBCs membranes for water in humans and more than 30 animal species. In Figure 57a, diagram of the values of the permeability P of RBCs for water in humans and some animal species at temperatures of 25 and 37 °C is shown. The authors put forward an interesting hypothesis about the relationship between the permeability of the RBCs membrane and the physical activity of the species. As shown by the results of NMR experiments, the permeability of RBCs to water is higher in species with greater physical activity, with a higher metabolic rate, or a higher average blood circulation rate.
In order to establish the mechanism of permeability, the permeability of erythrocytes was measured when protein channels in membranes are blocked. The inhibition of the permeability of the RBCs membrane for water molecules caused by insertion of p-Chloromercuribenzene sulfonate (PCMBS), which blocks the transfer of water through water protein channels, was calculated. A diagram showing the percentage of inhibition of permeability for different types of RBC is shown in Figure 58. The main result of the work is the correlation of the RBCs water permeability with the molecular basis of the process. It has been shown that the permeability of erythrocytes for water, inhibited by PCMBS, corresponds to the water transported through the membrane by water protein channels. A low percentage of inhibition indicates the absence of water protein channels in the RBCs membrane. This is observed for chicken red blood cells, for example. For human RBCs ~50% of the membrane’s permeability to water is due to water protein channels. Thus, the percentage inhibition of permeability is related to the number of blocked water protein channels and differs from species to species.
Despite the elegance of the paramagnetic doping NMR method for estimating the water permeability of RBCs, the presence of paramagnetic ions can affect the actual metabolic process, that was discussed in the works [161,162,163]. In addition, the obtained values of the permeability of RBCs are calculated from the indirect parameters of molecular diffusion.
From this point of view, pulsed field gradient NMR (PFG NMR) is preferable. This technique allows a direct quantitative estimate of molecular diffusion without a destructive effect on biological systems [160,164,165,166,167]. The originality of the method lies in the ability to measure the partial self-diffusion coefficients and the relative amounts (populations) of diffusant molecules in heterogeneous systems, which include biological cells [82,83,168,169,170].
Below we have discussed the self-diffusion of water molecules, lateral diffusion and self-diffusion of pentasubstituted fullerene derivatives C60 (with adducts of 3-mercaptopropanesulfonic acid C60[S(CH2)3SO3Na]5H (1), mercaptopropionic acid C60[S(CH2)2COOK]5H (2), proline C60[N(CH2)3CHCOOK]5Cl (3)) in a RBCs suspensions investigated by PFG NMR [171,172,173].
To measure the self-diffusion coefficients of diffusant molecules in a suspension of RBCs, the “stimulated echo” pulsed sequence was applied. The details of PFG NMR technique are given in supplementary.

3.2.2. Self-Diffusion of Water Molecules in a RBCs Suspension

Figure 59 shows the proton spectra of a suspension of RBCs at different values of pulsed field gradient amplitudes.
Figure 60a shows the diffusion decays of protons of water molecules in RBCs suspension at different diffusion times td [26]. This diffusion decays are decomposed into three exponential components approximated well (in four orders of magnitudes) by Equation (A1), as it is illustrated in Figure 60b. In Figure 60b capture the partial self-diffusion coefficients Ds1, Ds2, Ds3 and populations p1, p2, p3 are indicated.
Partial self-diffusion coefficients and populations dependences on the diffusion time td are given in Figure 61. With td increasing self-diffusion coefficients Ds1 and Ds2 are decreased, while a Ds3 is diffusion time independent. The Ds3 value closes to the self-diffusion coefficient of bulk water. The population p1 increases but p2, p3 decrease with increasing of td.
These three types of water in the RBCs suspension are intracellular, extracellular, and bulk water. Similar conclusions were obtained for water in other biological systems (chlorella and yeast) [82,83,170]. In the following discussion, the main attention is paid to intracellular water, which behavior is characterized by Ds1 and p1.
The Ds1(td) decreasing with increasing of diffusion time indicates the limited nature of water self-diffusion in erythrocytes, which corresponds to region of restricted diffusion in Figure A2 [174,175].
Figure 62a shows the dependences of the self-diffusion coefficient of intracellular water molecules Ds1 on the diffusion time at temperatures of 5, 20, and 35 °C. As it can be seen from the figure, the dependence Ds1(td) is flatter compared to td−1. The effective self-diffusion coefficients calculated from Equation A6 are shown in Figure 62b. Parameters D0 and Dp were: at 10 °C D0 = 2.7∙10−9 m2/s, Dp = 2.5∙10−12 m2/s; at 25 °C D0 = 2.7∙10−9 m2/s, Dp = 4.5∙10−12 m2/s; at 35 °C D0 = 2.7∙10−9 m2/s, Dp = 6.0∙10−12 m2/s.
If the values of D0, Dp, and the pore size a are known the value of the permeability P of the pore wall can be calculated using the Equation (A5) [176], where P is the permeability of the pore wall, a is the pore size. This equation was used to assess the permeability of the cell walls of the chlorella [83], yeast [82], and roots of corn [175].
For erythrocytes, the diffusion restriction size a, calculated from Equation (A4), does not depend on temperature, and its value was 2.1 μm. The value of the permeability P of the cell membrane of mouse erythrocytes, calculated in accordance with Equation (A5), varies from 0.3·10−5 m/s to 0.5·10−5 m/s with an increasing temperature from 5 to 35 °C
The hindered water self-diffusion coefficient Dp, which characterized permeability was calculated from Equation (A5). The temperature dependence of Dp is Arrhenius type (Figure 63).
As it follows from Equation (A5), under the condition D0 >> Dp, this dependence actually reflects the dependence of the cell wall permeability P on the temperature t. The activation energy of self-diffusion was ED = 24.1 ± 1.9 kJ/mol.
Due to the cell wall permeability, a molecular exchange between intracellular and extracellular water takes place. To estimate the residence time of water molecule inside the cell and the corresponding exchange constant k, a two-component model of molecular exchange is used, which implies an exponential distribution function of the residence time. As it follows from this model, the relative part of water molecules in the cell (population p1) dependence on diffusion time is (10) [43,176,177]:
p 1 t d   =   p 1 0   e x p   t d τ
where τ is the water molecules residence time in the cell.
Another reason of population p1 reducing is a spin–lattice relaxation. Therefore, the p1(td) is a biexpontial shape [178]:
p 1   =   p f   e x p t d τ + p s   e x p t d T 1
where τ is the lifetime of the molecule; T1 is a spin–lattice relaxation time, which is ≈600 ms; ps, pf are the weights of a slow and fast components, respectively. Figure 64 shows the dependence p1(td) and its decomposition into components in accordance with Equation (11). The value of the lifetime τ of the molecule inside the cell was 20 ± 2 ms at 30 °C.

3.2.3. Lateral Diffusion in the RBCc Membrane

To study the self-diffusion of lipids in the RBCs membrane, we analyzed the diffusion decay of 1 H spin–echo signal of (CH2) n and CH3 lipid groups in the range of 0–3.7 ppm. Diffusion decay was approximated by the sum of two exponentials in accordance with Equation (A1). The main part of the diffusion decay belongs to the tail, which is characterized by the self-diffusion coefficient DL, which is varied from 3∙10−12 m2/s to 10−11 m2/s depending on temperature and diffusion time.
Figure 65 shows the diffusion decays of 1 H spin–echo signals for different regions of integration. The relative part (population pL) of the smaller self-diffusion coefficient DL increases with the distance of the integration region of the spin–echo signal from the position of the water signal (Figure 65). This allows us to conclude that the self-diffusion coefficient DL characterizes the lateral diffusion of lipids. The components with large values of the self-diffusion coefficients, which determine the initial part of the diffusion decay, belong to the wings of the water signal.
The temperature dependence of DL obtained at a diffusion time td = 10 ms in Arrhenius coordinates is shown in Figure 66. The activation energy of lateral diffusion of lipid molecules ED was 25 ± 2.9 kJ/mol. This value is in a good agreement with the activation energy of lipid diffusion in model systems (bilayer phospholipid membranes) ED = 25 kJ/mol [179].
The lateral self-diffusion coefficient DL is decreased with increasing of diffusion time td. The dependence DL (td) is shown in Figure 67a. The effective self-diffusion coefficients DLeff, calculated from Equation (A5) is shown in Figure 67b. The estimated restriction size aL is about 1.4 μm.

3.2.4. Self-Diffusion of Fullerene C60 Derivative Molecules in a RBCs Suspension

The study of self-diffusion of fullerene C60 derivative molecules in a suspension of erythrocytes was carried out [172,173]. Water-soluble fullerene derivatives (WSFD) are promising materials for biomedical research and further pharmacological applications [180,181,182]. We have registered and analyzed the diffusion decays of molecules of the following pentasubstituted derivatives of fullerene C60 in a RBCs suspension: with adducts of 3-mercaptopropanesulfonic acid C60[S(CH2)3SO3Na]5H (1), mercaptopropionic acid C60[S(CH2)2COOK]5H (2), proline C60[N(CH2)3CHCOOK]5Cl (3).
Figure 68 and Figure 69 show examples of the spin-echo 1H NMR spectrum evolution in the pulsed field gradient experiment of RBCs suspension with and without addition of compound 1. Figure 68 shows the echo spectra at a small value of the gradient pulse amplitude g. In these spectra, for all systems, intense signals of rapidly moving water molecules are observed; Figure 68b,c also show signals of molecules of fullerene derivatives in the aqueous phase. At large gradient amplitude, the NMR signals of rapidly moving water molecules and WSFD in the aqueous phase are suppressed, which is clearly observed in Figure 69b for an aqueous solution of WSDF. However, in Figure 69c, signals of slowly moving components of the RBCs-WSFD system are observed. This fact gives opportunity to record the tail in the diffusion decay of the WSFD and to estimate the corresponding self-diffusion coefficient. Similar results were also obtained for other systems. This observation indicates a relatively slow movement of WSFD molecules in a RBCs suspension compared to their diffusion in an aqueous solution. Therefore, we can assume that the translational mobility of WSFD molecules strongly depend on their interaction with RBCs.
In Figure 70, an example of diffusion decay of compound 2 molecules in RBCs suspension at diffusion time td 10 and 300 ms is shown. The diffusion decay of compound 2 molecules in an aqueous solution is shown for comparison in insertion.
In aqueous solutions, diffusion decays of WSFD molecules are biexponential, with self-diffusion coefficients Ds1w and Ds2w. At the same time, an additional third component appears in the RBCs suspension. Thus, the mobility of WSFD molecules in aqueous solutions is characterized by two self-diffusion coefficients Ds1w and Ds2w, but in RBCs suspension—three diffusion components are observed Ds1s, Ds2s, and Ds3s (Table 12).
As it can be seen from Table 1, the self-diffusion coefficients Ds1w, Ds2w, and Ds1s, Ds2s are close to each other in aqueous solutions and in a RBCs suspension. The larger value of the coefficients Ds1s in the suspension compared to aqueous solutions is apparently due to the partial overlap of the signals of protons of water molecules and the recorded signal, which, in turn, leads to an increase in the component Ds1s. The presence of biexponentiality of diffusion decay in aqueous solutions of fullerene C60 derivatives is attributed to the formation of associates. In this case, the highest self-diffusion coefficient Ds1w corresponds to isolated molecules, and Ds2w to associated WSFD molecules.
As opposed to aqueous solutions in a RBCs suspension, the mobility of WSFD molecules is characterized by three coefficients of self-diffusion. The smallest coefficient Ds3s practically coincides with the coefficient of lateral diffusion of lipids in the erythrocyte membrane (≈7∙10−12 m2/s). The other two self-diffusion coefficients are close to the SDC of WSFD molecules in bulk aqueous solutions. Therefore, in RBCs suspension molecules of fullerene C60 derivatives are in an isolated and associated form in the aqueous phase, and are its also associated with the RBC membrane.
Due to the presence of permeability of RBCs membranes, there is an exchange between the WSFD molecules bound to the membrane and located in the aqueous phase. Quantitative information on metabolic processes can be obtained from an analysis of the dependence of the population p3 (the relative fraction of the WSFD molecules in RBCs) on the diffusion time td. Self-diffusion coefficients Ds3s, population p3(0), and lifetimes τ of molecules of C60 fullerene derivatives in erythrocytes are presented in Table 13.
Thus, using the PFG NMR method, the ability of molecules of fullerene C60 derivatives to penetrate into biological membranes was revealed. Partial coefficients of self-diffusion of WSFD molecules in RBCs were determined. It was found that in a RBCs suspension WSFD molecules are in the form of isolated and associated molecules in the aqueous phase or are associated with the cell membrane. The time of exchange of WSFD molecules between the RBC membrane and an aqueous solution was estimated.

4. Summaries

The interaction of proteins with other macromolecules or small partner molecules plays an important role in most biological processes. NMR spectroscopy techniques in combination with pulsed magnetic field gradient NMR have a unique ability to extract information about these interactions and are used more often to investigate protein systems of increasing complexity, including proteins with an internally disordered structure. A characteristic feature of such proteins is that sample an ensemble of rapidly interconverting alternative conformations ranging from random coils to more structured conformations with secondary structure and residual tertiary structure elements. It should be recognized that the amount of experimental material concerning the transport properties of proteins in solutions and mixtures remains insufficient. As shown in a review of experimental results of studying the self-diffusion of polymers and globular proteins in solutions, NMR diffusion (PFG NMR) is one of the most effective methods for studying the structural and dynamic properties of molecules in condensed systems.
1 H pulsed field gradient NMR technique was shown to be a versatile tool for monitoring the penetration of water and biologically active species (fullerene derivatives) in RBC as well as their adsorption on the cellular membrane. The analysis of the spin echo attenuation curves during the pulsed field gradient 1H NMR experiment provides an opportunity to extract the partial self-diffusion coefficients and estimate the fractions of the water-soluble fullerene derivative molecules or clusters persistent in the aqueous phase or bound to the blood cells in the mouse erythrocyte suspension. It was revealed that fullerene derivative molecules most probably get fixed on the RBC surface (or inside the membrane) due to the fact that their self-diffusion coefficient is equal to the lateral diffusion coefficient of the cell lipids and close to fullerene derivative sorbed on RBC ghost and liposome self-diffusion coefficients. The average characteristic time reflecting the exchange of the fullerene molecules between the RBC-bound state and aqueous solution was estimated. Thus, the pulsed field gradient NMR can provide essential information projected on the in-vivo behavior of potential drug candidates after their intravenous administration in animals and humans.

Author Contributions

Writing—original draft preparation of the first part, A.V.C., N.A.S.; writing—original draft preparation of the second part, I.A.A.; D.L.M.; writing—review and editing of the first part, V.I.V.; writing—review and editing of the second part, V.D.S. All authors have read and agreed to the published version of the manuscript.

Funding

NMR measurements were performed using equipment of the Multi-User Analytical Center of the Institute of Problems of Chemical Physics RAS and Science Center in Chernogolovka RAS with the support of State Assignment of the Institute of Problems of Chemical Physics RAS (state registration No 0089-2019-0010/AAAA-A19-119071190044-3).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Pulsed Field Gradient NMR Technique

In our review we will discuss the results of ion-exchange membranes and biological systems investigation by hetero nuclear high resolution NMR, NMR relaxation and pulsed field gradient NMR. High resolution and spin relaxation NMR techniques are well known and describe elsewhere. Pulsed field gradient NMR is not so widely applied; therefore, in our paper we will give the brief introduction in this technique. Pulsed field gradient NMR is attractive for studying diffusion because it is suitable for direct measurement of the diffusion coefficients and the relative fractions of the diffusant in different phases of heterogeneous systems. As a rule, the diffusion coefficients are measured by stimulated echo sequences [42,43].
The stimulated spin echo signal involves three radio- frequency 90° (RF) pulses (Figure A1). Magnetic field gradients are formed between the first and the second and after the third RF pulse. Under these conditions, the translational displacement of a species is accompanied by irreversible dephasing of the magnetization vector, which results in a decrease in the spin echo signal amplitude. Conventionally, the dependence of the echo signal intensity A on the magnetic field gradient amplitude g called the “diffusion decay” A(g) is analyzed. For slow exchange of molecules occurs between the phases, the diffusion time td = ∆−δ/3 much more compared to the diffusant lifetime in i-th phase diffusion decay is approximated by Equation (A1).
Figure A1. Stimulated echo pulse sequence with the magnetic field gradient pulses.
Figure A1. Stimulated echo pulse sequence with the magnetic field gradient pulses.
Membranes 11 00385 g0a1
Here, τ is the time interval between the first and second RF pulses, τ1 is the time interval between the second and the third ones, ∆ is the interval between the gradient pulses, δ is duration of the equivalent rectangular magnetic field gradient pulses, g is the amplitude of the magnetic field gradient pulse.
A g = A 2 τ , τ 1 , g A 2 τ , τ 1 , 0 = i = 1 m p i   e x p   ( γ 2 g 2 δ 2 t d D s i )
p i = p i   e x p   2 τ T 2 i τ 1 T 1 i / i = 1 m p i   e x p   2 τ T 2 i τ i T 1 i
i = 1 m p i = 1
where γ is the gyromagnetic ratio, td is ∆−δ/3, Dsi is the diffusion coefficient in the i-th domain (phase), pi is the relative fraction of molecules in the i-th domain (population of the i-th phase), T1i, T2i are the times of longitudinal and transverse relaxation times of nuclei in the i-th phase.
For the molecules undergoing unhindered isotropic Brownian motion, the evolution of spin echo signal is described by the following exponential Equation.
A 2 τ , τ 1 , g = A 2 τ , τ 1 , 0   exp   γ 2 g 2 δ 2 t d D s
where γ is gyromagnetic ratio, td=Δ-δ/3 is the diffusion time, and Ds is the self-diffusion coefficient, τ, τ1 and g are shown in Figure A1; A (2τ, τ1, 0) is expressed by the Equation:
A 2 τ , τ 1 , 0 = A 0 2   e x p   2 τ T 2 τ 1 T 1
where A(0) is the signal intensity after the first radio frequency (RF) pulse (Figure A1). T1 and T2 are the spin–lattice and spin–spin relaxation times, respectively. During measurement of echo signal evolution, τ and τ1 are fixed, and only the dependence of A on g is analyzed.
By analyzing the A(g) dependence under slow exchange conditions, it is possible to determine the partial diffusion coefficients and the relative proportions of the diffusant molecules in various nano- and microvolumes of membrane systems.
The procedure for calculation of diffusion coefficients in complex polymer electrolytes was described in several publications [48,50]. The modern NMR spectroscopic instrumentation allows the diffusion coefficient to be measured on various spatial scales starting from several tenths of micrometer; the measurable diffusion coefficients vary from 10−14 to 10−8 m2/s.
Pulsed field gradient NMR technique application requires sufficiently long (several milliseconds and more) spin-spin relaxation times. This complicates the measurements of diffusion coefficients of ions with quadrupole nuclei (e.g., 23 Na, 133 Cs). The majority of studies were carried out on the nuclei characterized by sufficiently long transverse relaxation times (1 H, 7 Li, 19 F).
Modern NMR spectrometers give possibility to record a Fourier transform spin echo spectra, therefore, we can measure partial self-diffusion coefficients of different molecules or molecule fragments. Note that, as in the case of relaxation NMR experiments, applying of the Fourier transform of the stimulated echo, and thus acquiring the high-resolution NMR spectrum for investigated molecules significantly simplifies interpretation of the diffusion data particularly in terms of corresponding data to different lines in spectrum therefore, to the different molecules.
The uniqueness of the pulsed field gradient NMR method lies in the possibility of studying the geometry and permeability of porous media and biological cells. Its idea is based on analysis of how the diffusion decay curve depends on the diffusion time.
In Figure A2 idealized S-shape dependence of the self-diffusion coefficient Ds of diffusant molecules on the diffusion time td in porous systems with permeable walls is shown [82,83,174,175,183,184,185]. It is known from the literature that the dependence of the SDC of water molecules on the diffusion time in porous systems and biological cells is S-shaped, as shown in Figure A2.
Figure A2. Illustration of the dependence of the apparent diffusion coefficient D on diffusion time td for the case of restricted diffusion. D0 is the free diffusion coefficient, Dp is the restricted diffusion coefficient, a is the size of restriction region (pore or cell diameter), a = (6Dtd)1/2; P is the permeability of pore (cell) [174].
Figure A2. Illustration of the dependence of the apparent diffusion coefficient D on diffusion time td for the case of restricted diffusion. D0 is the free diffusion coefficient, Dp is the restricted diffusion coefficient, a is the size of restriction region (pore or cell diameter), a = (6Dtd)1/2; P is the permeability of pore (cell) [174].
Membranes 11 00385 g0a2
At short diffusion times td (region free diffusion), the diffusant molecules are not able to reach the pore walls so a self-diffusion coefficient does not depend on td and Ds = D0, where D0 is the self-diffusion coefficient of a bulk liquid.
In the region of restricted diffusion diffusant, molecules are collided with the pore wall and Ds depend on td, becoming a decreasing function of time.
At long diffusion time the diffusant molecules as a result of the permeation effect to diffuse into adjacent pores. Their movement is averaged over the entire pore space, the measured diffusion coefficient Ds = Dp also does not depend on td (region of hindered diffusion), and Dp value is less than D0.
If self-diffusion is observed in a completely isolated pore, then the mean square displacement of molecules is limited by the size of the pore. Therefore, with the increasing of diffusion time Ds decreases as td−1 and in the limit td→∞, Ds→0. From the Einstein Equation (A4) pore size a may be calculated.
D s t d = a 2 6 t d
The permeability may be estimated from Equation (A5) [174]:
1 D p = 1 D 0 + 1 p a
Biological cells are permeable therefore in restriction region Ds(td)∝td−n, where n< 1. For quantitative determination of the limitation size and permeability, for permeable pores the scaling approach which is described in several papers was applied [82,83].
In this case, the effective self-diffusion coefficient Dseff is calculated in which dependence on the diffusion time td will be∝td−1:
D s e f f t d = D s t d D p D 0 D s t d · D 0
where D0 is the self-diffusion coefficient at td→0 (not restriction region), Ds(td) is the experimentally obtained dependence, Dp is a hindered self-diffusion coefficient (td→∞).

References

  1. Hwang, S.-T.; Kemmermeyer, K. Membranes in Separations; John Wiley & Sons: Hoboken, NJ, USA, 1975; p. 654. [Google Scholar]
  2. Chalykh, A.E. Diffuziya v Polimernykh Sistemakh (Diffusion in Polymer Systems); Khimiya: Moscow, Russia, 1987; p. 276. [Google Scholar]
  3. Shaposhnik, V.A. Kinetika Elektrodializa (Kinetics of Electro- dialysis); Voronezh State University: Voronezh, Russia, 1989; p. 187. [Google Scholar]
  4. Nikolaev, N.I. Diffuziya v Membranakh (Diffusion in Membranes); Khimiya: Moscow, Russia, 1980; p. 221. [Google Scholar]
  5. Timashev, S.F. Physical Chemistry of Membrane Processes; Ellis Horwood Series in Physical Chemistry Series; Ellis Horwood: New York, NY, USA; London, UK; Toronto, ON, Canada; Sydney, Australia; Tokyo, Japan; Singapore, 1991; p. 246. [Google Scholar]
  6. Yaroslavtsev, A.B. Khimiya Tverdogo Tela (The Chemistry of Solid); Nauchnyi Mir: Moscow, Russia, 2009; p. 312. [Google Scholar]
  7. Zabolotskii, V.I.; Nikonenko, V.V. Perenos Ionov v Membranakh (Ion Transfer in Membranes); Nauka: Moscow, Russia, 1996; p. 342. [Google Scholar]
  8. Creekmore, R.W.; Reilley, C.N. Nuclear Magnetic Resonance Study of Ion-Exchange Resins. Anal. Chem. 1970, 42, 570–575. [Google Scholar] [CrossRef]
  9. Gough, T.E.; Sharma, H.D.; Subramanian, N. Proton magnetic resonance studies of ionic solvation in ion-exchange resins. Part I. Sulfonated cation-exchange resins. Can. J. Chem. 1970, 48, 917–923. [Google Scholar] [CrossRef]
  10. Bystrov, G.S.; Grigor’eva, G.A.; Nikovaev, N.I. Nuclear Magnetic Resonance Study of Ion-exchange Resin-Solvent Systems. Russ. Chem. Rev. 1976, 45, 823–846. [Google Scholar] [CrossRef]
  11. Mank, V.V.; Kurilenko, O.D. Issledovanie Mezhmolekulyarnykh Vzaimodeistvii v Ionoobmennykh Smolakh Metodom YaMR (Investigation of Intermolecular Interactions in Ion-Exchange Resins by NMR); Naukova Dumka: Kiev, Ukraine, 1976; p. 79. [Google Scholar]
  12. Mank, V.V.; Lebovka, N.I. Voda v Dispersnykh Sistemakh (Water in Dispersion Systems); Khimiya: Moscow, Russia, 1989; p. 203. [Google Scholar]
  13. Khutsishvili, V.G.; Bogachev, Y.S.; Volkov, V.I.; Serebryanskaya, A.I.; Shapet’ko, N.N.; Timashev, S.F.; Orman, M.L. Water state investigation in sulfocationite CU-2 phase by proton magnetic resonance technique. Russ. J. Phys. Chem. 1983, 57, 2524–2527. [Google Scholar]
  14. Saldadze, G.K.; Tagirova, R.I.; Volkov, V.I.; Chizhanov, S.A. The structure of one charge cation and water mass transfer in sulfocationexchangers on NMR data. Russ. J. Phys. Chem. 1993, 67, 1818–1823. [Google Scholar]
  15. Volkov, V.I.; Saldadze, G.K.; Tagirova, R.I.; Kropotov, L.V.; Khutsishvili, V.G.; Shapet’ko, N.N. Water state and diffusion mobility in ion exchange membrane MC-40 studied by NMR. Russ. J. Phys. Chem. 1989, 63, 1005–1009. [Google Scholar]
  16. Volkov, V.I.; Chernyak, A.V.; Golubenko, D.V.; Shevlyakova, N.V.; Tverskoy, V.A.; Yaroslavtsev, A.B. Mobility of cations and water molecules in sulfocation-exchange membranes based on polyethylene and sulfonated grafted polystyrene. Membr. Membr. Technol. 2020, 1, 54–62. [Google Scholar] [CrossRef] [Green Version]
  17. Volkov, V.I.; Chernyak, A.V.; Golubenko, D.V.; Tverskoy, V.A.; Lochin, G.A.; Odjigaeva, E.S.; Yaroslavtsev, A.B. Hydration and diffusion of H+, Li+, Na+, Cs+ ions in cation-exchange membranes based on polyethylene and sulfonated-grafted polystyrene studied by NMR technique and ionic conductivity measurements. Membranes 2020, 10, 272. [Google Scholar] [CrossRef]
  18. Khutsishvili, V.G.; Bogachev, Y.S.; Volkov, V.I.; Tarasova, B.V.; Dreiman, N.A.; Shapetko, N.N.; Timashev, S.F. Water state investigation in perfluorinated sulfo cation exchange membranes by proton magnetic resonance technique. Russ. J. Phys. Chem. 1984, 58, 2633. [Google Scholar]
  19. Volkov, V.I.; Sidorenkova, E.A.; Timashev, S.F.; Lakeev, S.G. State and diffusive mobility of water in perfluorinated sulfocationite membranes according to proton magnetic resonance data. Russ. J. Phys. Chem. 1993, 67, 914–918. [Google Scholar]
  20. Skirda, V.D.; Volkov, V.I. Pulsed field gradient NMR for the molecular system physical-chemistry processes investigations. Russ. J. Phys. Chem. 1999, 73, 323–342. [Google Scholar]
  21. Volkov, V.I.; Volkov, E.V.; Sanginov, E.A.; Pavlov, A.A.; Timofeev, S.V.; Safronova, E.Y.; Stenina, I.A.; Yaroslavtsev, A.B. Diffusion mobility of alkali metals in perfluorinated sulfocationic and carboxylic membranes as probed by 1H, 7Li, 23Na, and 133Cs NMR spectroscopy. Russ. J. Inorg. Chem. 2010, 55, 318–324. [Google Scholar] [CrossRef]
  22. Bunce, N.J.; Sondheimer, S.J.; Fyle, C.A. Proton NMR Method for the Quantitative Determination of the Water Content of the Polymeric Fluorosulfonic Acid Nafion-H. Macromolecules 1986, 19, 333–339. [Google Scholar] [CrossRef]
  23. Ye, G.; Janzen, N.; Goward, G.R. Solid-State NMR Study of Two Classic Proton Conducting Polymers: Nafion and Sulfonated Poly(ether ether ketone)s. Macromolecules 2006, 39, 3283–3290. [Google Scholar] [CrossRef]
  24. Shestakov, S.L.; Pavlov, A.A.; Maksimychev, A.V.; Chernyak, A.V.; Volkov, V.I.; Timofeev, S.V. A NMR study of the hydration of sulfo and carboxyl groups in perfluorinated cation exchange membranes. Russ. J. Phys. Chem. B 2010, 4, 1005–1013. [Google Scholar] [CrossRef]
  25. Chernyak, A.V.; Vasiliev, S.G.; Avilova, I.A.; Volkov, V.I. Hydration and Water Molecules Mobility in Acid Form of Nafion Membrane Studied by 1H NMR Techniques. Appl. Magn. Reson. 2019, 5, 677–693. [Google Scholar] [CrossRef]
  26. Hammer, R.; Schönhoff, M.; Hansen, M.R. Comprehensive Picture of Water Dynamics in Nafion Membranes at Different Levels of Hydration. J. Phys. Chem. B 2019, 123, 8313–8324. [Google Scholar] [CrossRef]
  27. Komoroski, R.A.; Mauritz, K.A. A sodium-23 nuclear magnetic resonance study of ionic mobility and contact ion pairing in a perfluorosulfonate ionomer. J. Am. Chem. Soc. 1978, 100, 7487–7489. [Google Scholar] [CrossRef]
  28. Komoroski, R.A. A Multinuclear Fourier Transform NMR Study of Perfluorosulfonate Ionomers. Adv. Chem. Ser. 1980, 187, 155–168. [Google Scholar]
  29. Volkov, V.I.; Sidorenkova, E.A.; Korochkova, S.A.; Novikov, N.A.; Sokol’skaya, I.B.; Timashev, S.F. The nature of selectivity of perfluorinated sulfocation exchanger membranes to alkiline metal ions on 7Li, 23Na, 133Cs NMR. Russ. J. Phys. Chem. 1994, 68, 309–316. [Google Scholar]
  30. Volkov, V.I.; Sidorenkova, E.A.; Korochkova, S.A.; Novikov, N.A.; Sokol’skaya, I.B.; Timashev, S.F. Effect of electrolyte sorbed by nonion-exchange mechanism on the state and diffusive mobility of water and alkali metal ions in perfluorinated sulfocationic membranes from NMR data. Russ. J. Phys. Chem. 1994, 68, 559–564. [Google Scholar]
  31. Nesterov, I.A.; Volkov, V.I.; Pukhov, K.K.; Timashev, S.F. Magnetic-relaxation of 7Li+ nuclei and dynamics of movements of lithium counter-ions and water-molecules in perfluorinated sulfocationite membranes. Russ. J. Chem. Phys. 1990, 10, 1155–1162. [Google Scholar]
  32. Tromp, R.H.; Van der Maarel, J.R.C.; De Bleijser, J.; Leyte, J.C. Counter-ion dynamics in crosslinked poly( styrene sulfonate) systems studied by NMR. Biophys. Chem. 1991, 41, 81–100. [Google Scholar] [CrossRef]
  33. Halle, B.; Bratko, D.; Piculell, L. Interpretetion of Counterion Spin Relaxation in Polyelectrolyte Solutions. II. Effects of Finite Polyion Lenght. Ber. Bunsenges. Phys. Chem. 1985, 89, 1254–1260. [Google Scholar] [CrossRef]
  34. Halle, B.; Wennerstrom, H.; Piculell, L. Interpretation of Counterion Spin Relaxation in Polyelectrolyte Solutions. J. Phys. Chem. 1984, 88, 2482–2494. [Google Scholar] [CrossRef]
  35. Boyle, N.G.; McBrierty, V.J.; Douglass, D.C. The behavior of water in Nafion membranes. Macromolecules 1983, 16, 75–80. [Google Scholar] [CrossRef]
  36. Volkov, V.I.; Nesterov, I.A.; Chichagov, A.V.; Muromtsev, V.I.; Timashev, S.F. NMR methods for studying ion and molecular transport in polymer electrolytes. Russ. J. Chem. Phys. 1985, 4, 644–650. [Google Scholar] [CrossRef]
  37. Volkov, V.I.; Nesterov, I.A.; Sundukov, V.I.; Kropotov, L.V.; Timashev, S.F. The diffusion transfer of water in perfluorinated sulfocation exchange membranes as studied by pulse NMR. Russ. J. Chem. Phys. 1989, 8, 209. [Google Scholar]
  38. Volkov, V.I.; Vasilyak, S.L.; Park, I.-W.; Kim, H.J.; Ju, H.; Volkov, E.V.; Choh, S.H. Water Behavior in Perfluorinated Ion-Exchange Membranes. Appl. Magn. Reson. 2003, 25, 43–53. [Google Scholar] [CrossRef]
  39. Volkov, V.I.; Volkov, E.V.; Timashev, S.F. The Mechanism for Ionic and Water Transport in Nafion Membranes from Resonance Data. In Magnetic Resonance in Colloid and Interface Science; Kluwer Academic Publishers: Amsterdam, The Netherlands, 2002; p. 267. [Google Scholar]
  40. Volkov, V.I.; Volkov, E.V.; Vasilyak, S.L.; Hong, Y.S.; Lee, C.H. The Ionic and Molecular Transport in Polymeric and Biological Membranes on Magnetic Resonance Data. In Fluid Transport in Nanoporous Materials; Springer: New York, NY, USA, 2006; p. 48. [Google Scholar]
  41. Volkov, V.I.; Eliseev, Y.G.; Kirsh, Y.E.; Fedotov, Y.A.; Timashev, S.F. Water and lithium self-diffusion in bi-sulfocontaining aromatic polyamides. Russ. J. Phys. Chem. 1992, 66, 1618–1622. [Google Scholar]
  42. Stejskal, E.O.; Tanner, J.E. Spin Diffusion Measurements: Spin Echoes in the Presence of a Time Dependent Field Gradient. J. Chem. Phys. 1965, 42, 288. [Google Scholar] [CrossRef] [Green Version]
  43. Maklakov, A.; Skirda, V.; Fatkullin, N. Selfdiffusion in Polymer Solutions and Melts (Samodiffisia v Rastvorakh i Rasplavakh Polymerov); Soloviev, U., Chalykh, A., Eds.; Kazan University Press: Kazan, USSR, 1987; p. 224. (In Russian) [Google Scholar]
  44. Saldadze, G.K.; Tagirova, R.I.; Volkov, V.I.; Chizhanov, S.A. Self-Diffusion Of Li+ Counter-ions And Water Molecules In CU-23 Macroporous Sulfocationites According To Pulse Nuclear-Magnetic-Resonance. Russ. J. Phys. Chem. 1993, 67, 1941–1943. [Google Scholar]
  45. Saldadze, G.K.; Volkov, V.I.; Tagirova, R.I.; Chizhanov, S.A. Water Diffusion in Heterogeneous Sulfocationite Systems According to Proton Magnetic-Resonance. Russ. J. Phys. Chem. 1993, 67, 773–776. [Google Scholar]
  46. Volkov, V.I.; Korochtkova, S.A.; Nesterov, I.A.; Kirsh, Y.E.; Timashev, S.F. Diffusion Mobility Of Water-Molecules In Cation-Exchange Membranes Based On Sulfonate-Containing Polyphenylenephthalamides. Russ. J. Phys. Chem. 1994, 68, 1310–1316. [Google Scholar]
  47. Volkov, V.I.; Korochtkova, S.A.; Timashev, S.F. Characteristics of Self-Diffusion of Aliphatic Monoatomic Alcohols in Perfluorinated Sulfocationite Membranes. Russ. J. Phys. Chem. 1995, 69, 1124–1129. [Google Scholar]
  48. Volkov, V.I.; Korotchkova, S.A.; Ohya, H.; Guo, Q. Self-diffusion of water-ethanol mixtures in polyacrylic acid-polysulfone composite membranes obtained by pulsed-field gradient nuclear magnetic resonance spectroscopy. J. Membr. Sci. 1995, 100, 273–286. [Google Scholar] [CrossRef]
  49. Sokolova, S.A.; Djakonova, O.V.; Kotov, V.V.; Hong, Y.S.; Volkov, V.I.; Lee, C.H. The self-diffusion of water and membrane structure in the new type of cation-exchange polyamide-acid membrane. Magn. Res. Imag. 2001, 19, 588–589. [Google Scholar] [CrossRef]
  50. Volkov, V.I.; Korotchkova, S.A.; Nesterov, I.A.; Ohya, H.; Guo, Q.; Huang, J.; Chen, J. The self-diffusion of water and ethanol in cellulose derivative membranes and particles with the pulsed field gradient NMR data. J. Membr. Sci. 1996, 110, 1–11. [Google Scholar] [CrossRef]
  51. Volkov, V.I.; Pavlov, A.A.; Sanginov, E.A. Ionic transport mechanism in cation-exchange membranes studied by NMR technique. Solid State Ion. 2011, 188, 124–128. [Google Scholar] [CrossRef]
  52. Volkov, V.I.; Popkov, Y.M.; Timashev, S.F.; Bessarabov, D.G.; Sanderson, R.D.; Twardowski, Z. Self-diffusion of water and fluorine ions in anion-exchange polymeric materials (membranes and resin) as determined by pulsed-field gradient nuclear magnetic resonance spectroscopy. J. Membr. Sci. 2000, 180, 1–13. [Google Scholar] [CrossRef]
  53. Kotov, V.V.; Dyakonova, O.V.; Sokolova, S.A.; Volkov, V.I. Struktura i elektrokhimicheskie svoistva kationoobmennykh membran na osnove chastichno imidizirovannoi poliamidokisloty [Structure and electrochemical properties of cation exchange membranes based on partially imidizated polyamide acid. Elektrokhimiya 2002, 38, 994–997. [Google Scholar]
  54. Volkov, V.I.; Kotov, V.V.; Netesova, G.A. The self-diffusion of water and saturated aliphatic alcohols in cation-exchange membranes. Russ. J. Phys. Chem. 2008, 82, 1184–1188. [Google Scholar] [CrossRef]
  55. Ponomarev, A.N.; Dobrovol’skii, Y.A.; Abdrashitov, E.F.; Bokun, V.C.; Sanginov, E.A.; Volkov, E.V.; Volkov, V.I. The new approach to the modification of the per fluorinated ion-ex changing membranes promising in hydrogen energy. Bull. RAS Energetika 2008, 3, 124–134. [Google Scholar]
  56. Volkov, V.I.; Dobrovolsky, Y.A.; Nurmiev, M.S.; Sanginov, E.A.; Volkov, E.V.; Pisareva, A.V. Charge and mass transport in the phenol-2,4-disulfonic acid-polyvinyl alcohol ion exchange membranes studied by pulsed field gradient NMR and impedance spectroscopy. Solid State Ion. 2008, 179, 148–153. [Google Scholar] [CrossRef]
  57. Volkov, V.I.; Ponomarev, A.N.; Yaroslavtsev, A.B.; Sanginov, E.A.; Pavlov, A.A. NMR investigation of proton conductive ion-exchange membranes transport behaviour. ISJAEE 2008, 2, 101–106. [Google Scholar]
  58. Voropaeva, E.Y.; Sanginov, E.A.; Volkov, V.I.; Pavlov, A.A.; Shalimov, A.S.; Stenina, I.A.; Yaroslavtsev, A.B. Transport Properties of MF-4SC Membranes Modified with Inorganic Dopants. Russ. J. Inorg. Chem. 2008, 53, 1536–1541. [Google Scholar] [CrossRef]
  59. Volkov, V.I.; Pavlov, A.A.; Fedotov, Y.A.; Marinin, A.A. Self-diffusion of water and alkaline cations in bisulfur-containing aromatic polyamides-water systems. Russ. J. Phys. Chem. 2010, 84, 1705–1711. [Google Scholar] [CrossRef]
  60. Volkov, V.I.; Volkov, E.V.; Timofeev, S.V.; Sanginov, E.A.; Pavlov, A.A.; Safronova, E.Y.; Stenina, I.A.; Yaroslavtsev, A.B. Water self-diffusion and ionic conductivity in perfluorinated sulfocationic membranes MF4SK. Russ. J. Inorg. Chem. 2010, 55, 355–360. [Google Scholar] [CrossRef]
  61. Safronova, E.Y.; Volkov, V.I.; Yaroslavtsev, A.B. Ion mobility and conductivity of hybrid ion-exchange membranes incorporating inorganic nanoparticles. Solid State Ion. 2011, 188, 129–131. [Google Scholar] [CrossRef]
  62. Volkov, V.I.; Volkov, E.V. Ionic and molecular Self-Diffusion in Ion-Exchange materials for fuel energetics studied by pulsed field gradient NMR. Appl. Magn. Reson. 2005, 29, 495–501. [Google Scholar] [CrossRef]
  63. Tsushima, S.; Teranishi, K.; Hirai, S. Water diffusion measurement in fuel-cell SPE membrane by NMR. Energy 2005, 30, 235–245. [Google Scholar] [CrossRef]
  64. Baglio, V.; Arico, A.S.; Antonucci, V.; Nicotera, I.; Oliviero, C.; Coppola, L.; Antonucci, P.L. An NMR spectroscopic study of water and methanol transport properties in DMFC composite membranes: Influence on the electrochemical behaviour. J. Power Sources 2006, 163, 52–55. [Google Scholar] [CrossRef]
  65. Edmondson, C.A.; Stallworth, P.E.; Chapman, M.E.; Fontanella, J.J.; Wintersgill, M.C.; Chung, S.H.; Greenbaum, S.G. Complex impedance studies of proton-conducting membranes. Solid State Ion. 2000, 135, 419–423. [Google Scholar] [CrossRef]
  66. Every, H.A.; Hickner, M.A.; McGrath, J.E.; Zawodzinski, T.A., Jr. An NMR study of methanol diffusion in polymer electrolyte fuel cell membranes. J. Membr. Sci. 2005, 250, 183–188. [Google Scholar] [CrossRef]
  67. Jayakody, J.R.P.; Khalfan, A.; Mananga, E.S.; Greenbaum, S.G.; Dang, T.D.; Mantz, R. NMR investigation of water and methanol transport in sulfonated polyareylenethioethersulfones for fuel cell applications. J. Power Sources 2006, 156, 195–199. [Google Scholar] [CrossRef]
  68. Freger, V.; Korin, E.; Wisniak, J.; Korngold, E.; Ise, M.; Kreuer, K.D. Diffusion of water and ethanol in ion-exchange membranes: Limits of the geometric approach. J. Membr. Sci. 1999, 160, 213–224. [Google Scholar] [CrossRef]
  69. Gong, X.; Bandis, A.; Tao, A.; Meresi, G.; Wang, Y.; Inglefield, P.T.; Jones, A.A.; Wen, W.-Y. Self-diffusion of water, ethanol and decafluropentane in perfluorosulfonate ionomer by pulse field gradient NMR. Polymer 2001, 42, 6485–6492. [Google Scholar] [CrossRef]
  70. Khalfan, A.N.; Sanchez, L.M.; Kodiweera, C.; Greenbaum, S.G.; Bai, Z.; Dang, T.D. Water and proton transport properties of hexafluorinated sulfonated poly(arylenethioethersulfone) copolymers for applications to proton exchange membrane fuel cells. J. Power Sources 2007, 173, 853–859. [Google Scholar] [CrossRef]
  71. Huang, Y.F.; Chuang, L.C.; Kannan, A.M.; Lin, C.W. Proton-conducting membranes with high selectivity from cross-linked poly (vinyl alcohol) and poly (vinyl pyrrolidone) for direct methanol fuel cell applications. J. Power Sources 2009, 186, 22–28. [Google Scholar] [CrossRef]
  72. Hallberg, F.; Vernersson, T.; Pettersson, E.T.; Dvinskikh, S.V.; Lindbergh, G.; Furo, I. Electrokinetic transport of water and methanol in Nafion membranes as observed by NMR spectroscopy. Electrochim. Acta 2010, 55, 3542–3549. [Google Scholar] [CrossRef]
  73. Kritskaya, D.A.; Abdrashidov, E.F.; Bokun, V.C.; Ponomarev, A.N.; Chernyak, A.V.; Vasil’ev, S.G.; Volkov, V.I. An NMR Study of Sorption-Diffusion Properties of MF-4SK-Carbon Composite Membranes in Aqueous Methanol Solutions. Pet. Chem. 2011, 1, 266–272. [Google Scholar] [CrossRef]
  74. Guillermo, A.; Gebel, G.; Mendil-Jakani, H.; Pinton, E. NMR and pulsed field gradient NMR approach of water sorption properties in Nafion at low temperature. J. Phys. Chem. B 2009, 113, 6710–6717. [Google Scholar] [CrossRef]
  75. Ma, Z.; Jiang, R.; Myers, M.E.; Thompson, E.L.; Gittleman, C.S. NMR studies of proton transport in fuel cell membranes at sub-freezing conditions. J. Mater. Chem. 2011, 21, 9302–9311. [Google Scholar] [CrossRef]
  76. Kusoglu, A.; Weber, A.Z. New Insights into Perfluorinated Sulfonic-Acid Ionomers. Chem. Rev. 2017, 117, 987–1104. [Google Scholar] [CrossRef]
  77. Zhao, Q.; Majsztrik, P.; Benziger, J. Diffusion and Interfacial Transport of Water in Nafion. J. Phys. Chem. B 2011, 115, 2717–2727. [Google Scholar] [CrossRef]
  78. Zawodzinski, T.A.; Derouin, C.; Radzinski, S.; Sherman, R.J.; Smith, V.T.; Springer, T.E.; Gottesfeld, S. Water Uptake by and Transport Through Nafion 117 Membranes. J. Electrochem. Soc. 1993, 140, 1041–1047. [Google Scholar] [CrossRef]
  79. Zawodzinski, T.A.; Springer, T.E.; Davey, J.; Jestel, R.; Lopez, C.; Valerio, J.; Gottesfeld, S. A Comparative Study of Water Uptake by and Transport through Ionomeric Fuel Cell Membranes. J. Electrochem. Soc. 1993, 140, 1981–1985. [Google Scholar] [CrossRef]
  80. Cappadonia, M.; Erning, J.W.; Niaki, S.M.S.; Stimming, U. Conductance of Nafion 117 membranes as a function of temperature and water content. Solid State Ion. 1995, 77, 65–69. [Google Scholar] [CrossRef]
  81. Nicotera, I.; Coppola, L.; Rossi, C.O.; Youssry, M.; Ranieri, G.A. NMR Investigation of the Dynamics of Confined Water in Nafion-Based Electrolyte Membranes at Subfreezing Temperatures. J. Phys. Chem. B 2009, 113, 13935–13941. [Google Scholar] [CrossRef]
  82. Suh, K.-J.; Hong, Y.-S.; Skirda, V.D.; Volkov, V.I.; Lee, C.-Y.; Lee, C.-H. Water self-diffusion behavior in yeast cells studied by pulsed field gradient NMR. Biophys. Chem. 2003, 104, 121–130. [Google Scholar] [CrossRef]
  83. Cho, C.-H.; Hong, Y.-S.; Kang, K.; Volkov, V.I.; Skirda, V.D.; Lee, C.-Y.J.; Lee, C.H. Water self-diffusion in Chlorella sp. studied by pulse field gradient NMR. Magn. Reson. Imaging 2003, 21, 1009–1017. [Google Scholar] [CrossRef]
  84. Ozerin, A.N.; Rebrov, A.V.; Yakunin, A.N.; Bogovtseva, L.P.; Timashev, S.F.; Bakeev, N.F. Structural changes in perfluorinated membranes in the processes of saponification and orientation stretching. Polym. Sci. USSR 1986, 28, 275–282. [Google Scholar] [CrossRef]
  85. Ozerin, A.N.; Rebrov, A.V.; Yakunin, A.N.; Bessonova, N.P.; Dreiman, N.A.; Sokolov, L.F.; Bakeev, N.F. Effect of low molecular weight fractions on the structure of the perfluorinated sulphocationite ion-exchange membrane. Polym. Sci. USSR 1986, 28, 2559–2565. [Google Scholar] [CrossRef]
  86. Eisenberg, A.; King, A.E. Ion Containing Polymers: Physical Properties and Structure; Academic Press: New York, NY, USA, 1977; p. 287. [Google Scholar]
  87. Yaroslavtsev, A.B. Perfluorinated ion-exchange membranes. Polym. Sci. Ser. A 2013, 55, 674–698. [Google Scholar] [CrossRef]
  88. Fujimura, M.; Hashimoto, T.; Kawai, H. Small-angle x-ray scattering study of perfluorinated ionomer membranes. 2. Models for ionic scattering maximum. Macromolecules 1982, 15, 136–144. [Google Scholar] [CrossRef]
  89. Gierke, T.D.; Munn, G.E.; Wilson, F.C. The morphology in nafion perfluorinated membrane products, as determined by wide-and small-angle x-ray studies. J. Polym. Sci. Polym. Phys. Ed. 1981, 19, 1687. [Google Scholar] [CrossRef]
  90. Hsu, W.Y.; Gierke, T.D. Ion transport and clustering in nafion perfluorinated membranes. J. Membr. Sci. 1983, 13, 307–326. [Google Scholar] [CrossRef]
  91. Mauritz, K.A.; Moore, R.B. State of Understanding of Nafion. Chem. Rev. 2004, 104, 4535–4586. [Google Scholar] [CrossRef]
  92. Mologin, D.A.; Khalatur, P.G.; Khokhlov, A.R. Structural Organization of Water-Containing Nafion: A Cellular-Automaton-Based Simulation. Macromol. Theor. Simul. 2001, 11, 587–607. [Google Scholar] [CrossRef]
  93. Chibirova, F.K.; Zakharin, D.S.; Sedov, V.E.; Timashev, S.F.; Popkov, Y.M.; Reiman, S.I. Characteristics Of The Perfluorinated Ion-Exchange Sulfocation Membrane-Structure According To Mossbauer-Spectroscopy Data. Russ. J. Chem. Phys. 1987, 6, 1137–1145. [Google Scholar]
  94. Volfkovich, Y.M.; Luzhin, V.K.; Vanyulin, A.N.; Shkolnikov, E.I.; Blinov, I.A. Application of the standard porometry method to investigation of the porous structure of ion-exchange membranes. Elektrokhimiya 1984, 20, 656. [Google Scholar]
  95. Rubatat, L.; Diat, O.; Gebel, G. New structural model for Nafion® membranes. In Laboratoire Léon Brillouin Scientific Report 2003–2004; Imprilith S.A.: Pringy, France, 2005; p. 218. [Google Scholar]
  96. Volkov, V.I.; Gladkikh, S.N.; Timashev, S.F.; Smirnov, N.I.; Muromtsev, V.I.; Andreeva, A.I.; Timofeev, S.V. Perfluorinated sulfocation exchange membrane structure studied by NMR relaxation 19F and paramagnetic probe techniques. Russ. J. Chem. Phys. 1983, 13, 49–53. [Google Scholar]
  97. Volkov, V.I.; Timashev, S.F. Magnetic-Resonance Method in The Study of Perfluorinated Ion-Exchange Membranes. Russ. J. Phys. Chem. 1989, 63, 209–220. [Google Scholar]
  98. Volkov, V.I.; Muromtsev, V.I.; Pukhov, K.K.; Timashev, S.F.; Balaev, G.A. Issledovanie structury perftorirovannyh kationitovyh membran metodom matrichnogo dvoinogo elektronno-yadernogo resonansa. Dokl. Akad. Nauk SSSR 1984, 276, 395. (In Russian) [Google Scholar]
  99. Volkov, V.I.; Pukhov, K.K.; Choh, S.H.; Park, I.-W.; Volkov, E.V. Matrix ENDOR of ion-exchange systems. Appl. Magn. Reson. 2003, 24, 177. [Google Scholar] [CrossRef]
  100. Vishnevskaya, G.P.; Molochnikov, L.S.; Safin, R.S. EPR v Ionitakh (ESR in Ionites); Nauka: Moscow, Russia, 1992; p. 166. [Google Scholar]
  101. Hindman, J.C. Nuclear Magnetic Resonance Effects in Aqueous Solutions of 1-1 Electrolytes. J. Chem. Phys. 1962, 36, 1000. [Google Scholar] [CrossRef]
  102. Knapp, P.S.; Waite, R.O.; Malinowski, E.R. NMR Studies of Aqueous Electrolyte Solutions. Hydration of Sodium Perchlorate, Hydrochloric Acid, and Perchloric Acid from Temperature Effects on Proton Shifts. J. Chem. Phys. 1968, 49, 5459–5463. [Google Scholar] [CrossRef]
  103. Creekmore, R.W.; Reilley, C.N. Nuclear Magnetic Resonance Determination of Hydration Numbers of Electrolytes in Concentrated Aqueous Solutions. J. Phys. Chem. 1969, 73, 1563–1568. [Google Scholar] [CrossRef]
  104. Malinowski, E.R.; Knapp, P.S.; Feuer, B. NMR Studies of Aqueous Electrolyte Solutions. I. Hydration Number of NaCI Determined from Temperature Effects on Proton Shift. J. Chem. Phys. 1966, 45, 4274–4279. [Google Scholar] [CrossRef]
  105. Volkov, V.I.; Marinin, A.A. NMR methods for studying ion and molecular transport in polymer electrolytes. Russ. Chem. Rev. 2013, 82, 248–272. [Google Scholar] [CrossRef]
  106. Maldonado, L.; Perrin, J.-C.; Dillet, J.; Lottin, O. Characterization of polymer electrolyte Nafion membranes: Influence of temperature, heat treatment and drying protocol on sorption and transport properties. J. Membr. Sci. 2012, 389, 43–56. [Google Scholar] [CrossRef]
  107. Iwamoto, R.; Oguro, K.; Sato, M.; Iseki, Y. Water in perfluorinated, sulfonic acid Nafion membranes. J. Phys. Chem. B 2002, 106, 6973–6979. [Google Scholar] [CrossRef]
  108. Batamack, P.; Fraissard, J. Proton NMR studies on concentrated aqueous silfuric acid solutions and Nafion-H. Catal. Lett. 1997, 49, 129–136. [Google Scholar] [CrossRef]
  109. Chernyak, A.V.; Gnezdilov, O.I.; Skirda, V.D.; Volkov, V.I. Hydration, self-diffusion and ionic conductivity of Li+, Na+ and Cs+ cations in Nafion membrane studied by NMR. Solid State Ion. 2021, 364, 115627. [Google Scholar]
  110. Marcus, Y. Thermodynamics of solvation of ions. Part 6.—The standard partial molar volumes of aqueous ions at 298.15 K. J. Chem. Soc. 1993, 89, 713–718. [Google Scholar] [CrossRef]
  111. Moster, A.L.; Mitchell, B.S. Hydration and Proton Conduction in Nafion/Ceramic Nanocomposite Membranes Produced by Solid-State Processing of Powders from Mechanical Attrition. Appl. Polym. Sci. 2009, 113, 243. [Google Scholar] [CrossRef]
  112. Thompson, E.L.; Capehart, T.W.; Fuller, T.J.; Jorne, J. Investigation of Low-Temperature Proton Transport in Nafion Using Direct Current Conductivity and Differential Scanning Calorimetry. J. Electrochem. Soc. 2006, 153, A2351–A2362. [Google Scholar] [CrossRef]
  113. Kononenko, N.; Nikonenko, V.; Grande, D.; Larchet, C.; Dammak, L.; Fomenko, M.; Volfkovich, Y. Porous structure of ion exchange membranes investigated by various techniques. Adv. Colloid Interface Sci. 2017, 246, 196–216. [Google Scholar] [CrossRef]
  114. Plazanet, M.; Sacchetti, F.; Petrillo, C.; Deme, B.; Bartolini, P.; Torre, R. Water in a Palymeric Electrolyte Membrane: Sorption/Desorption and Freezing phenomena. J. Membr. Sci. 2014, 453, 419. [Google Scholar] [CrossRef] [Green Version]
  115. Galitskaya, E.; Privalov, A.F.; Weigler, M.; Vogel, M.; Kashin, A.; Ryzhkin, M.; Sinitsyn, V. NMR diffusion studies of proton-exchange membranes in wide temperature range. J. Membr. Sci. 2020, 596, 117691. [Google Scholar] [CrossRef]
  116. Yaroslavtsev, A.B. Solid electrolytes: Main prospects of research and development. Russ. Chem. Rev. 2016, 85, 1255–1276. [Google Scholar] [CrossRef]
  117. Yaroslavtsev, A.B.; Stenina, I.A.; Golubenko, D.V. Membrane materials for energy production and storage. Pure Appl. Chem. 2020, 92, 1147–1157. [Google Scholar] [CrossRef]
  118. Lide, D.R. CRC Handbook of Chemistry and Physics; CRC Press: Boca Raton, FL, USA, 2004; p. 2660. [Google Scholar]
  119. Stenina, I.A.; Yaroslavtsev, A.B. Ionic Mobility in Ion-Exchange Membranes. Membranes 2021, 11, 198. [Google Scholar] [CrossRef] [PubMed]
  120. Alberts, B.; Johnson, A.; Lewis, J.; Raff, M.; Roberts, K.; Walter, P. Molecular Biology of the Cell, 5th ed.; Garland Science: New York, NY, USA, 2008; p. 1616. [Google Scholar]
  121. Van Den Bedem, H.; Fraser, J.S. Integrative, dynamic structural biology at atomic resolution—It’s about time. Nat. Methods 2015, 12, 307. [Google Scholar] [CrossRef] [Green Version]
  122. Mao, A.H.; Lyle, N.; Pappu, R.V. Describing sequence–ensemble relationships for intrinsically disordered proteins. Biochem. J. 2013, 449, 307–318. [Google Scholar] [CrossRef] [Green Version]
  123. Wright, P.E.; Dyson, H.J. Intrinsically disordered proteins in cellular signaling and regulation. Nat. Rev. Mol. Cell Biol. 2015, 16, 18–29. [Google Scholar] [CrossRef]
  124. Hong, M.; Zhang, Y.; Hu, F. Membrane protein structure and dynamics from NMR spectroscopy. Annu. Rev. Phys. Chem. 2012, 63, 1–24. [Google Scholar] [CrossRef] [Green Version]
  125. Buffy, J.J.; Waring, A.J.; Hong, M. Determination of peptide oligomerization in lipid bilayers using 19F spin diffusion NMR. J. Am. Chem. Soc. 2005, 127, 4477–4483. [Google Scholar] [CrossRef]
  126. Frei, J.N.; Broadhurst, R.W.; Bostock, M.J.; Solt, A.; Jones, A.J.; Gabriel, F.; Tandale, A.; Shrestha, B.; Nietlispach, D. Conformational plasticity of ligand-bound and ternary GPCR complexes studied by 19F NMR of the β 1-adrenergic receptor. Nat. Commun. 2020, 11, 1–14. [Google Scholar] [CrossRef]
  127. McFeeters, R.L.; Oswald, R.E. Structural mobility of the extracellular ligand-binding core of an ionotropic glutamate receptor. Analysis of NMR relaxation dynamics. Biochemistry 2002, 41, 10472–10481. [Google Scholar] [CrossRef]
  128. de Gennes, P.G. Scaling Concepts in Polymer Physics, 1st ed.; Cornell University Press: Ithaca, NY, USA, 1979; p. 324. [Google Scholar]
  129. Grosberg, A.Y.; Khokhlov, A.R.; Stanley, H.E.; Mallinckrodt, A.J.; McKay, S. Statistical physics of macromolecules. Comput. Phys. 1995, 9, 171–172. [Google Scholar] [CrossRef] [Green Version]
  130. Nesmelova, I.V.; Skirda, V.D.; Fedotov, V.D. Generalized concentration dependence of globular protein self-diffusion coefficients in aqueous solutions. Biopolymers 2002, 63, 132–140. [Google Scholar] [CrossRef]
  131. Brune, D.; Kim, S. Predicting protein diffusion coefficients. Proc. Natl. Acad. Sci. USA 1993, 90, 3835–3839. [Google Scholar] [CrossRef] [Green Version]
  132. Skirda, V.D.; Sundukov, V.I.; Maklakov, A.I.; Zgadzai, O.E.; Gafurov, I.R.; Vasiljev, G.I. On the generalized concentration and molecular mass dependencies of macromolecular self-diffusion in polymer solutions. Polymer 1988, 29, 1294–1300. [Google Scholar] [CrossRef]
  133. Tokuyama, M.; Oppenheim, I. Dynamics of hard-sphere suspensions. Phys. Rev. E 1994, 50, R16. [Google Scholar] [CrossRef] [Green Version]
  134. Lowe, C.P.; Masters, A.J. The viscoelastic response of Brownian suspensions. J. Chem. Phys. 1999, 111, 8708–8720. [Google Scholar] [CrossRef]
  135. Dyson, H.J.; Wright, P.E. Intrinsically unstructured proteins and their functions. Nat. Rev. Mol. Cell Biol. 2005, 6, 197–208. [Google Scholar] [CrossRef]
  136. Ward, J.J.; Sodhi, J.S.; McGuffin, L.J.; Buxton, B.F.; Jones, D.T. Prediction and Functional Analysis of Native Disorder in Proteins from the Three Kingdoms of Life. J. Mol. Biol. 2004, 337, 635–645. [Google Scholar] [CrossRef]
  137. Li, J.I.E.; Zhu, M.I.N.; Manning-Bog, A.B.; Di Monte, D.A.; Fink, A.L. Dopamine and L-dopa disaggregate amyloid fibrils: Implications for Parkinson’s and Alzheimer’s disease. FASEB J. 2004, 18, 962–964. [Google Scholar] [CrossRef]
  138. Lu, J.X.; Qiang, W.; Yau, W.M.; Schwieters, C.D.; Meredith, S.C.; Tycko, R. Molecular structure of β-amyloid fibrils in Alzheimer’s disease brain tissue. Cell 2013, 154, 1257–1268. [Google Scholar] [CrossRef] [Green Version]
  139. Melnikova, D.L.; Skirda, V.D.; Nesmelova, I.V. Effect of Intrinsic Disorder and Self-Association on the Translational Diffusion of Proteins: The Case of α-Casein. J. Phys. Chem. B 2017, 121, 2980–2988. [Google Scholar] [CrossRef]
  140. Gafurov, I.R.; Skirda, V.D.; Maklakov, A.I.; Perevezentseva, S.P.; Zimkin, Y.A. NMR Study of the Structure of Aqueous Gelatine Gels and the Process of Their Formation. Polym. Sci. USSR 1989, 31, 292–300. [Google Scholar] [CrossRef]
  141. Skirda, V.D.; Gafurov, I.R.; Maklakov, A.I.; Doroginitskii, M.M.; Flyaisher, G. Translational Mobility of Macromolecules in Networks. Vysokomol. Soedin. B 1988, 30, 313–314. [Google Scholar]
  142. Skirda, V.D.; Doroginizkij, M.M.; Sundukov, V.I.; Maklakov, A.I.; Fleischer, G.; Häusler, K.G.; Straube, E. Detection of Spatial Fluctuations of Segments in Swollen Polybutadiene Networks by Nuclear Magnetic Resonance Pulsed Field Gradient Technique. Makromol. Chem. 1988, 9, 603–607. [Google Scholar] [CrossRef]
  143. Marchesseau, S.; Mani, J.C.; Martineau, P.; Roquet, F.; Cuq, J.L.; Pugniere, M. Casein Interactions Studied by the Surface Plasmon Resonance Technique. J. Dairy Sci. 2002, 85, 2711–2721. [Google Scholar] [CrossRef]
  144. Horne, D.S. Casein Interactions: Casting Light on the Black Boxes, the Structure in Dairy Products. Int. Dairy J. 1998, 8, 171–177. [Google Scholar] [CrossRef]
  145. Euston, S.R.; Naser, M.A. Simulating the Equation of State of Model Globular Proteins Adsorbed at a Surface. Langmuir 2005, 21, 4227–4235. [Google Scholar] [CrossRef]
  146. Burchard, W. Filamentous Supramolecular Structures. Macromol. Symp. 2010, 295, 49–58. [Google Scholar] [CrossRef]
  147. Antonov, V.F. Membrannyi transport. Sorosovskii obrazovatel’nyi zhurnal 1997, 6, 14–20. (In Russian) [Google Scholar]
  148. Agre, P.; Preston, G.; Smith, B.; Jung, J.; Raina, S.; Moon, C.; Guggino, W.; Nielsen, S. Aquaporin CHIP: The Archetypal Molecular Water Channel. Am. J. Physiol. 1993, 265, F463–F476. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  149. Agre, P. The Aquaporin Water Channels. Proc. Am. Thorac. Soc. 2006, 3, 5–13. [Google Scholar] [CrossRef] [PubMed]
  150. Kozono, D.; Yasui, M.; King, L.S.; Agre, P. Aquaporin water channels: Atomic structure molecular dynamics meet clinical medicine. J. Clin. Investig. 2002, 109, 1395–1399. [Google Scholar] [CrossRef] [PubMed]
  151. Carbrey, J.M.; Agre, P. Discovery of the Aquaporins and Development of the Field. Aquaporins 2009, 190, 3–28. [Google Scholar]
  152. Outhred, R.; Conlon, T. The volume dependence of the erythrocyte water diffusion permeability. BBA—Biomembranes 1973, 318, 446–450. [Google Scholar] [CrossRef]
  153. Macey, R.I. Transport of water and urea in red blood cells. Am. J. Physiol. Cell Physiol. 1984, 246, C195–C203. [Google Scholar] [CrossRef]
  154. Boss, D.; Kühn, J.; Jourdain, P.; Depeursinge, C.; Magistretti, P.J.; Marqueta, P. Measurement of absolute cell volume, osmotic membrane water permeability, and refractive index of transmembrane water and solute flux by digital holographic microscopy. J. Biomed. Opt. 2013, 18, 036007-1–036007-9. [Google Scholar] [CrossRef]
  155. Zhurova, M.; Olivieri, A.; Holt, A.; Acker, J. A method to measure permeability of red blood cell membrane to water and solutes using intrinsic fluorescence. Clin. Chim. Acta 2014, 431, 103–110. [Google Scholar] [CrossRef]
  156. Benga, G. Comparative studies of water permeability of red blood cells from humans and over 30 animal species: An overview of 20 years of collaboration with Philip Kuchel. Eur. Biophys. J. 2013, 42, 33–46. [Google Scholar] [CrossRef]
  157. Conlon, T.; Outhred, T. Water diffusion permeability of erythrocytes using an NMR technique. Biochim. Biophys. Acta 1972, 288, 354–361. [Google Scholar] [CrossRef]
  158. Morariu, V.V.; Benga, G. Water Diffusion through Erythrocyte Membranes in Normal and Pathological Subjects: Nuclear Magnetic Resonance Investigations. In Membrane Processes; Benga, G., Ed.; Springer: New York, NY, USA, 1984; pp. 121–140. [Google Scholar]
  159. Himmelreich, U.; Chapman, B.E.; Kuchel, P.W. Membrane permeability of formate in human erythrocytes: NMR measurements. Eur. Biophys. J. 1999, 28, 158–165. [Google Scholar] [CrossRef]
  160. Chapman, B.E.; Kuchel, P.W. Diffusion-Ordered (DOSY) NMR of Solute Exchange Across the Human Erythrocyte Membrane: Including an application of PostWidder Laplace inversion. Diffus. Fundam. 2007, 4, 8.1–8.15. [Google Scholar]
  161. Gianolio, E.; Ferrauto, G.; Gregorio, E.D.; Aime, S. Re-evaluation of the water exchange lifetime value across red blood. Biochim. Biophys. Acta 2016, 1858, 627–631. [Google Scholar] [CrossRef]
  162. Herbst, M.-D.; Goldstein, J.-H. A review of water diffusion measurement by NMR in human red blood cells. Am. J. Physiol. 1989, 256, C1097–C1104. [Google Scholar] [CrossRef]
  163. Anisimov, A.; Suslov, M. Estimating the MRI Contrasting Agents Effect on Water Permeability of Plant Cell Membranes Using the 1H NMR Gradient Technique. Appl. Magn. Reson. 2021, 52, 235–246. [Google Scholar] [CrossRef]
  164. Andrasko, J. Water diffusion permeability of human erythrocytes studied by a pulsed gradient NMR technique. Biochim. Biophys. Acta 1976, 428, 304–311. [Google Scholar] [CrossRef]
  165. Waldeck, A.R.; Nouri-Sorkhabi, M.H.; Sullivan, D.R.; Kuchel, P.W. Effects of cholesterol on transmembrane water diffusion in human erythrocytes measured using pulsed field gradient NMR. Biophys. Chem. 1995, 55, 197–208. [Google Scholar] [CrossRef]
  166. Waldeck, A.R.; Kuchel, P.W.; Lennon, A.J.; Chapman, B.E. NMR diffusion measurements to characterise membrane transport and solute binding. Progr. Nucl. Magn. Reson. Spectr. 1997, 30, 39–68. [Google Scholar] [CrossRef]
  167. Larkin, T.J.; Kuchel, P.W. Erythrocyte orientational and cell volume effects on NMR q-space analysis: Simulations of restricted diffusion. Eur. Biophys. 2009, 39, 139–148. [Google Scholar] [CrossRef]
  168. Hong, K.S.; Yeom, G.-S.; Cho, J.-H.; Kim, E.-H.; Lee, C.; Lee, S.D.; Cheong, C. PFG NMR Study of Intra-cellular Drug Uptake in Xenopus laevis Oocyte. J. Korean Magn. Reson. Soc. 2005, 9, 29. [Google Scholar]
  169. Babsky, A.M.; Topper, S.; Zhang, H.; Gao, Y.; James, J.R.; Hekmatyar, S.K.; Bansal, N. Evaluation of extra-and intracellular apparent diffusion coefficient of sodium in rat skeletal muscle: Effects of prolonged ischemia. Magn. Reson. Med. 2008, 59, 485. [Google Scholar] [CrossRef]
  170. Avilova, I.A.; Vasil’ev, S.G.; Rimareva, L.V.; Serba, E.M.; Volkova, L.D.; Volkov, V.I. Water metabolism in cells of Saccharomyces cerevisiae of races Y-3137 and Y-3327, according to pulsed-field gradient NMR data. Russ. J. Phys. Chem. A 2015, 89, 710. [Google Scholar] [CrossRef]
  171. Avilova, I.A.; Smolina, A.V.; Kotelnikov, A.I.; Kotelnikova, R.A.; Loskutov, V.V.; Volkov, V.I. Self-diffusion of water and blood lipids in mouse erythrocytes. Appl. Magn. Reson. 2016, 47, 335–347. [Google Scholar] [CrossRef]
  172. Avilova, I.A.; Khakina, E.A.; Kraevaya, O.A.; Kotelnikov, A.I.; Kotelnikova, R.A.; Troshin, P.A.; Volkov, V.I. Self-diffusion of water-soluble fullerene derivatives in mouse erythrocytes. BBA—Biomembranes 2018, 1860, 1537–1543. [Google Scholar] [CrossRef] [PubMed]
  173. Avilova, I.A.; Soldatova, Y.V.; Kraevaya, O.A.; Zhilenkov, A.V.; Dolgikh, E.A.; Kotel’nikova, R.A.; Troshin, P.A.; Volkov, V.I. Self-Diffusion of Fullerene C60 Derivatives in Aqueous Solutions and Suspensions of Erythrocytes According to Pulsed Field Gradient NMR Data. Russ. J. Phys. Chem. A 2021, 95, 285–291. [Google Scholar] [CrossRef]
  174. Mitra, P.P.; Sen, P.N.; Schawartz, L.M. Short-time behavior of the diffusion coefficient as a geometrical probe of porous media. Phys. Rev. B 1993, 47, 8565. [Google Scholar] [CrossRef] [PubMed]
  175. Anisimov, A.V.; Sorokina, N.Y.; Dautova, N.R. Water diffusion in biological porous systems: A NMR approach. Magn. Reson. Imaging 1998, 16, 565–568. [Google Scholar] [CrossRef]
  176. Crick, F. Diffusion in embryogenesis. Nature 1970, 225, 420–422. [Google Scholar] [CrossRef]
  177. Karger, J.; Pheifer, H.; Heink, W. Principles and Application of Self-diffusion Measurements by Nuclear Magnetic Resonance. Adv. Magn. Reson. 1988, 12, 1–89. [Google Scholar]
  178. Hong, Y.S.; Kim, K.C.; Volkov, V.I.; Skirda, V.D.; Lee, C.-H. Structural and dynamic properties of polyoxyethylene sorbitan monooleate micelle in water dispersion studied by pulsed field gradient NMR. Appl. Magn. Reson. 2005, 29, 351. [Google Scholar] [CrossRef]
  179. Filippov, A.; Oradd, G.; Lindblom, G. The Effect of Cholesterol on the Lateral Diffusion of Phospholipids in Oriented Bilayers. Biophys. J. 2003, 84, 3079. [Google Scholar] [CrossRef] [Green Version]
  180. Beuerle, F.; Lebovitz, R.; Hirsch, A. Antioxidant Properties of Water-Soluble Fullerene Derivatives. In Medicinal Chemistry and Pharmacological Potential of Fullerenes and Carbon Nanotubes. Carbon Materials: Chemistry and Physics; Cataldo, F., Da Ros, T., Eds.; Springer: Dordrecht, The Netherlands, 2008; Volume 1, pp. 51–78. [Google Scholar]
  181. Khakina, E.A.; Kraevaya, O.A.; Popova, M.L.; Peregudov, A.S.; Troyanov, S.I.; Chernyak, A.V.; Martynenko, V.M.; Kulikov, A.V.; Schols, D.; Troshin, P.A. Synthesis of different types of alkoxy fullerene derivatives from chlorofullerene C60Cl6. Org. Biomol. Chem. 2017, 15, 773–777. [Google Scholar] [CrossRef]
  182. Hsieh, F.-Y.; Zhilenkov, A.V.; Voronov, I.I.; Khakina, E.A.; Mischenko, D.V.; Troshin, P.A.; Hsu, S.-H. Water-Soluble Fullerene Derivatives as Brain Medicine: Surface Chemistry Determines If They Are Neuroprotective and Antitumor. ACS Appl. Mater. Interfaces 2017, 9, 11482–11492. [Google Scholar] [CrossRef]
  183. Valiullin, R.; Skirda, V. Time dependent self-diffusion coefficient of molecules in porous media. J. Chem. Phys. 2001, 114, 452–458. [Google Scholar] [CrossRef]
  184. Valiullin, R.R.; Skirda, V.D.; Stapf, S.; Kimmich, R. Molecular exchange processes in partially filled porous glass as seen with NMR diffusometry. R. Phys. Rev. 1996, 55, 2664. [Google Scholar] [CrossRef] [Green Version]
  185. Mutina, A.R.; Skirda, V.D. Porous media characterization by PFG and IMFG NMR. J. Magn. Reson. 2007, 188, 122–128. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic representation of Gierke model. Reprinted with permission from [87]. Copyright 2013 Springer Nature.
Figure 1. Schematic representation of Gierke model. Reprinted with permission from [87]. Copyright 2013 Springer Nature.
Membranes 11 00385 g001
Figure 2. New insight of a Nafion membrane multiscale structure (adapted from [95]).
Figure 2. New insight of a Nafion membrane multiscale structure (adapted from [95]).
Membranes 11 00385 g002
Figure 3. Structure of the amorphous part of a perfluorinated sulfonate cation-exchange membrane [5]. (1) Polymer backbone; (2) hydrated counter-ions and functional groups at a low moisture content; (3) transport channels for ions and water molecules at a high moisture content; L1 = 4 nm according to small-angle X-Ray scattering data [84,85]; L2 = 10 nm according to Mössbauer spectroscopy [93]; l1 = l.2–1 nm according to ENDOR and relaxation NMR data [96,97,98,99]; l3 = 1.5 nm according to standard porosimetry [94] and ENDOR [98,99] methods.
Figure 3. Structure of the amorphous part of a perfluorinated sulfonate cation-exchange membrane [5]. (1) Polymer backbone; (2) hydrated counter-ions and functional groups at a low moisture content; (3) transport channels for ions and water molecules at a high moisture content; L1 = 4 nm according to small-angle X-Ray scattering data [84,85]; L2 = 10 nm according to Mössbauer spectroscopy [93]; l1 = l.2–1 nm according to ENDOR and relaxation NMR data [96,97,98,99]; l3 = 1.5 nm according to standard porosimetry [94] and ENDOR [98,99] methods.
Membranes 11 00385 g003
Figure 4. 1 H NMR spectra in acid ionic form of Nafion 117 membrane at different humidity: (1)RH = 95%; (2) RH = 78%; (3) RH = 64%; (4) RH = 58%; (5) RH = 32%; (6) RH = 10%; (7) RH = 0%. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Figure 4. 1 H NMR spectra in acid ionic form of Nafion 117 membrane at different humidity: (1)RH = 95%; (2) RH = 78%; (3) RH = 64%; (4) RH = 58%; (5) RH = 32%; (6) RH = 10%; (7) RH = 0%. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Membranes 11 00385 g004
Figure 5. 1 H NMR spectra in acid ionic form of Nafion membranes at different temperatures. Membrane samples were equilibrated with water vapor at 58% relative humidity. (1) +25 °C; (2) +10 °C; (3)0 °C; (4) −20 °C; (5) −40 °C; (6) −60 °C. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Figure 5. 1 H NMR spectra in acid ionic form of Nafion membranes at different temperatures. Membrane samples were equilibrated with water vapor at 58% relative humidity. (1) +25 °C; (2) +10 °C; (3)0 °C; (4) −20 °C; (5) −40 °C; (6) −60 °C. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Membranes 11 00385 g005
Figure 6. The example of 1H NMR spectra of Nafion 117 membrane in 95 volume percent CCl4 + 5 volume percent of C6H12: (a) Nafion 117 membrane equilibrated at 75% RH, (b) Nafion 117 membrane equilibrated at 10% RH; 1–signal of C6H12, 2–signal of water and hydrated H+ protons. Figures under signals are the area of NMR lines. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Figure 6. The example of 1H NMR spectra of Nafion 117 membrane in 95 volume percent CCl4 + 5 volume percent of C6H12: (a) Nafion 117 membrane equilibrated at 75% RH, (b) Nafion 117 membrane equilibrated at 10% RH; 1–signal of C6H12, 2–signal of water and hydrated H+ protons. Figures under signals are the area of NMR lines. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Membranes 11 00385 g006
Figure 7. Dependences on humidity of 1H chemical shift (relative to the bulk water signal) in Nafion 117 membranes. λg is the gravimetrically measured amount of water molecules per sulfonate group after membrane drying to constant weight at 110 °C or equilibrated with dry P2O5, h is hydration number of H+(H2O)h); 1—h = 2 (H5O2+) line is Equation 1.2 at δH2O = 4.3 ppm and δc = 13.8 ppm, circles—experiment; 2—h = 1 (H3O+) is experimental curve given for comparison. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Figure 7. Dependences on humidity of 1H chemical shift (relative to the bulk water signal) in Nafion 117 membranes. λg is the gravimetrically measured amount of water molecules per sulfonate group after membrane drying to constant weight at 110 °C or equilibrated with dry P2O5, h is hydration number of H+(H2O)h); 1—h = 2 (H5O2+) line is Equation 1.2 at δH2O = 4.3 ppm and δc = 13.8 ppm, circles—experiment; 2—h = 1 (H3O+) is experimental curve given for comparison. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Membranes 11 00385 g007
Figure 8. Chemical shift (δ) differences between water protons in the membrane and in bulk water as a function of the moisture content of MF-4SC (1–4) and F-4CF (5–8) membranes in various ionic forms. Membrane ionic form: (1, 5) H+, (2, 6) Li+, (3, 7) Na+, (4, 8) Cs+. Reprinted with permission from [20]. Copyright 2011 Elsevier.
Figure 8. Chemical shift (δ) differences between water protons in the membrane and in bulk water as a function of the moisture content of MF-4SC (1–4) and F-4CF (5–8) membranes in various ionic forms. Membrane ionic form: (1, 5) H+, (2, 6) Li+, (3, 7) Na+, (4, 8) Cs+. Reprinted with permission from [20]. Copyright 2011 Elsevier.
Membranes 11 00385 g008
Figure 9. Dependences of chemical shifts of water on the moisture content of MF-4SC membrane. Membrane ionic form: (1) Li+, (2) Na+, (3) K+, (4) Rb+, (5) Cs+, (6) Ba2+, (7) Ca2+, (8) Mg2+, (9) H+ (adapted from [19]).
Figure 9. Dependences of chemical shifts of water on the moisture content of MF-4SC membrane. Membrane ionic form: (1) Li+, (2) Na+, (3) K+, (4) Rb+, (5) Cs+, (6) Ba2+, (7) Ca2+, (8) Mg2+, (9) H+ (adapted from [19]).
Membranes 11 00385 g009
Figure 10. Illustration of hydration of functional groups in the salt forms of MF-4SC membranes by the example of its Li+-form. Reprinted with permission from [39]. Copyright 2002 Springer Nature.
Figure 10. Illustration of hydration of functional groups in the salt forms of MF-4SC membranes by the example of its Li+-form. Reprinted with permission from [39]. Copyright 2002 Springer Nature.
Membranes 11 00385 g010
Figure 11. Temperature dependences of 1 H NMR chemical shifts in acid ionic form of Nafion membrane at different water content (δH2O is bulk water chemical shift at 25 °C). (1) bulk H2O; (2) RH 98% (λ = 17.5 ± 0.4); (3) RH 95% (λ = 12 ± 0.4); (4) RH 78% (λ = 7.4 ± 0.4); (5) RH 64% (λ = 6.4 ± 0.4); (6) RH 58% (λ = 5.8 ± 0.4); (7) RH 32% (λ = 4.4 ± 0.4); (8) RH 10% (λ = 3.2 ± 0.4); (9) RH 0% (λ = 1.9 ± 0.4, drying on P2O5 or at 110 °C). Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Figure 11. Temperature dependences of 1 H NMR chemical shifts in acid ionic form of Nafion membrane at different water content (δH2O is bulk water chemical shift at 25 °C). (1) bulk H2O; (2) RH 98% (λ = 17.5 ± 0.4); (3) RH 95% (λ = 12 ± 0.4); (4) RH 78% (λ = 7.4 ± 0.4); (5) RH 64% (λ = 6.4 ± 0.4); (6) RH 58% (λ = 5.8 ± 0.4); (7) RH 32% (λ = 4.4 ± 0.4); (8) RH 10% (λ = 3.2 ± 0.4); (9) RH 0% (λ = 1.9 ± 0.4, drying on P2O5 or at 110 °C). Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Membranes 11 00385 g011
Figure 12. Temperature dependence of water proton chemical shifts in the Li+, Na+, and Cs+ ionic forms of Nafion membranes at various relative humidity where δH2O is the bulk water chemical shift δH2O = 4.30 ppm relatively TMS, at 20 °C: Li+ ionic form: (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7, (7) λ = 12; (8) Na+ ionic form, λ = 10; (9) Cs+ ionic form, λ = 4; (10) bulk water; λ is amount of water molecules per sulfonated group. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Figure 12. Temperature dependence of water proton chemical shifts in the Li+, Na+, and Cs+ ionic forms of Nafion membranes at various relative humidity where δH2O is the bulk water chemical shift δH2O = 4.30 ppm relatively TMS, at 20 °C: Li+ ionic form: (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7, (7) λ = 12; (8) Na+ ionic form, λ = 10; (9) Cs+ ionic form, λ = 4; (10) bulk water; λ is amount of water molecules per sulfonated group. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Membranes 11 00385 g012
Figure 13. NMR spectra of 7 Li (a), 2 3Na (b), and 133 Cs (c) nuclei in appropriate ionic form of MSC membrane at RH = 95% [17].
Figure 13. NMR spectra of 7 Li (a), 2 3Na (b), and 133 Cs (c) nuclei in appropriate ionic form of MSC membrane at RH = 95% [17].
Membranes 11 00385 g013
Figure 14. Dependences of chemical shifts in 7 Li NMR spectra on the moisture content of the Li+ forms of perfluorinated sulfonate cation-exchange membranes MF-4SC (1), Nafion 117 (2), and carboxylic membranes F-4CF (3). Reprinted with permission from [21]. Copyright 2010 Springer Nature.
Figure 14. Dependences of chemical shifts in 7 Li NMR spectra on the moisture content of the Li+ forms of perfluorinated sulfonate cation-exchange membranes MF-4SC (1), Nafion 117 (2), and carboxylic membranes F-4CF (3). Reprinted with permission from [21]. Copyright 2010 Springer Nature.
Membranes 11 00385 g014
Figure 15. Chemical shifts in 23 Na NMR (a) and 133 Cs NMR (b) spectra vs. the moisture content of perfluorinated sulfonate cation-exchange membranes MF-4SC (1) and carboxylic membranes F-CF (2). Reprinted with permission from [21]. Copyright 2010 Springer Nature.
Figure 15. Chemical shifts in 23 Na NMR (a) and 133 Cs NMR (b) spectra vs. the moisture content of perfluorinated sulfonate cation-exchange membranes MF-4SC (1) and carboxylic membranes F-CF (2). Reprinted with permission from [21]. Copyright 2010 Springer Nature.
Membranes 11 00385 g015aMembranes 11 00385 g015b
Figure 16. Relative parts of contact ion pairs vs. the moisture content for sulfonate cation-exchange (1–3) and carboxylic (4–6) membranes in various ionic forms. Membrane ionic form: (1, 4) Li+; (2, 5) Na+; (3, 6) Cs+. Reprinted with permission from [20]. Copyright 2011 Elsevier.
Figure 16. Relative parts of contact ion pairs vs. the moisture content for sulfonate cation-exchange (1–3) and carboxylic (4–6) membranes in various ionic forms. Membrane ionic form: (1, 4) Li+; (2, 5) Na+; (3, 6) Cs+. Reprinted with permission from [20]. Copyright 2011 Elsevier.
Membranes 11 00385 g016
Figure 17. Dependences of mobile water molecule amount in acidic form of Nafion 117 membrane on temperature at different water content: (1) λ = 3.2 ± 0.4; (2) λ = 4.4 ± 0.4; (3) λ = 5.8 ± 0.4; (4) λ = 6.4 ± 0.4; (5) λ = 7.4 ± 0.4; (6) λ = 12 ± 0.4. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Figure 17. Dependences of mobile water molecule amount in acidic form of Nafion 117 membrane on temperature at different water content: (1) λ = 3.2 ± 0.4; (2) λ = 4.4 ± 0.4; (3) λ = 5.8 ± 0.4; (4) λ = 6.4 ± 0.4; (5) λ = 7.4 ± 0.4; (6) λ = 12 ± 0.4. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Membranes 11 00385 g017
Figure 18. DSC thermograms in H+ ionic forms of perfluorinated sulfo cation exchange MF-4SC membrane, which is complete analog of Nafion. Dry–sample was dried to constant weight at 120 °C (λ is about 2 water molecules per sulfonate group), air dry–sample was dried to constant weight at room temperature (λ is about 5 water molecules per sulfonate group). Reprinted with permission from [38]. Copyright 2003 Springer Nature.
Figure 18. DSC thermograms in H+ ionic forms of perfluorinated sulfo cation exchange MF-4SC membrane, which is complete analog of Nafion. Dry–sample was dried to constant weight at 120 °C (λ is about 2 water molecules per sulfonate group), air dry–sample was dried to constant weight at room temperature (λ is about 5 water molecules per sulfonate group). Reprinted with permission from [38]. Copyright 2003 Springer Nature.
Membranes 11 00385 g018
Figure 19. Dependence of the amount of mobile water in the H+ form of MSC membranes at different moisture contents in the temperature range from 25 to −40 °C. The number of water molecules per sulfo group is λ = 5.1 (1), 7 (2), 8.4 (3), 12.5 (4). Reprinted with permission from [16]. Copyright 2020 Springer Nature.
Figure 19. Dependence of the amount of mobile water in the H+ form of MSC membranes at different moisture contents in the temperature range from 25 to −40 °C. The number of water molecules per sulfo group is λ = 5.1 (1), 7 (2), 8.4 (3), 12.5 (4). Reprinted with permission from [16]. Copyright 2020 Springer Nature.
Membranes 11 00385 g019
Figure 20. The dependences of mobile water molecules per sulfonated group λ (a) and mobile Li+ cations on temperature in lithium ionic form of Nafion 117 membrane (b); (a) the amount of mobile water molecules per sulfonate group: (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7. (b) a relative intensity of 7 Li+ signals: (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7.
Figure 20. The dependences of mobile water molecules per sulfonated group λ (a) and mobile Li+ cations on temperature in lithium ionic form of Nafion 117 membrane (b); (a) the amount of mobile water molecules per sulfonate group: (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7. (b) a relative intensity of 7 Li+ signals: (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7.
Membranes 11 00385 g020
Figure 21. Temperature dependences of 1 H spin–lattice R1 (T1−1) relaxation rates in Nafion 117 membrane with different water content (1) λ = 3.2; (2) λ = 4.4; (3) λ = 5.8; (4) λ = 7.4. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Figure 21. Temperature dependences of 1 H spin–lattice R1 (T1−1) relaxation rates in Nafion 117 membrane with different water content (1) λ = 3.2; (2) λ = 4.4; (3) λ = 5.8; (4) λ = 7.4. Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Membranes 11 00385 g021
Figure 22. Temperature-dependent longitudinal relaxation rates (R1) at apparent hydration levels of (a) λwater = 2, (b) λwater = 3, (c), λwater = 5, (d) λwater = 7, and (e) λwater = 8 for Nafion N117. Black squares represent the experimental data and the green and orange lines represent the BPP fit of the slow (mode I) and fast (mode II) motional modes, respectively; black lines show the superpositions of both contributions. The purple dashed line in (a) represents the best fit assuming a single BPP model Figure 4 from. Reprinted with permission from [26]. Copyright 2019 American Chemical Society.
Figure 22. Temperature-dependent longitudinal relaxation rates (R1) at apparent hydration levels of (a) λwater = 2, (b) λwater = 3, (c), λwater = 5, (d) λwater = 7, and (e) λwater = 8 for Nafion N117. Black squares represent the experimental data and the green and orange lines represent the BPP fit of the slow (mode I) and fast (mode II) motional modes, respectively; black lines show the superpositions of both contributions. The purple dashed line in (a) represents the best fit assuming a single BPP model Figure 4 from. Reprinted with permission from [26]. Copyright 2019 American Chemical Society.
Membranes 11 00385 g022
Figure 23. Temperature dependences of spin–lattice (R1—curve 2) and spin–spin (R2—curve 1) relaxation rates of 1 H water molecule nuclei in Li+ ionic form of MF-4SC membrane at water content λ = 20.5 water molecules per sulfonate group (adapted from [37,38]).
Figure 23. Temperature dependences of spin–lattice (R1—curve 2) and spin–spin (R2—curve 1) relaxation rates of 1 H water molecule nuclei in Li+ ionic form of MF-4SC membrane at water content λ = 20.5 water molecules per sulfonate group (adapted from [37,38]).
Membranes 11 00385 g023
Figure 24. Temperature dependences of spin–lattice (R1—curve 2) and spin–spin (R2—curve 1) relaxation rates of 7 Li lithium cation nuclei in Li+ ionic form of MF-4SC membrane at water content λ = 20.5 water molecules per sulfonate group (adapted from [31,38]).
Figure 24. Temperature dependences of spin–lattice (R1—curve 2) and spin–spin (R2—curve 1) relaxation rates of 7 Li lithium cation nuclei in Li+ ionic form of MF-4SC membrane at water content λ = 20.5 water molecules per sulfonate group (adapted from [31,38]).
Membranes 11 00385 g024
Figure 25. Temperature dependences of population pl (dashed line) and ph (solid line) for Li+ ionic form of MF-4SC membrane at λ = 20.5 (adapted from [37,38]).
Figure 25. Temperature dependences of population pl (dashed line) and ph (solid line) for Li+ ionic form of MF-4SC membrane at λ = 20.5 (adapted from [37,38]).
Membranes 11 00385 g025
Figure 26. Dependences of water self-diffusion coefficient and proton conductivity on humidity (λ), where λ is water molecule amount per sulfonate group in acid form of MF-4SC membrane [57].
Figure 26. Dependences of water self-diffusion coefficient and proton conductivity on humidity (λ), where λ is water molecule amount per sulfonate group in acid form of MF-4SC membrane [57].
Membranes 11 00385 g026
Figure 27. (a) Diffusion coefficients of water molecules dependence on the moisture content for sulfonate cation-exchange MF-4SC (1–4) and carboxylic F-4CF (5–8) perfluorinated membranes in various ionic forms. (b) Diffusion coefficients dependence of water molecules and hydrated H+ counterions in the acidic form of F-4CF membrane on the diffusion time, relative humidity is 95% (b). Membrane ionic forms: (1, 5) H+, (2, 6) Li+, (3, 7) Na+, (4, 8) Cs+. Reprinted with permission from [20]. Copyright 2011 Elsevier.
Figure 27. (a) Diffusion coefficients of water molecules dependence on the moisture content for sulfonate cation-exchange MF-4SC (1–4) and carboxylic F-4CF (5–8) perfluorinated membranes in various ionic forms. (b) Diffusion coefficients dependence of water molecules and hydrated H+ counterions in the acidic form of F-4CF membrane on the diffusion time, relative humidity is 95% (b). Membrane ionic forms: (1, 5) H+, (2, 6) Li+, (3, 7) Na+, (4, 8) Cs+. Reprinted with permission from [20]. Copyright 2011 Elsevier.
Membranes 11 00385 g027aMembranes 11 00385 g027b
Figure 28. Conductivity dependence on the moisture content for sulfonate cation-exchange MF-4SC (1–4) and carboxylic F-4CF (5–8) perfluorinated membranes in various ionic forms. Membrane ionic forms: (1, 5) H+, (2, 6) Li+, (3, 7) Na+, (4, 8) Cs+. Reprinted with permission from [20]. Copyright 2011 Elsevier.
Figure 28. Conductivity dependence on the moisture content for sulfonate cation-exchange MF-4SC (1–4) and carboxylic F-4CF (5–8) perfluorinated membranes in various ionic forms. Membrane ionic forms: (1, 5) H+, (2, 6) Li+, (3, 7) Na+, (4, 8) Cs+. Reprinted with permission from [20]. Copyright 2011 Elsevier.
Membranes 11 00385 g028
Figure 29. (a) The temperature dependences of average H+ and water-self diffusion coefficients for Nafion 117 membrane with different water content: curve (1) λ is 12; curve (2) λ is 7.4; curve (3) λ is 6.4; curve (4) λ is 5.8; curve (5) λ is 4.4; and curve (6) λ is 3.2; (b)water self-diffusion activation energies in high temperature region (squares—curve 1) and low temperature region (circles—curve 2). Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Figure 29. (a) The temperature dependences of average H+ and water-self diffusion coefficients for Nafion 117 membrane with different water content: curve (1) λ is 12; curve (2) λ is 7.4; curve (3) λ is 6.4; curve (4) λ is 5.8; curve (5) λ is 4.4; and curve (6) λ is 3.2; (b)water self-diffusion activation energies in high temperature region (squares—curve 1) and low temperature region (circles—curve 2). Reprinted with permission from [25]. Copyright 2019 Springer Nature.
Membranes 11 00385 g029
Figure 30. Temperature dependences of self-diffusion coefficients—curve 1 and calculated from Equation (7) protonic conductivity—curve 2, λ = 4.4 [25]. Curve 3—experimental temperature dependence of conductivity for λg = 3.7 [112].
Figure 30. Temperature dependences of self-diffusion coefficients—curve 1 and calculated from Equation (7) protonic conductivity—curve 2, λ = 4.4 [25]. Curve 3—experimental temperature dependence of conductivity for λg = 3.7 [112].
Membranes 11 00385 g030
Figure 31. Temperature dependence of the diffusion coefficients D measured by pulsed field gradient NMR at different water contents in acid form of MSC membrane: λ = 5.1 (1), 7 (2), 8.4 (3). Reprinted with permission from [16]. Copyright 2020 Springer Nature.
Figure 31. Temperature dependence of the diffusion coefficients D measured by pulsed field gradient NMR at different water contents in acid form of MSC membrane: λ = 5.1 (1), 7 (2), 8.4 (3). Reprinted with permission from [16]. Copyright 2020 Springer Nature.
Membranes 11 00385 g031
Figure 32. Diffusion coefficients of H+ cations calculated from the proton conductivity data according to Equation 1.7 (lower curve) and the average diffusion coefficients of water molecules and H+ cations measured by pulsed field gradient NMR (upper curve) in the H+ form of the MSC membrane at different moisture contents. RH = 100%, the membrane is in contact with water. Reprinted with permission from [16]. Copyright 2020 Springer Nature.
Figure 32. Diffusion coefficients of H+ cations calculated from the proton conductivity data according to Equation 1.7 (lower curve) and the average diffusion coefficients of water molecules and H+ cations measured by pulsed field gradient NMR (upper curve) in the H+ form of the MSC membrane at different moisture contents. RH = 100%, the membrane is in contact with water. Reprinted with permission from [16]. Copyright 2020 Springer Nature.
Membranes 11 00385 g032
Figure 33. (a) Spin–echo signal attenuation of water molecule 1H nuclei dependences on gradient pulsed amplitude A(g) (diffusion decay) in Li+ ionic form of Nafion membrane at different water content: (1) λ = 2.0, (2) λ = 4.0, (3) λ = 5.7, (4) λ = 7.4, (5) λ = 10.7. (b) Spin-echo signal attenuation of 7 Li nuclei dependences on gradient pulsed amplitude A(g) (diffusion decay) in Li+ ionic form of Nafion membrane at different water content: (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Figure 33. (a) Spin–echo signal attenuation of water molecule 1H nuclei dependences on gradient pulsed amplitude A(g) (diffusion decay) in Li+ ionic form of Nafion membrane at different water content: (1) λ = 2.0, (2) λ = 4.0, (3) λ = 5.7, (4) λ = 7.4, (5) λ = 10.7. (b) Spin-echo signal attenuation of 7 Li nuclei dependences on gradient pulsed amplitude A(g) (diffusion decay) in Li+ ionic form of Nafion membrane at different water content: (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Membranes 11 00385 g033
Figure 34. Self-diffusion coefficients of Li+ ions (curve 1) and water molecules (curve 2) dependences on water content λ in Nafion 117 membrane Li+ ionic form, where λ is water amount per sulfonate group. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Figure 34. Self-diffusion coefficients of Li+ ions (curve 1) and water molecules (curve 2) dependences on water content λ in Nafion 117 membrane Li+ ionic form, where λ is water amount per sulfonate group. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Membranes 11 00385 g034
Figure 35. Li+ cation self-diffusion coefficient temperature dependences for Nafion membrane Li+ ionic form at different water content λ, where λ is water amount per sulfonate group. (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Figure 35. Li+ cation self-diffusion coefficient temperature dependences for Nafion membrane Li+ ionic form at different water content λ, where λ is water amount per sulfonate group. (1) λ = 0.9, (2) λ = 2.0, (3) λ = 4.0, (4) λ = 5.7, (5) λ = 7.4, (6) λ = 10.7. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Membranes 11 00385 g035
Figure 36. Water molecule self-diffusion coefficient temperature dependences for Nafion membrane Li+ ionic form at different water content λ, where λ is water amount per sulfonate group. (1) λ = 2.0, (2) λ = 4.0, (3) λ = 5.7, (4) λ = 7.4, (5) λ = 10.7. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Figure 36. Water molecule self-diffusion coefficient temperature dependences for Nafion membrane Li+ ionic form at different water content λ, where λ is water amount per sulfonate group. (1) λ = 2.0, (2) λ = 4.0, (3) λ = 5.7, (4) λ = 7.4, (5) λ = 10.7. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Membranes 11 00385 g036
Figure 37. Diffusion decay of 23 Na nuclei in Nafion membrane Na+ ionic form, RH = 95%, t = 20 °C. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Figure 37. Diffusion decay of 23 Na nuclei in Nafion membrane Na+ ionic form, RH = 95%, t = 20 °C. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Membranes 11 00385 g037
Figure 38. Diffusion decay of 133 Cs nuclei in Nafion membrane Cs+ ionic form, RH = 95%, t = 20 °C. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Figure 38. Diffusion decay of 133 Cs nuclei in Nafion membrane Cs+ ionic form, RH = 95%, t = 20 °C. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Membranes 11 00385 g038
Figure 39. Temperature dependences of Li+ (curve 3), Na+ (curve 2), Cs+ (curve 1) self-diffusion coefficients in Li+, Na+, Cs+ Nafion 117 membrane ionic forms. RH = 98%. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Figure 39. Temperature dependences of Li+ (curve 3), Na+ (curve 2), Cs+ (curve 1) self-diffusion coefficients in Li+, Na+, Cs+ Nafion 117 membrane ionic forms. RH = 98%. Reprinted with permission from [109]. Copyright 2021 Elsevier.
Membranes 11 00385 g039
Figure 40. Diffusion decays of 7 Li (a), 23 Na (b), 133 Cs (c) nuclei NMR signals in appropriate ionic form of MSC membrane at RH = 95% and different temperatures (1) 20 °C, (2) 30 °C, (3) 40 °C, (4) 50 °C, (5) 60 °C, (6) 70 °C, (7) 80 °C [17].
Figure 40. Diffusion decays of 7 Li (a), 23 Na (b), 133 Cs (c) nuclei NMR signals in appropriate ionic form of MSC membrane at RH = 95% and different temperatures (1) 20 °C, (2) 30 °C, (3) 40 °C, (4) 50 °C, (5) 60 °C, (6) 70 °C, (7) 80 °C [17].
Membranes 11 00385 g040aMembranes 11 00385 g040b
Figure 41. Temperature dependences of Cs+, Na+, Li+ diffusion coefficients in appropriate ionic form of MSC membrane at RH = 95%: (1) Cs+ ionic form, Ea = 18.1 kJ/mol; (2) Na+ ionic form, Ea = 16.5 kJ/mol; (3) Li+ ionic form, Ea = 17.6 kJ/mol [17].
Figure 41. Temperature dependences of Cs+, Na+, Li+ diffusion coefficients in appropriate ionic form of MSC membrane at RH = 95%: (1) Cs+ ionic form, Ea = 18.1 kJ/mol; (2) Na+ ionic form, Ea = 16.5 kJ/mol; (3) Li+ ionic form, Ea = 17.6 kJ/mol [17].
Membranes 11 00385 g041
Figure 42. Temperature dependences of experimental σexp (1–4) and calculated σcalc (2′–4′) ionic conductivities in H+ (1), Li+ (2) and (2′), Na+ (3) and (3′), Cs+ (4) and (4′) ionic forms of MSC membrane at RH 95%.
Figure 42. Temperature dependences of experimental σexp (1–4) and calculated σcalc (2′–4′) ionic conductivities in H+ (1), Li+ (2) and (2′), Na+ (3) and (3′), Cs+ (4) and (4′) ionic forms of MSC membrane at RH 95%.
Membranes 11 00385 g042
Figure 43. 1 H chemical shift dependences on concentration of LiCl (1), NaCl (2), CsCl (3) aqueous solutions [17].
Figure 43. 1 H chemical shift dependences on concentration of LiCl (1), NaCl (2), CsCl (3) aqueous solutions [17].
Membranes 11 00385 g043
Figure 44. Water molecule and cation diffusion coefficient concentration dependences in lithium, sodium, cesium chloride aqueous solutions. (1) H2O in LiCl, (2) H2O in NaCl, (3) H2O in CsCl. (1′) Li+ in LiCl, (2′) Na+ in NaCl, (3′) Cs+ in CsCl [17].
Figure 44. Water molecule and cation diffusion coefficient concentration dependences in lithium, sodium, cesium chloride aqueous solutions. (1) H2O in LiCl, (2) H2O in NaCl, (3) H2O in CsCl. (1′) Li+ in LiCl, (2′) Na+ in NaCl, (3′) Cs+ in CsCl [17].
Membranes 11 00385 g044
Figure 45. 7 Li (1), 23 Na (2), 133 Cs (3) nuclear NMR chemical shift concentration dependences in lithium, sodium, cesium chloride aqueous solutions [17].
Figure 45. 7 Li (1), 23 Na (2), 133 Cs (3) nuclear NMR chemical shift concentration dependences in lithium, sodium, cesium chloride aqueous solutions [17].
Membranes 11 00385 g045
Figure 46. Self-diffusion coefficient dependences on humidity of water molecules (1) and lithium cations (2). Reprinted with permission from [39]. Copyright 2002 Springer Nature.
Figure 46. Self-diffusion coefficient dependences on humidity of water molecules (1) and lithium cations (2). Reprinted with permission from [39]. Copyright 2002 Springer Nature.
Membranes 11 00385 g046
Figure 47. Diffusion coefficients of water on the moisture content in the acidic form of membranes MF-4SC with the exchange capacity of 0.86 (1) and 0.34 mg-equiv.∙g−1 (2) [39,105]. Reprinted with permission from [105]. Copyright 2013 Turpion Limited.
Figure 47. Diffusion coefficients of water on the moisture content in the acidic form of membranes MF-4SC with the exchange capacity of 0.86 (1) and 0.34 mg-equiv.∙g−1 (2) [39,105]. Reprinted with permission from [105]. Copyright 2013 Turpion Limited.
Membranes 11 00385 g047
Figure 48. Chemical structure of sulfocontaining aromatic polyamides PA-1 and PA-2.
Figure 48. Chemical structure of sulfocontaining aromatic polyamides PA-1 and PA-2.
Membranes 11 00385 g048
Figure 49. Diffusion coefficients of lithium cations in the system lithium salt-disulfophthalic acid-water ((1) is terephthalic, (2) is isophthalic) as a function of the moisture content in macroporous sulfonate cation-exchange CU-23 (4) and perfluorinated sulfonate cation-exchange membrane MF-4SC (3) [59]. Exchange capacity/mg-equiv∙g−1: 0.86 for MF-4SC, 5 for macroporous cation exchanger CU-23, 2 for aromatic disulfo-containing polyamides. Reprinted with permission from [59]. Copyright 2010 Springer Nature.
Figure 49. Diffusion coefficients of lithium cations in the system lithium salt-disulfophthalic acid-water ((1) is terephthalic, (2) is isophthalic) as a function of the moisture content in macroporous sulfonate cation-exchange CU-23 (4) and perfluorinated sulfonate cation-exchange membrane MF-4SC (3) [59]. Exchange capacity/mg-equiv∙g−1: 0.86 for MF-4SC, 5 for macroporous cation exchanger CU-23, 2 for aromatic disulfo-containing polyamides. Reprinted with permission from [59]. Copyright 2010 Springer Nature.
Membranes 11 00385 g049
Figure 50. Effect of the diffusant content on the diffusion coefficients of water (1), methanol (2), ethanol (3), propanol (4) in the Li+-form of perfluorinated membrane MF-4SC m is the number of diffusant molecules per sulfonate group (adapted from [47]).
Figure 50. Effect of the diffusant content on the diffusion coefficients of water (1), methanol (2), ethanol (3), propanol (4) in the Li+-form of perfluorinated membrane MF-4SC m is the number of diffusant molecules per sulfonate group (adapted from [47]).
Membranes 11 00385 g050
Figure 51. Diffusion coefficients of methanol, ethanol, propanol, and butanol in bulk liquids (1) and MF-4SC membranes in various ionic forms. Membrane ionic form: (2) Li+, (3) H+, (4) Na+, (5) Cs+, (6) Rb+, (7) K+ (adapted from [47]).
Figure 51. Diffusion coefficients of methanol, ethanol, propanol, and butanol in bulk liquids (1) and MF-4SC membranes in various ionic forms. Membrane ionic form: (2) Li+, (3) H+, (4) Na+, (5) Cs+, (6) Rb+, (7) K+ (adapted from [47]).
Membranes 11 00385 g051
Figure 52. Water self-diffusion coefficients dependences on concentration (water content) in (1) bovine serum albumin, (2) gelatin, and (3) perfluorinated MS-4SC membrane in acidic form (adapted from [55]).
Figure 52. Water self-diffusion coefficients dependences on concentration (water content) in (1) bovine serum albumin, (2) gelatin, and (3) perfluorinated MS-4SC membrane in acidic form (adapted from [55]).
Membranes 11 00385 g052
Figure 53. Generalized concentration dependence on concentration of the self-diffusion coefficients of globular proteins molecules [130]. The figure shows the experimental values for BSA (pH = 4.8–5.2), myoglobin (pH = 6.8–7.2), lysozyme (pH = 2.9–3.0 and pH = 7.4–7.8), barstar (pH = 8.0–8.2). The solid line is the generalized self-diffusion coefficients concentration dependence of the flexible-chain polymers [132]. φ ^ –critical concentration, found as the intersection of the asymptote with zero slope ϕ−3, drawn to the experimental concentration dependences of the self-diffusion coefficients of proteins.
Figure 53. Generalized concentration dependence on concentration of the self-diffusion coefficients of globular proteins molecules [130]. The figure shows the experimental values for BSA (pH = 4.8–5.2), myoglobin (pH = 6.8–7.2), lysozyme (pH = 2.9–3.0 and pH = 7.4–7.8), barstar (pH = 8.0–8.2). The solid line is the generalized self-diffusion coefficients concentration dependence of the flexible-chain polymers [132]. φ ^ –critical concentration, found as the intersection of the asymptote with zero slope ϕ−3, drawn to the experimental concentration dependences of the self-diffusion coefficients of proteins.
Membranes 11 00385 g053
Figure 54. Concentration dependence of the diffusion coefficients of α-casein [139], globular proteins [130], and linear flexible polymers [132]. The master curves for globular proteins (blue experimental points) and flexible linear polymers (red solid line) are shown. The asymptotes with the slopes of ϕ−0, ϕ−3, and ϕ−12 are indicated. Solid magenta diamonds represent the concentration dependence of the diffusion coefficient of α-casein < D >/D0 normalized by φ ^ .
Figure 54. Concentration dependence of the diffusion coefficients of α-casein [139], globular proteins [130], and linear flexible polymers [132]. The master curves for globular proteins (blue experimental points) and flexible linear polymers (red solid line) are shown. The asymptotes with the slopes of ϕ−0, ϕ−3, and ϕ−12 are indicated. Solid magenta diamonds represent the concentration dependence of the diffusion coefficient of α-casein < D >/D0 normalized by φ ^ .
Membranes 11 00385 g054
Figure 55. Dependence of the diffusion attenuation of spin–echo signal in solutions of α-casein on diffusion time [139]. Same diffusion attenuations are plotted using different coordinates: [43,140] (a) as a function of (γδg) 2td and (b) as a funcTable 2. Diffusion attenuations are shown for 15% α-casein solution. Curves 1–3 correspond to diffusion attenuations collected at diffusion times 200, 400, and 800 ms, respectively. Curve 4 is a control experiment. The diffusion attenuation shown by curve 4 was collected at 200 ms after the completion of the experiments carried out at different values of td.
Figure 55. Dependence of the diffusion attenuation of spin–echo signal in solutions of α-casein on diffusion time [139]. Same diffusion attenuations are plotted using different coordinates: [43,140] (a) as a function of (γδg) 2td and (b) as a funcTable 2. Diffusion attenuations are shown for 15% α-casein solution. Curves 1–3 correspond to diffusion attenuations collected at diffusion times 200, 400, and 800 ms, respectively. Curve 4 is a control experiment. The diffusion attenuation shown by curve 4 was collected at 200 ms after the completion of the experiments carried out at different values of td.
Membranes 11 00385 g055
Figure 56. Schematic structure of the aquaporin channel. The fixed positive charge on the side chains of Arg-195 (R195) and His-180 (H180) create a barrier to cations, including the protonated water (H3O+). The hydrogen bonds of the water chain adjacent to Arg-195 are broken, thereby preventing the formation of proton conductivity [150]. Reprinted with permission from [150]. Copyright 2002 American Society for Clinical Investigation.
Figure 56. Schematic structure of the aquaporin channel. The fixed positive charge on the side chains of Arg-195 (R195) and His-180 (H180) create a barrier to cations, including the protonated water (H3O+). The hydrogen bonds of the water chain adjacent to Arg-195 are broken, thereby preventing the formation of proton conductivity [150]. Reprinted with permission from [150]. Copyright 2002 American Society for Clinical Investigation.
Membranes 11 00385 g056
Figure 57. Values of the permeability P of RBCs of humans and some laboratory animals and pets for water at 25 and 37 °C. Reprinted with permission from [156]. Copyright 2012 Springer Nature.
Figure 57. Values of the permeability P of RBCs of humans and some laboratory animals and pets for water at 25 and 37 °C. Reprinted with permission from [156]. Copyright 2012 Springer Nature.
Membranes 11 00385 g057
Figure 58. Inhibition of the permeability of RBCs in humans and some laboratory animals and pets for water by PCMBS. Reprinted with permission from [156]. Copyright 2012 Springer Nature.
Figure 58. Inhibition of the permeability of RBCs in humans and some laboratory animals and pets for water by PCMBS. Reprinted with permission from [156]. Copyright 2012 Springer Nature.
Membranes 11 00385 g058
Figure 59. The RBCs suspension 1 H spectra with suppression of water signal—pulsed field gradient amplitude g = 0.375 T/m (a) and g = 10.5 T/m (b). 0–1.70 ppm: (CH2)n and CH3 groups of lipids; 6–9 ppm: protein component (protons of aromatic rings and NH and OH groups. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Figure 59. The RBCs suspension 1 H spectra with suppression of water signal—pulsed field gradient amplitude g = 0.375 T/m (a) and g = 10.5 T/m (b). 0–1.70 ppm: (CH2)n and CH3 groups of lipids; 6–9 ppm: protein component (protons of aromatic rings and NH and OH groups. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Membranes 11 00385 g059
Figure 60. (a) The diffusion decays of water molecules in RBCs measured at various td (1—10 ms, 2—20 ms, 3—50 ms, 4—100 ms, 5—250 ms, 6—500 ms) and t = 35 °C; (b) the procedure of decomposition of the original diffusion decay on three exponential components according to Equation A1 is shown: Ds3 = 0.60·10−9 m2/s, p3′ = 0.947; Ds2 = 0.76·10−10 m2/s, p2′ = 0.045; m2/s; Ds1 = 1.37·10−11 m2/s, p1′ = 0.008 (1–original diffusion decay, 2–1st residual, 3–2nd residual). Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Figure 60. (a) The diffusion decays of water molecules in RBCs measured at various td (1—10 ms, 2—20 ms, 3—50 ms, 4—100 ms, 5—250 ms, 6—500 ms) and t = 35 °C; (b) the procedure of decomposition of the original diffusion decay on three exponential components according to Equation A1 is shown: Ds3 = 0.60·10−9 m2/s, p3′ = 0.947; Ds2 = 0.76·10−10 m2/s, p2′ = 0.045; m2/s; Ds1 = 1.37·10−11 m2/s, p1′ = 0.008 (1–original diffusion decay, 2–1st residual, 3–2nd residual). Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Membranes 11 00385 g060
Figure 61. The water self-diffusion coefficient Ds1, Ds2, and Ds3 dependences on the diffusion time td at 35 °C (a). The population p1, p2, and p3 dependences on the diffusion time td at 35 °C (b) Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Figure 61. The water self-diffusion coefficient Ds1, Ds2, and Ds3 dependences on the diffusion time td at 35 °C (a). The population p1, p2, and p3 dependences on the diffusion time td at 35 °C (b) Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Membranes 11 00385 g061aMembranes 11 00385 g061b
Figure 62. The experimental dependence of self-diffusion coefficients Ds(td) (a) T = 35 (1), 20 (2), and 5 (3) °C and effective self-diffusion coefficients Dseff(td) calculated from Equation (A5) (b) T = 35 (1), 25 (2), and 10 (3) °C on the diffusion time td for intracellular water molecules in the RBCs. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Figure 62. The experimental dependence of self-diffusion coefficients Ds(td) (a) T = 35 (1), 20 (2), and 5 (3) °C and effective self-diffusion coefficients Dseff(td) calculated from Equation (A5) (b) T = 35 (1), 25 (2), and 10 (3) °C on the diffusion time td for intracellular water molecules in the RBCs. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Membranes 11 00385 g062aMembranes 11 00385 g062b
Figure 63. A temperature dependence of the hindered water self-diffusion coefficient Dp. The solid straight line is the Arrhenius equation approximation. The activation energy is 24.1 ± 1.9 kJ/mol. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Figure 63. A temperature dependence of the hindered water self-diffusion coefficient Dp. The solid straight line is the Arrhenius equation approximation. The activation energy is 24.1 ± 1.9 kJ/mol. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Membranes 11 00385 g063
Figure 64. (1) the dependence of intracellular water population p1(td) on diffusion time td. (2) The dependence p1(td) after subtraction of spin–lattice relaxation part. Temperature is 35 °C. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Figure 64. (1) the dependence of intracellular water population p1(td) on diffusion time td. (2) The dependence p1(td) after subtraction of spin–lattice relaxation part. Temperature is 35 °C. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Membranes 11 00385 g064
Figure 65. Diffusion decays of the 1 H spin echo–signal from different regions of integration: (1) 0–3.7 ppm, (2) 0–3 ppm, (3) 0–2 ppm. Diffusion time td = 50 ms. Temperature is 35 °C.
Figure 65. Diffusion decays of the 1 H spin echo–signal from different regions of integration: (1) 0–3.7 ppm, (2) 0–3 ppm, (3) 0–2 ppm. Diffusion time td = 50 ms. Temperature is 35 °C.
Membranes 11 00385 g065
Figure 66. The temperature dependence of lipid self-diffusion coefficients DL. The straight line is the Arrhenius equation approximation. The activation energy is 25 ± 2.9 kJ/mol. The diffusion time td is 10 ms. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Figure 66. The temperature dependence of lipid self-diffusion coefficients DL. The straight line is the Arrhenius equation approximation. The activation energy is 25 ± 2.9 kJ/mol. The diffusion time td is 10 ms. Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Membranes 11 00385 g066
Figure 67. The dependence of the blood lipid self-diffusion coefficient DL(td) on diffusion time td at 20 °C (1) and 25 °C (2) (a). The dependence of effective blood lipid self-diffusion coefficient DeffL(td) on diffusion time td at 20 °C (1) and 25 °C (2), Dp = 1.9·10−12 m2/s (b). Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Figure 67. The dependence of the blood lipid self-diffusion coefficient DL(td) on diffusion time td at 20 °C (1) and 25 °C (2) (a). The dependence of effective blood lipid self-diffusion coefficient DeffL(td) on diffusion time td at 20 °C (1) and 25 °C (2), Dp = 1.9·10−12 m2/s (b). Reprinted with permission from [171]. Copyright 2016 Springer Nature.
Membranes 11 00385 g067
Figure 68. 1 H spin–echo NMR spectra at low gradient amplitude pulses g: (a) RBCs suspension, g = 50 G/cm (0–1.70 ppm: (CH2)n and CH3 groups of lipids; 6–9 ppm: protein component (protons of aromatic rings and NH and OH groups)); (b) fullerene derivative 1 in aqueous solution, g = 55 G/cm (1.5–2.5 ppm: 10H, m, -CH2-CH2-CH2; 2.6–3.7 ppm: 20H, m, -CH2-S and -CH2-SO3Na); (c) RBCs suspension with added fullerene derivative 1, g = 52 G/cm. The signals of the fullerene derivative components in suspension are indicated by vertical lines. Reprinted with permission from [172]. Copyright 2018 Elsevier.
Figure 68. 1 H spin–echo NMR spectra at low gradient amplitude pulses g: (a) RBCs suspension, g = 50 G/cm (0–1.70 ppm: (CH2)n and CH3 groups of lipids; 6–9 ppm: protein component (protons of aromatic rings and NH and OH groups)); (b) fullerene derivative 1 in aqueous solution, g = 55 G/cm (1.5–2.5 ppm: 10H, m, -CH2-CH2-CH2; 2.6–3.7 ppm: 20H, m, -CH2-S and -CH2-SO3Na); (c) RBCs suspension with added fullerene derivative 1, g = 52 G/cm. The signals of the fullerene derivative components in suspension are indicated by vertical lines. Reprinted with permission from [172]. Copyright 2018 Elsevier.
Membranes 11 00385 g068
Figure 69. 1 H spin–echo NMR spectra at high gradient amplitude pulses g = 700 G/cm: (a) RBCs suspension; (b) fullerene derivative 1 in aqueous solution; (c) RBCs suspension with added fullerene derivative 1. The signals of the fullerene derivative components in suspension are indicated by vertical lines. Reprinted with permission from [172]. Copyright 2018 Elsevier.
Figure 69. 1 H spin–echo NMR spectra at high gradient amplitude pulses g = 700 G/cm: (a) RBCs suspension; (b) fullerene derivative 1 in aqueous solution; (c) RBCs suspension with added fullerene derivative 1. The signals of the fullerene derivative components in suspension are indicated by vertical lines. Reprinted with permission from [172]. Copyright 2018 Elsevier.
Membranes 11 00385 g069
Figure 70. Diffusion decays at different diffusion times td (indicated in the figure) of molecules 2 in RBCs suspension. Insert: diffusion decay of compound 2 in aqueous solution [173]. Reprinted with permission from [173]. Copyright 2021 Springer Nature.
Figure 70. Diffusion decays at different diffusion times td (indicated in the figure) of molecules 2 in RBCs suspension. Insert: diffusion decay of compound 2 in aqueous solution [173]. Reprinted with permission from [173]. Copyright 2021 Springer Nature.
Membranes 11 00385 g070
Table 1. Calculated from 1 H and 7 Li relaxation data (Dcalc) and experimentally measured (Dexp) water molecule and lithium cation self-diffusion coefficients.
Table 1. Calculated from 1 H and 7 Li relaxation data (Dcalc) and experimentally measured (Dexp) water molecule and lithium cation self-diffusion coefficients.
Amount of Water Molecules per One Sulfonate GroupDcalcH2O
m2/s
DexpH2O
m2/s
DcalcLi+
m2/s
DexpLi+
m2/s
45∙10−124∙10−122∙10−121∙10−12
20.53∙10−102∙10−104∙10−113∙10−11
Table 2. Hydration numbers of H+, Li+, Na+ and Cs+ cations in perfluorinated sulfonate (MF-4SC) and carboxylic (F-4CF) cation exchange membranes.
Table 2. Hydration numbers of H+, Li+, Na+ and Cs+ cations in perfluorinated sulfonate (MF-4SC) and carboxylic (F-4CF) cation exchange membranes.
Type of MembraneMF-4SC
Ionic FormH+Li+Na+Cs+
h ± 0.52.33.83.51.7
Type of membraneF-4CF
Ionic formH+Li+Na+Cs+
h ± 0.51.32.82.01.0
Table 3. Hydration number h at different absolute amount of water molecules per ionic site SO3 (λ) in acidic form Nafion 117 membrane [25].
Table 3. Hydration number h at different absolute amount of water molecules per ionic site SO3 (λ) in acidic form Nafion 117 membrane [25].
λ, [H2O]/[SO3H]Hydration Number h
1.9 ± 0.41.4 ± 0.5
3.2 ± 0.42.4 ± 0.5
4.4 ± 0.43.0 ± 0.3
5.8 ± 0.43.5 ± 0.3
6.4 ± 0.44.1 ± 0.3
7.4 ± 0.43.4 ± 0.3
12.0 ± 0.43.9 ± 0.3
17.5 ± 0.44.5 ± 0.5
Table 4. Hydration numbers h of Li+ cation in lithium form of Nafion 117 membrane at different water contents λ [109].
Table 4. Hydration numbers h of Li+ cation in lithium form of Nafion 117 membrane at different water contents λ [109].
λ, [H2O/SO3]0.92.04.05.77.410.712
H0.6 ± 0.31.2 ± 0.52.1 ± 0.52.6 ± 0.52.9 ± 0.54.2 ± 1.05.0 ± 1.0
Table 5. Crystallography radii, Stokes–Einstein hydrodynamic ion radii, hydration numbers (h) of Li+, Na+, and Cs+ cations in appropriate MSC membrane ionic forms at RH = 95% and in equimolar aqueous salt chloride solutions [17].
Table 5. Crystallography radii, Stokes–Einstein hydrodynamic ion radii, hydration numbers (h) of Li+, Na+, and Cs+ cations in appropriate MSC membrane ionic forms at RH = 95% and in equimolar aqueous salt chloride solutions [17].
CationLi+Na+Cs+
Crystallography ionic radius, Å [30] 0.691.021.67
Stokes–Einstein hydrodynamic ionic radius, Å [31] 2.381.841.19
Stokes–Einstein hydrodynamic radius, estimated from ionic diffusion coefficient in chloride aqueous solution at infinite dilute concentration2.702.201.50
Total water uptake of membrane (λ)24.0021.0016.00
Water amount per membrane sulfonate group (λs)13.8010.308.10
Hydration number of cations (h) in membrane4.10 ± 1.005.00 ± 1.003.10 ± 1.00
Hydration number of cations (h) in aqueous solution [103,104]4.004.603.90
Table 6. Ionic conductivities calculated based on diffusion coefficients of water and Li+ cations and measured for MF-4SC and F-4CF membranes in different ionic forms at the relative humidity of 90% [30].
Table 6. Ionic conductivities calculated based on diffusion coefficients of water and Li+ cations and measured for MF-4SC and F-4CF membranes in different ionic forms at the relative humidity of 90% [30].
Ionic FormσMF-4SC/S∙cm−1σMF-4SC/S∙cm−1
ExperimentCalculationExperimentCalculation
H+2.8·10−24.3·10−21.5·10−71.5·10−6
Li+6.2·10−31.1·10−2 (see a)1.0·10−34.3·10−3
6.5·10−3 (see b)
Na+6.1·10−31.2·10−21.2·10−34.6·10−3
Cs+3.4·10−48.3·10−43.7·10−47.2·10−4
The calculations were carried out based on diffusion coefficients MF-4SC, F-4CF [20]: a water molecules and b lithium cations.
Table 7. Activation energies of Li+ and water molecule self-diffusion in Li+ ionic form of Nafion membrane Ea at different water content λ [109].
Table 7. Activation energies of Li+ and water molecule self-diffusion in Li+ ionic form of Nafion membrane Ea at different water content λ [109].
λ, [H2O/SO3]0.92.04.05.77.410.7
Ea Li+ cations self-diffusion, kJ/mol40.3 ± 2.038.4 ± 2.028.8 ± 2.025.0 ± 2.025 ± 2.027.8 ± 2.0
Ea water molecules self-diffusion, kJ/mol-28.8 ± 2.025.0 ± 2.023.0 ± 2.021.1 ± 2.019.2 ± 2.0
Table 8. Moisture content λ, amount water molecules per cation, self-diffusion coefficients Ds, at 20 °C, self-diffusion activation energies Ea of Li+, Na+, and Cs+ cations in Nafion membrane, MSC membrane at RH = 95% and chloride aqueous solutions [17,109].
Table 8. Moisture content λ, amount water molecules per cation, self-diffusion coefficients Ds, at 20 °C, self-diffusion activation energies Ea of Li+, Na+, and Cs+ cations in Nafion membrane, MSC membrane at RH = 95% and chloride aqueous solutions [17,109].
Membrane TypeCationMoisture Content λ, Amount Water Molecules Per CationCation Self-Diffusion Coefficient at 20 °C Ds, m2/sCation Self-Diffusion Activation Energy Ea, kJ/mol
NafionLi+12(1.5 ± 0.1)10−1020.5 ± 1.0
Na+10(2.1 ± 0.3)10−1019.3 ± 1.5
Cs+4(0.6 ± 0.2)10−1024.8 ± 1.5
MSC [17]Li+243.7∙10−1017.6
Na+214.4∙10−1018.1
Cs+168.3∙10−1016.5
Chloride aqueous solutionLi+24(8.2 ± 0.3)10−1017.1 ± 0.5
Na+21(1.1 ± 0.2)10−918.3 ± 0.6
Cs+16(1.7 ± 0.2)10−916.8 ± 0.6
Table 9. Water content λ, hydration number h, self-diffusion coefficients Ds, calculated conductivity σc, measured conductivity σe of Li+, Na+, Cs+ cations in Li+, Na+, Cs+ Nafion 117 membrane ionic forms, RH = 98%, t = 20 °C [109].
Table 9. Water content λ, hydration number h, self-diffusion coefficients Ds, calculated conductivity σc, measured conductivity σe of Li+, Na+, Cs+ cations in Li+, Na+, Cs+ Nafion 117 membrane ionic forms, RH = 98%, t = 20 °C [109].
Ionic FormWater Amount Per Sulfonated Group λHydration Number hCation Self-Diffusion Coefficient Ds, m2/s Calculated Ionic Conductivity, σc S/cmMeasured Ionic Conductivity, σe S/cm
Li+125.0 ± 1.0(1.5 ± 0.1)10−10(1.6 ± 0.1)10−2(1.3 ± 0.1)10−2
Na+106.0 ± 1.0(2.0 ± 0.3)10−10(2.0 ± 0.3)10−2(1.1 ± 0.1)10−2
Cs+41.0 ± 0.2(0.6 ± 0.2)10−10(6.0 ± 0.2)10−3(2.3 ± 0.3)10−3
Table 10. Experimental values of ionic conductivity at 25 °C and conductivity activation energies of H+, Li+, Na+, Cs+ cations in MSC membrane with different humidity. Ion exchange capacity is 2.5 mg-eq/g [17].
Table 10. Experimental values of ionic conductivity at 25 °C and conductivity activation energies of H+, Li+, Na+, Cs+ cations in MSC membrane with different humidity. Ion exchange capacity is 2.5 mg-eq/g [17].
RH, %95755832
Ionic FormEa, kJ/molσexp
mS/cm
Ea, kJ/molσexp
mS/cm
Ea, kJ/molσexp
mS/cm
Ea, kJ/molσexp
mS/cm
H4.37.8116.0123.0230.600
Li7.51.9300.5390.2520.008
Na26.03.4310.7400.2680.010
Cs17.04.6320.8370.3600.020
Table 11. Partial diffusion coefficients of water and ethanol and separation factor (α) of their mixture in membranes at 50 °C [48].
Table 11. Partial diffusion coefficients of water and ethanol and separation factor (α) of their mixture in membranes at 50 °C [48].
MembraneDH2O·10−10/m2∙s−1DEtOH·10−10/m2∙s−1α
MF-4SC6.01.50048
PA-10.60.15060
PA-240.00.10010,000
PAA—PSF1.30.026800
Table 12. Self-diffusion coefficients of WSFD molecules in aqueous solutions and RBCs suspension.
Table 12. Self-diffusion coefficients of WSFD molecules in aqueous solutions and RBCs suspension.
CompoundAqueous SolutionRBCs Suspension
Ds1w·1010, m2/sDs2w·1011, m2/sDs1s·1010, m2/sDs2s·1011, m2/sDs3s·1012, m2/s
14.1 ± 0.47.4 ± 0.75.5 ± 0.83.9 ± 0.65.5 ± 0.8
24.3 ± 0.87.5 ± 1.57.1 ± 1.44.4 ± 0.95.0 ± 1.0
31.2 ± 0.14.9 ± 0.58.0 ± 1.03.8 ± 0.66.0 ± 1.0
Table 13. The self-diffusion coefficients Ds3s, population p3, and lifetimes τ of WSFD molecules in RBCs [173].
Table 13. The self-diffusion coefficients Ds3s, population p3, and lifetimes τ of WSFD molecules in RBCs [173].
CompoundDs3s·1012, m2/sp3(0)τ, ms
15.5 ± 0.80.33440 ± 70
25.0 ± 1.00.13470 ± 70
36.0 ± 1.00.061200 ± 300
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Volkov, V.I.; Chernyak, A.V.; Avilova, I.A.; Slesarenko, N.A.; Melnikova, D.L.; Skirda, V.D. Molecular and Ionic Diffusion in Ion Exchange Membranes and Biological Systems (Cells and Proteins) Studied by NMR. Membranes 2021, 11, 385. https://doi.org/10.3390/membranes11060385

AMA Style

Volkov VI, Chernyak AV, Avilova IA, Slesarenko NA, Melnikova DL, Skirda VD. Molecular and Ionic Diffusion in Ion Exchange Membranes and Biological Systems (Cells and Proteins) Studied by NMR. Membranes. 2021; 11(6):385. https://doi.org/10.3390/membranes11060385

Chicago/Turabian Style

Volkov, Vitaliy I., Alexander V. Chernyak, Irina A. Avilova, Nikita A. Slesarenko, Daria L. Melnikova, and Vladimir D. Skirda. 2021. "Molecular and Ionic Diffusion in Ion Exchange Membranes and Biological Systems (Cells and Proteins) Studied by NMR" Membranes 11, no. 6: 385. https://doi.org/10.3390/membranes11060385

APA Style

Volkov, V. I., Chernyak, A. V., Avilova, I. A., Slesarenko, N. A., Melnikova, D. L., & Skirda, V. D. (2021). Molecular and Ionic Diffusion in Ion Exchange Membranes and Biological Systems (Cells and Proteins) Studied by NMR. Membranes, 11(6), 385. https://doi.org/10.3390/membranes11060385

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop