Проведен анализ регулярности/хаотичности орбит 45 шаровых скоплений (ШС) в центральной области Галактики радиусом 3.5 кпк, подверженных наибольшему воздействию со стороны вытянутого вращающегося бара. Использованы различные методы анализа, а именно: методы вычисления максимальных характеристических показателей Ляпунова, MEGNO (Mean Exponential Growth factor of Nearby Orbits), метод сечений Пуанкаре, частотный метод, основанный на вычислении фундаментальных частот, а также предложен новый метод, основанный на вычислении амплитудного спектра орбиты как функции времени и энтропии амплитудного спектра как меры хаотичности орбит. Обнаружена бимодальность в гистограмме распределения положительных показателей Ляпунова, вычисленных в классическом варианте, без перенормировки теневой орбиты, позволяющая реализовать вероятностный метод классификации ШС, что является также новым подходом. Для построения орбит шаровых скоплений использована модель гравитационного потенциала с баром в виде трехосного эллипсоида. Приняты следующие параметры бара: масса 10¹⁰M_⊙, длина большой полуоси 5 кпк, угол поворота оси бара 25°, скорость вращения 40 км с⁻¹ кпк⁻¹. Для формирования 6D-фазового пространства, требуемого для интегрирования орбит, использованы самые точные на сегодняшний день астрометрические данные со спутника Gaia (EDR3), а также новые уточненные средние расстояния до ШС. Проведена классификация ШС с регулярной и хаотической динамикой. Как показал анализ, шаровые скопления, обладающие малыми перицентрическими расстояниями и большими эксцентриситетами, в наибольшей степени подвержены влиянию бара и демонстрируют наибольшую хаотичность. Показано, что результаты классификации шаровых скоплений по характеру орбитальной динамики, полученные с использованием рассмотренных в работе различных методов анализа, хорошо коррелируют друг с другом.

The regularity/chaoticity of orbits of 45 globular clusters in the central region of the Galaxy with a radius of 3.5 kpc, which are subject to the greatest influence of the elongated rotating bar, is analyzed. Various methods of analysis are used, namely, the methods of calculating the maximum characteristic Lyapunov exponents (MCLE), MEGNO (Mean Exponential Growth factor of Nearby Orbits), the Poincare´ section method, the frequency method based on calculating fundamental frequencies, and a new method is proposed based on calculating the orbit amplitude spectrum as a function of time and calculating the entropy of the amplitude spectrum as a measure of orbital chaos. Bimodality is found in the histogram of the distribution of positive Lyapunov exponents calculated in the classical version, without renormalizing the shadow orbit, which allows implementing a probabilistic method for GC classification, which is also a new approach. To construct the orbits of globular clusters, we used the gravitational potential model with a bar in the formof a triaxial ellipsoid. The following bar parameters were adopted: mass 10¹⁰M_⊙, length of the semi-major axis 5 kpc, angle of rotation of the bar axis 25°, rotation velocity 40 km s⁻¹ kpc⁻¹. To form the 6D phase space required for integrating the orbits, we used the most accurate astrometric data to date from the Gaia satellite (EDR3), as well as new refined average distances to globular clusters. Globular clusters with regular and chaotic dynamics were classified. As the analysis showed, globular clusters with small pericentric distances and large eccentricities are most susceptible to the influence of the bar and demonstrate the greatest chaos. It is shown that the results of the classification of globular clusters by the nature of their orbital dynamics, obtained using the various analysis methods considered in the work, correlate well with each other.