Coarse-grained molecular dynamics simulation on strain-hardening and fracture behaviors of polycarbonate: Effect of polydispersity and spatial distribution

Coarse-grained molecular dynamics simulation of polycarbonate is utilized to investigate the relationship between molecular structure (i.e., polydispersity and molecular spatial distribution) and strain-hardening and fracture behavior of polycarbonate. We find that strain-hardening modulus and chain extensibility, which are the constitutive parameters of the Eindhoven Glassy Polymer model are highly affected by spatial distribution but are insensitive to polydispersity. This is attributed to the higher rate of nonaffine deformation in the structure with a high radius of gyration. On the other hand, maximum stress at fracture is highly influenced by both spatial distribution and polydispersity due to the ability to sustain entanglements at larger strain. We suggest the phenomenological expression of maximum stress as a function of the radius of gyration, the number of entanglements, and polydispersity.


Introduction
Polycarbonate (PC) is a versatile amorphous thermoplastic material widely used in various applications ranging from automotive components to electronic devices.This is owing to its well-known properties, such as high-impact resistance, thermal stability, and optical transparency [1].These excellent mechanical responses of polycarbonate are attributed to its molecular structure [2].The molecular structure of polycarbonate is a system of macromolecular chains consisting of the arrangement of spatially distributed monomers restricted by topological constraints or molecular entanglements by neighboring chains.This molecular configuration plays a crucial role in determining the material's stress-strain behavior when subjected to uniaxial tension [3].Understanding the influence of molecular structure on mechanical response is essential to develop the material's performance.
The selection of the constitutive model and parameters used in numerical simulations, the finite element method (FEM), holds a critical role in describing the mechanical or stress-strain behaviors of PC.It is widely acknowledged that the molecular structure of polymers profoundly influences their mechanical properties [2].Previous research has focused on experimental data fitting rather than insights from molecular dynamics simulations to obtain the constitutive parameters [3][4][5][6].As a result, the intricate relationship between molecular structure and stress-strain behavior, as well as the underlying mechanisms, remains inadequately explained in existing studies.
The deformation behavior of polycarbonate can be categorized into three regions: pre-yield, post-yield, and fracture [4,7].The pre-yield behavior is purely elastic and consistent across all polycarbonate samples [7].In the post-yield region, the constitutive Eindhoven Glassy Polymer (EGP) model is employed to describe deformation behavior in this study [4,8,9].The EGP model characterizes strain-hardening 398 Page 2 of 13 behavior through two parameters: the strain-hardening modulus ( G r ) and chain extensibility ( ) [4].Within the strain- hardening region, the rate of local plastic activity is directly proportional to the strain-hardening modulus [10,11], and this local plastic activity can be quantified in MD simulations using nonaffine particle displacement.However, the EGP model has limitations in characterizing the fracture behavior of polycarbonate, as the maximum stress is influenced by various molecular parameters, including molecular weight, polydispersity index (PDI), the radius of gyration ( R g ), and molecular entanglement [12,13].The examination of the influence exerted by the polydispersity index and the radius of gyration on both the parameters of the EGP model and the maximum stress, therefore, represents a prospective avenue for gaining valuable insights into the enhancement of the mechanical properties of polycarbonate in the future.
While the molecular structure and behavior are difficult to observe and examine in conventional experimentation, a molecular dynamics (MD) simulation provides insight into atomic interaction and position over time during the application or testing.There are two techniques of MD: all-atom molecular dynamics (AAMD) and coarse-grained molecular dynamics (CGMD) simulations.In CGMD simulation, a group of atoms is treated as a single coarse-grained particle [14].The simplification reduces the computation time, allowing for a more extensive system and/or a more extended simulation.The principal advantage of computational simulation over experimental investigation is that it allows direct quantification of the molecular structure that affects macroscopic response.
The mechanical response of polycarbonate has been the subject of investigation in multiple recent studies, using both AAMD [10,[15][16][17][18] and CGMD [7,12,13,[19][20][21][22][23][24][25] simulations of tensile loading studies to gain insight into the behavior of this material.However, these MD studies investigated only either monodisperse models [7,12,19,20,22,25] or a typical realistic distribution of molecular weight model [10,15,16].These MD models were not appropriate for the investigation of the effect of molecular weight distribution on strain-hardening and fracture behaviors and the evolution of molecular entanglement of PC.Nevertheless, experimental studies of polycarbonate have indicated that the degree of molecular weight distribution measured by polydispersity index (PDI), which is defined as the ratio of weight-average molecular weight ( M w ) and number-average molecular weight ( M n ), i.e., PDI = M w /M n , can be a significant factor characterizing the mechanical performance of polycarbonate [8,[26][27][28].Additionally, recent studies revealed that the spatial distribution of each polymeric chain, which is quantified by the radius of gyration and segment length, also affects the deformation behaviors of polycarbonate [10,12,13,29].Therefore, it is urged to investigate how PDI and radius of gyration affect the strain-hardening and fracture behaviors of polycarbonate.
This study employs CGMD to investigate the influence of the polydispersity index and spatial distribution on the molecular behavior of glassy polycarbonate when subjected to uniaxial tensile deformation.Through the construction of a series of coarse-grained models characterized by the polydispersity index and radius of gyration, the research examines the effects of these two molecular parameters (PDI and R g ) on two constitutive EGP model parameters ( G r and ) in post-yielding strain-hardening region and maximum stress at fracture in detail.The strain-hardening modulus and chain extensibility of each structure are determined using the EGP model, and a mathematical expression is proposed to describe the maximum stress as a function of the polydispersity index, the radius of gyration, and molecular entanglement.The findings of this study contribute to the advancement of the understanding of the effect of the molecular structure of polycarbonate on strain-hardening and may inform future efforts to optimize these properties.

Coarse-grained particle model
The composition of a polycarbonate unit includes three different atomic arrangements: carbonate, phenylene, and isopropylidene, identified as particles A, B, and C, as shown in Fig. 1.The study employed Kubo et al.'s coarse-grained (CG) model [19], which reduces 33 atoms to four CG particles.This model decreases the simulation's degree of freedom by about 90% compared to the all-atom model, primarily by combining implicit hydrogen with the CG particle.As a consequence, this simulation can be up to 1000 times faster than the all-atom model.
The model defines potential energy into two sources: intermolecular interaction ( E inter ) and intramolecular inter- action ( E intra ).The intermolecular interaction originates from the non-bonding interaction determined by Morse's potential with a cut-off radius of 12 Å (1) where R ij is the displacement between a pair of non-bonded particles i-j; potential parameters D 0 , , and R 0 are the con- stants of Morse's potential determined for each combination of particles in an all-atom model.The intramolecular interaction originates from the bonding interaction and angle interaction and is expressed by where b and are the bond length and bond angle, respectively, within an individual molecular chain; the potential parameters K b l , K a l , b 0 , and 0 are constants determined for all combinations of bond length and bond angle types in an all-atom model.In the context of the physical properties, the parameters b 0 and 0 symbolize equilibrium bond lengths and bond angles, respectively.Recent research [19] has investigated the specific values of these parameters, which are presented in Table 1 for reference.All potential parameters are identical to the recent research on polycarbonate [19].In the CG representation, PC chains are stretched in a linear conformation during the deformation, and the dihedral term cannot be evaluated.Thus, the model does not include the dihedral term.Given that limited bond breaking is observed during polycarbonate fracture, it is notable that the model does not consider bond breaking [16].

Coarse-grained particle arrangement and radius of gyration control
The arrangement of CG particles within polycarbonate follows an equilibrium positioning strategy [13], wherein each particle is positioned in optimal equilibrium with the preceding particle along the chain (molecule).Noting that the first particle of the chain is placed in a random position in three-dimensional space.In the context of this study, the placement of CG particles is constrained by the equilibrium bond lengths ( b 0 ) and bond angles ( 0 ) on a two-dimensional plane, with the torsional angle ( ) maintained at zero.This constraint is upheld until a specified number of monomers, referred to as the segment length, is attained.Subsequently, a random torsional angle is assigned to the particle.A schematic illustration of the polycarbonate configuration generation process is presented in Fig. 2. Notably, for (2) a chosen segment length of three monomers, the random value selection applies to instead of zero as depicted in Fig. 2, while remains fixed at 0 or 76 degrees.Upon increasing the segment length, a corresponding increase in the radius of gyration is observed.Manipulating the segment length enables precise control over the radius of gyration within the polycarbonate structure.Figure 3 shows an example of a configuration generated with a segment length of 16 monomers.The following section will show the correlation between segment length and the resulting radius of gyration.

Preparation of simulation cells
Polycarbonate simulation cells were prepared with nine different M w ranging from 32.8 to 50.95 kg/mol, with an identical M n = 32.8kg/mol; therefore, there were nine different polydispersities of polycarbonate systems prepared for the investigation ranging from 1.00 to 1.56.The details and molecular distribution of the simulation cell are shown in Table 2.The study considered three different sizes of molecular chains: 32-monomer, 128-monomer, and 224-monomer, denoted as short, medium, and long chains, respectively.The monodisperse (PDI = 1.00) polycarbonate system is composed of 32 medium chains.To increase M w , four medium chains were split into two short and two long chains equally.The total molecular chains and M n of the system remain the same at 32 chains and 32.8 kg/mol, respectively, but M w gradually increases, and the length of the cubic structure is around 11 nm, see Fig. 4a.The study also addresses an issue related to simulation cell size, which is further elaborated in Appendix A of the research.The artificially polydispersed polycarbonate systems were used to study the effect of the presence of various sizes of molecular chains within the same polymer system on its molecular behavior and mechanical performance.
The spatial distribution of the polymeric chain is quantitatively described by the chain radius of the gyration where N p , m i , and Δr i are the number of particles, the mass of particles, the relative position of particle i with its center of mass of each molecular chain, respectively.The root-mean-square of the chain radius of gyration determines the spatial distribution of the polymer system and is denoted as R g .
To quantify the spatial distribution of mass-distributed structure, the root-mean-square radius of gyration normalized by molecular weight is used to describe the mechanical performance of the polycarbonate and is expressed as where N c and M j are the numbers of molecular chains and mass of molecule j, respectively.
To examine the combinatory effect of spatial distribution and polydispersity on molecular behavior, each case of polydispersity was created with eight different radii of gyration by giving eight different segment lengths: 2, 4, 6, 8, 12, 16, 24, and 32 monomers.Noting that the segment length is directly proportional to the radius of gyration, see [13] for more details.Ten initial structures of each segment length case are reproduced with different random seeds to cope with stochastic deviation.In total, 720 models of polycarbonate structure are investigated in this study, calculated by the number of polydispersities × the number of segment lengths × the number of random seeds ( 9 × 8 × 10 = 720 structures).
The initial structures of the polycarbonate system were prepared according to the above mass distributions and segment lengths.The structures were equilibrated at 1000 K for 0.1 ns using the isothermal-isobaric ensemble (NPT), then cooled down to 300 K at a rate of 1000 K/ns.Finally, the structures were under NPT ensemble at 300 K for 1 ns for the equilibration.The time step of the calculation is 1 fs in all simulations.A three-dimensional periodic boundary condition was implemented.All simulations were performed using a Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [30,31] and were visualized using Open Visualization Tool (OVITO) [32].
Apart from R n g , molecular entanglement plays an important role in determining polycarbonate's deformation and fracture behavior [7,12,13].The average number of Fig. 2 Schematic of polycarbonate configuration.Carbonate, phenylene, and isopropylidene particles are labeled blue, green, and red circles, respectively.All bond lengths and bond angles ( ) are constrained with equilibrium bond lengths and bond angles; see Table 1.The angles of 180 degrees are not labeled in the figure.Each monomer is connected on the same twodimensional plane, i.e. torsional angle ( ) = 0 degree (only shown in the third monomer) Fig. 3 Polycarbonate molecule with a segment length of 16 monomers (all types of particles are colored red) entanglement points ( Z e ) per molecule was evaluated using the Z1 code [33][34][35][36][37].The averaged molecular properties ( R n g and Z e ) among ten structures of each polydispersity and seg- ment length are shown in Table 3.The molecular parameters of each segment length across different polydispersities fall within the same range.R n g significantly increases as segment length, while Z e shows only a slight decrease.

Representative molecular structures
For explanation, in the result ("Results" section), the true stress-strain curves and molecular structure development of only four representative structures are described.The combination of each of the lowest spatial distribution (segment length of 2 monomers) of monodisperse (PDI = 1.00) and polydisperse (PDI = 1.56) and the highest spatial distribution (segment length of 32 monomers) of monodisperse and polydisperse were investigated, which were denoted as ML, MH, PL, and PH, respectively.The molecular properties of the representatives are shown in Table 4.

Deformation analysis setup
Monotonic loading uniaxial tension simulations under the NPT ensemble were performed on all the simulation cells by induced external strain at a constant true strain rate, d ∕dt = 2.398 × 10 −6 fs −1 , on the loading direction from true strain = 0 to 2.398 for 10 6 timesteps of 1 fs.In the other two lateral directions, the pressure is controlled at 1 atm during the deformation, see Fig. 4b.The simulations are carried out under an isothermal temperature of 300 K.The study adopted a true strain rate rather than an engineering strain rate because, from the experimentation, the strain-hardening of glassy polymer is closely related to a constant true strain rate [38].

Calculation and concept of nonaffine particle displacement
Nonaffine particle displacement ( r na i ) is calculated as the deviation between the current position and the affine position from the reference of the particle i.The affine position is the product of the stretch ratio and the position of the particle at the reference frame, which is at zero strain in this study.For the structure at a certain stretch ratio ( ), nonaffine particle displacement compared to the initial structure can be described as: where r na i , r i , and r i (0) are nonaffine particle displacement, the position, and the initial position of particle i, respectively. (5) To understand the connection between the nonaffine particle displacement and local plastic activity, we assume that a typical structure undergoes a bulk stretch with a ratio of 2.0.The microscopic distance between two particles must be doubled in length if they are purely affine transformed.However, the intermolecular and intramolecular interactions counteract and restrain purely affine transformation.The more nonaffine particle displacement, the more stress is required to move that particle further and relates to effective chain stiffness [39].Nonaffine particle displacement is detectable to a minimal extent within the elastic region, but becomes largely measurable in the plastic region, thereby describing local plastic deformation through nonaffinity.To describe the nonaffine particle displacement of the bulk structure, the root-mean-square of all particles in the structure is used and expressed as:

Eindhoven Glassy Polymer model
The post-yield behavior of polycarbonate undergoes strainhardening, described by the increase in stress after yielding at large strain, which is attributed to the entanglement network of polymeric molecules during the deformation.
The EGP model successfully captures the strain-hardening response of polycarbonate in the experiment by employing the Edwards-Vigils description in the molecular view [4,40].The simplified expression of true strain-hardening stress ( σH ) in isochoric uniaxial tension can be written as where is the stretch ratio defined as the ratio between the current length and the initial length.G r and fitting constant parameters of EGP, representing strain-hardening modulus and chain extensibility.The expression is evaluated by assuming the mobility of entanglement is chosen to be zero.σH is defined only in the post-yield region; σH = σ − σY , where σ is an observable true stress and σY is a yield stress.It is noteworthy that, in the case of is chosen to be zero, the expression of strain-hardening stress (Eq.6) reduces to a well-establish neo-Hookean model [29], showing a linear relationship between strain-hardening stress σH and the poly- nomial of stretch ratio 2 − −1 , i.e., σH = G r  2 −  −1 .
The strain-hardening modulus G r is the slope of the relation- ship.Hence, the nonlinearity is solely attributed to the chain extensibility.

Stress-strain relationship and strain-hardening & fracture behaviors
Figure 5 depicts the true stress-strain ( σ-ε ) relationships of four representative structures, namely ML, MH, PL, and PH.The magnification of a region with a relatively small strain from 0.0 to 0.6 is displayed in Fig. 5b, revealing that polydispersity and the normalized radius of gyration have a negligible effect on the elastic behavior and yield stress. ( Due to unclear yielding behavior, we adopt the scheme introduced in our previous work to evaluate yield stress [7].Nonetheless, in the post-yielding region, it is evident that the stress increases with a higher strain-hardening rate in the structure with a higher normalized radius of gyration, whereas the stress slightly increases in the structure with a lower normalized radius of gyration.The variation in polydispersity does not affect the deformation behavior in the low-strain region.The significant difference in the strain-hardening rate between the structures with a high and low normalized radius of gyration highlights the influence of chain distribution on the post-yielding strain-hardening behavior.In particular, the higher strain-hardening rate in the highly dispersed structures suggests that the chains undergo a more significant degree of nonaffine stretching, leading to higher resistance to deformation.Figure 6 displays the nonaffine particle displacement of four representative structures during the deformation.Observation of nonaffine particle displacement during the deformation reveals that nonaffine particle displacement of the structures with a higher normalized radius of gyration increases at a higher rate than those with a lower normalized radius of gyration.Since more stress is required to translate nonaffine particles [10], stresses of the cases MH and PH increase significantly. After the strain-hardening region, all the structures attain maximum stress, namely fracture.The σ-ε relationships of the cases ML and PL are similar, and so are the maximum stress values at around 2500 MPa.On the other hand, among the higher normalized radius of gyration structures (MH and PH), the stress-strain relationships are parallel from the beginning to around true strain = 1.10.Then, the case of MH separates and fractures soon at 11,000 MPa, but the case of PH continues strain-hardening and approaches maximum stress at 23,000 MPa.The polydispersity has a significant effect on the maximum stress of the high normalized radius of gyration polycarbonate system.The increment is attributed to the number of entanglements in the polymeric system.It is noted here that the observable maximum stress is remarkably high, especially in the case of PH.The value of 23,000 MPa is comparable with the theoretical strength of the ideal chemical bonds [7].Therefore, the bond-breaking may be effective in the actual situation, but for simplicity, we did not simulate the effect.
Figure 7 reveals the development of the number of entanglements per monomer during the deformation, = Z∕N m , where Z and N m are the total numbers of entanglements and monomers in the structure, respectively.In the cases of MH and PH, both display a consistent pattern of , from a similar starting value at around 0.3 entanglements per monomer and gradually increase between the range of true strain from 0.0 to 1.1, indicating a strong correlation between increments of stress and entanglements.From the observation around ε = 1.2 , an abrupt increase of is observed, which suggests the formation of The more entanglements are formed, the more degree of strain-hardening arises.In the case of PH, increases upon 0.5 entanglements per monomer, but in the case of MH reaches a maximum of 0.4 entanglements per monomer.Even with the same initial amount of entanglement, the development of under deformation is affected by polydispersity, particularly in the cases with a large normalized radius of gyration.

Effect of molecular structure on strain-hardening modulus and chain extensibility
From the observation of the true stress-strain relationship, the strain-hardening stress is obtained by subtracting yield stress from true stress.Figure 8 reveals the relationship between strain-hardening stress σH and the polynomial of stretch ratio 2 − −1 from the simulation data and the EGP model fitting, ranging from post-yielding to pre-fracture region.The fitting parameters ( G r and ) are reported in Table 4. Figure 8a shows the cases of ML and PL, in which the relationship is essentially linear because, from the fitting, chain extensibilities are approximately zero, see Table 4. Figure 8b reveals the nonlinear response of σH as a function 2 − −1 in two cases, MH and PH.The results suggest that the structure with a higher initial radius of gyration exhibits a higher strain-hardening modulus and chain extensibility.
The EGP model fitting is applied to all 720 structures in this research.Figure 9 shows the relationship between fitting parameters and the initial normalized radius of gyration.Figure 9a indicates that G r has a parabolic relationship with the R n g and is independent of polydispersity.While only a small variation of G r is observed in lower R n g , G r variates more widely in higher R n g .Due to the sensitivity of G r to the initial number of entanglements, the observed wide variation is attributed to the differences in the number of entanglements within the polymeric system, which were not investigated in this study.Taking into account the direct proportionality between the strain-hardening modulus and the true strain rate [41], the values of G r resulting from the simulation notably exceed those determined through experimentation.Figure 9b shows that the chain extensibility cannot be observed in the structure with lower R n g (less than 10∕ 2 √ kg∕mol ).The chain extensibility increases as R n g increases and is saturated at around 0.20.No structure has higher than 0.25.The values obtained from the model fitting are comparable with the values from the experimentation fitting, i.e., 0.23 [4].√ kg∕mol ) of different PDIs exhibit the same range below 5000 MPa.As R n g increases, the maximum stress of the higher PDI structure increases at a higher rate than that of the relatively lower PDI structure.The monodisperse structure exhibits maximum stresses that are confined to approximately 10000 MPa, while the polydisperse structure displays a significant increase in maximum stresses, reaching up to 25000 MPa.
According to our previous investigation, the phenomenological relationship between maximum stress, the normalized radius of gyration, and the average entanglements per molecule ( Z e ) follows the power law relationship [13].In addition to the previous model, the PDI is added to the relationship can be written as with four fitting parameters: 0 = 5.264 MPa, 1 = 1.774 , 2 = 0.536 , and 3 = 0.947 .Figure 10b shows the rela- tionship between the values from simulation and power law model.This phenomenological expression suggests that not only R n g and Z e have a significant effect on maxi- mum stress, but PDI also describes the fracture of polycarbonate.The expression only encourages the direction for further investigation of the connection between a macroscopic maximum stress and molecular parameters in the future.

Conclusion
In this study, we aim to investigate the effect of initial spatial distribution and polydispersity on the strain-hardening and fracture behavior of polycarbonate.We obtain a better understanding of deformation behavior by performing a series of CGMD uniaxial tensile deformations of polycarbonate with various initial radii of gyration and polydispersities.The structure with a larger radius of gyration exhibits a higher strain-hardening modulus and chain extensibility, while PDI has no effect on strain-hardening response.However, PDI has a positive correlation with maximum stress, which is attributed to the ability to maintain the entanglements of long chains in more significant strains.Lastly, the phenomenological relationship between molecular parameters and maximum stress is suggested to convey and inspire future investigation in this field.

Appendix A. Validation of simulation cell size
The validation of the selected simulation cell size (approximately 11 nm in length) employed in the present study is addressed through a comparative analysis involving structures of varying dimensions: 113 Å and 906 Å.The methodology applied for these structures is consistent with the approach described in "Coarse-grained particle arrangement and radius of gyration control" section, and they both undergo the same deformation setup as described in "Deformation analysis setup" section.Notably, the length of the simulation cell of the latter structure surpasses the radius of gyration by more than five times as recommended in the former literature [42].The true stress-true strain relationship of these two structures subjected to the deformation is depicted in Fig. 11.Remarkably, significant differences in both the maximum true stress and the strain-hardening behavior between the two structures are absent.Although the larger simulation cell size results in less stress fluctuation during deformation, it culminates in a comparable maximum true stress value at around 12000 MPa.This phenomenon can be attributed to the influence of the microscale void present in the structure characterized by the larger radius of gyration.The cell size plays an important role in constraining the void's dimensions during the deformation simulation; excessively diminutive cell sizes may yield erroneous outcomes.Conversely, when dealing with structures possessing large radii of gyration, the voids manifest a finely distributed nature in a small size [13].Consequently, the impact of the simulation cell size on the overall behavior becomes marginal in investigating maximum stress and strain-hardening behavior.

Appendix B. Validation of cooling rate
During the equilibration process, where artificial deterministic particle positions transition into a more realistic stochastic distribution, there are three distinct sub-processes.Firstly, an isobaric-isothermal (NPT) ensemble is maintained for 0.1 ns at 1000 K. Subsequently, the system experiences cooling within the NPT ensemble, transitioning from 1000 K to 300 K over a duration of 0.7 ns, resulting in a cooling rate of 1000 K/ns.Lastly, the system is further equilibrated in an NPT ensemble at 300 K for an additional 1 ns.
The validation of the selected cooling rate (1000 K/ns) employed in the present study is addressed through a comparative analysis involving structures of two cooling rates during equilibrium process.The structures are cooled down from 1000 K to 300 K with 1.4 ns and 0.7 ns or cooling rates of 500 K/ns and 1000 K/ns, respectively.The former cooling rate is performed to consider whether, or not, the cooling rate affects nonaffine displacement parameter during the deformation.The latter cooling rate is values applied throughout the study.The two cooling rates are subjected to the structures of low (L) and high (H) radii of gyration.Figure 12 shows the development of nonaffine particle displacement during the deformation of four structures.The result implies that there is no difference in nonaffine particle displacement due to changes in cooling rate.

Fig. 1
Fig. 1 Schematic of a polycarbonate (PC) molecular structure and coarse-grained CG representation.Particles A, B, and C represent carbonate, phenylene, and isopropylidene, respectively

Fig. 4 a
Fig. 4 a Simulation cell before deformation -each color represents an individual molecule.b Schematic of uniaxial tensile deformation

Fig. 7 Fig. 8
Fig. 7 Number of entanglements per ( ) as a function of true strain ( ε ) from for the case of MH and PH

Fig. 9 Fig. 10 a
Fig. 9 Fitting parameters of as a function of the normalized radius of gyration ( R n g ): a Strain-hardening modulus ( G r ) and b Chain extensibility ( ).Colors represent different polydispersities (PDI)

Fig. 11
Fig. 11 True stress-true strain curves of structure with two simulation cell sizes of 113 Å (black) and 906 Å (red)

Fig. 12
Fig.12 Nonaffine particle displacement ( r na ) as a function of true strain ( ε ) of structure with different radii of gyration and cooling rate.The structures of low (L) and high (H) radii of gyration have values of 47 Å and 170 Å, respectively.The cooling rates are 500 K/ns (dotted) and 1000 K/ns (line)

Table 3
Average molecular properties of molecular structures