From a Dinuclear System to Close Binary Cosmic Objects
From a Dinuclear System to Close Binary Cosmic Objects
Adamian, G. G.;Antonenko, N. V.;Lenske, H.;Sargsyan, V. V.
2023-04-20 00:00:00
Review 1, 1,2 3 1,4 G. G. Adamian *, N. V. Antonenko , H. Lenske and V. V. Sargsyan Joint Institute for Nuclear Research, 141980 Dubna, Russia; antonenk@theor.jinr.ru (N.V.A.); sargsyan@theor.jinr.ru (V.V.S.) Bogoliubov Laboratory of Theoretical Physics of Joint Institute of Nuclear Research, Tomsk Polytechnic University, 634050 Tomsk, Russia Institut für Theoretische Physik, Justus-Liebig-Universität, D-35392 Giessen, Germany; horst.lenske@physik.uni-giessen.de A. Alikhanyan National Science Laboratory (YerPhI), Yerevan 0036, Armenia * Correspondence: adamian@theor.jinr.ru Abstract: Applying the ideas from microscopic objects to macroscopic stellar and galactic systems, the evolution of compact di-stars and di-galaxies is studied in the mass asymmetry coordinate. The formation of stable binary systems is analyzed. The role of symmetrization of an initially asymmetric binary system is revealed in the transformation of gravitational energy into internal energy of stars or galaxies accompanied by the release of a huge amount of energy. For the contact binary stars, the change of the orbital period is explained by evolution to symmetry in mass asymmetry coordinates. The matter transfer in binary black holes is studied. The conditions for the merger of black holes in a binary system are analyzed regarding the radiation of gravitational waves. Using the model based on the Regge-like laws, the Darwin instability effect in binary systems is discussed. New analytical formulas are derived for the period of orbital rotation and the relative distance between the components of a binary system. The impossibility of the appearance of a binary cosmic object from a single cosmic object is revealed. Keywords: close binary stars; close binary galaxies; mass asymmetry; fission PACS: 95.30.-k; 26.90.+n Citation: Adamian, G.G.; Antonenko, N.V.; Lenske, H.; Sargsyan, V.V. From a Dinuclear System to Close Binary Cosmic 1. Introduction Objects. Astronomy 2023, 2, 58–89. https://doi.org/10.3390/ Di-star systems are widespread objects: about a quarter of all stars belong to the di-star astronomy2020006 family. Compact or close di-stars (binaries) with separation distances of only a few stellar diameters are of interest for stellar evolution, for example, for a merger process of stars. Academic Editor: Ignatios The distinctive feature of close binary systems is matter transfer (and ejection) [1–7], which Antoniadis is obviously impossible between two well-separated stars. A spectacular recent case is Received: 26 July 2022 KIC 9832227, which attracted a lot of interest from the media when this binary star was Revised: 11 November 2022 predicted [8] to merge in 2022, which would lead to the formation of a red nova. Note that Accepted: 29 March 2023 this has not happened yet. The luminous red novae have only recently been identified as Published: 20 April 2023 a separate class of stellar transients [9]. They show up by relatively long outbursts with spectral distributions centered in the red, ranging in luminosities as intermediate between classical novae and supernovae. The luminous red nova V1309 Sco has played a key role. The Optical Gravitational Lensing Experiment (OGLE) revealed a shape of the light curve Copyright: © 2023 by the authors. characteristic of a contact di-star with an exponentially decreasing orbital period [10]. These Licensee MDPI, Basel, Switzerland. data well confirmed the previous prediction [11] that merging contact binary stars lead to This article is an open access article luminous red novae. distributed under the terms and The merger process is represented as common-envelope evolution [12,13]. Despite the conditions of the Creative Commons short duration of this process, it has crucial consequences for stellar evolution. The binary Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ components may survive with a reduced orbital separation, leading in some cases to 4.0/). compact binaries consisting of white dwarfs, neutron stars, and black holes. Some of these Astronomy 2023, 2, 58–89. https://doi.org/10.3390/astronomy2020006 https://www.mdpi.com/journal/astronomy Astronomy 2023, 2 59 systems can merge with the gravitational wave emission, as recently observed with the LIGO and Virgo interferometers [14]. Alternatively, the binary object can merge into a single one having exotic properties. There is a general opinion that the process of merger is an insufficiently studied problem of stellar evolution [8,15,16]. However, there is some consensus on the mechanisms of the important stages of contact binary evolution. Initial separations in a binary are large, with orbital periods of days to months. Interaction with the third star and the tidal friction may reduce the binary period to a day or two, as found in the dynamical simulations [17] and observations [18]. Both magnetic braking and internal structure evolution may bring the binary into contact. If this happens, the stellar radii and Roche geometry prevent further evolution of the orbital period, ending in a stellar merger. As found in ref. [11], red novae with a wide range of luminosities can result from contact di-stars with various initial masses. As is commonly believed, contact binaries, for example, W UMa, end their evolution by merging into a single star [7]. The dissipation of orbital energy in the initial phase of merger led to the appearance of the glowing red nova V1309 Sco observed in 2008. However, the details of the matter transfer between the parts of the contact binary object or the mechanism triggering the merger are still unclear. The work [16] favors a gradual mass transfer from a less massive (with a smaller radius) secondary star to a primary one driven by the energy received from the primary. The contact is supported by the aforementioned magnetic braking or internal structure evolution. This results in an abrupt end when the mass ratio leads to the Darwin instability to start a merger which for a red nova seems to be the triggering mechanism for the outburst [10,19]. The Darwin instability occurs when the mass ratio becomes small enough, and the heavy star can no longer keep the light star synchronously rotating via tidal interaction. The orbital angular momentum transferred from the intrinsic spin changes the orbit more than the spin, which leads to a runaway. As found, for massive primary stars of the main sequence, this happens at a mass ratio of 0.09 [20]. There is also another scenario [21]: after contact of stars, there is a brief but intense mass transfer in a di-star system, changing the originally more massive star into a less massive one. This process may oscillate until a stable contact configuration is eventually achieved. A dynamic mass transfer without the Darwin instability was also investigated [22], and mergers triggered by a tidal runaway based on a non-equilibrium response to tidal dissipation [23] were considered. Thus, more detailed observational and theoretical studies are required. Overcontact, contact, and near-contact binaries, forming di-star compounds with the average distances between the stars close to the sum of their radii [4–7] are of great interest in the study of stellar evolution. Information about the evolution of compact binaries is necessary to understand the processes observed in isolated stars. Compact binary stars are a good laboratory for a wide range of astrophysical phenomena, such as mass transfer between stars. The observations of their evolution verify our understanding of the inner structure of stars. The problem of the origin of binary stars or binary galaxies is still unclear [4–7,24,25]. As shown in ref. [26], there is no dissociative equilibrium between single and binary stars in the galaxy. The number of binary stars is many orders of magnitude larger than expected for dissociative equilibrium. So the origin of binary stars is not related to the capture of one star by another into a bound orbit. In addition, there is no sharp difference between close and wide binary stars, and the angular momenta relative to their centers of gravity are extremely large. Angular momenta are in the interval between values close to the angular moments of stars with extremely high rotational velocity (close binary stars) and values exceeding these values by thousands of times (wide binary stars). These facts refute the assumption about the origin of binary stars due to the fission of individual stars. According to the law of conservation, the angular momentum of a binary star does not exceed the angular momentum of individual stars with extremely high rotational speeds. These conclusions are valid only if we ignore external influences on binary star evolution [26–28]. Astronomy 2023, 2 60 Since mass transfer is an important observable for close di-stars and di-galaxies, it is meaningful to study the evolution of the system in the mass asymmetry coordinate h = ( M M )/( M + M ), where M (i = 1, 2) are the masses of the components of 1 2 1 2 i the binary system at fixed total mass M = M + M and orbital angular momentum L 1 2 of the system. Classical Newtonian mechanics can be used to explore the evolution of close binary stars and galaxies in their center-of-mass coordinate system by analyzing the total potential energy as a function of h [29–35]. The limits of the formation and evolution of binary systems are of interest. The methods used were tested for similar processes in nuclear systems where mass asymmetry is an important collective coordinate governing the fusion of two nuclei [36,37]. A nuclear molecule or a dinuclear system consists of two individual touching nuclei. There are two main collective degrees of freedom in a dinuclear system which govern its dynamics: (i) the relative motion between the clusters leading to quasistationary states in the internuclear potential and to the decay into two fragments which is called quasifission since no compound system is first formed, (ii) the transfer of nucleons or light particles between two clusters of the dinuclear system leading to evolution in mass and charge asymmetries. There is a structural forbiddenness for the motion of the nuclei to smaller internuclear distances during the fusion process. Fusion of heavy nuclei in the internuclear distance R is impossible and can occur due to the transfer of nucleons, i.e., by a motion in h [36,37]. Nuclear dynamics is certainly different from the gravitational interactions of di-stars and di-galaxies. Nuclear reactions are governed by short-range strong interactions, onto which minor contributions of long-range (repulsive) Coulomb and centrifugal forces are superimposed. The dinuclear approach is a key tool for describing the fusion of two heavy nuclei. In the approaching phase and after fusion, there is a mass loss by the emission of protons, neutrons, and light clusters like Alpha particles. Once a critical distance and mass ratio are reached, fusion occurs. A highly excited compound nucleus is formed in thermal equilibrium at temperatures of the order of one to a few MeV, corresponding to 10 K, cooled down rapidly by ejection of nucleons, nuclear clusters, and g quanta. Hence, dinuclear dynamics covers essentially the same spectrum of phenomena as expected for di-stars and di-galaxies. Thus, it is worth exploring to what extent the method from the femtoscale of microscopic objects is applicable to macroscopic binary galactic and stellar systems [29–35]. 2. Theoretical Approach The differential of the total energy of a binary stellar or galactic system is expressed as a function of relative distance r, conjugate canonical momentum p, and mass asymmetry coordinate h as ¶E ¶E ¶E ¶E dE(r, p, h) = dt + dr + dp + dh. (1) ¶t ¶r ¶p ¶h As we consider the binary system as a closed system, the conservation of the total energy results in ¶E dh = 0 (2) ¶h dt where ¶E/¶h = 0 is the general solution. In the center of the mass system, the total energy of the binary system is a sum of radial and orbital parts of kinetic energies and potential energy. As seen below, we attach the orbital kinetic energy part to the interaction V between two components of the binary system. In this case, the expression of the total energy of the di-star system reads as E = + U, (3) 2m Astronomy 2023, 2 61 where p is the radial component of the momentum p and m = m(h) = M M / M r 1 2 = M(1 h )/4 is the reduced mass. The total potential energy of a binary stellar or galactic system, U = U + U + V, (4) 1 2 is the sum of the potential energies G M U = w , (5) k k 2R of its components (k = 1, 2), and the energy V of their interaction. The radiation energy is much smaller than the absolute values of U and V to be disregarded. In Equation (5), G, 1,2 w , M , and R are, respectively, the gravitational constant, the dimensionless structural k k k factor, the mass, and the radius of the component. In general, the value of w is determined by the density profile of a stellar or a galactic object. By employing the relation known from observations, we express the radius of the object in terms of its mass as [29,30] R = M where n and g are the constants. So 2 n Ggw M U = . (6) Since two objects rotate around the common center of mass, the star–star interaction potential contains, together with the gravitational energy of interaction V , the kinetic energy of orbital rotation V : V(r) = V + V = V + , (7) G R G 2mr where L is the orbital angular momentum of the binary system, which is conserved during the conservative mass transfer. At r r = R + R and r r , t 1 2 t G M M V (r) = (8) and G M M r 1 2 V (r) = 3 , (9) 2r t r respectively [38]. Here, r is the touching distance. From the conditions ¶V /¶rj = 0 t r=r 2 2 and ¶ V /¶r j > 0, we find the equilibrium relative distance between two objects r=r corresponding to the minimum of V: r = (10) Gm M at r r and m t 1/4 2 3 L r r = (11) Gm M at r < r . Finally, one can derive the expression for the object–object interaction potential m t Astronomy 2023, 2 62 G M M G 1 2 3 3 V(r ) = = w M M (12) m V 1 2 2r 2 at r r or m t " # 2 2 G M M 2r G g M M 2g 1 2 1 2 V(r ) = 3 = 3 (13) n n 2 2 4 4 n n 2r 2 M + M t r w M M ( M + M ) 1 2 t V 1 2 1 2 at r < r . Here, m t w = M m r mi M M M 2 1i 2i and r and m = m(h ) = = (1 h ) are, respectively, the distance between its mi i i M 4 i components and the reduced mass of the initial binary system. Deriving Equations (12) and (13), we employ the known relation r = r . (14) mi Using Equations (6), (12) and (13), we obtain the final expression for the total potential energy (4) of the binary system 2 n 2 n 3 3 U = g[w M + w M ] + w M M (15) 2 V 1 2 1 2 1 at r r and m t " #! G M M 2g 1 2 2 n 2 n U = g w M + w M + 3 (16) 1 2 n n 2 4 4 n n 2 2 M + M w M M ( M + M ) 1 2 V 1 2 1 2 at r < r . m t For the binary systems considered, the velocity v (r ) c, where c is the speed of r m light; therefore, the relativistic effects can be ignored. Since G M /r c , the gravitational field is rather weak, and the equations of Newtonian mechanics can be used instead of the equations of general relativity. Note that exotic binary systems with neutron stars, white dwarfs, and black holes are not considered in this section. 2.1. Binary Stars In order to calculate the factor w , we make use of the single-star model from ref. [6], because it describes well the observed relations between the temperature, radius, mass, and luminosity of stars; the mass distribution of stars; magnetic fields of stars; spectra of seismic vibrations of the Sun; and other features of stars. Employing the dimensionless structure factor 1/4 w = 1.644 (17) from [6], the observed radius-mass relation (n = 2/3) 2/3 1 M 2/3 R = M = R , g M M M and the relation between the star masses M = (1 + h) and M = (1 h) in the binary 1 2 2 2 system and the mass asymmetry coordinate h, we find from Equations (15) and (16) that G M 13/12 13/12 2 3 U = a[(1 + h) + (1 h) ] + b[1 h ] (18) 2R Astronomy 2023, 2 63 at r r and m t G M 13/12 13/12 U = a[(1 + h) + (1 h) ] 2R 1 h g + b 3 (19) 2/3 2/3 2 2/3 2/3 1/2 (1 + h) + (1 h) [1 h ][(1 + h) + (1 h) ] at r < r , where m t 13/12 a = 1.644 , 2M 5 3 G M R G M R M b = = , 2 2 2 64L M 2L 2M i i 4/3 b = , 2M 1/3 10/3 2 L M g = . 1/2 11/3 GR M Here, we assume that the orbital angular momentum L and total mass M are conserved during the evolution of a di-star in the mass asymmetry coordinate h. The orbital angular momentum L is calculated by using the observed masses M of stars and period P of i k,i orb,i their orbital rotation [39–43]. As seen from Equation (18), the stability of a binary stellar system depends on the orbital- 1/3 1/2 2 2 rotation period P or on the value of L = m (G Mr ) = m G M P /(2p) and orb,i i i mi i orb,i total mass M. Employing Equation (18), we can study the evolution of a binary stellar system in h. The extremes of the potential energy are determined from the numerical solution of the equation ¶U G M 13 1/12 1/12 2 2 = a[(1 + h) (1 h) ] 6bh[1 h ] = 0. (20) ¶h 2R 12 As follows, h = h = 0 (symmetric binary system) is the root of Equation (20) and at h = 0 the potential has a minimum if a < a = b cr or 1/2 7/8 128.5p M P < orb,i 2 3 G M 2M (1 h ) and a maximum if a > a . The minimum at h = 0 lies symmetrically relative to the cr two barriers at h = h . Expanding Equation (20) up to the third-order terms in h and solving it, the barrier positions are obtained as 1/2 864 b 22464a 1/2 h = 2 . 864 b + 3289a So at a < a the potential energy as a function of h has two symmetric maxima at h = h cr b and the minimum at h = h = 0. The fusion of two stars with jh j < h occurs only i b by overcoming the barrier at h = +h or h = h . With decreasing ratio a/b, the value b b of B = U(h ) U(h ) increases and the symmetric di-star system becomes more stable. h b i The evolution of two stars with 0 < jh j < h to the symmetric di-star configuration is i b energetically favorable. Thus, an initially asymmetric binary system with jhj = jh j < h is i b driven to mass symmetry, implying a flow of mass towards equilibrium and an increase in internal energy of stars by the amount DU = U(h ) U(h = 0) (Figure 1a). At a a , i cr Astronomy 2023, 2 64 h = h = 0 and the inverse U-type potential has a maximum at h = 0. In such a system, the fusion of stars (one star “swallows” another) is the only mode of motion in h transforming the di-star into a mono-star with the release of energy E = U(h ) U(h = 1) (Figure 1b). (a) (b) = =0 m b -1 - 0 1 -1 0 1 b i b Figure 1. The schematic driving potential energies of the star–star systems at a < a (a), and cr a > a (b). The arrows show the corresponding initial binary stars. The notation used in the text cr is indicated. 1 1/2 1/2 If b a, then h ! 2 0.71. In this case, the condition 0 < h < 2 means b b 1/2 2 that an asymmetric system with the mass ratio M / M > (1 + 2 ) 6 moves to more 1 2 asymmetric configurations and the relative distances between its components increase in accordance with Equation (11). These unstable binary stars with jhj > h are unlikely to live long. Indeed, close binary stars with a high mass ratio are very rare objects (Figure 2). Note that this constraint on the mass ratio is independent of the relation used between the mass and the radius of the object. 0.9 0.8 Binary stars Binary galaxies 0.8 0.7 0.7 0.6 0.6 0.5 0.5 0.4 0.4 0.3 0.3 0.2 0.2 0.1 0.0 0.1 11 12 13 1.5 2.0 2.5 3.0 3.5 10 10 10 M/M M/M Sun Sun Figure 2. Mass asymmetry h versus the total mass M/ M = ( M + M )/ M (M = M ) for i 1i 2i Sun the di-stars and di-galaxies from Tables 1–4. Astronomy 2023, 2 65 Table 1. Calculated (n = 2/3) values of jh j, DU = U(h ) U(h = 0), B = U(h ) U(h ), i i h b i B = U(h ) U(h = 1), and data from observations on M / M , M / M , P [8,42,43] for the close f b 1 2 orb binary stars (BS) listed. P DU B B M M orb h f 1 2 BS jh j M M (Days) (J) (J) (J) 41 40 39 AH Aur 1.68 0.28 0.71 0.4941 10 2 10 8 10 41 40 40 AP Aur 2.05 0.50 0.61 0.5694 10 3 10 10 40 40 39 DN Aur 1.44 0.30 0.66 0.6169 8 10 10 9 10 39 40 40 AW Vir 1.11 0.84 0.14 0.3540 4 10 4 10 10 41 39 39 AW UMa 1.38 0.14 0.82 0.4387 10 7 10 5 10 41 40 40 HV UMa 2.84 0.54 0.68 0.7108 2 10 3 10 2 10 40 40 39 KIC9832227 1.40 0.32 0.63 0.4583 7 10 2 10 9 10 40 40 39 HV Aqr 1.31 0.19 0.75 0.3734 9 10 10 6 10 40 40 39 QX And 1.23 0.29 0.62 0.4122 6 10 10 8 10 41 40 40 RR Cen 2.09 0.45 0.65 0.6060 10 2 10 10 39 40 40 EM Lac 1.06 0.67 0.23 0.3891 7 10 3 10 10 40 40 39 GW Cep 1.06 0.39 0.46 0.3188 3 10 2 10 8 10 40 40 40 V700 Cyg 0.92 0.60 0.21 0.3400 6 10 3 10 10 40 40 39 V870 Ara 1.34 0.11 0.85 0.3997 7 10 5 10 4 10 Table 2. The calculated (n = 2/5) values of jh j, DU = U(h ) U(h = 0), B = U(h ) U(h ), and i i b i B = U(h = h ) U(h = 1) and observational data (ordinal numbers and morphological types Ty f b of galaxies according to Hubble’s classification, total orbital mass M/ M of pairs in Sun’s mass units, projection of the linear distance X between the components, and linear diameters A (1i) = 2R and 25 1i A (2i) = 2R of the components) of close binary elliptical galaxies (BG) from the catalog of isolated 25 2i galaxy pairs [24]. X A (1i) A (2i) DU B B M 25 25 h f BG Ty jh j (kpc) (kpc) (kpc) (J) (J) (J) 11 50 50 51 194 E-E 2.72 10 27.0 36.9 29.3 0.28 2 10 2 10 3 10 11 51 50 51 279 E-E 3.76 10 20.6 27.2 17.5 0.50 3 10 2 10 6 10 11 48 50 52 399 E-E 5.84 10 28.2 27.1 26.5 0.03 3 10 2 10 2 10 13 50 52 54 501 E-E 1.05 10 38.3 36.0 35.7 0.01 10 3 10 4 10 12 50 48 52 554 E-E 1.05 10 57.6 52.2 42.7 0.25 2 10 2 10 3 10 11 49 50 51 577 E-E 2.27 10 21.3 29.0 25.5 0.16 6 10 2 10 3 10 Table 3. The same as in Table 2 but for other close binary elliptic-spiral galaxies. X A (1i) A (2i) DU B B M 25 25 h BG Ty jh j (kpc) (kpc) (kpc) (J) (J) (J) 11 49 49 51 144 E-Sa 1.43 10 17.1 18.3 16.0 0.17 3 10 3 10 2 10 11 49 49 51 254 E-Sb 5.16 10 47.7 48.9 43.3 0.15 6 10 8 10 8 10 10 51 48 50 331 SO-E 8.15 10 17.4 23.4 10.1 0.78 2 10 3 10 2 10 11 50 48 51 552 Sa-E 5.01 10 38.8 39.5 26.2 0.47 3 10 6 10 9 10 Astronomy 2023, 2 66 Table 4. The same as in Table 2 but for close binary spiral galaxies. X A (1i) A (2i) DU B B M 25 25 h f BG Ty jh j (kpc) (kpc) (kpc) (J) (J) (J) 11 50 50 52 1 Sb-Sb 3.70 10 9.3 12.1 11.8 0.28 2 10 6 10 10 11 50 49 51 105 Sb-Sb 1.89 10 23.5 28.8 21.2 0.37 2 10 5 10 2 10 11 50 50 50 165 Sa-Sb 2.23 10 10.0 12.9 9.9 0.32 5 10 2 10 5 10 11 49 48 50 201 Sa-Sb 1.05 10 32.5 34.8 22.8 0.48 9 10 10 4 10 11 51 206 Sb-Sb 1.28 10 16.8 13.2 9.0 0.45 2 10 12 51 50 53 237 Sb-Sb 4.20 10 46.3 42.9 35.7 0.23 5 10 4 10 8 10 11 51 243 Sb-SO 3.02 10 37.6 28.1 21.7 0.31 5 10 12 52 272 Sc-Sc 1.11 10 83.5 32.7 23.6 0.39 8 10 11 50 47 50 297 Sc-Sc 1.27 10 29.4 29.9 18.7 0.53 2 10 10 7 10 11 52 439 Sa-SO 4.95 10 29.4 23.4 13.0 0.63 10 11 50 50 52 474 Sa-Sb 4.77 10 21.7 27.1 23.3 0.19 3 10 6 10 10 11 51 47 52 490 Sa-Sa 4.90 10 30.0 29.6 19.6 0.47 2 10 6 10 10 12 52 50 53 507 Sb-Sa 1.60 10 32.2 35.5 25.1 0.41 10 8 10 10 11 48 49 51 516 Sc-Sb 2.59 10 20.0 21.1 20.3 0.05 4 10 9 10 4 10 11 50 50 52 524 Sb-Sc 4.83 10 15.9 17.4 14.6 0.22 4 10 3 10 2 10 10 50 49 50 539 SO-Sa 9.68 10 17.8 25.8 17.0 0.48 2 10 3 10 4 10 2.2. Binary Galaxies Employing w = 1 for the dimensionless structural factor, the radius-mass relation (n = 2/5) [24] 2/5 2/5 1 M 2/5 5/2 5/2 R = M = R + R , 1i 2i g M 2/5 10 1 observed for the galaxies of large mass (M > 10 M ), where M and R = M 1i,2i 1i,2i 1i,2i g 1i,2i are, respectively, the masses and radii of the components of the binary system before transferring the mass, and the relation between the coordinate h and the galaxy masses M and M in the binary system, we write Equations (15) and (16) as 8/5 8/5 2 3 U = a[(1 + h) + (1 h) ] b[1 h ] (21) at r r and m t 3(1 h ) 8/5 8/5 U = a (1 + h) + (1 h) + 2/5 2/5 (1 + h) + (1 h) + (22) 2 2/5 2/5 3 [1 h ][(1 + h) + (1 h) ] at r < r , where m t G M a = , 2/5 5/2 5/2 13/5 2 R + R 1i 2i 5/2 5/2 G M R + R 1i 2i b = , 5 5 128r R R mi 1i 2i 36/5 4 2 2 Gm r i mi b = . 6/5 5/2 5/2 M R + R 1i 2i In order to calculate a and b, the observed values of M, R , R , and r = X are used, 1i 2i mi where X is the projection of the linear distance between the components of the binary Astronomy 2023, 2 67 galaxy from ref. [24]. The potential energy has an extremum at h = h = 0; it is the minimum if a < a = b cr or 22/5 5/2 5/2 25 R + R 1i 2i r < mi 5 5 37/5 2 R R 1i 2i and the maximum if a > a . As readily seen, the extreme points of the potential depend cr only on r and R . Note that for touching binary systems (r R + R ) there is a mi 1i,2i mi 1i 2i symmetric minimum because the condition a < b or 22/5 5/2 5/2 25 R + R 1i 2i > 1 11/5 5/2 5/2 16 4R R 1i 2i ¶U always holds. Expanding the equation = 0 to the third-order terms in h inclusive and ¶h solving it, we find the barrier positions at 2 3 1/2 1/2 37/5 5 5 2 r R R 1875b 600a 6 mi 7 1/2 1/2 1i 2i h = 2 2 41 5 (23) 22/5 1875b + 28a 5/2 5/2 25 R + R 1i 2i 1/2 (a < a ). As a result, h < 2 0.71. So, in a strongly asymmetric binary system with cr 1/2 2 the mass ratio M / M > (1 + 2 ) 6, the galaxies should fly apart. Indeed, the binary 1i 2i galaxies with a high mass ratio are rare objects [24] (Figure 2). 3. Mass Transfer in Close Binary Cosmic Systems as a Source of Energy in the Universe 3.1. Close Binary Stars As assumed, the values of L and M are conserved during the binary star evolution in h. The orbital angular momentum L is found using the observed star masses M and i 1i,2i period P of orbital rotation at h = h . The binary stars differ in values of M and L and, orb,i i i accordingly, have potential energies of a different kind. One can express the potential energy (18) in units of 13/12 13/12 2 3 u = u + u + xv = (1 + h) + (1 h) + x(1 h ) , 1 2 3 47/12 G M R where x = . As seen from Figure 3, at relatively large x > 0.025 the u 2M 3.288L as a function of h has two asymmetric maxima and the minimum at h = 0. The decrease in x causes the change of the shape of potential energy in h: h approaches h = 0 and b m the height of barrier B in h decreases. At x = 0.025, there is an inverse U-type symmetric potential with a maximum at h = 0, and the di-star system becomes unstable with respect to the mass transfer coordinate. The asymmetrization of the system becomes energetically 47/12 2 favorable. The parameter x M /L depends on M and L . For example, the value of x decreases with increasing L or decreasing M. i Astronomy 2023, 2 68 -2.00 0.00 x=0.025 -0.02 -2.02 -0.04 -0.06 -2.04 0.05 -0.08 -0.10 -2.06 -0.12 0.1 -2.08 -0.14 -0.16 -2.10 0.2 -0.18 -0.20 -2.12 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -2.02 x=0.025 -2.04 -2.06 0.05 -2.08 -2.10 0.1 -2.12 -2.14 -2.16 -2.18 0.2 -2.20 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 3. The calculated dimensionless u, u + u , and xv vs. h at indicated x (see text). 1 2 In Figures 4–6, the potential energies (driving potentials) U(h) versus h are presented for the close di-star systems. For the systems considered, a < a or cr 47/24 3 1/2 L < [10.1GR M ] , 2M and the potential energies have symmetric barriers at h = h and a minimum at h = h = 0. As seen in Figure 4, the barrier in h appears as a result of the interplay between the gravitational energy U + U of the stars and the star–star interaction V. Both energies 1 2 have a different trend as a function of h: U + U decreases and V increases with changing 1 2 h from h = 0 to 1. One should note that the driving potentials U(h) for the di-star systems look similar to the driving potentials for the microscopic dinuclear systems [36,37]. The same conclusion was derived in ref. [29] with the h-independent (constant) structural factors w . The evolution of the di-star system depends on its initial mass asymmetry h = h . If the original di-star is asymmetric and jh j < jh j, then it evolves in h to the global i b minimum at h = 0 (a symmetric di-star system). The matter of a heavy star moves to an adjacent light star enforcing the symmetrization of a di-star without additional driving energy. The symmetrization of an asymmetric binary star is accompanied by the decrease of potential energy U or the transformation of the potential energy into the internal energy of 39 41 stars. A huge amount of energy DU 10 10 J is released during the symmetrization (see Figures 5 and 6 and Table 1). So the symmetric di-star is created at large excitation energy. Note that for the binary stars considered, the energy of a single star (jhj = 1) is larger than the energy of a symmetric binary (h = 0). So the merger of stars in a di-star is energetically unfavorable. u + u 1 2 x v Astronomy 2023, 2 69 -1.14 0.0 Cr B -1.15 -0.5 -1.16 -1.17 Cr B -1.0 -1.18 -1.5 -1.19 -1.20 -2.0 -1.21 -1.22 -2.5 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.16 Cr B -1.17 -1.18 -1.19 -1.20 -1.21 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 4. Calculated gravitational energies U + U of the stars; the energy of interaction V between 1 2 them; the total potential energy U versus the mass-asymmetry coordinate h for the close binary star a Cr B. The arrow indicates the value of h for the original binary system. When the symmetrization of a binary occurs, the distance between the stars decreases because of the increase in the gravitational interaction. After a symmetric di-star is formed, it performs zero-point oscillations near h = 0, which cause the radiation of gravitational waves. The transfer of matter in the symmetric di-stars is a mass asymmetry fluctuation near h = 0. If there is a transmission (fluctuation) in one direction, then there is a transmission in the opposite direction, as proposed in ref. [21]. There are indications of that effect [7]. Ifjh j > h or h = 0, the di-star system is unstable with respect to the asymmetrization. i b b The matter is transferred from the light star to the heavy one even without additional exter- nal energy. We only know one close binary system a Cr B (M = 2.58M , M = 0.92M , 1 2 w = 1.30, w = 1.68, b/a = 0.039), for which jh j = 0.47 > h = 0.33 (Figure 4). 1 2 i b Since the fusion barriers B in h are quite large in Table 1 for the systems withjh j < h , h i there is a strong hindrance to the formation of a single star from the di-star system due to thermal diffusion in h. The existence of a barrier in h may be the reason why very asymmetric close di-stars with jh j > h are rarely observed. Indeed, there can not be a i b stable di-star with a very light star, jh j > h . i b A spectacular recent case is KIC 9832227, which was predicted [8] to be merged in 2022, enlightening the sky as a red nova. For this system (h = 0.63, h = 0.84), we conclude i b that a fast merger is excluded (see Figure 5). Instead, the di-star is driven towards mass symmetry. Matter is transferring from a heavy star to a light one, and the relative distance between two stars and the period of the orbital rotation decreases. A huge amount of energy DU 10 J is released during the symmetrization. As seen in Figure 5 and Table 1, the di-stars KIC 9832227 and RR Cen (h = 0.65, h = 0.85) have almost the same h , h , and i b i b potential energy dependencies on h. So the observation of the RR Cen di-star for a possible merger is also desirable. Observational data in ref. [44] refute the prediction for the 2022 red nova merger [8]. U + U (10 J) 1 2 U (10 J) V (10 J) Astronomy 2023, 2 70 -5.4 -7.4 KIC 9832227 RR Cen -5.6 -7.6 -7.8 -5.8 -8.0 -6.0 -8.2 -6.2 -8.4 -8.6 -6.4 -8.8 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -4.6 -4.6 GW Cep -4.8 AW UMa -5.0 -4.7 -5.2 -4.8 -5.4 -4.9 -5.6 -5.0 -5.8 -5.1 -6.0 -6.2 -5.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -5.50 -4.80 V700 Cyg -5.55 EM Lac -4.85 -5.60 -4.90 -5.65 -4.95 -5.70 -5.75 -5.00 -5.80 -5.05 -5.85 -5.10 -5.90 -5.15 -5.95 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -4.6 -6.4 V870 Ara AH Aur -4.8 -6.6 -5.0 -6.8 -5.2 -7.0 -7.2 -5.4 -7.4 -5.6 -7.6 -5.8 -7.8 -6.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 5. Driving potentials U(h) calculated for the close binary stars indicated. The arrows indicate the values of h for the respective initial binary systems. 41 41 41 41 U (10 J) U (10 J) U (10 J) U (10 J) U (10 J) 41 41 41 U (10 J) U (10 J) U (10 J) Astronomy 2023, 2 71 -4.8 -6.3 AW Vir -4.9 QX And -6.4 -5.0 -5.1 -6.5 -5.2 -6.6 -5.3 -5.4 -6.7 -5.5 -6.8 -5.6 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -5.4 VY Cet -4.4 BO Ari -5.5 -4.6 -5.6 -4.8 -5.7 -5.0 -5.8 -5.2 -5.9 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -8.4 -5.6 AP Aur DN Aur -8.6 -8.8 -5.8 -9.0 -6.0 -9.2 -9.4 -6.2 -9.6 -9.8 -6.4 -10.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.15 -4.8 HV Aqr HV UMa -1.20 -5.0 -1.25 -5.2 -1.30 -5.4 -1.35 -5.6 -1.40 -5.8 -1.45 -6.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 6. The same as in Figure 5, but for other close binary stars. For the KIC 9832227, we studied the change of potential energy with the variations of the total mass M and orbital angular momentum L (Figure 7). We take effectively into consideration the losses of mass and angular momentum of the binary star during its evolution. As seen, simultaneous losses of M and L (D M DL ) by 20% or 50% weakly i i influence the features of potential energy. The decrease of M at fixed L affects stronger (the depth of the minimum in U(h) decreases) but the evolution to the global symmetric minimum is still energetically favorable. So the related losses of M and L almost do not influence the symmetrization of the system. 41 41 41 41 U ( 10 J) U (10 J) U (10 J) U (10 J) U (10 J) 41 41 41 U (10 J) U (10 J) U (10 J) Astronomy 2023, 2 72 -4.3 KIC 9832227 -4.3 KIC 9832227 0.8M, 1L -4.4 0.8M, 0.8L -4.4 -4.5 -4.6 -4.5 -4.7 -4.6 -4.8 -4.9 -4.7 -5.0 -4.8 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 KIC 9832227 -2.6 KIC 9832227 -2.56 0.5M, 1L 0.5M, 0.5L -2.7 -2.58 -2.8 -2.60 -2.62 -2.9 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 7. Calculated total potential energies U vs h for the binary star KIC 9832227. The notation (0.8 M, 1 L), (0.5 M, 1 L), (0.8 M, 0.8 L), and (0.5 M, 0.5 L) means that the calculations are performed with the losses of the total mass M and orbital angular momentum L = L by (20%, 0%), (50%, 0%), (20%, 20%), and (50%, 50%), respectively. The arrows show the corresponding initial h for binary stars. 3.2. Close Binary Galaxies In the calculations, the observational data from the catalog of isolated galaxy pairs [24] are employed: the linear diameters of galaxies, A (1i, 2i) = 2R ; the projection of 25 1i,2i the linear distance X between the components of binary galaxies; and the total orbital mass M of pairs. In the close binary galaxies considered, the average relative distance r = X between the galaxies is commensurate with the sum of the radii of these galax- mi ies [24]. Matter transfer between the galaxies in a binary system is closely related to 5/2 5/2 5/2 their radii and relative distance. If M = MR / R + R at n = 2/5, then 1i,2i 1i,2i 1i 2i 5/2 5/2 5/2 5/2 h = R R / R + R . For the systems shown in ref. [24], a < a = b; i cr 1i 2i 1i 2i 8 therefore, the potential energy has a symmetric barrier at h = h and a minimum at h = h = 0. There is a similarity of the driving potentials U(h) for macroscopic binary galaxies and stars to microscopic dinuclear systems [29,36,37]. Di-galaxies differ by r and R , and, correspondingly, by the potential energy mi 1i,2i dependencies. In Figures 8 and 9, the driving potentials U(h) are presented versus h for the close elliptic, spiral, and elliptic-spiral di-galaxy systems. For systems shown, except binary galaxy 272 [24], a < a = b and, thus, there are potential barriers at h = h and cr a minimum at h = h = 0. The barrier in h appears from the interplay between the total gravitational energy U + U of the galaxies and the galaxy–galaxy interaction potential 1 2 V. Note that the driving potentials U(h) for the di-galaxy systems look like the driving potentials for the close binary stars and the microscopic dinuclear systems [29,30,36,37]. 41 41 U (10 J) U (10 J) 41 41 U (10 J) U (10 J) Astronomy 2023, 2 73 -1.15 -3.3 -1.20 -3.4 -1.25 -3.5 -3.6 -1.30 -3.7 -1.35 -3.8 -1.40 -3.9 -4.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.5 -6.4 -6.6 -1.6 -6.8 -7.0 -1.7 -7.2 -7.4 -1.8 -7.6 -7.8 -1.9 -8.0 -8.2 -2.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.15 -1.00 -1.20 -1.05 -1.25 -1.30 -1.10 -1.35 -1.15 -1.40 -1.45 -1.20 -1.50 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 8. Calculated total potential energies U vs. h for the indicated (numbers) close binary elliptic galaxies at n = 2/5. The numbers are serial numbers of binary galaxies in the catalog of isolated pairs of galaxies [24]. The arrows show the corresponding initial h for binary galaxies. An important question is to what extent the calculated results depend on the choice of the parameter n. The dependencies of the potentials calculated with n = 2/3 and n = 2/5 are qualitatively the same, except that the ratio U(h = h )/U(h = 0) increases with n. So all binary systems considered evolve to the symmetric configurations independently of n. The evolution path of a close binary galaxy depends on its initial mass asymmetry h = h . If the original di-galaxy is asymmetric with jh j < h , then its evolution in h to the i i b global minimum at h = 0 is energetically favorable, that is, to form a symmetric di-galaxy system. The matter of a heavy galaxy can move freely to a light one enforcing symmetriza- tion of di-galaxy without additional energy. This symmetrization leads to the decrease in potential energy U, thus transforming the potential energy into the internal kinetic energy. For example, for the close elliptic binary galaxies 194 (h = 0.28), 279 (h = 0.50), and i i 554 (h = 0.25) [24], the internal energies of galaxies increase during symmetrization by the 50 51 50 amount DU = U(h ) U(h = 0) = 2 10 , 3 10 , and 2 10 J, respectively (see Figures 8 and 9, and Tables 2–4). As found, most of the close binary galaxies are asymmetric and their symmetrization leads to the release of a large amount of energy (see Tables 2–4). 53 52 52 U (10 J) U (10 J) U (10 J) 52 55 52 U (10 J) U (10 J) U (10 J) Astronomy 2023, 2 74 -2.8 -6.0 -6.2 -3.0 -6.4 -6.6 -3.2 -6.8 -7.0 -3.4 -7.2 -7.4 -3.6 -7.6 -7.8 -3.8 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -2.0 -2.2 -4.0 -2.4 -4.2 -2.6 -2.8 -4.4 -3.0 -3.2 -4.6 -3.4 -3.6 -4.8 -3.8 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -4.0 -7.4 -7.6 -4.2 -7.8 -4.4 -8.0 -8.2 -4.6 -8.4 -4.8 -8.6 -8.8 -5.0 -9.0 -9.2 -5.2 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -2.2 -2.0 -2.3 -2.1 -2.4 -2.5 -2.2 -2.6 -2.3 -2.7 -2.4 -2.8 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 9. The same as in Figure 8 but for the close binary elliptic-spiral and spiral galaxies indicated. If jh j > h or h = 0, the di-galaxy system is unstable to evolve towards a more i b b asymmetric system; thus, one galaxy “swallows” another. The matter flows freely from a light galaxy to a heavy one. This process leads to energy release. Representative examples of this mode of evolution are three close spiral binary galaxies, namely the galaxies 206 (jh j = 0.45, h = 0.34), 243 (jh j = 0.31, h = 0.07), and 439 (jh j = 0.63, h = 0.56) [24], i b i b i b respectively, for which jh j > h (Table 4). i b The binary galaxy withjh j < h can merge only by overcoming the barrier at h = +h i b b or h = h . With decreasing ratio a/b, the value of this barrier B = U(h ) U(h ) b h b i 51 52 51 51 U (10 J) U (10 J) U (10 J) U (10 J) 51 51 52 52 U (10 J) U (10 J) U (10 J) U (10 J) Astronomy 2023, 2 75 increases. Since the values of B are quite large for the systems with jh j < h in Tables 2–4, h i asymmetrization of the di-star system due to thermal diffusion in h is strongly suppressed. Indeed, very asymmetric close binary galaxies with jh j > h are rarely observed. This imposes restrictions on the asymmetric configurations with jhj > h of the di-galaxy systems. There can not be a stable di-galactical system with a very light galaxy. If the value of a becomes larger than a , the minimum in U(h) disappears, the inverted cr U-shaped potential appears with a maximum at h = 0, and the di-galaxy asymmetrization breaks the system finally apart. Hence, with increasing mass asymmetry, two constituents of the binary galaxy move away from each other to form two separate galaxies. The spiral binary galaxy 272 (see Table 4) is a good candidate to trace such evolution. As found in our calculations, the double galaxy undergos either symmetrization or asymmetrization, depending on the original mass asymmetry. Asymmetrization may cause a merger process, but if the critical mass asymmetry is reached, the system is driven apart and forms two isolated galaxies with a changed mass ratio. This scenario can be considered an incomplete merger. So the source of expansion of a binary galaxy is the transfer of matter from a lighter component to a heavier one. A necessary and sufficient condition for this is the fulfillment of the inequality 22/5 5/2 5/2 0.15 R + R 1i 2i r > . mi 5 5 R R 1i 2i The mechanism presented can be generalized for multiple galaxies, groups of galaxies, and galaxy associations. 4. Origin Of Orbital Period Change in Contact Binary Stars As found, the orbital period of W-type overcontact binary GW Cep (q = M / M = 0.37, 1i 2i h = 0.46) decreases with time [39–41]. For the overcontact binaries VY Cet (q = 0.67, h = 0.20) and V700 Cyg (q = 0.65, h = 0.21), cyclic oscillations were found to be i i superimposed on an increase in the secular period. This effect was explained either by strong external perturbation, i.e., by a close-by third object, or by the magnetic activity cycles. For the EM Lac (q = 0.63, h = 0.23) and AW Vir (q = 0.76, h = 0.14) binaries, i i the periods demonstrate secular increase. As concluded, the period variations in a W UMa-type binary star correlate with the mass ratio q and the mass M of the primary 1i component [39–41]. The lower the mass ratio q in binaries, the shorter the period. 2pmr 2p m Using the relation between L and r , we obtain the period P = = of orb w L orb orbital rotation with frequency w . At r > r and r r as orb m t m t 1/2 > m P = 2p (24) orb G M and 1/2 < t P = 2p , (25) orb G M respectively. As seen, at initial jh j < jh j and r > r (r r ), the system moves towards i b m t m t the symmetric configuration and, correspondingly, h decreases, r decreases (r increases), m t > < and finally, P decreases (P increases) [31,32]. orb orb For the KIC 9832227 system (jh j = 0.63), matter is transferred from a heavy star to a light one, the relative distance r > r between two stars and the period P of the orbital m t orb rotation decrease. The evolution in h pushes the system to the touching configuration (r = r ) at some critical mass asymmetry jhj = jh j 0.45 (Figure 10). Further evolution m t t in h leads to a configuration with partial overlap (r < r ) of stars. So at jhj jh j the m t t 3/2 < < period P slightly increases because P r and r increases with decreasing h. Thus, orb orb t Astronomy 2023, 2 76 if the system reaches the point h = h and partial overlap, the period abruptly changes. A similar period behavior is shown in Figure 10 for the other contact binaries considered. 2.2 1.75 KIC 9832227 1.70 GW Cep 2.0 1.65 1.8 1.60 1.6 1.55 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 2.0 2.4 RR Cen 1.9 2.3 QX And 1.8 2.2 1.7 2.1 1.6 2.0 1.5 1.9 1.4 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 2.0 2.0 AW UMa V870 Ara 1.8 1.8 1.6 1.6 1.4 1.4 1.2 1.2 1.0 1.0 0.8 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 10. The calculated relative r = R (solid lines) and touching r = R = R + R (dashed m m t t 1 2 lines) distances between the components of di-stars indicated in the units of the Sun radius R = R sun vs. h. For the binary GW Cep (jh j = 0.46), jh j > jh j, the system moves towards the i i t symmetry and the orbital period decreases. For almost symmetric EM Lac (jh j = 0.23) and AW Vir (jh j = 0.14) binaries, jh j < jh j and their periods increase. As found, the i i variations in the period of W UMa-type binary star correlate with the evolution towards a global symmetric minimum. At the low mass ratio q or the large h , the binaries usually demonstrate a decreasing period because jh j > jh j (r > r ), while the periods in systems i t m t with larger q (small h ) increase because jh j < jh j (r < r ). t m t i i 5. Stability of Macroscopic Binary Systems The Regge theory turned out to be very influential in the development of elementary- particle physics [45–49]. As shown in the application of Regge’s ideas to astrophysics, the 4/3 spins of planets and stars are well described by the Regge-like law for a sphere (S M ), R /R , R /R R /R , R /R R /R , R /R m sun t sun m sun t sun m sun t sun R /R , R /R R /R , R /R R /R , R /R m sun t sun m sun t sun m sun t sun Astronomy 2023, 2 77 while the spins of galaxies and clusters of galaxies obey the Regge-like law for a disk 3/2 (S M ) [25,34,50–53]. Unlike semi-phenomenological approaches, these expressions contain only fundamental constants and do not depend on any established empirical quantities [25,34,50–53]. In Refs. [25,34,50–53], a cosmic analog of the Chew–Frautschi plot with two important cosmological Eddington and Chandrasekhar points was also constructed. The Darwin instability effect can be studied in a binary star or galaxy by using the well-known model of Refs. [25,34,50–53] based on the Regge theory. The total angular mo- mentum J of a binary system is the sum of the orbital angular momentum L and the spins tot S (k = 1, 2) of individual components: J = L + S + S . (26) tot 1 2 The values of J and S are expressed using the Regge-like law for stars and planets (l = 3) tot i or galaxies (l = 2): (1+l)/l J = h ¯ (27) tot and (1+l)/l S = h ¯ , (28) where h ¯ , m , M (k =1,2), and M = M + M are the Planck constant, masses of proton k 1 2 and astrophysical objects (planets, stars, or galaxies), and the total mass of the system, respectively. The maximum (L and S + S are antiparallel) and minimum (L and S + S 1 2 1 2 are parallel) orbital angular momenta are L = J + S + S (29) max tot 1 2 and L = J S S , (30) min tot 1 2 respectively. Using the coordinate h instead of masses M and M [33] and Equations (28)–(30), 1 2 we derive (1+l)/l (1+l)/l S + S (1 + h) + (1 h) 1 2 = , (31) (1+l)/l (1+l)/l (1+l)/l 2 (1 + h) (1 h) min (1+l)/l (1+l)/l S + S (1 + h) + (1 h) 1 2 = . (32) (1+l)/l (1+l)/l (1+l)/l max 2 + (1 + h) + (1 h) At h = 0 (the symmetric binary system), we have S + S 1 1 2 = > 1 (33) 1/l L 2 1 min and S + S 1 1 1 2 = > . 1/l L 2 + 1 3 max For the symmetric ( h = 0) binary star (planet) and binary galaxy, (S + S )/L 0.44 1 2 max and 0.41, respectively. At h = 1, we have S + S 1 2 ! ¥ min Astronomy 2023, 2 78 and S + S 1 1 2 = . L 2 max As follows from these expressions, for very asymmetric binaries, the ratios (S + S )/L 2 max,min are almost independent of l. According to ref. [20], the Darwin instability occurs when the binary mass ratio is very small [q = M / M < 0.1] or the mass asymmetry is very 2 1 large [h = (1 q)/(1 + q) > 0.82]. As seen in Figure 11, the ratios (S + S )/L and 1 2 max (S + S )/L continuously increase with jhj. Since these ratios are larger than 1/3, all 1 2 min possible binary stars (planets) or binary galaxies, independently of their mass asymmetry h, have the Darwin instability (S + S L) [19] and, thus, should merge. However, 1 2 the observations do not confirm this, which probably means that there is no the Darwin instability effect in such binary systems and, accordingly, the mechanism of merger has a different origin and should be revealed. Since the spins of planets, stars, galaxies, and clusters of galaxies are well described by the Regge theory [25,34,50–53], we can be sure of the correctness of this conclusion. (b) (a) 0.50 0.49 0.48 0.47 0.46 0.45 0.44 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 2.8 (c) 2.4 2.0 1.6 1.2 0.8 0.4 0.0 -0.4 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 Figure 11. The calculated ratios (S + S )/L (a), (S + S )/L (b), (S S )/L [solid line] (c), 1 2 max 1 2 min 1 2 2 and (S S )/L [dashed line] (c) as functions of mass asymmetry at l = 3. In the cases of antiparallel 2 1 1 spins (c), L = J + S S and L = J S + S . 1 tot 1 2 2 tot 1 2 Note that in the cases of antiparallel spins with L = J + S S and tot 2 1 1 L = J S + S (Figure 11), the ratiosjS S j/L andjS S j/L are larger than 1/3 2 tot 1 2 2 1 1 1 2 2 in the asymmetric binaries with jhj 0.5. As seen in Figure 12, the dependencies of L , L , L , and L on h are different. max min 1 2 The matter transfer can increase or decrease the orbital angular momentum. For example, at h ! 0 the binary system has smaller L = L . The observations of L versus h may be max useful to distinguish the difference between the orientations of orbital and spin momenta. (S +S )/L 1 2 max (S -S )/L , (S -S )/L 1 2 1 2 1 2 (S +S )/L 1 2 min Astronomy 2023, 2 79 2.00 (b) (a) 0.20 1.96 0.16 1.92 0.12 1.88 0.08 1.84 0.04 1.80 1.76 0.00 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 2.0 (c) 1.8 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 0.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 Figure 12. The calculated ratios L / J (a), L / J (b), L / J [solid line] (c), and L / J max tot tot tot tot min 2 1 [dashed line] (c) as functions of mass asymmetry at l = 3. Employing Equations (26)–(28) and the results of refs. [29,31,32], we obtain new analytical formulas for the relative distance between the components of the binary " # (1+l)/l (1+l)/l (1+l)/l 2 2 ML h ¯ M M M M 1 2 r = = + e + e 1 2 2 2 2 2 G M M G M M m m m p p p 1 2 1 2 h i (1+l)/l (1+l)/l (1+l)/l 2 (2 l)/l 2 + e (1 + h) + e (1 h) 1 2 2h ¯ M = (34) 3 2 2m 2 Gm p [1 h ] at r r and m t ! " # 1/2 1/4 (1+l)/l (1+l)/l (1+l)/l h i 3/4 h ¯ Mg M M M 1 2 l l r = M + M + e + e m 1 2 1 2 2 2 m m m G M M p p p 1 2 1/4 3n 3 3 (2+[3n 1]l)/l 2h ¯ m g p M n n 3/4 = (1 + h) + (1 h) [ ] G 2m h i 1/2 (1+l)/l (1+l)/l (1+l)/l 2 + e (1 + h) + e (1 h) 1 2 (35) 1/2 [1 h ] at r < r . One can similarly derive the formulas for the orbital rotation period m t " # 1/2 (1+l)/l (1+l)/l (1+l)/l 3 3 r 2ph ¯ M M M M m 1 P = 2p = + e + e orb 1 2 3 3 G M m m m G M M p p p 1 2 h i (1+l)/l (1+l)/l (1+l)/l (3 2l)/l 3 2 + e (1 + h) + e (1 h) 1 2 4ph ¯ M = (36) 5 3 2m 2 G m p p [1 h ] L /J max tot L /J , L /J 1 tot 2 tot L /J min tot Astronomy 2023, 2 80 at r r and m t 1/2 1/2 3 3 n n r g M + M t 1 2 P = 2p = 2p orb G M G M 1/2 3 3n 1 (3n 1)/2 g m n n 3/2 = 2p [(1 + h) + (1 h) ] (37) 2G 2m at r < r [31,32]. The values e = e = 1 and e = e = 1 correspond to the cases m t 2 2 1 1 of antiparallel and parallel, respectively, orbital and spin angular momenta. In the cases of antiparallel spins, e = e = 1 and e = e = 1. The observation data result in 1 2 1 2 the relationship R = M /g between the radius and mass of the star, where n = and k 3 2/3 g = M /R [6] and the galaxy, where the value of n depending on mass is in the interval 2 2 , [24,30]. As seen in Figure 13, at r > r , m t 5 3 h i (1+l)/l (1+l)/l (1+l)/l 2 + e (1 + h) + e (1 h) 1 2 d d r = R = m m [1 h ] or r decreases with jhj and, finally, h i (1+l)/l (1+l)/l (1+l)/l 2 + e (1 + h) + e (1 h) 1 2 P = orb [1 h ] or P decreases. At r > r (Figure 13), the dependence of P as a function of mass asym- m t orb orb metry has similar behavior at antiparallel and parallel orbital and spin angular momenta. In the case of antiparallel spins with L = J + S S (L = J S + S ) and M M , tot 2 2 tot 2 2 1 1 1 1 the values of r and P decrease (increase ) with h decreasing from 1 to 0 (Figure 14). At m orb r r (Figure 15), the value of m t d n n 3/2 P = (1 + h) + (1 h) [ ] orb or P increases with decreasingjhj and does not depend on mutual orientations of orbital orb momentum and spins. In contrast, the distance h i 1/2 (1+l)/l (1+l)/l (1+l)/l 2 + e (1 + h) + e (1 h) 1 2 3/4 d n n R = [(1 + h) + (1 h) ] 1/2 [1 h ] or r depends on the orientations of orbital momentum and spins. As seen in Figure 15, the dependencies of R on h at different spins have various behaviors. m Astronomy 2023, 2 81 (b) (a) 0.40 0.36 0.32 0.28 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 0.19 (c) (d) 0.18 0.17 0.16 0.15 0.14 100 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 d d Figure 13. The calculated dimensionless relative distances r = R (a,b) and orbital rotation m m periods P (c,d) as functions of h at r > r (r = R , r = R ) and l = 3. The plots (a,c) and m t m m t t orb (b,d) correspond to the systems with the parallel and antiparallel orbital angular momentum and spin. (b) (a) 600 R >R +R R >R +R m 1 2 m 1 2 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 d d Figure 14. The calculated dimensionless relative distance r = R (a) and orbital rotation period m m P (b) as functions of h in the cases of antiparallel spins with L = J + S S (dashed line) and 1 tot 1 2 orb L = J S + S (solid line), l = 3, and r > r (r = R , r = R ). tot m t m m t t 2 1 2 d d P R orb m d d P R orb orb m Astronomy 2023, 2 82 (b) (a) 1.216 1.214 1.212 1.210 1.208 1.206 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 5.0 (c) (d) 2.8 4.5 4.0 2.7 3.5 2.6 3.0 2.5 2.5 2.4 2.0 2.3 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 d d d Figure 15. The calculated dimensionless r = R (a–c) and P (d) as functions of h at r r m t m m orb (r = R , r = R ), l = 3, and n = 2/3. The plots (a,b) correspond to the system with the parallel and m m t t antiparallel orbital angular momentum and spin, respectively. In the plot (c), the cases of antiparallel spins with L = J + S S (dashed line) and L = J S + S (solid line) are presented. Note tot tot 1 1 2 2 1 2 that the value of P does not depend on the orientations of the orbital momentum and spins. orb 6. Possibility of Formation of Binary Cosmic System From Single Cosmic Object 6.1. Close Binary Stars For transition from a single star with mass M and potential energy U(h = 1) to an asymmetric or almost symmetric binary star with jhj < jh j, the system should overcome the height barrier B = U(jhj = h ) U(jhj = 1) at h = h . The values of B are too f b b f large (see Figures 5 and 6 and Table 1) to allow the single star to be transformed into a binary star with jhj < jh j due to thermal fluctuations [35]. Thus, the fission in h (matter transfer from a heavier star to a lighter one) cannot be the origin of binary stars. Our conclusion is consistent with that of ref. [26,28]. If the fission hypothesis mentioned is not valid, then the only possible conclusion can be the assumption of a common origin of the components of a binary star from a pre-stellar state of matter [28]. Thus, some stars likely arise during the formation of stellar groups in the form of pairs, triplets, quartets, etc. The observed excess of binary stars in the stellar associations suggests that these binary stars transform into ordinary binary stars after leaving the associations [28]. 6.2. Close Binary Galaxies In the galactic systems from Tables 2–4, and Figures 8 and 9, the fission barriers B in h are too high to form almost symmetric binary and asymmetric galaxies withjhj < jh j from the single galaxy via thermal fluctuations of h [35]. So we arrive at the same conclusion as in the case of fission of single stars. 7. On Evolution of Compact Binary Black Holes The general view is that compact binary systems consisting of white dwarfs, neutron stars, and black holes eventually merge as a result of gravitational wave radiation, as recently observed with the LIGO and Virgo interferometers [14]. The process of merging is still insufficiently studied. We are going to add a few aspects highlighting the role of d d R R m m orb Astronomy 2023, 2 83 matter transfer in binary black holes (BBH) [34], which are certainly a special class of objects. The presence of the event horizon restrains the flow of mass between two cores of BBH. At a speculative level, quantum tunneling as the Hawking radiation [54] may play a faint, most likely insignificant role due to not yet received firm confirmation. However, a BBH is usually embedded in a cloud of remnant matter left over from the progenitor stars. Therefore, a BBH should always be considered together with the surrounding accretion disk. The material of that disk may serve as the matter transfer between two cores by the so-called sloshing effect, as illustrated by the hydrodynamical calculations in ref. [55]. In order to understand the possibility of matter exchange in a BBH system, we consider the potential energy of a BBH as a function of h. The possibility of matter transfer based on the potential energy of a BBH system should always be understood together with its accretion disk. In addition to the BBH’s own extremely strong gravitational fields, the potential energy contains the interaction potential between two black holes. Being well aware of the strong gravitational fields, we use Newtonian mechanics because, at a distance, post-Newtonian effects can alter the binary potential by no more than 25% [55], which will not affect the overall BBH properties. General relativity becomes of course, essential at separations close to touching. The total potential energy U = U + U + V + V + V (38) 1 2 1R 2R of a BBH is given by the sum of the potential U (k = 1, 2) and rotational energies V k kR (k = 1, 2) of two nonzero spin black holes, and the black hole–black hole interaction potential V. The energy of the black hole “k” is G M U = , (39) where G, M , and R are the gravitational constant, mass, and radius of the black hole, k k respectively. The rotation energy of a black hole is S W M c k k k V = = , (40) kR 2 2 where S = M R c and W = c/R are the spin and rotational frequency of the black hole k k k k k “k”, respectively. The radius of the event horizon (distance from the gravitating mass M at which the particle velocity is equal to c) [25,34,50–53] G M R = (41) is derived from the energy conservation law for a particle with mass m: G M m/R + mc = 0, k k 2 2 2 where mc is the sum of the kinetic mc /2 and rotation mc /2 energies. The derivation of the radius R is based on classical mechanics and the Newtonian law of gravity. In ref. [56], the expression R = 2G M /c was also obtained for a non-rotating black hole. k k Using Equations (39) and (41), we obtain U = M c . (42) k k Due to the rotation of two black holes around the common center of mass, the black hole–black hole interaction potential V(r) contains, together with the gravitational potential V , the kinetic energy of orbital rotation V : G R 2 2 L mv V(r) = V + V = V + = V + , (43) G R G G 2mr 2 Astronomy 2023, 2 84 1/2 where L, v = (G M[2/r 1/r ]) and r are the orbital angular momentum of a BBH, r m m the speed and the semi-major axis of an elliptical relative orbit, respectively. At r r , the black hole–black hole interaction potential is defined as in Equation (8). From the 2 2 conditions ¶V /¶rj = 0 and ¶ V /¶r j > 0, we find the relative equilibrium distance r=r r=r m m between two black holes corresponding to the minimum of V (see Equation (10)). The total angular momentum J = G M /c of a BBH [25,34,50–53] is assumed to be conserved during the conservative matter transfer and the orbital angular momentum of a BBH is 2 2 2 G( M M M ) 1 2 L = J S S = . (44) 1 2 Here the orbital momentum and spins are parallel and the value of L is minimal. Employing Equations (10) and (44), we derive the expression 4G M r = . (45) As seen, the larger M, the larger r . Since r /(R + R ) = 4, the separation of the m m 1 2 components increases at M ! M . From Equations (8), (43) and (45) we obtain the 1 2 simple formula 2 2 G M M M M c mc 2 2 1 1 V(r ) = = = (46) 2r 8M 8 for the interaction potential. The value of V(r ) depends only on the reduced mass m and the velocity of light. So using Equations (40), (42) and (46), we derive the final expression 2 2 Mc M M Mc m 1 2 U = 1 + = 1 + (47) 2 4M 2 4M for the total potential energy of a BBH. 1/2 For the BBH considered, v (r ) = (G M/r ) = c/2 and with good accuracy, one r m m can disregard relativistic effects and use the Newtonian law of gravity. Employing the coordinate h instead of masses M and M , we rewrite expressions (44) 1 2 and (47) for the orbital angular momentum G M L = 1 h (48) 2c and the total potential energy 2 2 Mc 1 h U = 1 + . (49) 2 16 Since the solution of the equation ¶U Mc = h = 0 ¶h 16 leads to h = h = 0 and 2 2 ¶ U Mc j = > 0, h=h ¶h 16 the potential landscape has a global minimum at h = h = 0 for an arbitrary total mass M of a BBH (Figure 16). So the transfer of matter between the black holes in the BBH is energetically favorable and can occur to reach a global minimum. The initial asymmetric system is easily driven to the symmetric BBH. This conclusion does not depend on the Astronomy 2023, 2 85 choice of parameters. The losses of the total mass and orbital angular momentum do not affect symmetrization of a BBH. The transfer of matter between two black holes becomes possible because they interact with their own extreme gravitational fields. However, the evolution of a BBH in h depends also on the parameter of inertia in this coordinate. Since the surfaces of two black holes are spaced R + R < r from each other, the parameter 2 m of inertia in h is expected to be very large and, accordingly, prohibits symmetrization of a BBH. -0.500 -0.505 -0.510 -0.515 -0.520 -0.525 -0.530 -0.535 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 16. Dependence of potential energy of the BBH (49) on h. Asymmetrization (the transfer of matter from a lighter component to a heavy one) of a BBH considered is energetically unfavorable, and, correspondingly, the merger channel in h is strongly suppressed. Thus, the question of the mechanism of the merger of two black holes and the origin of gravitational waves remains open [14]. At antiparallel orbital momentum and spins, the value of 2 2 2 G( M + M + M ) 1 2 L = J + S + S = (50) 1 2 is maximal. For a BBH with this L, 2 2 3 Mc (1 h ) U = 1 + (51) 2 2 2 16(3 + h ) 1/2 2 2 and r /(R + R ) 36, v(r ) = (G M/r ) c/6, h = 0, ¶ U/¶h j > 0, and m m m m h=h 1 2 m all conclusions given above are also valid in this case. Comparing Figures 16 and 17, we can draw the same conclusions in the case when black holes in a BBH spin in the opposite directions. U/(Mc ) Astronomy 2023, 2 86 -0.5000 -0.5005 -0.5010 -0.5015 -0.5020 -0.5025 -0.5030 -0.5035 -0.5040 -1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 Figure 17. Dependence of potential energy of the BBH (51) on h. 8. Summary The evolution of isolated binary stellar and galactic systems was considered in the framework of the approach originally formulated for the study of nuclear fusion reactions, thus combining scales that differ by many orders of magnitude. Despite the differences in details, the common aspects of macroscopic objects (di-star, di-galaxy, BBH) and mi- croscopic dinuclear systems prevail. Both types of systems evolve along well-defined trajectories in the classical phase space. The conserved quantity is the total energy, that is, the Hamilton function, of a binary system. Using general arguments, we have shown that the conservation of energy is sufficient to fix the trajectory of a binary system in the landscape of potential energy determined by the masses of the objects and their interaction. Exploiting the stationarity of the total energy, the stability conditions were derived and investigated as a function of mass asymmetry. We emphasized that the interpretation of mass asymmetry as a collective coordinate has been successfully used to describe nuclear reactions, the cluster structure of nuclei, the decay of dinuclear systems, and the fusion of two heavy nuclei [36,37]. Here we have shown that this collective degree of freedom plays a comparably important role in space objects. In close di-star or di-galaxy systems, the coordinate h can govern the merger and symmetrization (due to the matter transfer) processes. An interesting aspect is that once h is determined, for example, by observation, it is possible to draw a conclusion about the stellar or galactic structure. The new theoretical interpretation of a di-star or di-galaxy system is based on the fact that, after the formation, the lifetime of a binary system is long enough to reach equilibrium conditions in the mass asymmetry coordinate. Therefore, we could conclude that the system will be included in the sample of all binary and single configurations with the probability depending on the potential energy of this configuration. For the systems considered, a < a and there are potential barriers at h = h and minimum at h = h = 0. So cr m two distinct evolution scenarios arise: if jh j < h , the system is driven to the symmetric i b configuration (towards a global minimum of the potential landscape). However, ifjh j > h , i b the binary system evolves towards the mono-object system. All the considered asymmetric close binary systems, with the exception of the di-star a Cr B and di-galaxies 206, 243, 272, 439, satisfy the condition jh j < h and a symmetrization process occurs in these systems. The loss of the total mass and orbital angular momentum has little effect on the process of symmetrization. The merger of the binary stars, including the KIC 9832227, and binary galaxies considered is energetically unfavorable. The formation of a single object from the binary system (jh j < h ) by thermal diffusion in the mass asymmetry coordinate is i b strongly suppressed. Symmetrization of a binary system (jh j < h ) due to the matter transfer is one of the i b important sources of conversion of gravitational energy into other types of energy in the universe. For example, in the cases of compact binary stars and compact binary galaxies, U/(Mc ) Astronomy 2023, 2 87 39 41 48 52 the released energies are about 10 10 and 10 10 J, respectively. Symmetrization of a binary system will lead to close objects with equal masses, temperatures, luminosities, and radii, which are observable quantities. During the process of mass-symmetrization, the channels of binary decay in the relative coordinate are closed. The central result 1/2 of the approach reviewed is that stable binary systems exist only at jhj < h = 2 1/2 2 (M / M < (1 + 2 ) ' 6) because otherwise, the stars are getting closer to merging. 1 2 Thus, di-stars or di-galaxies with jhj > h can not exist for long enough. Indeed, binary systems with a large ratio M / M are very rare objects in the universe. 1 2 Asymmetrization (the transfer of matter from a lighter component to a heavy one) is equivalent to incomplete merging. Asymmetrization is also the origin of the expansion of a binary galaxy. The separation of components from each other is represented as an analog of the expansion of the universe within a binary system. The condition under which the asymmetrization process is realized depends mainly on the relative distance between the components of binary systems and their linear dimensions. For contact di-stars, changes in the orbital period can be well explained by evolution in h towards symmetry. We predicted that the decrease and increase in orbital periods are associated, respectively, with the non-overlapping (jh j > jh j, r R + R ) and i t m 1 2 overlapping (jh j < jh j, r < R + R ) stages of a binary star during its symmetrization. t m i 1 2 Thus, observing the change of periods allows us to distinguish between these two stages of a di-star. Based on calculations of potential energy, we have demonstrated that the mass asym- metry coordinate (matter transfer) is useful for analyzing the transformation of a single star (galaxy) into a binary star (galaxy). Mass ejection from a single star could create a binary star in the mass asymmetry coordinate. However, the barrier preventing symmetrization in h is quite high. So the formation of asymmetric and almost symmetric di-star or di-galaxy systems from the respective mono-star or mono-galaxy by thermal overcoming of barriers in the driving potential (the diffusion process in h) is hardly possible. Binary stars (galaxies) with a mass ratio M / M > 6 evolve to more asymmetric configurations and are rather 1 2 unstable with respect to the decay in the relative distance. Thus, matter ejection from a single star is either captured back or escapes, which is consistent with the findings of refs. [26,28]. Thus, it is possible to assume a common origin of the components of a binary star (galaxy) from the pre-stellar (pre-galactic) state of matter. Based on the Regge-like laws [25,34,50–53], we demonstrated that all possible binary stars (planets) or binary galaxies, regardless of their mass asymmetry, satisfy the Darwin instability condition (S + S L), which contradicts observations. This output is not 1 2 sensitive to model parameters. Therefore, we should look for another mechanism that triggers the merger of binary contact components. Employing the Regge-like laws, we derived new analytical formulas for the relative distance and orbital rotation period of the binary system, which depend on the total mass M, mass asymmetry h, and the fundamental constants G, h ¯ , and m . Employing Newton’s law of gravity and considering the BBH potential energy as a function of h, we have shown the possibility of transferring matter between black holes in a BBH. The evolution of an asymmetric BBH to a symmetric one (h = h = 0) is energetically favorable. Although black holes have their own strong gravitational fields, the transfer of matter in BBH occurs due to the energy of interaction between two black holes. A BBH does not send any signals during its evolution. Perhaps the result of this evolution can be indirectly observed. Symmetrization of a BBH leads to a decrease of U, thus converting the potential energy into internal kinetic energy. For example, for the BBH 4M + 2M (h = 0.33) and 36M + 29M (h = 0.11), the internal energies of black holes increase during symmetriza- 2 2 47 tion by the amount DU = U(h ) U(h = 0) = Mc h 10 J. So the BBH is also the source of thermal energy. The transfer of matter from a lighter component to a heavier one ( the merger of black holes in a BBH) is not an energetically advantageous process, and, thus, the question of the origin of gravitational waves remains open. Astronomy 2023, 2 88 Within the presented model, we can also perform dynamic calculations of the evolution of a binary system in h. For example, we can calculate the relaxation time (symmetrization) and the asymmetrization time. In addition to the total potential energy, it is also necessary to calculate the mass parameter and the coefficient of friction for this coordinate. However, this extension of the model is the subject of future research. Acknowledgments: We are grateful to Alexis Diaz-Torres and Avas Khugaev for fruitful discussions. Author Contributions: All authors contributed equally to this investigation. All authors have read and agreed to the published version of the manuscript. Funding: This work was partially supported by the Alexander von Humboldt-Stiftung (Bonn). Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: Not applicable. Conflicts of Interest: The authors declare no conflict of interest. References 1. Kopal, Z. Close Binary Systems; Shapman and Hall Ltd.: London, UK, 1978. 2. Shore, S.N.; Livio, M.; van den Heuvel, E.P.J. Interacting Binaries; Springer: Berlin/Heidelberg, Germany, 1994. 3. Hilditch, R.W. An Introduction to Close Binary Stars; Cambridge University Press: Cambridge, UK, 2001. 4. Boyarchuk, A.A. Mass Transfer in Close Binary Stars; Taylor and Francis: London, UK; New York, NY, USA, 2002. 5. Eggleton, P.P. Evolutionary Processes in Binary and Multiple Stars; Cambridge University Press: Cambridge, UK, 2006. 6. Vasiliev, B.V. Astrophysics and Astronomical Measurement Data. Available online: http://astro07.narod.ru (accessed on 1 January 2022 ). 7. Cherepashchuk, A.M. Close Binary Stars, Volumes I and II; Fizmatlit: Moscow, Russia, 2013. (In Russian) 8. Molnar, L.A.; Van Noord, D.M.; Kinemuchi, K.; Smolinski, J.P.; Alexander, C.E.; Cook, E.M.; Jang, B.; Kobulnicky, H.A.; Spedden, C.J.; Steenwyk, S.D. Prediction of a red nova outburst in kic 9832227. arXiv 2017, arXiv:1704.05502. 9. Kulkarni, S.R.; Ofek, E.O.; Rau, A.; Cenko, S.B.; Soderberg, A.M.; Fox, D.B.; Gal-Yam, A.; Capak, P.L.; Moon, D.S.; Li, W.; et al. An unusually brilliant transient in the galaxy M85. Nature 2007, 447, 458. [CrossRef] [PubMed] 10. Tylenda, R.; Hajduk, M.; Kaminski, ´ T.; Udalski, A.; Soszynski, ´ I.; Szyman ´ ski, M.K.; Kubiak, M.; Pietrzynski, ´ G.; Poleski, R.; Ulaczyk, K. V1309 Scorpii: merger of a contact binary. Astron. Astrophys. 2011, 528, A114. [CrossRef] 11. Soker, N.; Tylenda, R. Main-Sequence Stellar Eruption Model for V838 Monocerotis. APJL 2003, 582, L105. [CrossRef] 12. Sana, H.; de Mink, S.E.; de Koter, A.; Langer, N.; Evans, C.J.; Gieles, M.; Gosset, E.; Izzard, R.G.; Le Bouquin, J.-B.; Schneider, F.R.N.; et al. Binary Interaction Dominates the Evolution of Massive Stars. Science 2012, 337, 444. [CrossRef] 13. Kochanek, C.S.; Adams, S.M.; Belczynski, K. Stellar mergers are common. MNRAS 2014 , 443, 1319. [CrossRef] 14. Abbott, B.P.; LIGO Scientific and Virgo Collaboration. Erratum: Tests of general relativity with gw150914. Phys. Rev. Lett. 2016, 116, 241102. [CrossRef] [PubMed] 15. Pejcha, O.; Metzger, B.D.; Tyles, J.G.; Tomida, K. Pre-explosion Spiral Mass Loss of a Binary Star Merger. Astrophys. J. 2017, 850, 59. [CrossRef] 16. Eggleton, P.P. Formation and Evolution of Contact Binaries. Available online: https://koreascience.kr/article/JAKO20121855965 6013.pdf (accessed on 1 January 2022). 17. Fabrycky, D.; Tremaine, S. Shrinking Binary and Planetary Orbits by Kozai Cycles with Tidal Friction. Astrophys. J. 2007, 669, 1298. [CrossRef] 18. Tokovinin, A.; Thomas, S.; Sterzik, M.; Udry, S. Tertiary companions to close spectroscopic binaries. Astron. Astrophys. 2006, 450, 681. [CrossRef] 19. Darwin, G.H. The determination of the secular effects of tidal friction by a graphical method. Proc. R. Soc. 1879, 29, 168. 20. Rasio, F.A. The Minimum Mass Ratio of W Ursae Majoris Binaries. APJL 1995, 444, L41. [CrossRef] 21. Ste ¸ pien, ´ K. Evolution of the progenitor binary of V1309 Scorpii before merger. Astron. Astrophys. 2011, 531, A18. [CrossRef] 22. D’Souza , M.C.R.; Motl, P.M.; Tohline, J.E.; Frank, J. Numerical Simulations of the Onset and Stability of Dynamical Mass Transfer in Binaries. Astrophys. J. 2006, 643, 381. [CrossRef] 23. Koenigsberger, G.; Moreno, E. Stability of Macroscopic Binary Systems. Rev. Mex. Astron. Astrofis. 2016, 52, 113. 24. Karachentsev, I.D. Binary Galaxies; Nauka: Moscow, Russia, 1987. (In Russian) 25. Muradian, R.M. The Regge law for heavenly bodies. Phys. Part. Nucl. 1997, 28, 471. [CrossRef] 26. Ambartsumian, V.A. On the origin of double stars. Vistas Astron. 1937, 14, 207. [CrossRef] 27. Ambartsumian, V.A.; Mirzoyan, L.V. An observational approach to the early stages of stellar evolution. Natl. Acad. Sci. Armen. 1947, 84, 317–330 Astronomy 2023, 2 89 28. Ambartsumian, V.A. Vistas in Astronomy; Beer, A., Ed.; Elsevier: Amsterdam, The Netherland, 1956; Volume 2. 29. Sargsyan, V.V.; Lenske, H.; Adamian, G.G.; Antonenko, N.V. From dinuclear systems to close binary stars: Application to mass transfer. Int. J. Mod. Phys. E 2018, 27, 1850063. [CrossRef] 30. Sargsyan, V.V.; Lenske, H.; Adamian, G.G.; Antonenko, N.V. Close Binary Galaxies: Application to Source of Energy and Expansion in Universe. Int. J. Mod. Phys. E 2019, 28, 1950031. [CrossRef] 31. Sargsyan, V.V.; Lenske, H.; Adamian, G.G.; Antonenko, N.V. Origin of the orbital period change in contact binary stars. Int. J. Mod. Phys. E 2019, 28, 1950044. [CrossRef] 32. Sargsyan, V.V.; Lenske, H.; Adamian, G.G.; Antonenko, N.V. From dinuclear systems to close binary stars and galaxies. Phys. Atom. Nucl. 2020, 83, 60. 33. Adamian, G.G.; Antonenko, N.V.; Lenske, H.; Sargsyan, V.V. Stability of Macroscopic Binary Systems. Commun. Theor. Phys. 2019, 71, 1335. [CrossRef] 34. Adamian, G.G.; Antonenko, N.V.; Lenske, H.; Sargsyan, V.V. On the evolution of compact binary black holes. Int. J. Mod. Phys. E 2020, 29, 2050094. [CrossRef] 35. Adamian, G.G.; Antonenko, N.V.; Lenske, H.; Sargsyan, V.V. On the possibility of formation of binary cosmic systems from the single cosmic objects. Int. J. Mod. Phys. E 2022, 31, 2250071. [CrossRef] 36. Adamian, G.G.; Antonenko, N.V.; Scheid, W. Clustering effects within the dinuclear model. in book. In Clusters in Nuclei; Lecture Notes in Physics; Beck, C., Ed.; Springer: Berlin, Germany, 2012; Volume 2. 37. Adamian, G.G.; Antonenko, N.V.; Zubov, A.S. Dinuclear systems in complete fusion reactions. Phys. Part. Nucl. 2014, 45, 848. [CrossRef] 38. Devries, R.M.; Clover, M.R. Coulomb potentials in heavy-ion interactions. Nucl. Phys. A 1975, 243, 528. [CrossRef] 39. Qian, S.-B. A possible relation between the period change and the mass ratio for W-type contact binaries. Mon. Not. R. Astron. Soc. 2001, 328, 635–644. [CrossRef] 40. Qian, S.-B. Are overcontact binaries undergoing thermal relaxation oscillation with variable angular momentum loss? Mon. Not. R. Astron. Soc. 2003, 342, 1260. 41. Yang, Y.-G.; Qian, S.-B. Deep, low mass ratio overcontact binary systems. XIV. A statistical analysis of 46 sample binaries. Astron. J. 2015, 150, 69. 42. Yakut, K.; Eggleton, P.P. Evolution of close binary systems. Astrophys. J. 2005, 629, 1055. [CrossRef] 43. Gazeas, K.; Stepien, K. Angular momentum and mass evolution of contact binaries. Mon. Not. R. Astron. Soc. 2008, 390, 1577. 44. Socia, Q.J.; Welsh, W.F.; Short, D.R.; Orosz, J.A.; Angione, R.J.; Windmiller, G.; Caldwell, D.A.; Batalha, N.M. KIC 9832227: Using Vulcan Data to negate the 2022 Red Nova merger prediction. Astrophys. J. Lett. 2018, 864, L32. [CrossRef] 45. Regge, T. Introduction to complex orbital momenta. Nuovo Cim. 1959, 14, 951. [CrossRef] 46. Collins, P.D.B. An Introduction to Regge Theory and High Energy Physics; Cambringe University Press: Cambringe, UK, 1977. 47. Barone, V.; Predazzi, E. High-Energy Particle Diffraction; Springer: Berlin/Heidelberg, Germany, 2002. 48. Chew, G.F.P.; Frautschi, S.C. Principle of Equivalence for all Strongly Interacting Particles within the S-Matrix Framework. Phys. Rev. Lett. 1961, 7, 394. [CrossRef] 49. Gribov, V.N. Partial waves with complex orbital angular momenta and the asymptotic behavior of the scattering amplitude. Sov. Phys. JETP 1961, 14, 1395. 50. Muradian, R.M. The Primeval Hadron—Origin of Stars Galaxies and Astronomical Universe. Astrophys. Space Sci. 1980, 69, 339. [CrossRef] 51. Muradian, R.M. On the origin of the rotation of galaxies in Ambartsumian’s cosmogony. Astrofiz 1975, 11, 237. (In Russian) 52. Muradian, R.M. Cosmic numbers and the rotation of the metagalaxy. Astrofiz 1977, 13, 63. (In Russian) 53. Muradian, R.M. Origin of magnetic fields and superdense cosmogony. Astrofiz 1978, 14, 439. (In Russian) 54. Hawking, S.W. Black hole explosions? Nature 1974, 248, 30. [CrossRef] 55. Bowen, D.B.; Campanelli, M.; Krolik, J.H.; Mewes, V.; Noble, S.C. Electromagnetic Emission from Supermassive Binary Black Holes Approaching Merger. Astrophys. J. 2017, 838, 42. 56. Michell, J., VII. On the means of discovering the distance, magnitude, &c. of the fixed stars, in consequence of the diminution of the velocity of their light, in case such a diminution should be found to take place in any of them, and such other data should be procured from observations, as would be farther necessary for that purpose. Phylosophical Trans. R. Soc. 1784, 74, 35. Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png
Astronomy
Multidisciplinary Digital Publishing Institute
http://www.deepdyve.com/lp/multidisciplinary-digital-publishing-institute/from-a-dinuclear-system-to-close-binary-cosmic-objects-vczy1xLBiz