Phantom chain simulations for fracture of star polymer networks with various strand densities

Yuichi Masubuchi,Takato Ishida,Yusuke Koide,Takashi Uneyama

2024-06-14

Abstract:Despite many attempts, the relation between fracture and structure of polymer networks is yet to be clarified. For this problem, a recent study for phantom chain simulations [Macromolecules, 56, 9359 (2023)] has demonstrated that the fracture characteristics obtained for polymer networks with various node functionalities and conversion ratios lie on master curves if they are plotted against cycle rank. In this study, we extended the simulation to the effect of prepolymer concentration on the relationships between cycle rank and fracture characteristics within the concentration range of 1<c/c^*<8, concerning the overlapping concentration c^*. We created networks from sols of star-branched phantom bead-spring chains via end-linking reaction between different chains through Brownian dynamics simulations with varying the number of branching arms f from 1 to 8, and the conversion ratio {\phi}_c from 0.6 to 0.95. For the resultant networks, cycle rank {\xi} was consistent with the mean-field theory. The networks were uniaxially stretched with energy minimization until break to obtain modulus G, strain at break {\epsilon}_b, stress at break {\sigma}_b, and work for fracture W_b. With the branch point density \u{psion}_br, G/\u{psion}_br, {\epsilon}_b, {\sigma}_b/\u{psion}_br, and W_b/\u{psion}_br of the data for various f and {\phi}_c draw master curves if plotted against {\xi}. The master curves depend on c; as c increases, all the mechanical characteristics monotonically increase. If we plot {\sigma}_b/\u{psion}_br and W_b/\u{psion}_br against G/\u{psion}_br, the data for various f and {\phi}_c lie on master curves but depending on c. Consequently, the fracture characteristics are not solely described by modulus for the examined energy-minimized phantom chain networks.

Soft Condensed Matter

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to explore the influence of prepolymer concentration on the fracture behavior of star - polymer networks. Specifically, by simulating the fracture characteristics of star - polymer networks with different branch densities, the researchers aim to understand the relationship between network structure and fracture performance. The paper mentions that although previous studies have shown that the cycle rank (\(\xi\)) can be a key parameter for describing the fracture characteristics of polymer networks, most of these studies have focused on the influence of node functionality and conversion rate. This paper further extends these studies by considering the influence of prepolymer concentration.

### Main research questions:
1. **Influence of prepolymer concentration on the fracture characteristics of star - polymer networks**: By changing the concentration of prepolymer (\(c/c^*\) ranging from 1 to 8), the researchers explore its influence on network fracture characteristics (such as modulus \(G\), fracture strain \(\epsilon_b\), fracture stress \(\sigma_b\) and fracture work \(W_b\)).
2. **Relationship between cycle rank and fracture characteristics**: By simulating star - polymer networks under different prepolymer concentrations, the researchers verify whether the cycle rank \(\xi\) can still be a key parameter for describing fracture characteristics and explore whether these relationships depend on the prepolymer concentration.

### Research methods:
- **Model and simulation**: Use the phantom chain model for Brownian dynamics simulation to create star - polymer networks. Generate networks through end - linking reactions and conduct uniaxial tensile experiments after energy minimization.
- **Parameter variation**: The researchers change the functionality \(f\) of star - polymers (from 1 to 8), the conversion rate \(\phi\) (from 0.6 to 0.95) and the concentration \(\rho\) of prepolymer (from 2 to 16).
- **Data analysis**: Analyze the stress - strain curves during the stretching process, extract the modulus \(G\), fracture strain \(\epsilon_b\), fracture stress \(\sigma_b\) and fracture work \(W_b\), and correlate them with the cycle rank \(\xi\).

### Main findings:
- **Modulus and fracture characteristics increase with increasing prepolymer concentration**: As the prepolymer concentration increases, the modulus \(G\) of the network and fracture characteristics (such as fracture strain \(\epsilon_b\), fracture stress \(\sigma_b\) and fracture work \(W_b\)) all increase monotonically.
- **Role of cycle rank \(\xi\)**: For different functionality \(f\) and conversion rate \(\phi\), when the cycle rank \(\xi\) is taken as the abscissa, the modulus \(G/\nu_{br}\), fracture strain \(\epsilon_b\), fracture stress \(\sigma_b/\nu_{br}\) and fracture work \(W_b/\nu_{br}\) all fall on the master curves. However, these master curves depend on the prepolymer concentration \(\rho\).
- **Relationship between fracture characteristics and modulus**: Fracture characteristics are not only determined by the modulus, but also influenced by the prepolymer concentration. Even with the same modulus, networks with different prepolymer concentrations exhibit different fracture characteristics.

### Conclusion:
The research results show that the prepolymer concentration has a significant influence on the fracture characteristics of star - polymer networks, and the cycle rank \(\xi\) is still a key parameter for describing these characteristics. However, fracture characteristics are not only determined by the modulus, but also influenced by the prepolymer concentration. These findings provide new insights for understanding the fracture mechanism of polymer networks.