Received:2019-07-30Online:2019-11-1
Fund supported: |
Abstract
Keywords:
PDF (598KB)MetadataMetricsRelated articlesExportEndNote|Ris|BibtexFavorite
Cite this article
G. G. Adamian, N. V. Antonenko, H. Lenske, V. V. Sargsyan. Stability of Macroscopic Binary Systems*. [J], 2019, 71(11): 1335-1340 doi:10.1088/0253-6102/71/11/1335
It is now commonly believed that the contact binaries, for example, the W UMa binaries, end their evolution by merging into a single star.[1] The dissipation of the orbital energy in the initial violent phase of merger resulted in the luminous red nova V1309 Sco observed in 2008. The luminous red novae have recently been identified as a distinct class of stellar transients.[2] They are characterized by relatively long outbursts with spectral distributions centered in the red, ranging in luminosities between classical novae and supernovae. The observations in the Optical Gravitational Lensing Experiment have revealed a shape of the light curve characteristic for contact binary system with an exponentially decreasing orbital period[3] and confirmed the earlier conjecture of Ref. [4] that the luminous red novae arise from merging contact binary stars. The spectacular case is KIC 9832227 which was predicted[5] to be merge in 2022, enlightening the sky as a red nova. The compact binaries composed of white dwarfs, neutron stars, and black holes eventually merge through gravitational wave emission.[6]
The details of specific mechanism, which triggers the merger of the contact binary components, are still controversial.[7-12] There is assumption[11] on the gradual mass transfer from less massive (and hence smaller radius) secondary star to the primary one (heavier star). This mass transfer is driven by the structural change in the secondary star caused by the energy received from the primary star. The contact is sustained by the magnetic braking or thermonuclear evolution. When the mass ratio is extreme enough for the Darwin instability, a merger starts that triggers the outburst in red novae.[3, 13] The Darwin instability happens when the spin angular momentum of the system is more than one third of the orbital angular momentum. This instability plays a role once the mass ratio became small enough that the companion star can no longer keep the primary star synchronously rotating via the tidal interaction. For most of the primary massive stars, this occurs at the mass ratio $q=M_2/M_1< 0.1$.[14] There is alternative scenario:[15] At first contact in a binary system, a brief, but intense, mass transfer sets in changing originally more massive star into less massive one. This process may oscillate until eventually a stable contact configuration is reached. Also the dynamic mass transfer without the Darwin instability[16] and mergers triggered by a tidal runaway based on a non-equilibrium response to tidal dissipation[17] have been investigated. Thus, there is a large quest for detailed observational and theoretical investigations.[15-17]
The Regge-theory proved to be very influential in the development of elementary-particle physics. The application of Regge ideas to astrophysics has shown that the spins of planets and stars are well described by the Regge-like law for a sphere ($S\sim M^{4/3}$), while the spins of galaxies and clusters of galaxies obey the Regge-like law for a disk ($S\sim M^{3/2}$).[18-20] In contrast to earlier semi-phenomenological approaches these expressions contain only fundamental constants as the parameters and are independent of any fitted empirical quantities.[18-20] In Refs. [18--20], the cosmic analog of the Chew-Frautschi plot with two important cosmological Eddington and Chandrasekhar points on it has been also constructed.
The aim of the present article is to study the Darwin instability effect in the binary star or galaxy by using the well-known model of Refs. [18--20] based on the Regge-theory. The total angular momentum ${\pmb J}_{\rm tot}$ of binary system is the sum of orbital angular momentum ${\pmb L}$ and the spins ${\pmb S}_k$ ($k=1, 2$) of the individual components:
The $J_{\rm tot}$ and $S_k$ are expressed using the Regge-like law for stars and planets ($n=3$) or galaxies ($n=2$):
and
where $\hbar$, $m_p$, $M_k$ ($k=1,2$), and $M=M_1+M_2$ are the Planck constant, masses of proton and astrophysical objects (planets, stars or galaxies), and the total mass of system, respectively. The maximum (the antiparallel orbital and spin angular momenta) and minimum (the parallel orbital and spin angular momenta) orbital angular momenta are
and
respectively. Using the mass asymmetry (mass transfer) coordinate $\eta=(M_1-M_2)/M$ instead of masses $M_1=\frac{M}{2}(1+\eta)$ and $M_2=\frac{M}{2}(1-\eta)$[21-23] and Eqs. (3)--(5), we derive
At $\eta=0$ (the symmetric binary system), we have
$$\frac{S_1+S_2}{L_{\rm min}}=\frac{1}{2^{1/n}-1}>1$$
and
$$\frac{S_1+S_2}{L_{\rm max}}=\frac{1}{2^{1/n}+1}>\frac{1}{3}\,.$$
For the symmetric binary star (planet) and binary galaxy with $\eta=0$, $(S_1+S_2)/L_{\rm max}\approx$ 0.44 and 0.41, respectively. At $\eta=1$, we have
$$\frac{S_1+S_2}{L_{\rm min}}\to\infty$$
and
$$\frac{S_1+S_2}{L_{\rm max}}=\frac{1}{2}.$$
As follows from last two expressions, for very asymmetric binaries, the ratios $(S_1+S_2)/L_{\rm max,min}$ are almost independent of the value of $n$. According to Ref. [14], the Darwin instability can occur when the binary mass ratio is very small [$q=M_2/M_1< 0.1$] or the mass asymmetry is very large [$\eta=(1-q)/(1+q)>0.82$]. As seen in Fig. 1, the ratios $(S_1+S_2)/L_{\rm max}$ and $(S_1+S_2)/L_{\rm min}$ continuously increases with $\eta$ increasing from 0 to 1. Because their absolute values are larger than 1/3, all possible binary stars (planets) or binary galaxies, independently of their mass asymmetry $\eta$, should have the Darwin instability ($S_1 + S_2 \ge \frac{1}{3}L$)[13] and, correspondingly, should merge. However, the observations do not support this conclusion which probably means that there is no Darwin instability effect in such binary systems and, correspondingly, the mechanism of merger has other origin. Thus, it is necessary to search another mechanism that triggers the merging of the contact binary components.
Fig. 1
New window|Download| PPT slideFig. 1The calculated ratios $(S_1+S_2)/L_{\rm max}$ (a), $(S_1+S_2)/L_{\rm min}$ (b), $(S_1-S_2)/L_{2}$ [solid line] (c), and $(S_2-S_1)/L_{1}$ [dashed line] (c) as functions of mass asymmetry at $n=3$. In the cases of antiparallel spins (c), $L_1=J_{\rm tot}+S_1-S_2$ and $L_2=J_{\rm tot}-S_1+S_2$.
Note that in the cases of antiparallel spins with $L_1=J_{\rm tot}+S_1-S_2$ and $L_2=J_{\rm tot}-S_1+S_2$ (Fig. 1), the ratios $|S_2-S_1|/L_{1}$ and $|S_1-S_2|/L_{2}$ are larger than $\frac{1}{3}$ for the asymmetric binaries with $|\eta|\ge 0.5$.
As seen in Fig. 2, the dependencies of $L_{\rm max}$, $L_{\rm min}$, $L_{1}$, and $L_{2}$ on mass asymmetry have different behavior. The evolution of system in mass asymmetry (mass transfer) can increase or decrease the orbital angular momentum. For example, at $\eta\to 0$ the binary system has smaller $L=L_{\rm max}$. The observations of the dependence of $L$ on $\eta$ may be useful to distinguish the difference between the orientations of orbital and spin angular momenta.
Fig. 2
New window|Download| PPT slideFig. 2The calculated ratios $ L_{\rm max}/J_{\rm tot}$ (a), $ L_{\rm min}/J_{\rm tot}$ (b), $L_{2}/J_{\rm tot}$ [solid line] (c), and $L_{1}/J_{\rm tot}$ [dashed line] (c) as functions of mass asymmetry at $n=3$.
Employing Eqs. (1)--(3) and results of Refs. [21, 23], we obtain new analytical formulas for the relative distance between the components of the binary
at $R_{m}>R_{t}=R_1+R_2$ ($R_{k}$, ${k}=1,2$, are the radii of binary components and $R_{t}$ is the touching distance) and
at $R_{ m}\le R_{t}$.One can similarly derive the formulas for the orbital rotation period
at $R_{m}>R_{t}$ and
at $R_{m}<R_{t}$.[21, 23] Here, $G$ is the gravitational constant. The values $\epsilon_1=\epsilon_2=1$ and $\epsilon_1=\epsilon_2=-1$ correspond to the cases of antiparallel and parallel orbital and spin angular momenta. The values $\epsilon_1=-\epsilon_2=1$ and $\epsilon_1=-\epsilon_2=-1$ correspond to the cases of antiparallel spins. The observational data result in the relationship $R_{k}=gM_{k}^{l}$ between the radius and mass of the star, where the constants $l=\frac{2}{3}$ and $g=R_{\odot}/M_\odot^{2/3}$ ($M_{\odot}$ and $R_{\odot}$ are mass and radius of the Sun)[24] and the galaxy, where the constant $l$ depending on mass is in the interval $\left[\frac{2}{5},\frac{2}{3}\right]$.[22, 25] As seen in Fig. 3, at $R_{m}>R_{t}$,
$$R_{m}^{d}=\frac{\left[2^{(1+n)/n}+\epsilon_1 (1+\eta)^{(1+n)/n}+\epsilon_2(1-\eta)^{(1+n)/n}\right]^2}{\left[1-\eta^2\right]^2}$$
Fig. 3
New window|Download| PPT slideFig. 3The calculated dimensionless relative distances $R_{m}^{d}$ [(a), (b)] and orbital rotation periods $ P_{\rm orb}^{d}$ [(c), (d)] as functions of mass asymmetry at $R_{m}>R_{t}=R_1+R_2$ and $n=3$. The plots [(a), (c)] and [(b), (d)] correspond to the systems with the parallel and antiparallel orbital and spin angular momenta, respectively.
Or $R_{m}$ decreases with decreasing $|\eta|$ and, finally, the
$$P_{\rm orb}^{d}=\frac{\left[2^{(1+n)/n}+\epsilon_1 (1+\eta)^{(1+n)/n}+\epsilon_2 (1-\eta)^{(1+n)/n}\right]^3}{\left[1-\eta^2\right]^3}$$
or $P_{\rm orb}$ decreases. At $R_{m}>R_{t}$ (Fig. 3), the dependence of $P_{\rm orb}$ as a function of mass asymmetry has similar behavior in the cases when the orbital and spin angular momenta are antiparallel and parallel. In the case of antiparallel spins with $L_1=J_{\rm tot}+S_1-S_2$ ($L_2=J_{\rm tot}-S_1+S_2$) and $M_1\ge M_2$, the values of $R_{m}$ and $P_{\rm orb}$ decrease (increase ) with $\eta$ going from 1 to 0 (Fig. 4). At $R_{m}\le R_{t}$ (Fig. 5), the value of
$$P_{\rm orb}^{d}=\left[(1+\eta)^l+(1-\eta)^l\right]^{3/2}$$
Fig. 4
New window|Download| PPT slideFig. 4The calculated dimensionless relative distance $R_{m}^{d}$ (a) and orbital rotation period $ P_{\rm orb}^{d}$ (b) as functions of mass asymmetry in the cases of antiparallel spins with $L_1=J_{\rm tot}+S_1-S_2$ (dashed line) and $L_2=J_{\rm tot}-S_1+S_2$ (solid line), $n=3$, and $R_{m}>R_{t}=R_1+R_2$.
Fig. 5
New window|Download| PPT slideFig. 5The calculated dimensionless $R_{m}^{d}$ [(a), (b), (c)] and $ P_{\rm orb}^{d}$ [(d)] as functions of mass asymmetry at $R_{m}\le R_{t}=R_1+R_2$, $n=3$, and $l=2/3$. The plots (a) and (b) correspond to the system with the parallel and antiparallel orbital and spin angular momenta, respectively. In the plot (c), the cases of antiparallel spins with $L_1=J_{\rm tot}+S_1-S_2$ (dashed line) and $L_2=J_{\rm tot}-S_1+S_2$ (solid line) are presented. Note that the value of $ P_{\rm orb}^{d}$ does not depend on the orientations of orbital momentum and spins.
or $P_{\rm orb}$ increases with decreasing $|\eta|$ and does not depend on orientations of orbital momentum and spins. In contrast, the distance
$$\begin{eqnarray*} R_{\rm m}^{d}= \left[(1+\eta)^l+(1-\eta)^l\right]^{3/4}\\ \times\frac{\left[2^{(1+n)/n}+\epsilon_1 (1+\eta)^{(1+n)/n}+\epsilon_2 (1-\eta)^{(1+n)/n}\right]^{1/2}}{\left[1-\eta^2\right]^{1/2}} \end{eqnarray*} $$
or $R_{m}$ depends on the orientations of orbital momentum and spins. As seen in Fig. 5, the dependencies $R_{m}^{d}$ at different orientations on mass asymmetry have various behavior.
In conclusion, within the model[18-20] based on the Regge-like laws, we have shown that all possible binary stars (planets) or binary galaxies, independently of their mass asymmetry $\eta$, satisfy the Darwin instability condition ($S_1 + S_2 \ge \frac{1}{3}L$) which contradicts to the observations. This conclusion is not sensitive to the parameters of model. Therefore, one should search for other mechanism that triggers the merger of the contact binary components.
Employing the Regge-like laws, we have derived the new analytical formulas for the relative distance and orbital rotation period of the binary system, which depend on the total mass $M$, mass asymmetry $\eta$, and the fundamental constants $G$, $\hbar$, $m_p$. We have predicted that decreasing and increasing periods as functions of mass asymmetry (mass transfer) are related, respectively, with the non-overlapping ($R_{m} > R_1+R_2$) and overlapping ($R_{m} \le R_1+R_2$) stage of the binary object evolution.
Reference By original order
By published year
By cited within times
By Impact factor
[Cited within: 1]
[Cited within: 1]
[Cited within: 2]
[Cited within: 1]
[Cited within: 1]
[Cited within: 1]
[Cited within: 1]
[Cited within: 1]
[Cited within: 1]
[Cited within: 2]
[Cited within: 2]
[Cited within: 2]
[Cited within: 1]
[Cited within: 2]
[Cited within: 5]
[Cited within: 5]
[Cited within: 3]
[Cited within: 1]
[Cited within: 3]
[Cited within: 1]
[Cited within: 1]