Phantom-Chain Simulations for the Effect of Node Functionality on the Fracture of Star-Polymer Networks

Yuichi Masubuchi,Yuya Doi,Takato Ishida,Naoyuki Sakumichi,Takamasa Sakai,Koichi Mayumi,Kotaro Sato,Takashi Uneyama

DOI: https://doi.org/10.1021/acs.macromol.3c01291

2023-07-19

Abstract:The influence of node functionality (f) on the fracture of polymer networks remains unclear. While many studies have focused on multi-functional nodes with f>4, recent research suggests that networks with f=3 exhibit superior fracture properties compared to those with f=4. To clarify this discrepancy, we conducted phantom chain simulations for star-polymer networks varying f between 3 and 8. Our simulations utilized equimolar binary mixtures of star branch prepolymers with a uniform arm length. We employed a Brownian dynamics scheme to equilibrate sols and induce gelation through end-linking reactions. We prevented the formation of odd-order loops owing to the binary reaction and second-order loops algorithmically. We stored network structures at various conversion ratios ({\phi}_c) and minimized energy to reduce computation costs induced by structural relaxation. We subjected the networks to stretching until fracture to determine stress and strain at break and work for fracture, {\epsilon}_b, {\sigma}_b, and W_b. These fracture characteristics are highly dependent on {\phi}_c for networks with small f but relatively insensitive for those with large f. Thus, the networks with small f exhibit greater fracture properties than those with large f at high {\phi}_c, whereas the opposite relationship occurs at low {\phi}_c. We analyzed {\epsilon}_b, {\sigma}_b, and W_b concerning the cycle rank {\xi} and the broken strand fraction {\phi}_bb. We found {\epsilon}_b, {\sigma}_b/{\phi}_bb, and W_b/{\phi}_bb monotonically decrease with increasing {\xi}, and the data for various f and {\phi}_c superpose with each other to draw master curves. These results imply that the mechanical superiority of the networks with small f comes from their smaller {\xi} that gives higher {\epsilon}_b, {\sigma}_b/{\phi}_bb, and W_b/{\phi}_bb than the networks with large f.

Soft Condensed Matter

What problem does this paper attempt to address?

The problem that this paper attempts to solve is that the influence of the node functionality (\( f \)) on the fracture properties of polymer networks is still unclear. Although many studies have focused on multi - functional nodes (\( f>4 \)), recent research has shown that a network with a functionality of 3 exhibits better fracture properties than a network with a functionality of 4. To clarify this difference, the authors carried out phantom - chain simulations of star - polymer networks and studied the influence of the change in functionality between 3 and 8 on the network fracture properties.

Specifically, the paper explores this problem through the following methods:
1. **Simulation method**: Use the Brownian dynamics scheme to balance the sol and induce the gelation reaction while preventing the formation of odd - order rings and secondary rings.
2. **Structure storage and energy minimization**: Store the network structure at different conversion rates (\( \varphi_c \)) and perform energy minimization to reduce the computational cost.
3. **Tensile experiment**: Stretch the network until it breaks, and measure the stress (\( \sigma_b \)) at break, the strain (\( \varepsilon_b \)) and the fracture energy (\( W_b \)).
4. **Data analysis**: Analyze the relationship between the fracture characteristics and the cycle rank (\( \xi \)) and the fraction of broken chains (\( \varphi_b \)). It is found that the fracture characteristics decrease monotonically with the increase of \( \xi \), and the data of different \( f \) and \( \varphi_c \) can be overlapped into a master curve.

Through these methods, the paper aims to reveal the specific influence mechanism of the node functionality on the fracture properties of polymer networks, especially why the networks with smaller functionality exhibit better fracture properties at high conversion rates.