(Non)Resonance Bonds in Molecular Dynamics Simulations: A Case Study Concerning C60 Fullerenes

Jacek Siódmiak

DOI: https://doi.org/10.3390/e26030214

IF: 2.738

2024-02-29

Entropy

Abstract:In the case of certain chemical compounds, especially organic ones, electrons can be delocalized between different atoms within the molecule. These resulting bonds, known as resonance bonds, pose a challenge not only in theoretical descriptions of the studied system but also present difficulties in simulating such systems using molecular dynamics methods. In computer simulations of such systems, it is often common practice to use fractional bonds as an averaged value across equivalent structures, known as a resonance hybrid. This paper presents the results of the analysis of five forms of C60 fullerene polymorphs: one with all bonds being resonance, three with all bonds being integer (singles and doubles in different configurations), one with the majority of bonds being integer (singles and doubles), and ten bonds (within two opposite pentagons) valued at one and a half. The analysis involved the Shannon entropy value for bond length distributions and the eigenfrequency of intrinsic vibrations (first vibrational mode), reflecting the stiffness of the entire structure. The maps of the electrostatic potential distribution around the investigated structures are presented and the dipole moment was estimated. Introducing asymmetry in bond redistribution by incorporating mixed bonds (integer and partial), in contrast to variants with equivalent bonds, resulted in a significant change in the examined observables.

physics, multidisciplinary

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to deal with resonant bonds in molecular dynamics simulations, especially for different forms of C60 fullerenes. Specifically, the author analyzed five different forms of C60 fullerenes: one form in which all bonds are resonant bonds, three forms in which all bonds are integer bonds (different configurations of single and double bonds), and one form in which most bonds are integer bonds but ten bonds (located within two opposite pentagons) have a value of 1.5. The study explored the influence of different bond types on the structural stability of fullerenes by analyzing the Shannon entropy values of the bond - length distributions of these structures, the characteristic frequencies of the inherent vibrations (the first vibration mode), and the electric potential distribution maps. The introduction of asymmetric bond redistribution (i.e., a mixture of integer and partial bonds) led to significant changes in the observables under investigation compared to variants of equivalent bonds.