During the carbonization process of raw polymer precursors, graphitic structures can evolve. The presence of these graphitic structures affects mechanical properties of the carbonized carbon fibers. To gain a better understanding of the chemistry behind the evolution of these structures, we performed atomistic scale simulations using the ReaxFF reactive force field. Three different polymers were considered as a precursor: idealized ladder PAN (polyacrylonitrile), a proposed oxidized PAN and PBO (poly(p-phenylene-2,6-benzobisoxazole)). We determined the underlying molecular details of polymers conversion into a carbon fiber structure. Since these are C/H/O/N-based polymers, we first developed an improved force field for C/H/O/N chemistry based on the Density Functional Theory (DFT) data with a particular focus on N2 formation kinetics and its interactions with polymer-associated radicals formed during the carbonization process. Then, using this improved force field, we performed atomistic scale simulations of the initial stage of the carbonization process for the considered polymers. Based on our simulation data we determined the molecular pathways for the formation of low-molecular weight gas-species, all-carbon rings crucial for further graphitic structures evolution and possible alignment of the evolved all-carbon 6-membered rings clusters.

Irradiation of polymer films by a CO 2 infrared laser under ambient conditions converts the polymer into porous graphene or laser-induced graphene (LIG). Here, we simulate the formation of LIG from five different commercially available polymers using reactive molecular dynamics. We determined that the molecular structure of the parent polymer has a significant effect on the final graphitic structure. CO is liberated during the initial part of the LIG formation process when the polymer is converted into an amorphous structure, while H 2 is evolved steadily as the amorphous structure is converted to an ordered graphitic structure. The LIG structure has out-of-plane undulations and bends due to a significant number of 5-and 7-member carbon rings present throughout the structure. We find that the simulated molecular structure compares well with recent experimental observations from the literature. We also demonstrate that the yield of LIG is higher in inert conditions, compared to environments with oxygen. Polybenzimidazole-derived LIG has the highest surface area and yield among the five polymers examined. These findings provide knowledge of LIG formation mechanisms that can be leveraged for bulk LIG applications such as sensors, electrocatalysts, microfluidics, and targeted heating for welding polymers.

The superlative strength-to-weight ratio of carbon fibers (CFs) can substantially reduce vehicle weight and improve energy efficiency. However, most CFs are derived from costly polyacrylonitrile (PAN), which limits their widespread adoption in the automotive industry. Extensive efforts to produce CFs from low cost, alternative precursor materials have failed to yield a commercially viable product. Here, we revisit PAN to study its conversion chemistry and microstructure evolution, which might provide clues for the design of low-cost CFs. We demonstrate that a small amount of graphene can minimize porosity/defects and reinforce PAN-based CFs. Our experimental results show that 0.075 weight % graphene-reinforced PAN/graphene composite CFs exhibits 225% increase in strength and 184% enhancement in Young’s modulus compared to PAN CFs. Atomistic ReaxFF and large-scale molecular dynamics simulations jointly elucidate the ability of graphene to modify the microstructure by promoting favorable edge chemistry and polymer chain alignment.

Using Monte Carlo simulations we study two-dimensional prey-predator systems. Measuring the variance of densities of prey and predators on the triangular lattice and on the lattice with eight neighbours, we conclude that temporal oscillations of these densities vanish in the thermodynamic limit. This result suggests that such oscillations do not exist in two-dimensional models, at least when driven by local dynamics. Depending on the control parameter, the model could be either in an active or in an absorbing phase, which are separated by the critical point. The critical behaviour of this model is studied using the dynamical Monte Carlo method. This model has two dynamically nonsymmetric absorbing states. In principle both absorbing states can be used for the analysis of the critical point. However, dynamical simulations which start from the unstable absorbing state suffer from metastable-like effects, which sometimes renders the method inefficient.

No abstract