Revealing the High-Modulus Mechanism of Polyimide Films Prepared with 3,4′-ODA

To prepare PIs (polyimides) with desirable thermal and mechanical properties is highly demanded due to their widespread applications in flexible optoelectronic devices and printed circuit boards. Here, the PI films of BPDA/4,4′-ODA, BPDA/3,4′-ODA, PMDA/4,4′-ODA, PMDA/3,4′-ODA systems were prepared, and it was found that the PIs with 3,4′-ODA always exhibit a high modulus compared with the PIs with 4,4′-ODA. To disclose the mechanism of high-modulus PI films with 3,4′-ODA, amorphous PI models and uniaxial drawing PI models were established and calculated based on MD simulation. The PI structural deformations at different length scales, i.e., molecular chain cluster scale and repeat unit scale, under the same stress were detailed and analyzed, including the variation of chain conformation, bond length, bond angle, internal rotation energy, and torsion angle. The results indicate that PIs with 3,4-ODA have higher internal rotation energy and smaller deformation with the same stress, consistent with the high modulus.


Introduction
Polyimide (PI), originating from the polycondensation of diamine and dianhydride, is a kind of polymer possessing prominent mechanical properties, good temperature stability and a low coefficient of thermal expansion. Among its various structures, the PI based on the pyromellitic dianhydride (PMDA)/4,4 -oxidianiline (4,4 -ODA) system is widely used in microelectronics, engineering, and aerospace fields because of its simple synthesis process, low cost, and high performance [1][2][3][4]. For applying to different situations, polyimide could be prepared for diverse productions by varying the processing technique, including fibers, films, resin, and so on. During the age of electronic information, one significant application of polyimide films is the robust substrate for flexible optoelectronic devices and printed circuits [5], in which the PI films require both a high modulus to resist deformation and excellent thermal stability to withstand high temperatures during processing [6,7].
The traditional PI film is difficult to dissolve in organic solvents, so it is mainly prepared by a two-step method, that is, aromatic diamine and dianhydride are condensed and polymerized into the precursor of poly (amic acid) (PAA) in polar aprotic solvents, such as N-methyl-2-pyrrolidone (NMP), dimethylformamide (DMF), and dimethylacetamide (DMAc), and then dehydrated to form PI [8]. Therefore, the properties of PI depend on the type of diamine and dianhydride [9,10]. The different chemical structures of the two monomers could produce a great influence on the molecular chain conformation, condensed structure, and mechanical performance of PI [5,11,12].
In the process of designing and developing new materials, molecular simulation technology has been gradually established and widely used in biology, materials, and other fields [13,14]. Molecular dynamics (MD) simulations and Monte Carlo (MC) [15][16][17] techniques have been demonstrated to be highly consistent with experimental results. This

Preparation of PI Films
The molar ratio of PMDA and ODA was 1:1, ODA was dissolved in DMAc, and PMDA was added step by step. The PAA solution of 15 wt% PMDA/ODA system with appropriate viscosity was obtained by stirring at 0 • C for 6 h. After casting the PAA solution, the PAA solution was placed on the ultra-clean table for 48 h to eliminate foaming and ensure a uniform thickness of the film. After putting the film into a blast oven, the atmosphere was heated up to 300 • C according to a certain program, held for 2 h and slowly cooled.
The preparation process of the other three systems of PI films was the same as above.

Characterization
Fourier Transform Infrared Spectrometer (Nexus 670, Nicolet Company, Madison, WI, USA) measurements were performed with a scanning wavenumber ranging from 4000 cm −1 to 400 cm −1 at a scanning number of 16 in an ambient atmosphere. The FT-IR spectra from molecular simulation were calculated based on an optimized geometry at DFT/B3LYP/6-31G(d) level.
Dynamic thermomechanical analysis was employed on a DMA Q800 (TA Instruments, New Castle, DE, USA) at a temperature from 50 to 450 • C under an N 2 atmosphere. The heating rate and load frequency were set as 5 • C/min and 1 Hz, respectively.
The modulus was obtained by a universal tensile testing machine. The film samples were cut to 4 mm × 20 mm and tested at 25 • C.

MD Simulations
In this work, four kinds of polyimide were selected for molecular simulations. These polyimides were derived from two kinds of dianhydride, i.e., PMDA, BPDA, and two isomeric kinds of diamine, i.e., 4,4 -ODA, 3,4 -ODA. The corresponding molecular models were constructed using Material Studio (MS) 8.0.
A molecular chain containing 20 repeat units [30][31][32] was constructed by the Homopolymer tool in MS. Cells with ten molecular chains were constructed for each system by the Amorphous Cell module in MS. The initial density of each box is 0.1 g/cc. After the geometry optimization of models, a model with minimum energy was selected for dynamic equilibrium with the COMPASS force field [16,33].
The dynamic equilibrium was divided into four parts [23,24,34]: (i) Compress the cell by a 40 ps MD in the NPT ensemble (fixed number of particles, pressure, and temperature) at high pressure (0.5 GPa, approximate to the experimental density). (ii) Anneal the cell by subsequent MD runs in the NVT ensemble (fixed number of particles, volume, and temperature) at 598 K and 298 K for about 40 ps. (iii) Relax the cell by a 500 ps NPT MD and 0.0001 GPa to check the course of the density fluctuations. (iv) Equilibrate the cell by a 1000 ps NVT MD run in order to further improve the equilibration.
The Nose thermostat (Q ratio of 1) and the Berendsen barostat (decay constant of 0.1 ps) were used in the MD simulation.
The tensile model was constructed by applying fixed stress to the amorphous model to make the model produce strain and ensure that the strain of the model was within the elastic strain range.
In order to analyze the conformation changes of molecular chains during uniaxial drafting, radial distribution function (RDF), mean-square displacement (MSD) length distribution, and torsion distribution were estimated using Forcite Analysis tools.
According to the covalent and non-covalent interactions, the possibility of finding another atom within a certain distance of one atom is calculated by using RDF, which can be used to calculate the bond length and other ordered distances [35]. The definition of RDF is per Equation (1): where n B is the number of B atoms at a distance of r in a shell of thickness ∆r from atom A, N B is the total number of B atoms in the system, and V is the total volume of the system. MSD analysis was performed to determine the mode of displacement of particles followed over time [36]. The MSD is defined as Equation (2): where r i (0) and r i (t) are the initial and final positions of molecules (mass centers of particle i) over the time interval t, N is the atom number, and brackets denote the ensemble average. Fractional free volume (FFV) was calculated using the "Atom volume and surface" tool under the condition of applying a Connolly probe with 0 nm radius. The free volume is identified as the sum of static voids created by chain packing and transient gaps generated by thermal-induced chain rearrangement, which provides the diffusing molecules with a low-resistance path for transport. Internal free space and effective arrangement of molecular chains can be calculated by fractional free volume. FFV is given by Equation (3): In this equation, V and V W represent the polymer volume and the van der Waals volume, respectively.
Glass transition temperature (T g ) is calculated by the density change of the model from high temperature to low temperature, step by step, at 25 K, and the step cooling method is closer to reality than the continuous cooling method [37], as T g of stepwise cooling is closer to the experiment than continuous cooling [38].
The torsion angle is the dihedral angle of two conjugated planes, which can also be used to reflect the internal rotation of the polymer [39]. It is an important indicator of polymer conformation. Ψ1 is defined as the dihedral angle (C-N-C-C), and Ψ2 and Ψ3 are defined as the dihedral angle besides the ether linkage (C-O-C-C, C-C-O-C). The statistical torsion angle is 180 steps with 2 • per step.

Model Validation Mechanical Properties of Polyimide and FT-IR
To reveal the high-modulus mechanism of PIs with 3,4 -ODA in the present research, four types of polyimide models were obtained by binary condensation polymerization of PMDA, BPDA, 4,4 -ODA, and 3,4 -ODA were constructed for molecular modeling using Material Studio. The model validity and equilibrium check was performed by stating their density, energy, and gyration radius fluctuations during the long MD runs [40,41], as shown in Figure 1, by taking PMDA/4,4 -ODA systems as an example. The small density and energy fluctuations during the long MD runs can be clearly observed, and the gyration radius leveled off, demonstrating the stability of our molecular models. As shown in Figure 2, the elastic modulus E was experimentally measured through the tensile test, while the modulus was simulated through the MD simulation. According to the results, the experimental modulus of the synthesized PIs, in descending order, is PMDA/3,4′-ODA (3.00 GPa) > BPDA/3,4′-ODA (2.80 GPa) > BPDA/4,4′-ODA (2.60 GPa) > PMDA/4,4′-ODA (2.10 GPa), indicating that the PIs with 3,4′-ODA exhibit a high modulus  As shown in Figure 2, the elastic modulus E was experimentally measured through the tensile test, while the modulus was simulated through the MD simulation. According to the results, the experimental modulus of the synthesized PIs, in descending order, is PMDA/3,4 -ODA (3.  As shown in Figure 2, the elastic modulus E was experimentally measured through the tensile test, while the modulus was simulated through the MD simulation. According to the results, the experimental modulus of the synthesized PIs, in descending order, is PMDA/3,4′-ODA (3.00 GPa) > BPDA/3,4′-ODA (2.80 GPa) > BPDA/4,4′-ODA (2.60 GPa) > PMDA/4,4′-ODA (2.10 GPa), indicating that the PIs with 3,4′-ODA exhibit a high modulus compared with the PIs with 4,4′-ODA. Similarly, the simulated modulus in descending order is PMDA/3,4′-ODA (4.26 GPa) > BPDA/3,4′-ODA (3.19 GPa) > BPDA/4,4′-ODA (2.74 GPa) > PMDA/4,4′-ODA (2.07 GPa), proving the rationality of our molecular models. Thus, the constructed models were used for further exploration of the properties and chain behaviors of the polyimides.  The chemical structure of the PI films characterized by FT-IR is shown in Figure 3a. We found three characteristic absorption bands of polyimide at 1773 cm −1 , 1709 cm −1 , and 1368 cm −1 , which are attributed to the C-O asymmetrical stretching of imide groups, C-O symmetrical stretching of imide groups, and C-N stretching of the imide ring. The characteristic absorption bands at 1660 cm −1 and 1550 cm −1 of the precursor PAA did not appear at all, indicating that our films are completely cyclized into polyimide.
Interestingly, the peaks of the benzene ring of BPDA/3,4 -ODA and PMDA/3,4 -ODA around 1500 cm −1 are not single peaks. The FT-IR of these systems was simulated through Gaussian simulation, and we show the results taking the BPDA systems as examples. Figure 3b is the FT-IR simulation curve of BPDA/3,4 -ODA, and Figure 3c is the FT-IR simulation curve of BPDA/4,4 -ODA. The oscillator strength of the two peaks in Figure 3b is 351 (left) and 450 (right), and the peak intensity in Figure 3c is 912 cm −1 . Figure 3d is the vibration trajectory of the functional groups belonging to BPDA/4,4 -ODA. Figure 3e,f is the vibration trajectory of the functional groups belonging to BPDA/3,4 -ODA when the wavenumber is larger and smaller, respectively. It can be observed that the higher wavenumber is the para-substituted peak of the benzene ring, and the lower wavenumber is the meta-substituted peak of the benzene ring. It indicates that the synthesis of polyimide material was successful.
Polymers 2021, 13, 3175 6 of 16 351 (left) and 450 (right), and the peak intensity in Figure 3c is 912 cm −1 . Figure 3d is the vibration trajectory of the functional groups belonging to BPDA/4,4′-ODA. Figure 3e,f is the vibration trajectory of the functional groups belonging to BPDA/3,4′-ODA when the wavenumber is larger and smaller, respectively. It can be observed that the higher wavenumber is the para-substituted peak of the benzene ring, and the lower wavenumber is the meta-substituted peak of the benzene ring. It indicates that the synthesis of polyimide material was successful.

Analysis of Molecular Chain Conformation
In order to compare the deformation resistance and the changes of the internal microstructure of the material under the same stress, the molecular chain was extracted from the cell after 1000 ps NVT MD equilibrium of each PI system as the initial conformation of the PI molecular chain of the system. Then again, the molecular chain was extracted from the cell after applying 8 GPa stress in the same direction of each system as the stretching conformation of the PI molecular chain of the system, as shown in Figure 4.
It can be observed that the molecular chains of BPDA/3,4 -ODA in Figure 4a and BPDA/4,4 -ODA systems in Figure 4b are folded. The molecular chains of PMDA/3,4 -ODA in Figure 4c systems are also curled without stretching, but some segments are straight chains, and it is easier to form a chain orientation. After stretching, the proportion of straight chains is higher. The molecular chains of PMDA/4,4 -ODA systems in Figure 4d are spring-like spirals; it is difficult to have a relatively regular structure in some parts, and it is difficult to form a chain orientation and a molecular chain orientation. After careful analysis of the molecular chain conformation before and after stretching, it is observed that the molecular chains of the BPDA/3,4 -ODA system and the BPDA/4,4 -ODA system have a certain probability of straight parts in the short range. However, due to the asymmetry of the BPDA monomer, it is difficult to form a real regular molecular chain. However, because of the symmetry of the PMDA monomer, the PMDA/3,4 -ODA system can form a ladder-like regular structure, which leads to short-range order and easy to form a regular molecular chain. The PMDA/4,4 -ODA system is difficult to form an ordered structure because the bond angle of the C-O-C chemical bond naturally forms a relatively curled structure. The order of the molecular chain will also affect the rigidity of the polymer and also can explain the PMDA/3,4 -ODA and PMDA/4,4 -ODA system, although it also contains a PMDA monomer, which has a similar chemical structure; however, the PMDA/3,4 -ODA system can obtain higher modulus, the PMDA/4,4 -ODA system modulus is relatively minimum, and the BPDA/3,4 -ODA and BPDA/4,4 -ODA modulus are in the middle of the two.

Analysis of Molecular Chain Conformation
In order to compare the deformation resistance and the changes of the internal microstructure of the material under the same stress, the molecular chain was extracted from the cell after 1000 ps NVT MD equilibrium of each PI system as the initial conformation of the PI molecular chain of the system. Then again, the molecular chain was extracted from the cell after applying 8 GPa stress in the same direction of each system as the stretching conformation of the PI molecular chain of the system, as shown in Figure 4. It can be observed that the molecular chains of BPDA/3,4′-ODA in Figure 4a′ and BPDA/4,4′-ODA systems in Figure 4b′ are folded. The molecular chains of PMDA/3,4′-ODA in Figure 4c′ systems are also curled without stretching, but some segments are straight chains, and it is easier to form a chain orientation. After stretching, the proportion of straight chains is higher. The molecular chains of PMDA/4,4′-ODA systems in Figure  4d′ are spring-like spirals; it is difficult to have a relatively regular structure in some parts, and it is difficult to form a chain orientation and a molecular chain orientation. After careful analysis of the molecular chain conformation before and after stretching, it is observed that the molecular chains of the BPDA/3,4′-ODA system and the BPDA/4,4′-ODA system have a certain probability of straight parts in the short range. However, due to the asymmetry of the BPDA monomer, it is difficult to form a real regular molecular chain. However, because of the symmetry of the PMDA monomer, the PMDA/3,4′-ODA system can form a ladder-like regular structure, which leads to short-range order and easy to form a regular molecular chain. The PMDA/4,4′-ODA system is difficult to form an ordered structure because the bond angle of the C-O-C chemical bond naturally forms a relatively curled structure. The order of the molecular chain will also affect the rigidity of the polymer and also can explain the PMDA/3,4′-ODA and PMDA/4,4′-ODA system, although it also contains a PMDA monomer, which has a similar chemical structure; however, the The calculated mean-square end-to-end distance <Ree 2 > and the mean-square radius of gyration <S 2 > of each system before and after stretching are listed in Table 1. We also calculated the <Ree 2 > of each chain fully straightened without destroying its molecular structure, and the results are shown in Table 1. The <Ree 2 >/<S 2 > of the PMDA/4,4 -ODA system is the largest but has the lowest elastic modulus, indicating that the modulus is not directly related to <Ree 2 >/<S 2 > and molecular chain rigidity. By calculating <S 2 > -stretched/<S 2 >, we found that the change range of the cluster size of the PMDA/3,4 -ODA system was much smaller than that of the PMDA/4,4 -ODA system under the same strain. This indicates that the deformation resistance of the PMDA/3,4 -ODA system is significantly better than that of PMDA/4,4 -ODA in the length scale of wire cluster size. The <S 2 > -stretched/<S 2 > of BPDA/3,4 -ODA system is slightly smaller than that of BPDA/4,4 -ODA system, which is also the same as the modulus law of these systems. After analyzing the deformation resistance of each system at the scale of the molecular chain line, it was necessary to analyze the deformation resistance of each system at a smaller scale. The length of each repeating unit was statistically analyzed by MS before and after stretching. The distribution of repeating unit lengths was plotted into probability curves and noise reduction, and the result is shown in Figure 5. Through analysis of the curve, we found that the length of the repeat unit in the 3,4 -ODA system has little change from Figure 5a. There are slight peaks at 11.5 Å, 13.5 Å, and 14.5 Å, and the curve decreases at 11 Å. This fully illustrates that the 3,4 -ODA system has good deformation resistance under 8 GPa stress. The length of amorphous repeat units is largely located near 10 Å, and a large number of repeat units are still located near 10 Å under strain. The length of repeating units in the 4,4 -ODA system has changed greatly, from a large distribution near 10 Å to a concentrated distribution near 12 Å-14 Å and a strong peak at 14 Å. This shows that in the system under the action of stress, deformation occurred in a large number of repeating units, showing a very poor ability to resist deformation. Because the BPDA monomer itself is relatively long, the average length of the repeating unit of the system containing the BPDA monomer is longer than that of the system containing the PMDA monomer. For BPDA/3,4 -ODA, the curve moves to the right as a whole, forming new peaks at 14 Å, 16 Å, and 18 Å, but the overall distribution does not change significantly, which also reflects good deformation resistance. For the BPDA/4,4 -ODA system, it can be observed from the curve that the length of the repeat unit increases greatly between 16 Å and 19 Å, and the curve shifts to the right as a whole. It can be observed from the comparison that the system containing 3,4 -ODA has strong deformation resistance, and the system containing 4,4 -ODA has poor deformation resistance. Among them, the PMDA/3,4 -ODA system has the strongest deformation resistance, and the PMDA/4,4 -ODA system has the worst resistance. This result is consistent with the test results in the experiment, reflecting the rationality of the model. More importantly, a method for comparing the modulus of polymer is established, i.e., a method for qualitative analysis of the PI modulus and the variation of the conformation of molecular chains and repeated units under strain in the range of elastic deformation.
The deformation resistance at the length scale of the two repeat units is shown in Figure 6. For the PMDA system, it can still be observed that the deformation resistance of the PMDA/3,4 -ODA system is better than that of the PMDA/4,4 -ODA system. However, the strength of the deformation resistance of BPDA system cannot be clearly compared under more than two repeat unit lengths, so it does not continue to expand the statistical scale.
From the above analysis, it can be concluded that the deformation resistance of the PMDA/3,4 -ODA system is significantly better than that of the PMDA/4,4 -ODA system based on the structural variation analysis from one repeating unit, two repeating units, and molecular chain scale. Similarly, the deformation resistance of the BPDA/3,4 -ODA system is slightly better than that of BPDA/4,4 -ODA. The deformation resistance at the length scale of the two repeat units is shown in Figure 6. For the PMDA system, it can still be observed that the deformation resistance of the PMDA/3,4′-ODA system is better than that of the PMDA/4,4′-ODA system. However, the strength of the deformation resistance of BPDA system cannot be clearly compared under more than two repeat unit lengths, so it does not continue to expand the statistical scale.

Mechanism Analysis
For material deformation, there are four influencing factors: intermolecular forces, the change of covalent bond length, the change of internal rotation angle, and the change of bond angle. Because the deformation of each tensile model is in the range of elastic deformation, no molecular chain slip occurs, so the intermolecular force is not discussed.
To discuss the changes of the covalent bond length before and after stretching, the covalent bond changes in the PI of each system were analyzed by all-atom RDF. This method is to count the probability of atoms appearing at a certain distance from an atom, which might be used to reflect the bond length. The results are shown in Figure 7. The peaks at 1.1 Å and 1.2 Å are the bond lengths of C-H and C=O, respectively. The peak at 1.4 Å is very strong, and the composition is more complex, consisting of C-N on the imide ring, C-C on the benzene ring, and C-O in two kinds of ODA. Through comparison, it is found that the RDF curve changes before and after stretching are very small, indicating that the change in PI covalent bond length before and after stretching is negligible, which is not the main reason for the material modulus. From the above analysis, it can be concluded that the deformation resistance of the PMDA/3,4′-ODA system is significantly better than that of the PMDA/4,4′-ODA system based on the structural variation analysis from one repeating unit, two repeating units, and molecular chain scale. Similarly, the deformation resistance of the BPDA/3,4′-ODA system is slightly better than that of BPDA/4,4′-ODA.

Mechanism Analysis
For material deformation, there are four influencing factors: intermolecular forces, the change of covalent bond length, the change of internal rotation angle, and the change of bond angle. Because the deformation of each tensile model is in the range of elastic deformation, no molecular chain slip occurs, so the intermolecular force is not discussed.
To discuss the changes of the covalent bond length before and after stretching, the covalent bond changes in the PI of each system were analyzed by all-atom RDF. This method is to count the probability of atoms appearing at a certain distance from an atom, which might be used to reflect the bond length. The results are shown in Figure 7. The peaks at 1.1 Å and 1.2 Å are the bond lengths of C-H and C=O, respectively. The peak at 1.4 Å is very strong, and the composition is more complex, consisting of C-N on the imide ring, C-C on the benzene ring, and C-O in two kinds of ODA. Through comparison, it is found that the RDF curve changes before and after stretching are very small, indicating that the change in PI covalent bond length before and after stretching is negligible, which is not the main reason for the material modulus.  The difference in the chemical structure affects the conformation of the repeating unit. The internal rotation angle energy barrier of the repeating unit is obtained by MS. The rotation barriers of the internal rotation angle Ψ1, Ψ2, and Ψ3 of these systems are shown in Figure 8. By comparison, it is observed that the Ψ1′s energy barrier of each system has little fluctuation and is much smaller than Ψ2 and Ψ3, which is roughly distributed between 7 kcal/mol and 10 kcal/mol. Thus, the modulus provided by Ψ1 is almost the same. The PMDA/4,4′-ODA system has the smallest Ψ2′s and Ψ3′s rotation barriers of 108 kcal/mol, and the PMDA/3,4′-ODA system has the largest Ψ2′s and Ψ3′s rotation barriers of 212 kcal/mol. The order of Ψ2′s and Ψ3′s rotation barriers from big to small is PMDA/3,4′-ODA>BPDA/3,4′-ODA>BPDA/4,4′-ODA>PMDA/4,4′-ODA; it is consistent with the order of the modulus, thus indicating that the internal rotation angle rotation barrier is the key factor affecting the modulus of the PI film. This indicates that the mod- The difference in the chemical structure affects the conformation of the repeating unit. The internal rotation angle energy barrier of the repeating unit is obtained by MS. The rotation barriers of the internal rotation angle Ψ1, Ψ2, and Ψ3 of these systems are shown in Figure 8. By comparison, it is observed that the Ψ1 s energy barrier of each system has little fluctuation and is much smaller than Ψ2 and Ψ3, which is roughly distributed between 7 kcal/mol and 10 kcal/mol. Thus, the modulus provided by Ψ1 is almost the same. The PMDA/4,4 -ODA system has the smallest Ψ2 s and Ψ3 s rotation barriers of 108 kcal/mol, and the PMDA/3,4 -ODA system has the largest Ψ2 s and Ψ3 s rotation barriers of 212 kcal/mol. The order of Ψ2 s and Ψ3 s rotation barriers from big to small is PMDA/3,4 -ODA>BPDA/3,4 -ODA>BPDA/4,4 -ODA>PMDA/4,4 -ODA; it is consistent with the order of the modulus, thus indicating that the internal rotation angle rotation barrier is the key factor affecting the modulus of the PI film. This indicates that the modulus provided by the internal rotation of ether bonds in each system is quite different, which is a key factor affecting the mechanical properties of each system. Through MS, the torsion angle changes of Ψ3 in each system before and after stretching were counted, the curve of torsion angle distribution and probability was drawn, and the noise was reduced. The result is shown in Figure 9. Compared with the PMDA system and BPDA system, it is obvious that the torsional angle of the system containing 3,4′-ODA changes less, and the torsional angle of the system containing 4,4′-ODA changes more. During stretching, the torsional angles of PMDA/3,4′-ODA system at −145° and 145° decreased slightly, and the peak values moved to −158° and 158°. The torsional angles at −35° and 35° increased correspondingly, but more torsional angles remained near the original angle. For the PMDA/4,4′-ODA system, the peaks at −70° and 100° decreased significantly, and this conformation disappeared in large quantities and transformed into the conformations with torsion angles of −155° and 55°. This reflects that the torsional angle of the PMDA/3,4′-ODA system has changed in the deformation process. However, due to its high torsional energy barrier, it needs to absorb a large amount of energy to complete the torsion. The total number of changed conformations is significantly less than that of PMDA/4,4′-ODA, which is shown as a high modulus in macroscopic properties. The law of the BPDA system is similar to that of the PMDA system. For the BPDA/3,4′-ODA system, the peaks at −157° and 158° slightly weaken and move to −153° and 145°, and the Through MS, the torsion angle changes of Ψ3 in each system before and after stretching were counted, the curve of torsion angle distribution and probability was drawn, and the noise was reduced. The result is shown in Figure 9. Compared with the PMDA system and BPDA system, it is obvious that the torsional angle of the system containing 3,4 -ODA changes less, and the torsional angle of the system containing 4,4 -ODA changes more. During stretching, the torsional angles of PMDA/3,4 -ODA system at −145 • and 145 • decreased slightly, and the peak values moved to −158 • and 158 • . The torsional angles at −35 • and 35 • increased correspondingly, but more torsional angles remained near the original angle. For the PMDA/4,4 -ODA system, the peaks at −70 • and 100 • decreased significantly, and this conformation disappeared in large quantities and transformed into the conformations with torsion angles of −155 • and 55 • . This reflects that the torsional angle of the PMDA/3,4 -ODA system has changed in the deformation process. However, due to its high torsional energy barrier, it needs to absorb a large amount of energy to complete the torsion. The total number of changed conformations is significantly less than that of PMDA/4,4 -ODA, which is shown as a high modulus in macroscopic properties. The law of the BPDA system is similar to that of the PMDA system. For the BPDA/3,4 -ODA system, the peaks at −157 • and 158 • slightly weaken and move to −153 • and 145 • , and the peaks near −47 • and 47 • slightly increase. For the BPDA/4,4 -ODA system, the peaks near −95 • and 70 • on the curve decrease significantly, and become multiple peaks at −161 • , −109 • , −37 • , 19 • , 63 • , and 146 • , indicating that more conformations have changed. Compared with BPDA/3,4 -ODA, the conformational changes of the BPDA/3,4 -ODA system are more and more obvious, and the macroscopic performance is that the modulus is lower than that of the BPDA/3,4 -ODA system. Compared with BPDA/3,4′-ODA, the conformational changes of the BPDA/3,4′-ODA system are more and more obvious, and the macroscopic performance is that the modulus is lower than that of the BPDA/3,4′-ODA system. There is a wide range of C-O-C bond angles. In order to study the change of bond angle before and after stretching, the C-O-C bond angle before and after stretching was counted, as shown in Figure 10. It can be observed that during the stretching process, the bond angle of each system increases slightly, and the peak in the figure shifts about 1.5° to the right. However, this change is far from enough to double the deformation of the polymer. There is a wide range of C-O-C bond angles. In order to study the change of bond angle before and after stretching, the C-O-C bond angle before and after stretching was counted, as shown in Figure 10. It can be observed that during the stretching process, the bond angle of each system increases slightly, and the peak in the figure shifts about 1.5 • to the right. However, this change is far from enough to double the deformation of the polymer.

Thermal Performance Analysis
For polyimide film, its thermal property is an important indicator for engineering applications. The Tgs of these systems were tested through the DMA test. The curves of the DMA are shown in Figure 11A, and the test results of Tgs are listed in Table

Thermal Performance Analysis
For polyimide film, its thermal property is an important indicator for engineering applications. The Tgs of these systems were tested through the DMA test. The curves of the DMA are shown in Figure 11A, and the test results of Tgs are listed in Table 2  To calculate Tgs, we performed 500 ps NPT and NVT MD simulations at different temperatures, then the density-temperature relationship of polyimides obtained at different temperatures was depicted in Figure 11B, and the Tgs can be obtained from the inflection point of the fitting curves. It can be seen from Table 2 that the calculated values are relatively close to the test results, which further illustrates the rationality of the models. Tgs can usually reflect the rigidity and mobility of the molecular chains. The results

Thermal Performance Analysis
For polyimide film, its thermal property is an important indicator for engineering applications. The Tgs of these systems were tested through the DMA test. The curves of the DMA are shown in Figure 11A, and the test results of Tgs are listed in Table 2  To calculate Tgs, we performed 500 ps NPT and NVT MD simulations at different temperatures, then the density-temperature relationship of polyimides obtained at different temperatures was depicted in Figure 11B, and the Tgs can be obtained from the inflection point of the fitting curves. It can be seen from Table 2 that the calculated values are relatively close to the test results, which further illustrates the rationality of the models. Tgs can usually reflect the rigidity and mobility of the molecular chains. The results To calculate Tgs, we performed 500 ps NPT and NVT MD simulations at different temperatures, then the density-temperature relationship of polyimides obtained at different temperatures was depicted in Figure 11B, and the Tgs can be obtained from the inflection point of the fitting curves. It can be seen from Table 2 that the calculated values are relatively close to the test results, which further illustrates the rationality of the models. Tgs can usually reflect the rigidity and mobility of the molecular chains. The results showed that the orders of calculated Tgs are PMDA/4,4 -ODA > PMDA/3,4 -ODA > BPDA/4,4 -ODA ≈ BPDA/3,4 -ODA. The FFV of each system calculated by MS is shown in Figure 12a. The density, specific volume, van der Waals volume, and FFV of each system are listed in Table 2. It can be observed that the FFV from large to small is BPDA/4,4 -ODA > BPDA/3,4 -ODA > PMDA/3,4 -ODA > PMDA/4,4 -ODA. This law and Tg from low to high has been, according to the free volume theory, that the greater the FFV, the more cavities can provide the motion of the chain, so there will be lower Tg. The mean-square displacement (MSD) of each system was calculated to analyze the dynamic motion ability of the molecular chain. The MSD curves are shown in Figure 12b. The higher the MSD curve is, the better the motion ability of the polymer molecular chain. The MD simulation of the PI molecular chain at 298 K for 1000 ps was carried out. By comparison, we found that the stronger the molecular chain movement ability, the lower the Tg.  The FFV of each system calculated by MS is shown in Figure 12a. The density, specific volume, van der Waals volume, and FFV of each system are listed in Table 2. It can be observed that the FFV from large to small is BPDA/4,4′-ODA > BPDA/3,4′-ODA > PMDA/3,4′-ODA > PMDA/4,4′-ODA. This law and Tg from low to high has been, according to the free volume theory, that the greater the FFV, the more cavities can provide the motion of the chain, so there will be lower Tg. The mean-square displacement (MSD) of each system was calculated to analyze the dynamic motion ability of the molecular chain. The MSD curves are shown in Figure 12b. The higher the MSD curve is, the better the motion ability of the polymer molecular chain. The MD simulation of the PI molecular chain at 298 K for 1000 ps was carried out. By comparison, we found that the stronger the molecular chain movement ability, the lower the Tg.

Conclusions
The purpose of the current research was to determine factors that have the greatest impact on modulus after chemical structure changes. Analysis of the change of covalent bond length, bond angle, and torsion angle clearly indicates that the rotation in molecular chains is a key factor in providing modulus. What is more, by intuitively comparing the deformation resistance of different PIs systems at different length scales under the same strain, i.e., one repeating unit, two repeating units, and the whole molecular chain length scale, we have developed an efficient method in the present work to qualitatively predict the order of the modulus of the polymers. Both the experimental test and MD simulation show that the PI of PMDA/3,4′-ODA has high Tg and good deformation resistance, namely a high modulus. This PI itself has short-range rigid-straight units, and the molecular chain is straighter after stretching, which can be more neatly arranged, accounting for its better mechanical properties. Further research should be undertaken to exert the

Conclusions
The purpose of the current research was to determine factors that have the greatest impact on modulus after chemical structure changes. Analysis of the change of covalent bond length, bond angle, and torsion angle clearly indicates that the rotation in molecular chains is a key factor in providing modulus. What is more, by intuitively comparing the deformation resistance of different PIs systems at different length scales under the same strain, i.e., one repeating unit, two repeating units, and the whole molecular chain length scale, we have developed an efficient method in the present work to qualitatively predict the order of the modulus of the polymers. Both the experimental test and MD simulation show that the PI of PMDA/3,4 -ODA has high Tg and good deformation resistance, namely a high modulus. This PI itself has short-range rigid-straight units, and the molecular chain is straighter after stretching, which can be more neatly arranged, accounting for its better mechanical properties. Further research should be undertaken to exert the potential of this structure to prepare polyimide materials with long-range ordered structures such as polyimide fibers.