Modeling C540-C20 fullerene collisions

Collisions of C20 and C540 fullerenes are studied in a wide range of velocities by means of classical molecular dynamics. The simulations show that the collision scenario strongly depends on the collision velocity of the fullerenes. At low collision energies, the fullerenes are repelled by the van der Waals forces, and after bouncing off a part of the kinetic energy of their translational motion is converted into the energy of cage vibrations. At higher collision energies, the fullerenes overcome the potential barrier of Pauli repulsion with the formation of the new chemical bonds, and a significant change in the geometry of the molecules is observed. At very high collision energies, structure of the molecules is destroyed with the formation of new carbon clusters. Analysis of collisions for fullerenes of significantly different size revealed some new features associated with their asymmetric energy exchange.Published versio

Graphene nanoribbon as an elastic damper

Heterostructures composed of dissimilar two-dimensional nanomaterials can have nontrivial physical and mechanical properties which are potentially useful in many applications. Interestingly, in some cases, it is possible to create heterostructures composed of weakly and strongly stretched domains with the same chemical composition, as has been demonstrated for some polymer chains, DNA, and intermetallic nanowires supporting this effect of two-phase stretching. These materials, at relatively strong tension forces, split into domains with smaller and larger tensile strains. Within this region, average strain increases at constant tensile force due to the growth of the domain with the larger strain, at the expense of the domain with smaller strain. Here, the two-phase stretching phenomenon is described for graphene nanoribbons with the help of molecular dynamics simulations. This unprecedented feature of graphene that is revealed in our study is related to the peculiarities of nucleation and the motion of the domain walls separating the domains of different elastic strain. It turns out that the loading-unloading curves exhibit a hysteresis-like behavior due to the energy dissipation during the domain wall nucleation and motion. Here, we put forward the idea of implementing graphene nanoribbons as elastic dampers, efficiently converting mechanical strain energy into heat during cyclic loading-unloading through elastic extension where domains with larger and smaller strains coexist. Furthermore, in the regime of two-phase stretching, graphene nanoribbon is a heterostructure for which the fraction of domains with larger and smaller strain, and consequently its physical and mechanical properties, can be tuned in a controllable manner by applying elastic strain and/or heat

Transverse discrete breathers in unstrained graphene

Discrete breathers (DB) are spatially localized vibrational modes of large amplitude in defect-free nonlinear lattices. The search for DBs in graphene is of high importance, taking into account that this one atom thick layer of carbon is promising for a number of applications. There exist several reports on successful excitation of DBs in graphene, based on molecular dynamics and ab initio simulations. In a recent work by Hizhnyakov with co-authors the possibility to excite a DB with atoms oscillating normal to the graphene sheet has been reported. In the present study we use a systematic approach for finding initial conditions to excite transverse DBs in graphene. The approach is based on the analysis of the frequency-amplitude dependence for a delocalized, short-wavelength vibrational mode. This mode is a symmetry-dictated exact solution to the dynamic equations of the atomic motion, regardless the mode amplitude and regardless the type of interatomic potentials used in the simulations. It is demonstrated that if the AIREBO potential is used, the mode frequency increases with the amplitude bifurcating from the upper edge of the phonon spectrum for out-of-plane phonons. Then a bell-shaped function is superimposed on this delocalized mode to obtain a spatially localized vibrational mode, i.e., a DB. Placing the center of the bell-shaped function at different positions with respect to the lattice sites, three different DBs are found. Typically, the degree of spatial localization of DBs increases with the DB amplitude, but the transverse DBs in graphene reported here demonstrate the opposite trend. The results are compared to those obtained with the use of the Savin interatomic potential and no transverse DBs are found in this case. The results of this study contribute to a better understanding of the nonlinear dynamics of graphene and they call for the ab initio simulations to verify which of the two potentials used in this study is more precise