Molecular-Level Insight of CP52/NBR Damping Composites through a Combination of Molecular Dynamics Simulation and Experimental Method

To enhance the damping properties of nitrile butadiene rubber (NBR), the elastomer used was blended with chlorinated paraffin 52 (CP52) to prepare NBR/CP52 composites. The results showed that CP52 could significantly enhance the damping properties of NBR and shift the glass transition temperature (Tg) to lower temperatures. Molecular dynamics models of the CP52/NBR system were established, and the damping properties of the CP52-reinforced NBR were investigated using molecular dynamics (MD) simulations. Through the combination of MD simulations and the experimental results, the essential mechanism of the enhanced damping properties of the NBR was methodically expatiated and was ascribed to the Cl-CP-H····NC-NBR (type I) and CP-Cl····H-NBR-CN (type II) analogous hydrogen bonds formed between NBR and CP52. The higher the CP52 content, the higher the analogous hydrogen bond concentration, and the better the damping properties of the CP52/NBR composites. The experimental results were very consistent with the MD simulation results, meaning that the combination method can provide a new means to optimize the design of damping materials and broaden the application range of small polar molecules in the damping modification of polar rubber materials.


Introduction
Currently, vibration and noise problems are becoming increasingly prominent, and the application of damping materials is one of the most direct and effective methods for reducing vibration and noise [1]. Nitrile butadiene rubber (NBR) has excellent mechanical strength, damping properties, promising oil and heat resistance, and is widely used in damping materials [2,3]. However, due to the influence of the glass transition temperature (T g ), the single-component NBR's high damping temperature range is lower than room temperature, and the effective damping temperature range is generally narrow and cannot necessarily meet the needs of different engineering fields and application environments [4,5].
In general, the dynamic mechanical behavior of polymers can be adjusted by conventional means such as the additions of plasticizers and fillers, mechanical blending of various polymers, and copolymerization. However, the existing damping modification methods and research methods have certain limitations. For example, both blending and copolymerization can broaden the effective damping temperature range of rubber-based damping materials to a certain extent, but the blending materials have problems such as easy phase separation of different components and unstable mechanical properties. Copolymerization also has problems such as complicated preparation process, high batch production cost, and serious environmental pollution. NBR (4155) with an acrylonitrile mass fraction of 41% was provided by Nnatex Industry Co., Ltd. (Zhenjiang, China). CP52 was provided by the Damao Chemical Reagent Factory. (Tianjin, China). magnesium oxide (MgO), stearic acid (SA), DM (2,2dithiodibenzothiazole) and TMTD (tetramethyl thiuram disulfide) was provided by yuanhao Co., Ltd. (Hebi, China). All of the other raw materials are commercially available industrial products.

Preparation of CP52/NBR Composites
First, the NBR was kneaded with a mixing roller for 5 min, and CP52 was added to the kneaded NBR, which was then kneaded at room temperature for 5 min. Then, the rubber additives were added to the above samples, including 5 phr of MgO, 1 phr of SA, 1.5 phr of accelerator DM, 0.5 phr of accelerator TMTD, and 0.5 phr of sulfur. Evenly mixing for 10 min. Finally, the samples were hot-pressed and vulcanized for 15 min at 12 MPa and 150 • C and then cooled naturally to room temperature to obtain the specimens.

Characterization
Dynamic mechanical analysis (DMA, DMA1, METTLER TOLEDO, Zurich, Switzerland) was used to test the dynamic mechanical properties in the tensile mode. The test frequency was 10 Hz, the test temperature range was −60 • C~80 • C, and the heating rate was 3 • C/min. Fourier-transform infrared spectroscopy (FTIR, Nicolet IS20, Thermo Fisher Scientific, Madison, WI, USA) was used to test the samples pre-pressed into films using the ATR module and total reflection mode. The scanning wave number ranged from 400~4000 cm −1 .

Model and Simulation Details
The Materials Studio (MS) 2019 software was used to reveal the effect of CP52 on the NBR composites from the microscopic aspect, providing theoretical guidance for the experimental study of damping materials.
The MD simulation was carried out by using the Forcite and Amorphous Cell modules. During the simulation, an Andersen thermostat was used for temperature control, a Berendsen barostat controlled the pressure, and the condensed phase optimized molecular potentials for atomistic simulation studies (COMPASS) force field was adopted. Figure 1 shows the construction procedure for the CP52/NBR composite amorphous cells. The acrylonitrile repeat units ( Figure 1a) and butadiene repeat units ( Figure 1b) were first randomly copolymerized as an NBR polymer chain (Figure 1c) with an acrylonitrile mass fraction of 41%. Then, the NBR polymer chains and CP52 small molecules ( Figure 1d) were built in a periodic boundary cell ( Figure 1e); according to the mass ratios, each cell consisted of 4 NBR chains with 50 repeat units and different CP52 molecules and are listed in Table 1. Subsequently, the cell's energy was minimized by using the Smart Minimizer method to relax the state of the minimal potential energy. After the completion of energy minimization, in order to eliminate the internal stress, the cell was annealed from the low temperature of 200 K to the upper temperature of 400 K with 100 cycles to prevent the system from being trapped in a local high-energy minimum. Finally, the dynamics simulation was performed using an NVT ensemble (constant atomic number, volume, and temperature) at 298 K for 500 ps and 1000 ps of NPT (constant number of particles, pressure, and temperature). The simulation was carried out at 0.1 MPa to relax the polymer structure further. Then, the relevant physical parameters required for the system to reach perfect equilibrium were calculated. frequency was 10 Hz, the test temperature range was -60 °C~80 °C, and the heating rate was 3 °C/min. Fourier-transform infrared spectroscopy (FTIR, Nicolet IS20, Thermo Fisher Scientific, Madison, USA) was used to test the samples pre-pressed into films using the ATR module and total reflection mode. The scanning wave number ranged from 400~4000 cm −1 .

Model and Simulation Details
The Materials Studio (MS) 2019 software was used to reveal the effect of CP52 on the NBR composites from the microscopic aspect, providing theoretical guidance for the experimental study of damping materials.
The MD simulation was carried out by using the Forcite and Amorphous Cell modules. During the simulation, an Andersen thermostat was used for temperature control, a Berendsen barostat controlled the pressure, and the condensed phase optimized molecular potentials for atomistic simulation studies (COMPASS) force field was adopted. Figure 1 shows the construction procedure for the CP52/NBR composite amorphous cells. The acrylonitrile repeat units ( Figure 1a) and butadiene repeat units ( Figure 1b) were first randomly copolymerized as an NBR polymer chain (Figure 1c) with an acrylonitrile mass fraction of 41%. Then, the NBR polymer chains and CP52 small molecules ( Figure  1d) were built in a periodic boundary cell ( Figure 1e); according to the mass ratios, each cell consisted of 4 NBR chains with 50 repeat units and different CP52 molecules and are listed in Table 1. Subsequently, the cell's energy was minimized by using the Smart Minimizer method to relax the state of the minimal potential energy. After the completion of energy minimization, in order to eliminate the internal stress, the cell was annealed from the low temperature of 200 K to the upper temperature of 400 K with 100 cycles to prevent the system from being trapped in a local high-energy minimum. Finally, the dynamics simulation was performed using an NVT ensemble (constant atomic number, volume, and temperature) at 298 K for 500 ps and 1000 ps of NPT (constant number of particles, pressure, and temperature). The simulation was carried out at 0.1 MPa to relax the polymer structure further. Then, the relevant physical parameters required for the system to reach perfect equilibrium were calculated.   The reliability of the model was tested in relation to energy and density. When the energy fluctuation in the system was less than 10% of the total energy and the density fluctuation was within 3%, the system was considered stable. The simulated density of the system and the density of the real sample are listed in Table 2. Compared with the experimental data, the simulated density deviates from the density measured in the real sample by no more than 3%. It shows that the model is close to the actual situation and can reflect the real performance of the material. The relevant data calculated by this model have reference value and can be used to analyze the damping mechanisms of CP52/NBR materials. The effect of CP52 on the damping properties of NBR is shown in Table 3 and Figure 2. With the increase in CP52 content, the loss of the tangent (tanδ) and the peak area (TA, tanδ > 0.3) of the CP52/NBR composites increased, the effective damping temperature range (∆T, tanδ > 0.3) broadened, and T g shifted to a lower temperature. When the added amount of CP52 is 100 phr, the ∆T of the composite is 46.8 • C, and the tanδ max is as high as 1.78, which is 62% higher than NBR without CP52. It shows that CP52 can significantly enhance the damping properties of NBR.

FTIR Analysis of CP52/NBR Material
The position and intensity of the absorption peaks on the FTIR spectra can reflect the state of the functional groups inside the material and the intermolecular forces. Figure 3a is the FTIR spectra of NBR and CP52 and the CP52/NBR materials, wherein (1) and (2) are the characteristic absorption peaks of -CN and -Cl. It can be found that, with increased CP52 content, the peaks are red-shifted [11], and the stretching vibration peak of -CN gradually changes from 2236 cm −1 to 2234 cm −1 . The characteristic absorption peak of -Cl also moves to a lower frequency and continues to increase. The reason is that the strong induction effect leads to the averaging electron cloud distribution of the polar groups, causing the chemical bond between the polar groups and the main molecular chain to elongate, decreasing the wave number of the characteristic peak [12,13].

FTIR Analysis of CP52/NBR Material
The position and intensity of the absorption peaks on the FTIR spectra can reflect the state of the functional groups inside the material and the intermolecular forces. Figure 3a is the FTIR spectra of NBR and CP52 and the CP52/NBR materials, wherein (1) and (2) are the characteristic absorption peaks of -CN and -Cl. It can be found that, with increased CP52 content, the peaks are red-shifted [11], and the stretching vibration peak of -CN gradually changes from 2236 cm −1 to 2234 cm −1 . The characteristic absorption peak of -Cl also moves to a lower frequency and continues to increase. The reason is that the strong induction effect leads to the averaging electron cloud distribution of the polar groups, causing the chemical bond between the polar groups and the main molecular chain to elongate, decreasing the wave number of the characteristic peak [12,13].  Therefore, it can be speculated that, based on the crosslinking network established by vulcanization, CP52 and NBR molecular chains can also form a network structure through intermolecular interactions. As shown in Figure 4, the thick black line represents the NBR molecular chain segment, the thin blue line represents the CP52 molecule, and the Cl-CP-H····NC-NBR and CP-Cl····H-NBR-CN intermolecular interactions are represented by imaginary lines as type I and type II, respectively. The interaction between these two molecules ensures the formation of a spatial network structure but is not as stable as a covalent bond. When subjected to alternating forces, such as vibration, the intermolecular interactions are continuously dissociated and generated; the mechanical energy is rapidly dissipated; thus, the damping performance of the CP52/NBR material improves.

FTIR Analysis of CP52/NBR Material
The position and intensity of the absorption peaks on the FTIR spectra can reflect the state of the functional groups inside the material and the intermolecular forces. Figure 3a is the FTIR spectra of NBR and CP52 and the CP52/NBR materials, wherein (1) and (2) are the characteristic absorption peaks of -CN and -Cl. It can be found that, with increased CP52 content, the peaks are red-shifted [11], and the stretching vibration peak of -CN gradually changes from 2236 cm −1 to 2234 cm −1 . The characteristic absorption peak of -Cl also moves to a lower frequency and continues to increase. The reason is that the strong induction effect leads to the averaging electron cloud distribution of the polar groups, causing the chemical bond between the polar groups and the main molecular chain to elongate, decreasing the wave number of the characteristic peak [12,13].  Therefore, it can be speculated that, based on the crosslinking network established by vulcanization, CP52 and NBR molecular chains can also form a network structure through intermolecular interactions. As shown in Figure 4, the thick black line represents the NBR molecular chain segment, the thin blue line represents the CP52 molecule, and the Cl-CP-H····NC-NBR and CP-Cl····H-NBR-CN intermolecular interactions are represented by imaginary lines as type I and type II, respectively. The interaction between these two molecules ensures the formation of a spatial network structure but is not as stable as a covalent bond. When subjected to alternating forces, such as vibration, the intermolecular interactions are continuously dissociated and generated; the mechanical energy is rapidly dissipated; thus, the damping performance of the CP52/NBR material improves. Therefore, it can be speculated that, based on the crosslinking network established by vulcanization, CP52 and NBR molecular chains can also form a network structure through intermolecular interactions. As shown in Figure 4, the thick black line represents the NBR molecular chain segment, the thin blue line represents the CP52 molecule, and the Cl-CP-H····NC-NBR and CP-Cl····H-NBR-CN intermolecular interactions are represented by imaginary lines as type I and type II, respectively. The interaction between these two molecules ensures the formation of a spatial network structure but is not as stable as a covalent bond. When subjected to alternating forces, such as vibration, the intermolecular interactions are continuously dissociated and generated; the mechanical energy is rapidly dissipated; thus, the damping performance of the CP52/NBR material improves.

Radial Distribution Function Analysis
The radial distribution function (RDF) refers to the radial distribution of a particle centered with respect to other particles, which can reflect the degree of order within the polymer molecules and intermolecular interactions [14]. As Figure 5 shows, (a) is the RDF spectra between the nitrile group on NBR and the para-hydrogen of the chlorine atom on CP52, while (b) is the RDF spectra between the chlorine atom on CP52 and the parahydrogen of the nitrile group on NBR; corresponding to the type I and type II intermolecular forces proposed in the previous section, respectively. It can be seen that the RDF spectra of both forces have a peak near 3.0 Å, which proves that these two interactions exist in the CP52/NBR system, which is consistent with the FTIR test results.

Radial Distribution Function Analysis
The radial distribution function (RDF) refers to the radial distribution of a particle centered with respect to other particles, which can reflect the degree of order within the polymer molecules and intermolecular interactions [14]. As Figure 5 shows, (a) is the RDF spectra between the nitrile group on NBR and the para-hydrogen of the chlorine atom on CP52, while (b) is the RDF spectra between the chlorine atom on CP52 and the para-hydrogen of the nitrile group on NBR; corresponding to the type I and type II intermolecular forces proposed in the previous section, respectively. It can be seen that the RDF spectra of both forces have a peak near 3.0 Å, which proves that these two interactions exist in the CP52/NBR system, which is consistent with the FTIR test results. According to reports [15,16], different intermolecular interaction forces form different probability densities. The type of intermolecular interaction can be determined by analyzing the morphology and peak position of the RDF spectrum. Generally, hydrogen bonds are distributed in the range of 2.6 Å to 3.1 Å [17]; strong van der Waals forces are between 3.1 Å and 5.0 Å [18], and weak van der Waals forces are more than 5.0 Å. The peaks of Cl-CP-H····NC-NBR and CP-Cl····H-NBR-CN are wide, near 2.9Å and 2.8Å, and the peak's shape is influenced by some hydrogen bonds and strong van der Waals forces.

Radial Distribution Function Analysis
The radial distribution function (RDF) refers to the radial distribution of a particle centered with respect to other particles, which can reflect the degree of order within the polymer molecules and intermolecular interactions [14]. As Figure 5 shows, (a) is the RDF spectra between the nitrile group on NBR and the para-hydrogen of the chlorine atom on CP52, while (b) is the RDF spectra between the chlorine atom on CP52 and the para-hydrogen of the nitrile group on NBR; corresponding to the type I and type II intermolecular forces proposed in the previous section, respectively. It can be seen that the RDF spectra of both forces have a peak near 3.0 Å, which proves that these two interactions exist in the CP52/NBR system, which is consistent with the FTIR test results. According to reports [15,16], different intermolecular interaction forces form different probability densities. The type of intermolecular interaction can be determined by analyzing the morphology and peak position of the RDF spectrum. Generally, hydrogen bonds are distributed in the range of 2.6 Å to 3.1 Å [17]; strong van der Waals forces are between 3.1 Å and 5.0 Å [18], and weak van der Waals forces are more than 5.0 Å. The peaks of Cl-CP-H····NC-NBR and CP-Cl····H-NBR-CN are wide, near 2.9Å and 2.8Å, and the peak's shape is influenced by some hydrogen bonds and strong van der Waals forces. According to reports [15,16], different intermolecular interaction forces form different probability densities. The type of intermolecular interaction can be determined by analyzing the morphology and peak position of the RDF spectrum. Generally, hydrogen bonds are distributed in the range of 2.6 Å to 3.1 Å [17]; strong van der Waals forces are between 3.1 Å and 5.0 Å [18], and weak van der Waals forces are more than 5.0 Å. The peaks of Cl-CP-H····NC-NBR and CP-Cl····H-NBR-CN are wide, near 2.9 Å and 2.8 Å, and the peak's shape is influenced by some hydrogen bonds and strong van der Waals forces. The strength of the interaction is close to that of hydrogen bonds, and the bonding mechanism is similar to that of hydrogen bonds. Therefore, their intermolecular interaction is defined as an analogous hydrogen bond, where Cl-CP-H····NC-NBR is a type I analogous hydrogen bond, and CP-Cl····H-NBR-CN is a type II analogous hydrogen bond. Their structures are shown in Figure 6.
The strength of the interaction is close to that of hydrogen bonds, and the bonding mechanism is similar to that of hydrogen bonds. Therefore, their intermolecular interaction is defined as an analogous hydrogen bond, where Cl-CP-H····NC-NBR is a type I analogous hydrogen bond, and CP-Cl····H-NBR-CN is a type II analogous hydrogen bond. Their structures are shown in Figure 6. The number of the two types of analogous hydrogen bonds in the periodic cells and the molar concentrations of the analogous hydrogen bonds in the last 50 conformations of each system were calculated. As Figure 7 shows, it can be seen that, with the increased CP52 content, the mass fraction of NBR decreased, the concentration of the type I analogous hydrogen bonds in the system increased, and the type II analogous hydrogen bonds tended to decrease. However, the total concentration of the analogous hydrogen bonds (CI+II) increased. It was found that the variation trend of CI+II is consistent with the measured tanδmax and TA. To further determine the relationship between the molar concentration of the analogous hydrogen bonds and the damping property of the material, the above data were analyzed using the linear regression method. The regression curve is shown in Figure 8. The concentration of the analogous hydrogen bonds in the CP52/NBR materials is closely related to the damping properties. The correlation coefficient between CI+II and tanδmax is 0.891, and the correlation coefficient with TA is 0.951. The number of the two types of analogous hydrogen bonds in the periodic cells and the molar concentrations of the analogous hydrogen bonds in the last 50 conformations of each system were calculated. As Figure 7 shows, it can be seen that, with the increased CP52 content, the mass fraction of NBR decreased, the concentration of the type I analogous hydrogen bonds in the system increased, and the type II analogous hydrogen bonds tended to decrease. However, the total concentration of the analogous hydrogen bonds (C I+II ) increased. The number of the two types of analogous hydrogen bonds in the per the molar concentrations of the analogous hydrogen bonds in the last 50 con each system were calculated. As Figure 7 shows, it can be seen that, with CP52 content, the mass fraction of NBR decreased, the concentration of the gous hydrogen bonds in the system increased, and the type II analogous hy tended to decrease. However, the total concentration of the analogous hy (CI+II) increased. It was found that the variation trend of CI+II is consistent with the me and TA. To further determine the relationship between the molar concen analogous hydrogen bonds and the damping property of the material, th were analyzed using the linear regression method. The regression curve is ure 8. The concentration of the analogous hydrogen bonds in the CP52/NB closely related to the damping properties. The correlation coefficient betw tanδmax is 0.891, and the correlation coefficient with TA is 0.951. It was found that the variation trend of C I+II is consistent with the measured tanδ max and TA. To further determine the relationship between the molar concentration of the analogous hydrogen bonds and the damping property of the material, the above data were analyzed using the linear regression method. The regression curve is shown in Figure 8. The concentration of the analogous hydrogen bonds in the CP52/NBR materials is closely related to the damping properties. The correlation coefficient between C I+II and tanδ max is 0.891, and the correlation coefficient with TA is 0.951. It is inferred that the analogous hydrogen bond is the main factor affecting the damping performance of the CP52/NBR material. The plasticity of CP52 makes the NBR molecular chain segment easier to move under the action of external forces. The chain segment is continuously dissociated and generated in the movement process, and the external energy is continuously consumed, resulting in a good damping effect. The higher the CP52 content, the higher the total concentration of analogous hydrogen bonds, so the better the damping performance.

Binding Energy Analysis
When the system contains two or more components, the components are combined through intermolecular interactions, and the energy consumed to separate the components is the binding energy (Ebinding). Ebinding is defined as the negative value of intermolecular interaction energy (Ebinding = −Einter = −(Etotal − ECP52 − ENBR)) [19]. Ebinding can characterize the strength of the intermolecular interaction. The larger Ebinding is, the stronger the interaction between the components in the system and the better the compatibility. The energyrelated parameters of the CP52/NBR system are listed in Table 4. All of the Ebinding values are positive and enlarged gradually with the increased CP52 content, indicating good compatibility between CP52 and NBR and a more potent intermolecular force; this is because the analogous hydrogen bonds are formed between CP52 and NBR. The higher the CP52 content, the higher the analogous hydrogen bond content in the system, Ebinding is more prominent, and more external energy is consumed in the dissociation process. Therefore, it can be determined that the material has better damping performance.

Fraction Free Volume and Cohesive Energy Density Analysis
The fraction free volume (FFV) represents the percentage of free volume in the system's total volume, which can reflect the tightness of the molecular accumulation in the system. Cohesive energy refers to the total energy that 1 mol molecules gather or the energy consumed to overcome intermolecular forces and transfer a molar of aggregated matter outside the range of the intermolecular forces [20]. Cohesive energy per unit volume It is inferred that the analogous hydrogen bond is the main factor affecting the damping performance of the CP52/NBR material. The plasticity of CP52 makes the NBR molecular chain segment easier to move under the action of external forces. The chain segment is continuously dissociated and generated in the movement process, and the external energy is continuously consumed, resulting in a good damping effect. The higher the CP52 content, the higher the total concentration of analogous hydrogen bonds, so the better the damping performance.

Binding Energy Analysis
When the system contains two or more components, the components are combined through intermolecular interactions, and the energy consumed to separate the components is the binding energy (E binding ). E binding is defined as the negative value of intermolecular interaction energy (E binding = −E inter = −(E total − E CP52 − E NBR )) [19]. E binding can characterize the strength of the intermolecular interaction. The larger E binding is, the stronger the interaction between the components in the system and the better the compatibility. The energy-related parameters of the CP52/NBR system are listed in Table 4. All of the E binding values are positive and enlarged gradually with the increased CP52 content, indicating good compatibility between CP52 and NBR and a more potent intermolecular force; this is because the analogous hydrogen bonds are formed between CP52 and NBR. The higher the CP52 content, the higher the analogous hydrogen bond content in the system, E binding is more prominent, and more external energy is consumed in the dissociation process. Therefore, it can be determined that the material has better damping performance.

Fraction Free Volume and Cohesive Energy Density Analysis
The fraction free volume (FFV) represents the percentage of free volume in the system's total volume, which can reflect the tightness of the molecular accumulation in the system. Cohesive energy refers to the total energy that 1 mol molecules gather or the energy consumed to overcome intermolecular forces and transfer a molar of aggregated matter outside the range of the intermolecular forces [20]. Cohesive energy per unit volume is called cohesive energy density (CED) [21]. Generally, the changing trend of FFV and CED is the opposite [22]. The higher CED is, the stronger the intermolecular force, the closer the molecular stacking in the system, the smaller the FFV, and the lower the T g of the material [23,24]. The calculation results of FFV and CED of the CP52/NBR system are shown in Figure 9. It can be seen that the changes in FFV and CED are not noticeable. With the increased CP52 content, FFV is always stable at around 38%, while CED fluctuates in the range of 330~380 J/cm −3 . The changing trend regarding FFV and CED is roughly the opposite. Therefore, FFV is not the main factor affecting the damping property and the T g of the CP52/NBR materials. Increasing the CP52 addition does not significantly change FFV. It is mainly due to the CP52 having a small molecular weight and filling in the macromolecular chains of NBR, making its chain segments move more easily, resulting in the glass transition region moving to a lower temperature. olymers 2023, 15,604 is called cohesive energy density (CED) [21]. Generally, the changing tren CED is the opposite [22]. The higher CED is, the stronger the intermolecu closer the molecular stacking in the system, the smaller the FFV, and the lo the material [23,24]. The calculation results of FFV and CED of the CP52/NB shown in Figure 9. It can be seen that the changes in FFV and CED are not no the increased CP52 content, FFV is always stable at around 38%, while CED the range of 330~380 J/cm −3 . The changing trend regarding FFV and CED i opposite. Therefore, FFV is not the main factor affecting the damping prope of the CP52/NBR materials. Increasing the CP52 addition does not signifi FFV. It is mainly due to the CP52 having a small molecular weight and fillin romolecular chains of NBR, making its chain segments move more easily, re glass transition region moving to a lower temperature.

Conclusions
In this article, CP52/NBR composites were prepared experimentally, an CP52 on the damping properties of NBR was studied. The molecular mo structed to calculate the intermolecular interaction parameters of the CP52 Then, the damping mechanism of the CP52/NBR composites was analyzed perimental results. The conclusions are as follows: (1) With increased CP52 content, the damping properties of the CP52/NB were greatly improved, and the characteristic absorption peaks of -CN infrared spectrum were red-shifted, indicating that intermolecular f formed between CP52 and NBR. (2) The RDF simulation results and FTIR test results show that there are hydrogen bonds in the CP52/NBR damping materials, Cl-CP-H····NC-Cl····H-NBR-CN. The linear regression analysis and Ebinding results showe gen bonding was the main factor affecting the damping properties of t materials.
(3) The FFV and CED of the CP52/NBR system did not change significa creased CP52 content, and the changing trend was roughly the oppo cated that FFV is not the main factor affecting the damping performa

Conclusions
In this article, CP52/NBR composites were prepared experimentally, and the effect of CP52 on the damping properties of NBR was studied. The molecular model was constructed to calculate the intermolecular interaction parameters of the CP52/NBR system. Then, the damping mechanism of the CP52/NBR composites was analyzed using the experimental results. The conclusions are as follows: (1) With increased CP52 content, the damping properties of the CP52/NBR composites were greatly improved, and the characteristic absorption peaks of -CN and -Cl on the infrared spectrum were red-shifted, indicating that intermolecular force could be formed between CP52 and NBR. (2) The RDF simulation results and FTIR test results show that there are two kinds of hydrogen bonds in the CP52/NBR damping materials, Cl-CP-H····NC-NBR and CP-Cl····H-NBR-CN. The linear regression analysis and E binding results showed that hydrogen bonding was the main factor affecting the damping properties of the CP52/NBR materials.
(3) The FFV and CED of the CP52/NBR system did not change significantly with increased CP52 content, and the changing trend was roughly the opposite. It is indicated that FFV is not the main factor affecting the damping performance and T g of the CP52/NBR system. The effect of CP52 on the glass transition region is mainly attributed to its plasticity.
On the basis of this research, MD simulation can be used to establish other polar small molecule and polar rubber system models, and experimental verification can be carried out to develop a series of damping materials with hydrogen bond network structure. Moreover, according to the principle of hydrogen bond-like, from the perspective of molecular structure and functional group, some polar small molecule substances were designed and synthesized for simulation and experiment, which further verified the 'hydrogen bond-like theory' and obtained excellent damping materials.
Funding: Independent Project of Naval University of Engineering, 425317T007.
Institutional Review Board Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available upon request from the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.