Quantitative Study on Solubility Parameters and Related Thermodynamic Parameters of PVA with Different Alcoholysis Degrees

In recent years, inverse gas chromatography (IGC) and molecular dynamics simulation methods have been used to characterize the solubility parameters and surface parameters of polymers, which can provide quantitative reference for the further study of the surface and interface compatibility of polymer components in the future. In this paper, the solubility parameters and surface parameters of two kinds of common alcoholysis, PVA88 and PVA99, are studied by using the IGC method. The accuracy of the solubility parameters obtained by the IGC experiment is verified by molecular dynamics simulation. On the basis of this, the influence of repeated units of polyvinyl alcohol (PVA) on solubility parameters is studied, so as to determine the appropriate chain length of the PVA for simulation verification calculation. The results show that the solubility parameters are not much different when the PVA chain length is 30 and above; the numerical trends of the solubility parameters of PVA88 and PVA99 at room temperature are the same as the results of molecular dynamics simulation; the dispersive surface energy γsd and the specific surface energy γssp are scattered with the temperature distribution and have a small dependence on temperature. On the whole, the surface energy of PVA99 with a higher alcoholysis degree is higher than that of PVA88 with a lower alcoholysis degree. The surface specific adsorption free energy (ΔGsp) indicates that both PVA88 and PVA99 are amphoteric meta-acid materials, and the acidity of PVA99 is stronger.


Introduction
Polyvinyl alcohol (PVA) is a kind of environmental friendly water-soluble polymer material [1] with good film-forming property, excellent performance, non-toxicity, and degradability [2][3][4]. It is a polymer obtained by the polymerization of vinyl acetate and hydrolysis of polyvinyl acetate [5,6], and it is widely used in the fields of food packaging, pharmaceutical film [7], biological tissue, and other composite materials [8][9][10][11]. As the application of PVA becomes more and more extensive, its good performance makes it a hot spot in the field of materials research.
The properties of PVA are mainly affected by its average polymerization degree and alcoholysis degree [12][13][14][15]. The properties of PVA are different because of the different degree of polymerization and hydrolysis (alcoholysis) of PVA monomer. There are a lot of hydroxyl groups in complete alcoholysis PVA [16], and there are residual acetate groups in incomplete alcoholysis PVA, which has a great influence on the properties of PVA [17]. Zhao Caixia et al. [18]. Studied the effect of alcoholysis degree on the formation and properties of polyvinyl alcohol hydrogel. The results showed that the temperature sensitivity and dynamic rheological properties of PVA hydrogels with different alcoholysis degrees were also different. Huang Xiuling et al. [19] prepared PVA and thymol slow-release composite films, and they found that different alcoholysis degrees of PVA have a great influence on thymol release.
The research on PVA as the matrix of composite materials not only studies the influence of PVA structure [20,21] but also studies the interface compatibility of PVA composites. The interface compatibility directly affects the functionality of the composites [22]; for example, Ji Jianchao et al. [23] studied the interface characteristics of polyvinyl alcohol (PVA) film and carbon powder for printing. By studying the influence of different alcoholysis degrees of PVA on the interface interaction of materials, it is found that the more appropriate PVA type makes the film and carbon powder have better compatibility and results in a better printing effect. Jawalkar et al. [24] believe that the mixed ratio of poly-L-lysine (PLL)/PVA has a great influence on the compatibility or incompatibility trend of composite materials, which leads to some limitations in the application of biotechnology. Therefore, they use molecular dynamics simulation to calculate the internal energy density of PLL/PVA with different mixing ratios and Flory-Huggins parameters, and they explain the relationship between the interface compatibility and molecular structure. Hailong Xu et al. [25] prepared the porous f-Ti 2 CT x /PVA foam by the freeze-drying process, and it possesses excellent EM absorption ability. This lightweight absorption-dominated electromagnetic interference (EMI) shielding materials minimize the twice pollution of the reflected electromagnetic (EM) wave and have a broad application prospect. Tatiana et al. [26] prepared composite films containing polyvinyl alcohol with different amounts of graphene oxide and modified the composite films with glycerin. The ductile properties of graphene oxide-filled polyvinyl alcohol composite films have increased significantly by using glycerin. However, all of the above studies are limited to one-sided macroscopic research and computer molecular scale research, and few researchers have provided a quantitative discussion on the compatibility parameters of PVA with different alcoholysis degrees.
In this paper, the solubility parameters and interaction parameters (important parameters to measure the interfacial compatibility) of two kinds of commonly used PVA with different alcoholysis degrees were studied, and we calculated the data of part of the surface energy to transparentize the physical parameters such as the surface of the PVA. The solubility parameters, thermodynamic parameters, and some surface properties of PVA were calculated indirectly by the inverse gas chromatography (IGC) method and computer simulation. At the same time, the accuracy of solubility parameters obtained by the IGC method is calculated and verified by using the Forcite module of computer simulation software (Materials Studio 8.0, San Diego, CA, USA). Materials Studio is a material calculation software that can help researchers construct and analyze structural models of materials and predict the relevance properties of materials.
This paper provides a quantitative reference for future research on the relationship between the surface characteristics and structure of PVA composites by using the IGC method and computer simulation to calculate the relevant parameters of PVA with different alcoholysis degrees. It can also help other researchers choose the appropriate type of PVA when studying the PVA composites.

Raw Materials and Reagents
The experimental drugs and instruments used in this study are listed in Table 1. Preparation of PVA Powder The particles of PVA2488 and PVA2499 were crushed by a pulverizer, and the powder of more than 180 mesh was selected and put into a wide-mouth bottle for standby. The difference of PVA2488 and PVA2499 is that they have different alcoholysis degrees; the alcoholysis degree of PVA2488 is 88, while that of PVA2499 is 99.

IGC Determination of PVA Solubility Parameters and Surface Properties
PVA2488 was dissolved in acetone, and then, 6201 pickled red vectors were added to the PVA solution (volume of the discs was 1-1.2 times that of acetone) at a ratio of 1:10 (w/w)). The dry mixture was filled in a stainless-steel column having a length of 1.0 m and an inner diameter of 1/8 inch. The column was connected to an Agilent 6890 N gas chromatogram, and the column was aged for 24 h with N 2 as the protection gas at 180 • C column temperature. After the aging, an IGC test could be carried out. A Hamilton syringe was used to inject 0.4 µL of probe solvent by hand, and three parallel injections were performed. The retention time was recorded. The solubility parameters and surface parameters of PVA were calculated according to the IGC theoretical formula. Test temperature: 110 • C, 120 • C, 130 • C, 140 • C, 150 • C.

Molecular Dynamics Simulation to Calculate the Solubility Parameters of PVA
PVA is a kind of copolymer obtained by the free radical polymerization of vinyl acetate, which is formed by the alcoholysis of the acetate group located along the PVAc chain. Since the reaction process is incomplete, it usually has two unit structures, as shown in Figure 1. According to the percentage of hydroxyl group in the product after alcoholysis, it is divided into partial alcoholysis (alcoholysis degree 80.0-98.5%), high alcoholysis (alcoholysis degree > 98.5%), and complete alcoholysis (alcoholysis degree 100.0%) [27,28].

Raw Materials and Reagents
The experimental drugs and instruments used in this study are listed in Table 1. The particles of PVA2488 and PVA2499 were crushed by a pulverizer, and the powder of more than 180 mesh was selected and put into a wide-mouth bottle for standby. The difference of PVA2488 and PVA2499 is that they have different alcoholysis degrees; the alcoholysis degree of PVA2488 is 88, while that of PVA2499 is 99.

IGC Determination of PVA Solubility Parameters and Surface Properties
PVA2488 was dissolved in acetone, and then, 6201 pickled red vectors were added to the PVA solution (volume of the discs was 1-1.2 times that of acetone) at a ratio of 1:10 (w/w)). The dry mixture was filled in a stainless-steel column having a length of 1.0 m and an inner diameter of 1/8 inch. The column was connected to an Agilent 6890 N gas chromatogram, and the column was aged for 24 h with N2 as the protection gas at 180 °C column temperature. After the aging, an IGC test could be carried out. A Hamilton syringe was used to inject 0.4 μL of probe solvent by hand, and three parallel injections were performed. The retention time was recorded. The solubility parameters and surface parameters of PVA were calculated according to the IGC theoretical formula. Test temperature: 110 °C, 120 °C, 130 °C, 140 °C, 150 °C.

Molecular Dynamics Simulation to Calculate the Solubility Parameters of PVA
PVA is a kind of copolymer obtained by the free radical polymerization of vinyl acetate, which is formed by the alcoholysis of the acetate group located along the PVAc chain. Since the reaction process is incomplete, it usually has two unit structures, as shown in Figure 1. According to the percentage of hydroxyl group in the product after alcoholysis, it is divided into partial alcoholysis (alcoholysis degree 80.0-98.5%), high alcoholysis (alcoholysis degree > 98.5%), and complete alcoholysis (alcoholysis degree 100.0%) [27,28].

Determine the PVA Chain Length
A PVA chain was generated by the visualizer tool of MS software, and a PVA repeat unit was constructed by the homopolymer tool. As a result of the limited computing resources, it is very important to determine the appropriate PVA chain length. Based on the experience of many related researchers, this paper establishes a fully alcoholysis PVA long chain with repetition units of 10, 20, 30, 40, 50, and 60; then, it optimizes its structure by the Forcite module. The optimized PVA chain was used to establish its amorphous structure through amorphous cells (ACs). Finally, the solubility parameters of PVA with different chain lengths were calculated by the Forcite module, so as to determine the appropriate PVA chain length for further molecular dynamics research.

Calculation of Solubility Parameters
Since PVA with an alcoholysis degree of 99 is close to complete alcoholysis, for the convenience of modeling, we use complete alcoholysis PVA (PVA 100 ) instead of PVA with an alcoholysis degree of 99 to study. PVA with an alcoholysis degree of 88 (PVA 88 ) was established by changing the occupancy ratio of two kinds of repeat units, and the single-chain models of PVA 88 and PVA 100 were iterated 5000 times by the steepest descent method and conjugate gradient method respectively to obtain the lowest energy. Then, the single-chain structure was annealed several times by the annealing program in the Forcite module, which was heated from 290 to 590 K at 50 K intervals and cooled back at 20 K intervals until the energy was no longer changed. Through this process, the energy balance of the system was accelerated. Then, the optimized PVA chain was used to establish its amorphous structure through amorphous cells (ACs). In order to avoid the ring-pricking phenomenon of a PVA long chain, the lower initial density of 0.6 g/cm −3 was set, and the smart algorithm was used to optimize its geometry for 10,000 iterations to obtain the equilibrium system. The AC amorphous equilibrium system was simulated by molecular dynamics (MD) [29][30][31]. First, the NVT ensemble was used to calculate the equilibrium of 100 ps, relax the structure, and release the tension of the system to achieve the equilibrium of geometric conformation; then, 400 ps was balanced under the NPT ensemble, and then, 100ps was run through NVT, and the trajectory file was analyzed to calculate the solubility parameters of two alcoholysis degrees at 298.15 K. (The force field selection for all molecular dynamics simulation processes: COMPASSII Version 1.2, temperature control method: nose, pressure control method: Andersen, step length 1fs, accuracy: 0.001 kcal/mol, cutoff distance: 12.5 A.)

Characteristic Retention Volume of Probe Solvent
Using the IGC method [32], the probe solvent retention volume V 0 g can be obtained by Formulas (1) and (2). The specific values are shown in Table 2. Among them, J is the non-ideal gas compression factor, F is the carrier velocity, ∆t is the difference between the retention time t r of the probe solvent and n-pentane t m , m is the mass of the fixed-phase PVA, and T is the column temperature.
P i and P o are the inlet and outlet pressure of the stainless-steel column.
The retention volumes of probe solvents on PVA2488 and PVA2499 at different temperatures are listed in Table 2. On the whole, the increase of temperature destroys and weakens the interaction between the probe solvent and PVA, so the interaction time between the probe solvent and PVA decreases, and the retention volume of the probe solvent and PVA decreases; in homologues, the larger the molecular weight of the probe solvent is, the stronger the interaction between the probe solvent and PVA, and the larger the retention volume. For different kinds of probe solvents, the retention time of some probe solvents (ethanol, 1-propanol, tetrahydrofuran, ethyl acetate, acetonitrile) on PVA was longer than other probe solvents at the same temperature. The larger retention volume of ethanol and 1-propanol may be due to the formation of intermolecular hydrogen bonds between small alcohols containing more hydroxyl groups and PVA; tetrahydrofuran (THF), ethyl acetate, and acetonitrile are the most polar solvents in common mobile phases, which may interact with PVA with strong polarity, resulting in the increase of retention volume.
Comparing the retention volumes of PVA2488 and PVA2499, it can be seen that the retention volumes of PVA2488 and alkane probe solvents are larger than those of PVA2499, but the retention volumes of PVA2499 and non-alkane probe solvents are larger than those of PVA2488, which may be due to the fact that PVA2499 has more hydroxyl groups and weaker interaction with alkane non-polar probe solvents than other solvents.

Thermodynamic Parameters of Probe Molecule
From the V 0 g of the probe solvent, the molar absorption enthalpy ∆H s l , the infinite dilution molar mixing enthalpy ∆H ∞ l , and the molar evaporation enthalpy ∆H v can be calculated [33,34]. The calculation Formulas (3)-(5) are as follows: where Ω ∞ 1 is the infinite dilution activity coefficient [35], which is obtained by Equation (6): where R is the universal gas constant, M 1 is the molar mass of the probe solvent, and P 0 1 , B 11 , and V 1 are the saturated vapor pressure at the temperature of the probe solvent, the second virial coefficient, and the molar volume, respectively. The specific values are shown in Tables 3 and 4. The values of ∆H s l , ∆H ∞ l , and ∆H v of PVA2488 and PVA2499 are presented in Table 3. The thermodynamic parameters were compared between PVA2488 and PVA2499 in the same kind of probes, and there was no significant rule. However, in the homologues, the molar evaporation enthalpy ∆H v values increased with the increasing carbon number, which indicated that the higher energy was required for the evaporation of the probe solvent from PVA, which was associated with stronger PVA adsorption.

12
The Flory-Huggins interaction parameter χ ∞ 12 reflects the change of interaction energy when PVA is mixed with the probe solvent [36], and the smaller the χ ∞ 12 value, the stronger the dissolving ability of the probe solvent to PVA, which can be calculated by Equation (7), and the results are shown in Table 4. Table 4 shows that the Flory-Huggins interaction parameter χ ∞ 12 between the probe solvent and PVA at 383.15-423.15 K is less than 0.5. The results show that these probes are good solvents for PVA2488. In PVA2499, the χ ∞ 12 values of butanone, p-xylene, o-xylene, ethylbenzene, ether, tetrahydrofuran, ethyl acetate, and acetonitrile are all less than 0.5, indicating that these probes are good solvents for PVA2499, and they provide comparable solvents. The lower χ ∞ 12 value of PVA2499 indicates that PVA2499 is more soluble in these solvents.

Solubility Parameters
The solubility parameter δ 1 of the probe solvent is expressed by Formula [37] (8): The solubility parameters δ 2 of PVA can be obtained by using δ 1 of the probe solvent and Formula (9).
where χ ∞ S is the entropy term of the Flory-Huggins interaction parameter. The calculated results of solubility parameter δ 2 of PVA at 383.15-423.15 K are shown in Table 5. These δ 2 values can be obtained by plotting δ 2 1 /RT − χ ∞ 12 V 1 and δ 1 (as shown in Figure 2), and the solubility parameter of PVA at 298.15 K at room temperature can be obtained by extrapolation (as shown in Figure 3).
where χ is the entropy term of the Flory-Huggins interaction parameter. The calculated results of solubility parameter δ of PVA at 383.15-423.15 K are shown in Table 5. These δ values can be obtained by plotting δ /RT − and δ (as shown in Figure 2), and the solubility parameter of PVA at 298.15 K at room temperature can be obtained by extrapolation (as shown in Figure 3).      Since the IGC method needs to vaporize the solvent, the solubility parameter at 298.15 K cannot be measured, so extrapolation is required. By measuring the solubility parameters at 383.15 K, 393.15 K, 403.15 K, 413.15 K, and 423.15 K, using the origin software for linear fitting, the linear relationship between the solubility parameters of PVA2488 and PVA2499 and the temperature was obtained, and the solubility parameter at 298.15 K was obtained.
It can be seen from Table 5 that the solubility parameters of PVA2488 and PVA2499 at room temperature are 23.15 and 25.62 (J/cm −3 ) 0.5 , respectively. PVA2499 has a higher solubility parameter, which is probably due to the fact that PVA2499 has more -OH Since the IGC method needs to vaporize the solvent, the solubility parameter at 298.15 K cannot be measured, so extrapolation is required. By measuring the solubility parameters at 383.15 K, 393.15 K, 403.15 K, 413.15 K, and 423.15 K, using the origin software for linear fitting, the linear relationship between the solubility parameters of PVA2488 and PVA2499 and the temperature was obtained, and the solubility parameter at 298.15 K was obtained.
It can be seen from Table 5 that the solubility parameters of PVA2488 and PVA2499 at room temperature are 23.15 and 25.62 (J/cm −3 ) 0.5 , respectively. PVA2499 has a higher solubility parameter, which is probably due to the fact that PVA2499 has more -OH groups, and compared with -COOCH 3 groups, -OH has a greater contribution to the solubility parameter delta, and PVA2499 has more hydrophilic hydroxyl groups, forcing a large number of hydrogen bonds between molecules, which makes the cohesive force of the polymer chain stronger, thus making the solubility parameter of PVA2499 larger.

Calculation of Solubility Parameters by Molecular Dynamics Simulation
The trend diagram of the PVA repeat unit numbers and solubility parameters is shown in Figure 4a. It can be seen from Figure 4 and Table 6 that the solubility parameter difference of PVA with a repeat unit number of 30 or more chain length is small, so the PVA chain with a repeat unit number of 30 or more can meet the requirements of calculation accuracy. As a result of the long chain length, the calculation speed is slow. In order to improve the calculation efficiency and make the repeat unit ratio of the PVA 88 model reasonable, the PVA chain with 50 repeat units is selected for molecular dynamics simulation. The random stereoscopic model of PVA 100 with a repeat unit n = 50 and the random stereoscopic model of PVA 88 with a repeat unit ratio m:n = 6:44 (see Figure 4b) were established by the copolymer tool. The amorphous structure system (see Figure 4c) was established by the model, and then, the molecular dynamics simulation calculation was carried out. The molecular dynamics calculation process is the same as that in Section 2.3.2. From the final NVT output stage, we can evaluate whether the system reaches equilibrium. The fluctuation of temperature and energy distribution in a small range can be considered as an equilibrium system (see Figure 4d), which can be used to calculate the solubility parameters. From the figure, we can see that the fluctuation of energy and temperature is very small, which indicates that the system has reached equilibrium.
was established by the model, and then, the molecular dynamics simulation calculation was carried out. The molecular dynamics calculation process is the same as that in Section 2.3.2. From the final NVT output stage, we can evaluate whether the system reaches equilibrium. The fluctuation of temperature and energy distribution in a small range can be considered as an equilibrium system (see Figure 4d), which can be used to calculate the solubility parameters. From the figure, we can see that the fluctuation of energy and temperature is very small, which indicates that the system has reached equilibrium.     Table 7 that the value of PVA100 is larger than that of PVA88, and the van der Waals effect has a greater contribution to the total solubility parameter of PVA100, which may be due to the relatively large molecular weight of PVA88, which makes the molecular volume larger and the value of PVA88 smaller.

Comparison of Solubility Parameters between IGC Method and Molecular Dynamics Simulation of PVA
The solubility parameters of PVA with different alcoholysis degrees were studied by the IGC method and molecular dynamics method. The solubility parameters of PVA2488 and PVA2499 measured by the IGC method were 23.74 and 25.62 (J/cm −3 ) 0.5 , respectively. The solubility parameters of PVA 88 and PVA 100 calculated by molecular dynamics simulation were 22.68 and 24.36 (J/cm −3 ) 0.5 , respectively.
The solubility parameters of PVA with different alcoholysis degrees calculated by the two methods are significantly similar; comparing the data measured by the IGC method and the PVA molecular dynamics method can prove the reliability of the IGC method, and it can prove that the PVA molecular dynamics method can also be used to calculate and verify the experimental results. Therefore, the PVA molecular dynamics method can also be used in other studies to verify or directly calculate the corresponding parameters, which is very helpful for some data that are difficult to measure.

Surface Properties
The surface energy γ s refers to the extra energy of the solid surface relative to the interior [38,39]. In this part, the surface energy of PVA with different alcoholysis degrees was calculated by IGC technology at the temperature of 383.15-423.15 K. The specific values are listed in Table 8. The surface energy dispersive component γ d s represents the interaction force between the surface of solid materials and non-polar molecules [40], which can be calculated according to the Dorris-Gray method [41]. See Formulas (10) and (11): γ CH 2 = 36.8 − 0.058t (11) where n is the Avogadro constant; ∆G CH 2 is the increment of adsorption free energy, which is the slope of the linear relationship between RT ln V 0 g and the number of carbon atoms, as shown in Figure 5; γ CH 2 is the solid surface energy composed of a methylene group, and a CH 2 is the cross-sectional area of the methylene group, which is 0.06 nm 2 .    15-423.15 K, the more specific the value that can be obtained from Table 8, and in conjunction with Figure 6, it can be seen directly that the dispersion component γ d s does not change significantly with the increase in temperature, indicating that the dispersion effect between the nonpolar probe and the PVA surface under the temperature gradient of 383.15-423.15 K is less affected by temperature, the γ d s value is relatively close and less sensitive to temperature. This is most likely due to the too high temperature, which makes a certain degree of deformation on the PVA surface, so that the surface dispersion force is not only affected by the nonpolar functional groups. It may also have some relationship with the surface roughness of PVA, but due to the complexity of the situation, a more regular γ d s cannot be obtained, so the dispersion component γ d s of PVA at room temperature cannot be obtained by extrapolation. Comparing the γ d s of PVA2488 and PVA2499, PVA2499 has a higher γ d s at the same temperature.  It can be seen from Figure 6 that the specific component of PVA and the total surface energy and dispersion component do not change much with temperature, but in general, the specific component and total surface energy of PVA2499 are larger than those of PVA2488. It may be due to the fact that PVA2499 has more exposed active hydroxyl groups than PVA2488 under the same degree of polymerization, which makes it have higher surface energy. It is likely that PVA2488, which has a lower degree of alcoholysis, contains more residual acetate and is a hydrophobic group, which reduces the tension of the probe molecule on its surface, so it has a lower surface energy than PVA2499. The surface-specific component γ sp s reflects the surface energy of the dispersive part between the filler and the probe, which includes polarity, acid-base, and hydrogen bonding. It can be calculated by the Good-van Oss method [42] combined with the specific adsorption free energy ∆G sp . The formula is as follows: The γ + l and γ − l [35] of dichloromethane are 124.58 and 0 mJ/m 2 , respectively, and the γ + l and γ − l of toluene are 16.23 and 0 mJ/m 2 , respectively. From this, the specific component and total surface energy of the filler can be calculated. The specific component γ sp s of the PVA in addition to the dispersion part is calculated using dichloromethane and toluene. The calculation results are shown in Table 8. In addition, the total surface energy of the PVA can be obtained by Formula (14).
It can be seen from Figure 6 that the specific component of PVA and the total surface energy and dispersion component do not change much with temperature, but in general, the specific component and total surface energy of PVA2499 are larger than those of PVA2488. It may be due to the fact that PVA2499 has more exposed active hydroxyl groups than PVA2488 under the same degree of polymerization, which makes it have higher surface energy. It is likely that PVA2488, which has a lower degree of alcoholysis, contains more residual acetate and is a hydrophobic group, which reduces the tension of the probe molecule on its surface, so it has a lower surface energy than PVA2499.

Surface Acidity and Alkalinity
In addition to the dispersion effect on the polar probe molecule, the surface of PVA also has an acid-base effect. According to the Sawyer method [43,44], Formula (15) can be used to obtain the surface-specific adsorption free energy ∆G sp of PVA at different temperatures. Figure 7 shows the calculation diagram of the acidic probe chloroform, the basic probe tetrahydrofuran, and the amphoteric probe ethyl acetate on the surface ∆G sp of the PVA calculation diagram; the straight line in the figure is the relationship between the RT ln V 0 g of n-alkane and the boiling point t b of the probe solvent. This straight line is a reference, and the distance between the unipolar probe and the amphipathic probe that deviates from this straight line is the ∆G sp of the probe on the PVA surface.
Polymers 2021, 13, x FOR PEER REVIEW 13 of 16  Table 9 lists the characteristic adsorption free energies ∆G of different probes on the surface of PVA2488 and PVA2499. Overall, tetrahydrofuran, trichloromethane, as well as ethyl acetate are all located above the straight line, illustrating the existence of both acidic and basic sites on the surface of PVA, which are amphoteric materials. Since the ∆G value of basic probe tetrahydrofuran on the surface of PVA is larger relative to the acidic probe trichloromethane, this indicates that the PVA surface has a stronger acidic effect; that is, PVA is a zwitterionic acid material; from the table, the characteristic adsorption free energy of the two PVAs can be compared. At the same temperature, the ∆G of trichloromethane on the surface of PVA2488 is larger and that of tetrahydrofuran on PVA2499 is larger, which also indicates that PVA2499 is relatively more acidic with respect to the PVA2488 surface.   Table 9 lists the characteristic adsorption free energies ∆G sp of different probes on the surface of PVA2488 and PVA2499. Overall, tetrahydrofuran, trichloromethane, as well as ethyl acetate are all located above the straight line, illustrating the existence of both acidic and basic sites on the surface of PVA, which are amphoteric materials. Since the ∆G sp value of basic probe tetrahydrofuran on the surface of PVA is larger relative to the acidic probe trichloromethane, this indicates that the PVA surface has a stronger acidic effect; that is, PVA is a zwitterionic acid material; from the table, the characteristic adsorption free energy of the two PVAs can be compared. At the same temperature, the ∆G sp of trichloromethane on the surface of PVA2488 is larger and that of tetrahydrofuran on PVA2499 is larger, which also indicates that PVA2499 is relatively more acidic with respect to the PVA2488 surface.

Conclusions
In order to choose the suitable PVA for the subsequent composites, in this paper, the solubility parameters of PVA with two alcoholysis degrees were studied by the IGC method, aiming to provide a quantitative reference for the subsequent interfacial compatibility study, and a PVA model corresponding to two alcoholysis degrees was built by molecular dynamics computer simulation and solubility parameter calculation, which verified the experimental calculation by the IGC method, and the results showed that the two interrater calculations had consistency. At the same time, the surface energy and surface energy components (dispersive surface energy component γ d s , specific surface energy component γ sp s ) are calculated with the help of the IGC method. The results show that the dispersion of the nonpolar probe and the PVA surface is affected by temperature at a temperature above 110 • C. The dispersion component γ d s of PVA at room temperature cannot be obtained by extrapolation, but the dispersion component contributes a lot to the surface energy as a whole, and compared to PVA2499, the surface energy of PVA2488 is larger. In addition, the surface acid-base properties of PVA with two alcoholysis degrees were investigated at a certain temperature, and the results showed that PVA was an amphoteric acid material, and the surface with a larger alcoholysis degree (PVA2499) had strong acidity.
This paper, through the study of PVA's surface properties with solubility parameters, as well as interaction parameters, and so on, can provide a valuable reference for the relationship between its structure and properties of composites, so as to make the corresponding composites have better interfacial bonding properties, for example, improving the performance of PVA hydrogel and PVA composite film.
In addition, the MD method used in this article can not only calculate the parameters of the surface properties of PVA but also simulate and calculate other data; for example, many studies have shown that the glass transition temperature (TG) of the material can be derived from the temperature-volume diagram obtained by MD. In the following research, we will continue to study the molecular dynamics in depth and explore the simulation calculation of other parameters such as TG by this method.