The swelling of a poly (methyl methacrylate) in supercritical carbon dioxide was studied by means of full atomistic classical molecular dynamics simulation. In order to characterize the polymer swelling, we calculated various properties related to the density, structure, and dynamics of polymer chains as a function of the simulation time, temperature, and pressure. In addition, we compared the properties of the macromolecular chains in supercritical CO₂ with the properties of the corresponding bulk system at the same temperature and atmospheric pressure. It was shown that diffusion of CO₂ molecules into the polymer led to a significant increase in the chain mobility and distances between them. Analysis of diffusion coefficients of CO₂ molecules inside and outside the poly(methyl methacrylate) sample has shown that carbon dioxide actively interacts with the functional groups of poly (methyl methacrylate). Joint analysis of the radial distribution functions obtained from classical molecular dynamics and of the averaging interatomic distances from Car-Parrinello molecular dynamics allows us to make a conclusion about the possibility of formation of weak hydrogen bonds between the carbon dioxide oxygen atom and the hydrogen atoms of the polymer methyl groups.

Classical molecular dynamics simulations of polyacrylamide (PAM) adsorption on cellulose nanocrystals (CNC) in a vacuum and a water environment are carried out to interpret the mechanism of the polymer interactions with CNC. The structural behavior of PAM is studied in terms of the radius of gyration, atom–atom radial distribution functions, and number of hydrogen bonds. The structural and dynamical characteristics of the polymer adsorption are investigated. It is established that in water the polymer macromolecules are mainly adsorbed in the form of a coil onto the CNC facets. It is found out that water and PAM sorption on CNC is a competitive process, and water weakens the interaction between the polymer and CNC.

Models of interaction between a poly(N-vinyl-2-pyrrolidone) macromolecule and a fragment of I_β-cellulose were built in a vacuum and water environment. The models were made to interpret the mechanism of interaction of the polymer and cellulose nanocrystals by the classical molecular dynamics method. The structural behavior of a poly(N-vinyl-2-pyrrolidone) macromolecule in water has been studied in terms of the radius of gyration, atom–atom radial distribution functions and number of hydrogen bonds. It was found that the polymer has a high affinity with the solvent and each monomer unit has on average 0.5 hydrogen bonds. The structural and energy characteristics of the polymer adsorption were investigated at different initial positions of the poly(N-vinyl-2-pyrrolidone) macromolecule relative to the cellulose fragment. It was observed that the polymer macromolecule was mainly adsorbed on the cellulose fragment in the globular form. Moreover, in the solvent the interaction of poly(N-vinyl-2-pyrrolidone) with the cellulose hydrophobic surface was stronger than that with the hydrophilic one. This study will show that the presence of water makes the interaction between the polymer and cellulose weaker than in a vacuum, and the polymer and cellulose mainly interact through their solvation shells.

Poly(3-hydroxybutyrate-co-3-hydroxyvalerate)/polycaprolactone (PHBV/PCL) polymer mixtures reinforced by cellulose nanocrystals (CNCs) have been obtained. To improve the CNC compatibility with the hydrophobic PHBV/PCL matrix, the CNC surface was modified by amphiphilic polymers, i.e., polyvinylpyrrolidone (PVP) and polyacrylamide (PAM). The polymer composites were characterized by FTIR, DSC, TG, XRD, microscopy, BET surface area, and tensile testing. The morphological, sorption, thermal, and mechanical properties of the obtained composites have been studied. It was found out that with an increase in the CNC content in the composites, the porosity of the films increased, which was reflected in an increase in their specific surface areas and water sorption. An analysis of the IR spectra confirms that hydrogen bonds can be formed between the CNC hydroxyl- and the –CO– groups of PCL and PHBV. The thermal decomposition of CNC in the PHBV/PCL/CNC composites starts at a much higher temperature than the decomposition of pure CNC. It was revealed that CNCs can either induce crystallization and the polymer crystallite growth or act as a compatibilizer of a mixture of the polymers causing their amorphization. The CNC addition significantly reduces the elongation and strength of the composites, but changes Young’s modulus insignificantly, i.e., the mechanical properties of the composites are retained under conditions of small linear deformations. A molecular-dynamics simulation of several systems, starting from simplest binary (solvent-polymer) and finishing with multi-component (CNC—polymer mixture—solvent) systems, has been made. It is concluded that the surface modification of CNCs with amphiphilic polymers makes it possible to obtain the CNC composites with hydrophobic polymer matrices.