The Dynamic State of a Pseudo-Crystalline Structure of B 42 Molecules

: In this paper, we have developed a mathematical model of the dynamics of molecular structures arising from an approach using B 42 molecules. The model is based on equations that determine the motion of the centers of mass of the molecules in question, and equations for the projections of the moments of quantities of motion of these particles on the axes associated with them. The thermal motions of boron atoms within each individual molecule are calculated using a bond-ordered potential. The intermolecular interactions of atoms have a potential that appears to coincide with that of free atoms. However, the interaction parameters differ. The obtained columns of molecular disks for B 42 are stable structures. We determine the characteristic vibration amplitudes and rotation frequency of the molecules in question.


Introduction
Boron has made a significant impact in our lives through its often unrecognized use in fertilizers, fungicides, soaps, detergents, and heat-resistant glassware [1]. Different applications have stimulated the development of boron science, namely, energy, organic synthesis, petroleum technology, nuclear physics, catalysts, chemistry, materials, and medicine [1,2].
Pham el al. [3] performed quantum-chemical calculations of boron structures by studying the geometry and stability of clusters B30, B36, B40, B42. In this paper, we investigate the structure that can be formed from B42 (C7v) and the dynamics of motion of each individual flat disk of B42. Disk B 42 is a molecule with a thickness of 1 atom. It should be noted that 2D materials or graphene analogs stand out among other materials for their unique properties, namely their excellent physical, chemical, optical, and biological qualities due to their uniform shape, high surface to volume ratio, and surface charge [4]. These properties are necessary for use in various fields, including medicine [5][6][7], electronics [8,9], mechanics and materials [10], thermoelectricity [11], and others.
Since 2D materials are of great interest, a large amount of work has also been also devoted to their synthesis [11,12]. In reference [13], a new two-dimensional all-carbon graphidine was obtained. The following are the properties of some 2D materials. Graphene doped with boron nitride was investigated [14] using density functional theory (DFT). The results show that the formation energy of these systems can increase or decrease with doping. Using DFT analysis, Yin et al. [15] showed that zigzag boron nitrides doped with Fe exhibit semiconducting properties with a direct forbidden band, while the primary nanoribbon from boron nitride is an insulator with an indirect wide energy gap. The same authors, using the molecular dynamics method, report that such iron-doped BN tapes can stably exist in real conditions. In reference [16], the thermal properties of hybrid nanosheets from boron graphene nitride with various states were studied by the molecular dynamics method. Reference [17] showed the possibility of combining a large volume of boron nitride plates with a small amount of polymer to manufacture plastic. Heat-conducting composite materials for thermoregulation of flexible electronic devices is a difficult task. Reference [18] managed to obtain a sample of polycrystalline cubic boron nitride. In reference [19], the authors study THz detectors based on 2D materials. The analysis performed in reference [20] shows the importance of understanding the influence of a tunable electronic structure on the mechanical response of materials made of carbon and boron nitride.
It is also important to consider other features related to boron. Chauhan et al. [21] discuss boronbased polymers, their capabilities, and applications. Issues regarding the impact of boron on the environment and human health are considered in reference [22]. Yuan et al. [23] obtained coralloid amorphous boron that crystallized upon electron-beam irradiation of the Cu-B-Mg alloy. The resulting boron nanospheres had a size of 5 to 20 nm and exist in an α-B structure. Gas adsorption on nanotubes and fulburenes is considered in reference [24]. Reference [25] reported that boron isotopes give boron arsenide single crystals an unusual and ultrahigh thermal conductivity. Reference [26] performed extensive optimization calculations for B80 isomers to find out which principles underlie the formation of large hollow frameworks from boron. Reference [27] reports on the method of growing single crystals of boron arsenide on a centimeter scale. The resulting crystals have unique physical properties. Their modeling showed that borofen can remain more resistant to external mechanical loads than graphene and BN sheets. The authors of [28] state that a polymer-crystalline electrolyte containing anion-trapping boron fragments improves the ionic conductivity of sodiumion batteries at room temperature. The results show that the adhesion energy for all five of these new allotropes is positive since all of these systems are stable; therefore, these structures can be synthesized in experiments. The possibility of using boron and its compounds for use in hydrogen production and storage systems, in lithium-ion batteries, supercapacitors, and thermal elements is discussed in reference [29].

Materials and Methods
In this paper, we consider the relative displacements of B42 molecular disks located in a structure that we call a column. To describe these displacements and thermal vibrations, we use a bondoriented intramolecular potential and a centrally symmetric potential for intermolecular interactions. Thus, we make the completely natural assumption that there are chemical bonds between the boron atoms inside the B42 molecule, and there are none between the molecules. The motion of disks as molecular formations obeys the equations of rotational motion and the equations of motion of their centers of mass.
The mesh structures of boron molecules are based on the triangles in Figure 1, as well as pentagons, hexagons, and heptagons. Triangles determine the local structure, and figures with a large number of angles determine the global structure since they are at the periphery of the surface structures under consideration. For example, a B30 molecule has a pentagon inside, with a side equal to one B-B bond. It has a pentagon on the periphery, too, but with a side equal to three covalent bonds. The same can be said with respect to B36 and B42 except that the basis of their designs is the hexagon and heptagon, respectively. Assuming that the molecules are deformed during the movement, we write the equations of the theorem on the angular momentum of the molecules for relative motion near their centers of mass: Here, , , L L L denote projections of the angular momentum vector, and , , p q r denote projections of the angular velocity vector. In this case, the projections of the angular momentum vector through the projections of the instantaneous angular velocity vector and the components of the inertia tensor are expressed as follows: Axial and centrifugal moments included in the last relations are expressed by the following formulas: , .
Here, m denotes the mass of any atom in the molecule, S1 denotes the number of these atoms in the described molecule. Equations (1)-(3) can easily be solved with respect to time derivatives of the projections of angular velocities on the moving axes. After that, they are integrated together with the Euler kinematic relations, which are also resolved with respect to the derivatives of the Euler angles.
The dynamic equations of rotational motion are initially presented in projections on the axis of the moving reference system associated with a single molecule. At the same time, it is convenient to represent the resulting motion in an absolute, motionless coordinate system. In this regard, we construct the absolute basis using the force characteristics: the projections of interatomic interaction forces and their moments at the first stage of calculations. Then, using the rotation matrix, which has components in the form of combinations of trigonometric functions of Euler angles, a transition is made to the projections of the moments of forces in a moving basis. These projections are then incorporated into the dynamic Euler equations. The latter equations are ordinary first-order differential equations with respect to the angular velocities of rotating molecules. The system of these equations will be closed if we include the kinematic Euler relations connecting the derivatives of Euler angles with the projections of the angular velocities and trigonometric functions of the rotation angles. Supplementing all these equations with initial data, we obtain the Cauchy problem for determining Euler angles as functions of time. The dynamic Euler equations are initially presented in projections on the axis of the moving reference system associated with a single B42 molecule. At the same time, it is convenient to represent the resulting motion in an absolute, motionless coordinate system. In this regard, we construct the absolute basis using the force characteristics: the projections of interatomic interaction forces and their moments at the first stage of calculations. Then, using the rotation matrix, which has components in the form of combinations of trigonometric functions of Euler angles, a transition is made to the projections of the moments of forces in a moving basis. These projections are then included in the dynamic Euler equations.
The magnitude of the force acting on a single atom of molecules B42 can be determined as follows: Here, , , k k k X Y Z denote projections of the resultant of all interatomic forces directed at the kth atom of the molecule in question; ki ρ denotes the distance between atoms inside an individual molecule; rkj denotes the distance between the boron atoms of various disks; S1 denotes the number of atoms in the described molecule, S2 denotes the number of atoms located close to the calculated molecule, Ubo denotes bond ordered force field (intramolecular potential), and vw U denotes the effective potential of interactions between atoms of various B42 disks.
The first group of terms on the right-hand sides of Equations (10)-(12) is responsible for the interaction inside the molecule, and the second group is responsible for cross interactions. The equations of motion of any atom written in the absolute reference frame will have the classical form: This theorem can be formulated as follows: The center of mass of the disk B42 moves as a material point having a mass equal to 42 B m under the action of a force equal to the sum of all cross "atomatomic" interactions in the system under consideration.
Here mB is the mass of the boron atom; v is the velocity of the B42 center of mass; Uvw(rij) is the potential of cross-actions; and grad is the gradient operator.
The right-hand side of the equation of the theorem under consideration does not include forces determined by the gradient of the bond-ordered force field. The latter is due to the fact that internal forces do not change the motion of the center of mass. The diatropic ring current present in the particle [1] causes plate B42 to be an elementary magnet. This allows the system of molecules to fit into columns. Similar processes occur in the blood of a person if red blood cells receive an excess charge for some reason. Conventional magnetic reactions resulting from the movement of electrons through the crystal structure of B42 should not be a reason for refusing to use the classical description of the dynamics of the aggregate of such particles. Magnetic forces of a classical nature enter the righthand sides of Equations (1)-(3). Differential Equations (1)-(3), (13), and (14) are integrated under the following initial conditions: Here, 0 0 0 , , k k k x y z denote the initial coordinates of boron atoms in a single molecule, 0 0 0 , , k k k u v w denote the initial velocity of the thermal motion of atoms. The Cauchy problem (1)-(18) was solved numerically using the fourth-order Runge-Kutta scheme of accuracy and a constant time step 6 10 t ns − ∆ = .

Algorithm Testing
We mention only three of the tests performed: 1. Inertia rotation of the free disk B42 in three different planes of the coordinate space. If there are no torques, then the initial speed remains constant throughout the movement. 2. Let the velocity vector of the test molecule representing the external action be in the plane of the disk B42 and not directed to its center of mass, which has a fixed position. Then, at the particle approaching site, the disk rotates at a certain angle, and, approaching according to the defined distance, the particle returns to its original position. 3. Two B42 plates located in the same plane but separated by a sufficient distance begin to move towards each other. Subsequently, the movement remains flat. The total energy of the system remains constant throughout the movement. The kinetic energy at the beginning of the test is equal to the kinetic energy of two molecules when they are removed from each other. All test examples showed an error no greater than 0.01%. Figures 2-9 show the dynamic state calculation data for columns containing two, three, and five B42 molecular disks.          -4 show the results of calculations of the dynamic state of two disks initially located in each other's potential wells. The calculations were performed without taking into account the influence of magnetic forces caused by the action of diatropic ring currents. Even in this case, the discs are held together, and their movements are symmetrical angular oscillations of large amplitude. As can be seen from Figure 2, the amplitude of these oscillations is about 170°. From Figure 3 it follows that the average oscillation frequency of the two disks is of the order 10 10 s −1 . With these oscillations, the trajectory of the node of the plate under consideration fills ¾ of the circle (Figure 4). Figures 5-7 show the results related to the "coin" column of three disks. As shown in Figure 5, the central disk in its angular oscillations makes almost a full revolution, and the peripheral disks make only half a revolution. This is also indicated by the trajectories shown in Figure 7. Figures 8 and 9 show the dynamic state of a column of five disks. As calculations show, with an increase in the number of disks, the symmetry of motion decreases. Now, the extreme disks make angular oscillations of the greatest amplitude (Figure 9). In this case, the average frequency of angular oscillations remained approximately the same. Along with torsional vibrations, the B42 molecules located in the "coin" column experience small vibrations in the heel and pitch angles (Figure 7). With a further increase in the number of disks in the column, the oscillation amplitude of the extreme disks is ¾ of the circle, and the internal disks have the amplitudes ¼ and ½. Therefore, the situation shown in Figure 9 is typical. Calculations show that linear stacking of the disks exhibits dynamic stability even if some of the internal disks participate in regular rotations. The latter is achieved by modifying the B42 molecule by introducing an ionizing group of atoms into it, as well as using an external electromagnetic field.

Discussion
Calculations show that planar molecules can be combined into "coin" columns. This is facilitated by both van der Waals forces and magnetic forces due to the presence of diatropic ring currents in B42 molecules. However, the resulting linear molecular structures are not static but are characterized by angular vibrations of large amplitude disks. An effective mathematical description of these oscillations is proposed. The angular velocities of the disks in the pseudo-crystalline structures under consideration are determined. Molecular crystals, such as plastic fullerite, have the ability to accumulate energy on internal, rotating degrees of freedom. B42. Molecular crystals that must be obtained from such columns must have outstanding hardness characteristics, since each molecule of such crystals is an elementary gyroscope and prevents access to it. In this study, an application of mechanical methods to describe the rotating dynamics of molecules in molecular crystals is performed.
The obtained calculated data can be a guideline when conducting measurements. It is also natural that the best verification of the correctness of the constructed mathematical model is a comparison with experimental data. In experiments, one can record the period of angular vibrations of B42 disks using nuclear magnetic resonance. Similar experiments were performed by Johnson et al. [30], for the case of rotation of fullerenes in the plastic phase of fullerite. In addition, in experiments, it is possible to determine the average distance between the B42 molecules in the columns and the binding energy between the parts of the column, which determines the strength of the considered molecular structures.