Phase Equilibrium and Interdiffusion in Poly(Vinyl Methyl Ether)-Water System

The phase state diagram of the poly(vinyl methyl ether)-water system in a wide concentration range was obtained by the optical interferometry method. It was shown that this system was characterized by a complicated phase equilibrium with two lower critical solution temperatures, one of which was located in the concentrated region at 21 °C, and the other one in the region of a dilute solution at 31 °C. In the framework of the Flory–Huggins theory, pair interaction parameters were calculated for different parts of the binodal curves, and an attempt was made to reverse simulate the diagram in different conditions. It was suggested that the unusual character of the diagram was associated with the formation of a complicated complex between PVME and water in the middle region of the compositions. Concentration profiles for different temperatures were constructed. For the first time for this system, the numerical values of the diffusion coefficients of poly(vinyl methyl ether) (PVME) into water and water in PVME were obtained. Concentration and temperature dependences of diffusion coefficients were constructed and analyzed. The kinetics of water sorption in PVME was plotted, the clustering integral was calculated, and the approximate number of molecules in a water cluster was estimated. It was shown that in the dilute solution region upon passing through the binodal curve, the interphase disappeared immediately, and the remaining fluctuation of the concentration decreased in size with time. The kinetics of this process was estimated from the change in the size of such a particle.


Introduction
The study of phase equilibria and translational mobility of the components of polymers aqueous solutions is currently a promising and important direction in the field of physical chemistry of macromolecular compounds [1,2]. First of all, this is associated with the creation of biodegradable polymeric materials, in particular, packaging, insulating films, tapes, casings, coatings, etc. [3,4]. Modern medicine increasingly relies on effective and targeted treatment using biocompatible transdermal polymer systems supplied with ingredients, including drugs, that regulate their phase state in the human body [5]. Thus, technologies for the synthesis and application of intelligent materials based on hydrogels that are sensitive to temperature and allow regulating blood flow in blood vessels are actively developed [6,7]. It is obvious that aqueous solutions of such polymers should have a lower critical solution temperature (LCST) located in the area of the human body's vital activity 20-36 • C. These conditions are met by poly(vinyl methyl ether) (PVME) and its copolymers.
Currently, aqueous solutions of PVME are traditionally used in these studies as model systems. Numerous works  show that the PVME-water system is characterized by a specific phase state diagram: two LCST coexist simultaneously in a binary system at the temperature-concentration field This article presented new experimental data on the phase equilibrium of the PVME-water system in a wide range of temperatures and concentrations. The thermodynamic parameters of mixing were also calculated, on the basis of which we would consider, in a new way, the features of the interaction of the components. For the first time, the coefficients of translational diffusion of PVME and water were obtained and generalized, as well as the results of simulation binodal curves in the framework of the quasi-chemical theory of solutions [1,31].

Materials and Methods
Poly (vinyl methyl ether) (PVME) and distilled water were used as the objects of the research. The original PVME (Basf, Ludwigshafen, Germany) was a 50% aqueous solution: M n = 45,000 g/mol, M w /M n = 1.73, T g = −31 • C. The refractive index of water is n 25 = 1.332. PVME films were formed from its 50% aqueous solution by casting onto the polyethylene terephthalate substrate at room temperature. After two days of drying at room temperature and 40% humidity, the films were annealed at 328 K in a thermo-cabinet XF050 (Moscow, Russia) and then for 4 h in a vacuum (VUP-5, Moscow, Russia) at a pressure of 10-5 mmHg. The remaining water content of the samples prepared in this way did not exceed 1.5 wt.%, according to thermal gravimetric analysis data (Netzsch TG209F1 Iris, Selb, Germany).
Optical interferometry [32]. The method was based on the principle of in situ registration of optical density distribution in the area of the conjugation of polymer and water and recording its change in time under isobaric-isothermal conditions of the process. The measurement method consisted of fixation of a PVME sample of 3 mm × 10 mm in size and about 120 µm thick between the diffusion cell glasses, the inner surfaces of which were covered with a layer of translucent metal (Ni+Cr) with a high reflection index. Using special fixtures, which were~5% thinner than the PVME film, the sample was put into optical contact with the surface of the plates. While a small wedge angle θ ≤ 5 • was established between the glasses. The polymer sample was placed in a flat capillary of the diffusion cell so that the interference strips oriented parallel to the edge of the wedge were to the longest side of the PVME sample. After assembly, the cell was thermoregulated at a set temperature for at least 30 min. Then, the capillary was filled with water.
All measurements were carried out on an ODA-2 IPCE diffusiometer (Moscow, Russia) [33] in the temperature range from 20 to 35 • C. The experiment was carried out in the heating-cooling mode with a step of 1 • C and thermostatting at each stage for at least 10-15 min. Special attention was paid to the kinetics of the phase particles' formation and dissolution. For this purpose, a targeted photographic survey of the phase decay zones at different temperatures was carried out.
A helium-neon laser (λ = 632.8 nm) was used as a light source. Methods of processing of the interferograms, interdiffusion zones, and phase diagrams construction did not differ from those described earlier [34,35].
Refractometry. Information about the temperature dependences of the refractive indices of polymers was of fundamental importance in interpreting the optical density distribution profiles in the interdiffusion zone. In the range from 20 to 60 • C, a large difference was observed in the values of the PVME and water refractive indices (∆n is about 0.137). An increment of the refractive index per interference strip equaled to~0.022, providing a sufficiently high sensitivity in identifying changes in the composition in the phase conjugation zone. The refractometry data of the studied samples were presented earlier [35].
Sorption measurements [32]. The experiments on the sorption of water vapor by the polymer were carried out by the desiccator method at 25 ± 1 • C. The polymer was preliminarily kept in a desiccator with calcined CaCl 2 granules at 0% humidity until a constant weight was established, then they were evacuated at a residual pressure of 10 −5 mmHg. We measured the kinetics of changes in the mass of samples placed in desiccators with different vapor activity. Dried samples with a thickness of about 100 µm were placed in desiccators with different humidity levels (p/p s 0.2, 0.4, 0.49, 0.67, 0.9, and 0.97). At regular intervals, the weight measurements of each sample were made on HR-200 A&D Co. LTD scales (Tokyo, Japan). The measurements were carried out until the sorption equilibrium was established, which, in the kinetic dependencies, was determined by the constancy of the sample mass over a period of 2-3 times longer than the time it was reached.

Conjugate Phase Zones
Optical interferometry allowed us to observe in situ the process of interdiffusion and phase decay at different temperatures in a wide range of compositions. Typical interferograms of diffusion zones obtained in the mode of a stepwise rise of temperature are shown in Figure 1. At T < 21 ± 0.5 • C, PVME is completely soluble in water. The interdiffusion zone at these temperatures corresponds to the region of smooth and continuous change in composition during the transition from PVME to water. As the temperature rises to T = 21 • C, the system undergoes phase decomposition in the region of the concentrated solution (I in Figure 1a). At this temperature, the interdiffusion zone is a superposition of three regions: region A-the dissolution of water in PVME, region B-the dissolution of PVME in water, region C-phase decay, within which are concentrated (localized) dispersions of PVME in water corresponding to the compositions of the coexisting phases of this area of the phase state diagram. With a further temperature increase to T = 31 ± 0.5 • C, the system undergoes a second phase decomposition in the region of dilute solutions of PVME in water ( Figure 1b). The interdiffusion zone, in this case, is a superposition of five regions: region A -the dissolution of water in PVME, region B-the dissolution of PVME in water, region of phase decomposition in the zone of concentrated (C) and dilute solutions (D), and also the intermediate region of solubility between them. Finally, at T = 34 ± 0.5 • C, zones I and II merge into a single heterophase region III (Figure 1c), and the system as a whole takes on a classical form characteristic of systems with amorphous separation [1,31]. established, which, in the kinetic dependencies, was determined by the constancy of the sample mass over a period of 2-3 times longer than the time it was reached.

Conjugate Phase Zones
Optical interferometry allowed us to observe in situ the process of interdiffusion and phase decay at different temperatures in a wide range of compositions. Typical interferograms of diffusion zones obtained in the mode of a stepwise rise of temperature are shown in Figure 1. At T < 21 ± 0.5 °C, PVME is completely soluble in water. The interdiffusion zone at these temperatures corresponds to the region of smooth and continuous change in composition during the transition from PVME to water. As the temperature rises to T = 21 °C, the system undergoes phase decomposition in the region of the concentrated solution (I in Figure 1a). At this temperature, the interdiffusion zone is a superposition of three regions: region A-the dissolution of water in PVME, region B-the dissolution of PVME in water, region C-phase decay, within which are concentrated (localized) dispersions of PVME in water corresponding to the compositions of the coexisting phases of this area of the phase state diagram. With a further temperature increase to T = 31 ± 0.5 °C, the system undergoes a second phase decomposition in the region of dilute solutions of PVME in water ( Figure  1b). The interdiffusion zone, in this case, is a superposition of five regions: region A -the dissolution of water in PVME, region B-the dissolution of PVME in water, region of phase decomposition in the zone of concentrated (C) and dilute solutions (D), and also the intermediate region of solubility between them. Finally, at T = 34 ± 0.5 °C, zones I and II merge into a single heterophase region III (Figure 1c), and the system as a whole takes on a classical form characteristic of systems with amorphous separation [1,31].
The analysis showed that even for a pseudo variant binary system in the presence of one variable parameter (temperature), the simultaneous existence of four phases is impossible. Thus, either the PVME-water system is not two-component (for example, due to the presence of complicated complexes between the components), or the influence of factors other than temperature and constant pressure should be taken into account, and then the system is a projection of a more complex volumetric diagram.  The analysis showed that even for a pseudo variant binary system in the presence of one variable parameter (temperature), the simultaneous existence of four phases is impossible. Thus, either the PVME-water system is not two-component (for example, due to the presence of complicated complexes Polymers 2020, 12, 2445 5 of 18 between the components), or the influence of factors other than temperature and constant pressure should be taken into account, and then the system is a projection of a more complex volumetric diagram. Figure 2 shows the concentration profiles of the compositions' distribution in accordance with the above interferograms. The gray areas of the concentration profiles correspond to the zones of the system phase decomposition, on both sides of which the compositions of the coexisting phases (ϕ 1 , ϕ 1 , ϕ 1 , ϕ 1 ) are established. It is shown that ϕ 1 , ϕ 1 , ϕ 1 , ϕ 1 under isothermal conditions are constant in time, but they change with temperature. Attention should be paid to the reproducibility of the values ϕ 1 , ϕ 1 , ϕ 1 , ϕ 1 in the modes of increasing and decreasing the temperature, which allows us to speak of the equilibrium of the obtained values of the coexisting phases. temperature corresponds to (a) 85 min, (b) 370 min, and (c) 397 min from the beginning of the experiment. Figure 2 shows the concentration profiles of the compositions' distribution in accordance with the above interferograms. The gray areas of the concentration profiles correspond to the zones of the system phase decomposition, on both sides of which the compositions of the coexisting phases (1, 1, 1, 1) are established. It is shown that 1, 1, 1, 1 under isothermal conditions are constant in time, but they change with temperature. Attention should be paid to the reproducibility of the values 1, 1, 1, 1 in the modes of increasing and decreasing the temperature, which allows us to speak of the equilibrium of the obtained values of the coexisting phases.
It is clearly seen that the concentration profile for the mutual solubility regions is characterized by smooth curves (to the left and to the right from I in Figure 2a). In the zone I of phase decay, a sharp change in concentration can be observed. For a temperature of 32 °C, the concentration profile consists of three sections of smooth curves corresponding to a concentrated solution, an intermediate region, and a diluted solution, and also two sections of a sharp change in concentration for the zones of phase decomposition I and II. For a temperature of 35 °C, the concentration profile of such a system has a classical form for a phase-separated system.  It is clearly seen that the concentration profile for the mutual solubility regions is characterized by smooth curves (to the left and to the right from I in Figure 2a). In the zone I of phase decay, a sharp change in concentration can be observed. For a temperature of 32 • C, the concentration profile consists of three sections of smooth curves corresponding to a concentrated solution, an intermediate region, and a diluted solution, and also two sections of a sharp change in concentration for the zones of phase decomposition I and II. For a temperature of 35 • C, the concentration profile of such a system has a classical form for a phase-separated system.

Phase State Diagram
The phase diagram of the PVME-water system was constructed from the compositions of the coexisting phases ( Figure 3a). It is characterized by the presence of two LCSTs (LCST 1 and LCST 2 ) and is in good agreement with the literature data [8]. Alekseev's straight-line diameter rule was found, according to which the midpoints of the LCST 1 and LCST 2 regions form straight lines (dashed lines in Figure 3) passing through the minimums of the delamination regions boundaries. It can be seen that in the concentrated solution region, LCST 1 is localized at 21 ± 0.5 • C and a water concentration of 0.32 vol.%, and in the diluted solution region, LCST 2 corresponds to 31 ± 0.5 • C and a water concentration of about 0.95 vol.%. Upon reaching 34 • C at a water concentration of approximatelỹ 0.6 vol.%, regions I and II merge, forming a single zone III of phase decomposition. It should be noted that the phase reversal region corresponds to the Alekseev diameter for the studied systems.

Phase State Diagram
The phase diagram of the PVME-water system was constructed from the compositions of the coexisting phases (Figure 3a). It is characterized by the presence of two LCSTs (LCST1 and LCST2) and is in good agreement with the literature data [8]. Alekseev's straight-line diameter rule was found, according to which the midpoints of the LCST1 and LCST2 regions form straight lines (dashed lines in Figure 3) passing through the minimums of the delamination regions boundaries. It can be seen that in the concentrated solution region, LCST1 is localized at 21 ± 0.5 °C and a water concentration of 0.32 vol.%, and in the diluted solution region, LCST2 corresponds to 31 ± 0.5 °C and a water concentration of about 0.95 vol.%. Upon reaching 34 °C at a water concentration of approximately 0.6 vol.%, regions I and II merge, forming a single zone III of phase decomposition. It should be noted that the phase reversal region corresponds to the Alekseev diameter for the studied systems. Figure 3b also shows microphotographs illustrating the phase structure of dispersions formed along isotherms in different zones of the phase diagram. If in the zone I, these are small formations from 10 to 50 μm, which retain their size for a long time, in zone II, the dispersed phase of PVME first precipitates in the form of small particles 10-20 μm, then the particle size rapidly increases to 400 μm, and inside these formations, a secondary phase decay occurs. Zone III is characterized by an abundance of fine phases up to 20 μm. We assume that the presence of large phase particles in this region is associated with the kinetics of their formation at the initial stage of phase decomposition.  Figure 3b also shows microphotographs illustrating the phase structure of dispersions formed along isotherms in different zones of the phase diagram. If in the zone I, these are small formations from 10 to 50 µm, which retain their size for a long time, in zone II, the dispersed phase of PVME first precipitates in the form of small particles~10-20 µm, then the particle size rapidly increases to 400 µm, and inside these formations, a secondary phase decay occurs. Zone III is characterized by an abundance of fine phases up to 20 µm. We assume that the presence of large phase particles in this region is associated with the kinetics of their formation at the initial stage of phase decomposition. In accordance with the classical Flory-Huggins theory [1,31], using the previously developed method [35], we calculated pair interaction parameters χ for different regions of the phase diagram using the equation: where ϕ 1 and ϕ 2 are volume concentrations of the first and second components, respectively, r 1 and r 2 are their polymerization degrees, ' and " refer to different coexisting phases, χ is the Flory-Huggins interaction parameter assuming the absence of its concentration dependence and is the average value for the pair interaction parameters χ 12 = χ + ∂χ ∂ϕ 1 ϕ 1 and χ 21 = χ + ∂χ ∂ϕ 2 ϕ 2 . Figure 4a shows that the temperature dependences of the Flory-Huggins parameter are linear. Dashed line 4 shows the calculated value χ cr = 1 for the PVME-water system.
The expansion of the parameter χ into the entropy (first term of the polynomial) and enthalpy (second term of the polynomial) components (χ S and χ H , respectively) shows that the contributions of χ H are close to each other, while χ S for all three bands is significantly different. Based on this, we assume that the thermodynamics of mixing of components in different zones is entropic in nature.
It can be seen that the values of the pair interaction parameter for zone II are close to χ cr , and the critical temperature, estimated from the point of intersection of the dependence χ − 1/T, coincides with the experimentally found value. The calculated values of LCST for zones I and III, estimated by the same method, differ significantly from the experimental ones. To explain this fact, the values of the pair interaction parameters χ 12 and χ 21 were determined. It was found that the values of the pair interaction parameters χ 12 and χ 21 in zones I and III differ significantly from each other, while for zone II, they completely coincide. This suggests that zones I and III are characterized by a concentration dependence of the Flory-Huggins parameter associated with the low thermal stability of the PVME·nH 2 O complex.
Nevertheless, based on the obtained equations for χ − 1/T, we tried to simulate fragments of the phase state diagram (Figure 4b) according to the equations for binodal and spinodal curves [1,31]: The analysis showed that for the diluted solution region, there is a good agreement with the experimental data for the PVME-water system; 6 in Figure 4b corresponds to binodal and spinodal curves. In the concentrated solution area (left side of the diagram), simulation in the framework of the same approach gave us a mixed result. Good agreement with the experimental data (curves 5 in Figure 4b) was obtained for the case of extrapolating the left side of the diagram to the entire concentration range under the assumption that one of the components is PVME, and the other is the PVME·nH 2 O complex with several water molecules from 2 to 6 hydrated for each PVME functional group. In this case, the experimental points (black circles) correlate well with the calculated binodal curve (curve 5). It should bear in mind that the calculated value of LCST 1 is~10 • C lower than that obtained experimentally.  Thus, the assessment of the phase state diagram and the thermodynamics of mixing of components in different temperature-concentration regions allowed us to assume that the phase  (4) shows the calculated value of χ cr , according to the Flory-Huggins theory. 5-phase equilibrium simulation for the PVME-PVME·nH 2 O system according to the χ − 1/T dependence for zone I, 6-phase equilibrium simulation for the PVME-water system according to the χ − 1/T dependence for zone II. Black points are experimental data of phase equilibrium.
Thus, the assessment of the phase state diagram and the thermodynamics of mixing of components in different temperature-concentration regions allowed us to assume that the phase equilibrium in the PVME-water system is due to the formation of the PVME·nH 2 O complex, and in general, the system should be considered not as a binary but as a three-component PVME-water-PVME·nH 2 O. Then, the left part of the diagram (to the left of the vertical dashed line) corresponds to the amorphous decomposition of the PVME-PVME·nH 2 O system, and the right side to the amorphous decomposition of PVME-water, when the PVME·nH 2 O complex is destroyed. Figure 5a shows the kinetic curves of interval sorption of water vapor at different humidity. It can be seen that the higher the vapor activity p/p s under isothermal conditions, the larger the slope of the initial portion of the curve, and the greater the amount of water sorbed by the polymer. In accordance with Rogers' classification [32], the kinetics of the establishment of an equilibrium state has a Fick character. It is known that in this case, the boundary concentrations of the sorbate quickly reach equilibrium values and remain unchanged during the entire process of establishing sorption equilibrium. Usually, such conditions are met when the rate of the structure rearrangement in the sorbent surface layer is higher than the rate of diffusion of sorbate molecules in the bulk of the polymer.

Sorption
Polymers 2020, 12, x FOR PEER REVIEW 9 of 19 equilibrium in the PVME-water system is due to the formation of the PVME·nH2O complex, and in general, the system should be considered not as a binary but as a three-component PVME-water-PVME·nH2O. Then, the left part of the diagram (to the left of the vertical dashed line) corresponds to the amorphous decomposition of the PVME-PVME·nH2O system, and the right side to the amorphous decomposition of PVME-water, when the PVME·nH2O complex is destroyed. Figure 5a shows the kinetic curves of interval sorption of water vapor at different humidity. It can be seen that the higher the vapor activity p/ps under isothermal conditions, the larger the slope of the initial portion of the curve, and the greater the amount of water sorbed by the polymer. In accordance with Rogers' classification [32], the kinetics of the establishment of an equilibrium state has a Fick character. It is known that in this case, the boundary concentrations of the sorbate quickly reach equilibrium values and remain unchanged during the entire process of establishing sorption equilibrium. Usually, such conditions are met when the rate of the structure rearrangement in the sorbent surface layer is higher than the rate of diffusion of sorbate molecules in the bulk of the polymer. For determining the effective diffusion coefficients, kinetic curves were constructed in relative

Sorption
where / is the amount of water absorbed at each moment of time t; / is equilibrium amount of sorbed water. The diffusion coefficient was determined from the initial rectilinear section of the integral curve in accordance with the formula [32]: where l is a sample thickness in cm; D is a diffusion coefficient (cm 2 /s); α is the angle of inclination of the straight line according to the dependencies The numerical values of the obtained effective diffusion coefficients are presented in Table 1. For determining the effective diffusion coefficients, kinetic curves were constructed in relative coordinates (x/m) (x/x max ) − t 1/2 , where x/m is the amount of water absorbed at each moment of time t; x/x max is equilibrium amount of sorbed water. The diffusion coefficient was determined from the initial rectilinear section of the integral curve in accordance with the formula [32]: where l is a sample thickness in cm; D is a diffusion coefficient (cm 2 /s); α is the angle of inclination of the straight line according to the dependencies (x/m) (x/x max ) − t 1/2 . The numerical values of the obtained effective diffusion coefficients are presented in Table 1. Water vapor sorption isotherms were plotted from the equilibrium parts of the kinetic curves (Figure 5b). The isotherm has a weakly expressed S-shaped character and belongs to the IV type according to Rogers classification [36]. This type of isotherm is traditionally described in the framework of the "double sorption" model [37]: The model assumes that some of the sorbed water molecules dissolve in PVME in accordance with the Flory-Huggins solutions theory (C FH ). Another smaller part of the sorbate is immobilized on the accessible active sites of macromolecules monomer units, which is described by the Langmuir monolayer sorption equation (C L ). Accordingly, the relationship between the relative water vapor pressure and the concentration of sorbate in PVME is described by the Flory-Huggins equation: where P P S is the relative vapor pressure, C FH is the bulk capacity of the sorbate in the sorbate + sorbent system, χ is the Flory-Huggins interaction parameter, r is the degree of polymerization of the polymer sorbent, and the Langmuir equation is as follows: where C ∞ is the sorption capacity of active centers in PVME, k is the energy characteristic. The constants of Equations (6) and (7) were determined by the numerical method of "gradient descent", minimizing the average deviation of the calculated isotherm from the experimentally obtained one.
In the initial region of the isotherm up to p/p s < 0.3, where the Langmuir model is predominantly realized, water molecules fill free active centers, which are -O-CH 3 groups located in vacancies of the excess free volume of the polymer. The share of such centers, estimated according to the ratio C ∞ /C max , is 15 mass.%.
In the range of relative humidity from 0.35 to 1.0, the sorption isotherms are well approximated by the Flory-Huggins equation with pair interaction parameters χ varying in the range from 2 to 3. These data correlate well with the χ data in Figure 4a. Additional information on the state of water molecules in PVME was obtained using the Starkweather equations [38]. Recall that the quantitative characteristic of the formation of associate-clusters of water molecules in the sorbent matrix is the clustering integral G 11 υ 1 (Figure 6a): where υ 1 , ϕ 1 , a 1 is the molar volume, volume composition, and activity of the sorbate, respectively. For an ideal solution, For the systems studied by us, the border value of the sorption capacity characterizing the start of cluster formation, P P S = 0.4 (marked with a dotted line in Figure 6). The average number of solvent molecules in a cluster can be obtained according to the Starkweather-Brown equation [38]:

Particles Dissolution Kinetics of Dispersed Phases
This study is of practical interest associated with the formation and subsequent destruction of functional phase particles, which were previously mentioned in [6,7]. The starting point of the model study was the formation of a gradient diffusion zone in the PVME-water system in the dilute region According to the sorption data at 25 • C and an increase in volume water concentration from 0.5 to 6.9 vol.%, N c rises from 0.5 to 6-7 (Figure 6b). In this case, the concentration of 6.9 vol.% actually corresponds to the binodal curve in the concentrated solution region. Thus, it can be assumed that the phase decomposition in the zone I is associated with a rapid increase in the number of water molecules in the cluster, which is formed around the PVME functional group. While the phase decay in zone II is caused by the decomposition of a similar cluster with an increase in temperature to 32 • C.

Particles Dissolution Kinetics of Dispersed Phases
This study is of practical interest associated with the formation and subsequent destruction of functional phase particles, which were previously mentioned in [6,7]. The starting point of the model study was the formation of a gradient diffusion zone in the PVME-water system in the dilute region near the binodal curve. If, after the formation of the diffusion zone, the temperature is raised to 32 • C (above the LCST 2 ), a phase decomposition occurs in a certain region of compositions (Figure 7). Within a few seconds (up to 100 s), dispersions of the aqueous phase in the PVME matrix and the PVME phase in water are formed. These dispersions are separated by a phase boundary corresponding to the position of the phase reversal region on the phase diagram in the zone II of phase decay. Under isothermal conditions, the rapid growth of microphases from 10 to 100 μm is observed over time. Simultaneously, their coalescence occurs according to Smolukhovsky [39] with the formation of macroparticles ranging in size from 300 to 400 μm. The process of the Oswald ripening is also possible. Similar effects were observed earlier on oligomeric systems in [40,41].
It is interesting to note that with a repeated increase in temperature from 30 to 33 °C, a secondary phase decomposition is observed, which is clearly identified in macrophases (Figure 8a). This process is completely equilibrium and is reproduced in heating-cooling cycles. Under isothermal conditions, the rapid growth of microphases from 10 to 100 µm is observed over time. Simultaneously, their coalescence occurs according to Smolukhovsky [39] with the formation of macroparticles ranging in size from 300 to 400 µm. The process of the Oswald ripening is also possible. Similar effects were observed earlier on oligomeric systems in [40,41].
It is interesting to note that with a repeated increase in temperature from 30 to 33 • C, a secondary phase decomposition is observed, which is clearly identified in macrophases (Figure 8a). This process is completely equilibrium and is reproduced in heating-cooling cycles.
To study the mechanism of the dissolution kinetics and diffusion of PVME macromolecules in aqueous solutions, the case of a sharp cooling of dispersions after the formation of disperse systems is of great importance. That is, the figurative point of the systems is outside the boundary curve. over time. Simultaneously, their coalescence occurs according to Smolukhovsky [39] with the formation of macroparticles ranging in size from 300 to 400 μm. The process of the Oswald ripening is also possible. Similar effects were observed earlier on oligomeric systems in [40,41].
It is interesting to note that with a repeated increase in temperature from 30 to 33 °C, a secondary phase decomposition is observed, which is clearly identified in macrophases (Figure 8a). This process is completely equilibrium and is reproduced in heating-cooling cycles. To study the mechanism of the dissolution kinetics and diffusion of PVME macromolecules in aqueous solutions, the case of a sharp cooling of dispersions after the formation of disperse systems is of great importance. That is, the figurative point of the systems is outside the boundary curve.
It was found that at the first stage of the gradient structure evolution, the particle dissolution proceeds mainly due to local mass transfer in the interphase region. At the second stage, the interphase boundary transforms into a diffusion one, the dimensions of which are continuously changing. At the third stage, concentration profiles-vertical interference strips-reach the "particle" center. At this stage, the "particle" can be considered as a heterophase fluctuation [42], and the It was found that at the first stage of the gradient structure evolution, the particle dissolution proceeds mainly due to local mass transfer in the interphase region. At the second stage, the interphase boundary transforms into a diffusion one, the dimensions of which are continuously changing. At the third stage, concentration profiles-vertical interference strips-reach the "particle" center. At this stage, the "particle" can be considered as a heterophase fluctuation [42], and the process of its dissolving can be considered in the framework of the classical diffusion theory of finite size "particle" dissolution ( Figure 8b). From an experimental point of view, the kinetics of this stage of the process can be estimated from the change in the amplitude of the concentration profile (Figure 8c) using the equation: where Q is the amount of a substance in a nanoparticle (droplet), D is the diffusion coefficient, and t is time [43]. The diffusion coefficient of PVME in a diluted solution was determined using the ratio logC = log Const − 1 2 Dt. The numerical value of the calculated effective diffusion coefficient is then 2.5×10 −7 cm 2 /s.

Diffusion
It has been established that the diffusion fronts motion in the region of concentrated and dilute solutions in coordinates X − t 1/2 under isobaric isothermal conditions of the process is linear. This indicates the diffusion mechanism of the spontaneous mixing of components in the region of true solutions (below the binodal curve). Here, X is the diffusion coordinate at the selected concentration, t is the observation time. The numerical values of the interdiffusion coefficients with a step equal to the concentration increment were calculated by the Matano-Boltzmann method [43,44]. Figure 9 shows the concentration dependences of the interdiffusion coefficients for temperatures of 25, 30, and 35 • C (curves 1, 2, and 3, respectively). For comparison, the figure shows a phase diagram (curve 4). It is clearly seen that for all temperatures in the concentrated solutions region (to the left of the binodal curve), a decrease in logD is observed as we approach the two-phase region boundary. It should be noted that the calculated values of the effective diffusion coefficient according to the sorption data completely coincide with curve 1, supplementing it at lower composition values. Such behavior of the interdiffusion coefficient is associated with the presence of a labile region inside the amorphous decomposition diagram, at the boundary of which the thermodynamic multiplier dlna dlnϕ in the equation D = D * dlna dlnϕ tends to zero [45]. The temperature dependence of the interdiffusion coefficients is plotted in the coordinates of the Arrhenius equation and is shown in Figure 10. It is clearly seen that in the investigated range of temperatures and compositions, it has a linear character both in the case of diffusion of water molecules into the PVME matrix (1 in Figure 10) and during diffusion of PVME into water (2 in Figure  10). The slope of such straight lines allowed us to calculate the diffusion activation energy, the numerical values of which correspond to 5.8-6.5 kJ/mol. In the dilute solutions region, diffusion has a more complex character since, in this case, the contribution of the thermodynamic correction to the value of the interdiffusion coefficient is determined by the distance of the figurative point from the binodal curve. The concentration region of mutual solubility at 35 • C is very small here (about 2 vol.% PVME) and logD = −5.55. This value is in good agreement with the self-diffusion coefficient of water [10]. For a temperature of 25 • C, the interdiffusion coefficient sharply decreases with an increase in the PVME content in the solution, but at medium values of the concentration (up to the binodal curve), it is monotonic. We assume that this is due to the presence of a complex, which hinders the translational mobility of molecules. The behavior of curve 2 in Figure 9 in the middle region of compositions is the most interesting: after the interdiffusion coefficient sharp decrease in the dilute solution region (similar to how it was at 25 • C), logD passes through a minimum and slightly increases. We associate this effect with the approach of the system to the second phase decay (zone II), from which it is separated by 1 • C. Consequently, it can be assumed that the existing complex begins to undergo changes (practically disintegrate), which, on the one hand, leads to an increase in the diffusion rate of components in this zone, but on the other hand, predetermines a rapid phase decomposition.
The temperature dependence of the interdiffusion coefficients is plotted in the coordinates of the Arrhenius equation and is shown in Figure 10. It is clearly seen that in the investigated range of temperatures and compositions, it has a linear character both in the case of diffusion of water molecules into the PVME matrix (1 in Figure 10) and during diffusion of PVME into water (2 in Figure 10).
The slope of such straight lines allowed us to calculate the diffusion activation energy, the numerical values of which correspond to 5.8-6.5 kJ/mol.  Figure 10. Temperature dependence of the interdiffusion coefficients into the PVME matrix (1) and into the water matrix (2).

Conclusions
We have shown that the phase equilibrium in aqueous solutions of PVME is a complex projection of a multi-component system that coexists at different temperatures. Thus, in the temperature range below 31 °C, the phase diagram corresponds to amorphous separation in the PVME-PVME·nH2O system, while above 31 °C, the PVME·nH2O complex is destroyed, and the phase diagram corresponds to the binary PVME-water system. At the same time, the intermediate regions of mutual solubility are obviously associated with the inertial effects of dehydration of functional groups of PVME and the redistribution of unbound water molecules. The kinetics of the formation and dissolution of phase particles is consistent with the Oswald mechanism, and the phase size in the diluted solution region can reach 400 μm. In this case, a secondary phase decay is observed within such macrophases.
Sorption studies have shown that in the concentrated solution region, the number of water molecules in a cluster formed around the functional group of PVME sharply increases. According to our calculations, up to 6-7 water molecules are combined in such a cluster, which is in good agreement with the literature data.
For the first time for the PVME-water system, we obtained numerical values for the diffusion coefficient. The logD values (cm 2 /s) vary in the range from −7 to −5.5, with an increase in temperature from 25 to 35 °C. A good correlation is observed between the data obtained both from sorption measurements and from direct contact of polymers by the optical interferometry method. The analysis of the concentration dependences of the diffusion coefficients showed that the diffusion rate decreases sharply when approaching the binodal curve and changes monotonically in the middle region of the compositions. A noticeable effect on diffusion is exerted by the decomposition of the PVME·nH2O complex in the middle region of the compositions-the diffusion of the components here slightly increases. Based on the temperature dependences, the diffusion activation energy was estimated, the numerical values of which correspond to 5.8-6.5 kJ/mol.

Conclusions
We have shown that the phase equilibrium in aqueous solutions of PVME is a complex projection of a multi-component system that coexists at different temperatures. Thus, in the temperature range below 31 • C, the phase diagram corresponds to amorphous separation in the PVME-PVME·nH 2 O system, while above 31 • C, the PVME·nH 2 O complex is destroyed, and the phase diagram corresponds to the binary PVME-water system. At the same time, the intermediate regions of mutual solubility are obviously associated with the inertial effects of dehydration of functional groups of PVME and the redistribution of unbound water molecules. The kinetics of the formation and dissolution of phase particles is consistent with the Oswald mechanism, and the phase size in the diluted solution region can reach 400 µm. In this case, a secondary phase decay is observed within such macrophases.
Sorption studies have shown that in the concentrated solution region, the number of water molecules in a cluster formed around the functional group of PVME sharply increases. According to our calculations, up to 6-7 water molecules are combined in such a cluster, which is in good agreement with the literature data.
For the first time for the PVME-water system, we obtained numerical values for the diffusion coefficient. The logD values (cm 2 /s) vary in the range from −7 to −5.5, with an increase in temperature from 25 to 35 • C. A good correlation is observed between the data obtained both from sorption measurements and from direct contact of polymers by the optical interferometry method. The analysis of the concentration dependences of the diffusion coefficients showed that the diffusion rate decreases sharply when approaching the binodal curve and changes monotonically in the middle region of the compositions. A noticeable effect on diffusion is exerted by the decomposition of the PVME·nH 2 O complex in the middle region of the compositions-the diffusion of the components here slightly