Effects of stellar wind,dynamical friction,and star mergers on the dynamical evolution of multiple stars

The dynamical evolution of small stellar groups composed of N=6 components was numerically simulated within the framework of a gravitational N-body problem. The effects of stellar mass loss in the form of stellar wind, dynamical friction against the interstellar medium, and star mergers on the dynamical evolution of the groups were investigated. A comparison with a purely gravitational N-body problem was made. The state distributions at the time of 300 initial system crossing times were analyzed. The parameters of the forming binary and stable triple systems as well as the escaping single and binary stars were studied. The star-merger and dynamical-friction effects are more pronounced in close systems, while the stellar wind effects are more pronounced in wide systems. Star-mergers and stellar wind slow down the dynamical evolution. These factors cause the mean and median semimajor axes of the final binaries as well as the semimajor axes of the internal and external binaries in stable triple systems to increase. Star mergers and dynamical friction in close systems decrease the fraction of binary systems with highly eccentric orbits and the mean component mass ratios for the final binaries and the internal and external binaries in stable triple systems. Star mergers and dynamical friction in close systems increase the fraction of stable triple systems with prograde motions. Dynamical friction in close systems can both increase and decrease the mean velocities of the escaping single stars, depending on the density of the interstellar medium and the mean velocity of the stars in the system.     

New orbits of wide visual double stars

Based on photographic and CCD observations with the Pulkovo 26-inch refractor, radial velocity measurements with the 1.5-m RTT-150 telescope (TUBITAK National Observatory, Turkey), and highly accurate observations published in the WDS catalog, we have obtained the orbits of ten wide visual double stars by the apparent motion parameter method. The orientation of the orbits in the Galactic coordinate system has been determined. For the outer pair of the multiple star HIP 12780 we have calculated a family of orbits with a minimum period P = 4634 yr. Two equivalent solutions with the same period have been obtained for the stars HIP 50 (P = 949 yr) and HIP 66195 (P = 3237 yr). We have unambiguously determined the orbits of six stars: HIP 12777 (P = 3327 yr), HIP 15058 (P = 420 yr), HIP 33287 (P = 1090 yr), HIP 48429 (P = 1066 yr), HIP 69751 (P = 957 yr), and HIP 73846 (P = 1348 yr). The orbit of HIP 55068 is orientated perpendicularly to the plane of the sky, P >1000 yr. The star HIP 48429 is suspected to have an invisible companion.     

Gravitational loss-cone instability: Monotonic and nonmonotonic angular momentum distributions of stars

Small perturbations of spherical star clusters around massive black holes are studied. The presence of a black hole gives rise to peculiar distributions that have no stars with low angular momenta (falling into the so-called “loss cone”). The stability of such a distribution has been found to depend significantly on whether it monotonically increases with angular momentum L (from the loss cone up to L = L _circ in circular orbits) or has a maximum at some intermediate L = L _*. In the case of spherical systems under consideration, the loss-cone instability is shown to be possible only for nonmonotonic distributions.     

A unified theory of the formation of galactic bars

We give arguments for a basically unified formation mechanism of slow (Lynden-Bell) and fast (common) galactic bars. This mechanism is based on an instability that is akin to the well-known instability of radial orbits and is produced by the mutual attraction and alignment of precessing stellar orbits (so far, only the formation of slow bars has been explained in this way). We present a general theory of the low-frequency modes in a disk that consists of orbits precessing at different angular velocities. The problem of determining these modes is reduced to integral equations of moderately complex structure. The characteristic pattern angular velocities Ω_p of the low-frequency modes are of the order of the mean orbital precession angular velocity \(\bar \Omega _{pr}\). Bar modes are also among the low-frequency modes; while \(\Omega _p \approx \bar \Omega _{pr}\) for slow bars, Ω_p for fast bars can appreciably exceed even the maximum orbital precession angular velocity in the disk Ω _pr ^max (however, it remains of the order of these precession angular velocities). The possibility of such an excess of Ω_p over Ω _pr ^max is associated with the effect of “repelling” orbits. The latter tend to move in a direction opposite to the direction in which they are pushed. We analyze the pattern of orbital precession in potentials typical of galactic disks. We note that the maximum radius of an “attracting” circular orbit r_c can serve as a reasonable estimate of the bar length l_b. Such an estimate is in good agreement with the available results of N-body simulations.     

Simulations of slow bars in anisotropic disk systems

The instability of anisotropic disk systems with elongated stellar orbits has been investigated. N-body generalized polytropic models of stellar disks have been constructed. They are shown to be unstable with respect to the bar formation at any degree of anisotropy. This result differs from the results of the studies of such models by other authors. The bar pattern speed and amplitude have been found. The initial distribution of precession rates and the adiabatic invariants of stellar orbits have been calculated. A bar is shown to be formed in such systems due to the radial orbit instability.     

Magnetic field and chemical composition of the peculiar star HD 10221

We analyzed the chemical composition of the chemically peculiar (CP) star HD 0221=43 Cas using spectra taken with the NES spectrograph of the 6-m telescope with a spectral resolution of 45 000. The Hβ line profile corresponds most closely to T_eff = 11 900 K and log g = 3.9. The rotational velocity is v_e sin i = 27 ± 2 km s^?1, and the microturbulence is ξ_t = 1 km s^?1. The results of our abundance determination by the method of synthetic spectra show that the star has chemical anomalies typical of SrCrEu stars, although its effective magnetic field is weak, B_e < 100 G. For silicon, we obtained an abundance distribution in atmospheric depth with a sharp jump of 1.5 dex at an optical depth of log τ₅₀₀₀ = ?0.3 and with silicon concentration in deep atmospheric layers. Similar distributions were found in the atmospheres of cooler stars with strong and weak magnetic fields. A comparison of the chemical peculiarities in HD 10221 with known CP stars with magnetic fields of various strengths leads us to conclude that a low rotational velocity rather than amagnetic field is the determining factor in the formation mechanism of chemical anomalies in the atmospheres of CP stars.     

Astrometric Study of Two Multiple Stars: ADS 2668 and ADS 8236

Based on CCD observations with the Pulkovo 26-inch refractor in 2003–2018, we have obtained the orbit of the visual double star ADS 2668 AB (P = 947 yr, a = 2.9″, e = 0.41, ω = 246°, Ω = 131°, i = 114°, T = 1456 yr) for the first time by the apparent motion parameter (AMP) method, which is consistent with the inner orbit of ADS 2668 Aa-Ab, and improved the orbit of ADS 8236 AB (P = 1996 yr, a = 4.69″, e = 0.39, ω = 201°, Ω = 166°, i = 110°, T = 1246 yr). The inner orbit of the photocenter of ADS 8236 with a period of 4.627 yr has been calculated from the residuals. This orbit of ADS 8236 Ba–Bb supplements the spectroscopic orbit by the elements specifying the orbital plane (i and Ω). In both cases, the planes of the inner and outer orbits are noncoplanar. The presence of an additional companion in the system ADS 2668 is discussed.     

Occultations of stars brighter than 15<Superscript>m</Superscript> by the largest trans-Neptunian objects in 2004–2014

We computed the occultations of stars brighter than 15^m by the largest trans-Neptunian objects (TNOs) for the next ten years. In our search, we used the following catalogs: Hipparcos; Tycho2 with the coordinates of 2838666 stars taken from UCAC2 (Herald 2003); and UCAC2 (Zacharias et al. 2003) with 16356096 stars between \(12\mathop .\limits^m 00\) and \(14\mathop .\limits^m 99\) north of ?45° declination. We predicted the occultations of stars by the seventeen largest numbered TNOs, the recently discovered 2004 DW, and four known binary Kuiper Belt objects. We selected 64 events at solar elongations of no less than 30°, including the extremely rare occultation of a \(6\mathop .\limits^m 5\) star by the double asteroid (66652) 1999 RZ₂₅₃ on October 4, 2007. Observations of these events by all available means are extremely important, since they can provide unique information about the sizes of TNOs and improve our knowledge of their orbits dramatically.     

3D stellar reddening map from 2MASS photometry: An improved version

An improved version of the 3D stellar reddening map in a space with a radius of 1200 pc around the Sun and within 600 pc of the Galactic midplane is presented. As in the previous 2010 and 2012 versions of the map, photometry with an accuracy better than 0.05 ^m in the J and Ks bands for more than 70 million stars from the 2MASS catalogue is used in the new version. However, the data reduction technique is considerably more complicated. As before, an analysis of the distribution of stars near the main-sequence turnoff on the (J ? Ks)?Ks diagram, where they form a distribution maximum, provides a basis for the method. The shift of this maximum, i.e., the mode (J ? Ks), along (J ? Ks) and Ks, given the spatial variations of the mean dereddened color (J ? Ks)₀ of these stars, is interpreted as a growth of the reddening with increasing distance. The main distinction of the new method is that instead of the fixed mean absolute magnitude, dereddened color, distance, and reddening for each cell, the individual values of these quantities are calculated for each star by iterations when solving the system of equations relating them. This has allowed one to increase the random accuracy of the map to 0.01 ^m and its spatial resolution to 20 pc in coordinates and distance and to 1° in longitude and latitude. Comparison with other reddening estimates for the same spatial cells and Gaia DR1 TGAS stars shows that the constructed map is one of the best maps for the space under consideration. Its systematic errors have been estimated to be σ(E(J ? Ks)) = 0.025 ^m , or σ(E(B ? V)) = 0.04 ^m . The main purpose of the map is to analyze the characteristics of Galactic structures, clouds, and cloud complexes. For this purpose, the reddening map within each spatial cell has also been computed by analyzing the reddening along each line of sight.     

Features of magnetic field structures in chemically peculiar stars. IV

We present results of modeling of the sample of magnetic stars. We have obtained such important for magnetic star physics parameters as the mean surface magnetic field B_s, the magnetic field at magnetic poles—B_p, the dipole inclination to the rotation equatorial plane α, and the distance to monopoles from the center of the star Δa. We present some information onmagnetic star physics that helps to understand the derived results better.     

Orbital stability of systems of closely-spaced planets,II: configurations with coorbital planets

We numerically investigate the stability of systems of 1 \({{\rm M}_{\oplus}}\) planets orbiting a solar-mass star. The systems studied have either 2 or 42 planets per occupied semimajor axis, for a total of 6, 10, 126, or 210 planets, and the planets were started on coplanar, circular orbits with the semimajor axes of the innermost planets at 1 AU. For systems with two planets per occupied orbit, the longitudinal initial locations of planets on a given orbit were separated by either 60° (Trojan planets) or 180°. With 42 planets per semimajor axis, initial longitudes were uniformly spaced. The ratio of the semimajor axes of consecutive coorbital groups in each system was approximately uniform. The instability time for a system was taken to be the first time at which the orbits of two planets with different initial orbital distances crossed. Simulations spanned virtual times of up to 1 × 10⁸, 5 × 10⁵, and 2 × 10⁵ years for the 6- and 10-planet, 126-planet, and 210-planet systems, respectively. Our results show that, for a given class of system (e.g., five pairs of Trojan planets orbiting in the same direction), the relationship between orbit crossing times and planetary spacing is well fit by the functional form log(t _c /t ₀) = b β + c, where t _c is the crossing time, t ₀ = 1 year, β is the separation in initial orbital semimajor axis (in terms of the mutual Hill radii of the planets), and b and c are fitting constants. The same functional form was observed in the previous studies of single planets on nested orbits (Smith and Lissauer 2009). Pairs of Trojan planets are more stable than pairs initially separated by 180°. Systems with retrograde planets (i.e., some planets orbiting in the opposite sense from others) can be packed substantially more closely than can systems with all planets orbiting in the same sense. To have the same characteristic lifetime, systems with 2 or 42 planets per orbit typically need to have about 1.5 or 2 times the orbital separation as orbits occupied by single planets, respectively.     

Steady state obliquity of a rigid body in the spin–orbit resonant problem: application to Mercury

In this paper, we consider the elliptic collinear solutions of the classical n-body problem, where the n bodies always stay on a straight line, and each of them moves on its own elliptic orbit with the same eccentricity. Such a motion is called an elliptic Euler–Moulton collinear solution. Here we prove that the corresponding linearized Hamiltonian system at such an elliptic Euler–Moulton collinear solution of n-bodies splits into \((n-1)\) independent linear Hamiltonian systems, the first one is the linearized Hamiltonian system of the Kepler 2-body problem at Kepler elliptic orbit, and each of the other \((n-2)\) systems is the essential part of the linearized Hamiltonian system at an elliptic Euler collinear solution of a 3-body problem whose mass parameter is modified. Then the linear stability of such a solution in the n-body problem is reduced to those of the corresponding elliptic Euler collinear solutions of the 3-body problems, which for example then can be further understood using numerical results of Martínez et al. on 3-body Euler solutions in 2004–2006. As an example, we carry out the detailed derivation of the linear stability for an elliptic Euler–Moulton solution of the 4-body problem with two small masses in the middle.     

Einstein–Maxwell Field Equation in Isotropic Coordinates: An Application to Modeling Superdense Star

We present a charged analogue of Pant et al. (2010, Astrophys. Space Sci., 330, 353) solution of the general relativistic field equations in isotropic coordinates by using simple form of electric intensity E that involve charge parameter K. Our solution is well behaved in all respects for all values of X lying in the range 0 <X≤ 0.11, K lying in the range 4 <K≤ 6.2 and Schwarzschild compactness parameter u lying in the range 0 <u≤ 0.247. Since our solution is well behaved for wide ranges of the parameters, we can model many different types of ultra-cold compact stars like quark stars and neutron stars. We have shown that corresponding to X = 0.077 and K = 6.13 for which u = 0.2051 and by assuming surface density ρ _b =4.6888×10¹⁴ g cm ^?3 the mass and radius are found to be 1.509M _⊙, 10.906 km respectively which match with the observed values of mass 1.51M _⊙ and radius 10.90 km of the quark star XTE J1739-217. The well behaved class of relativistic stellar models obtained in this work might have astrophysical significance in the study of more realistic internal structures of compact stars.     

Kinematic control of the inertiality of ICRS catalogs

We perform a kinematic analysis of the Hipparcos and TRC proper motions of stars by using a linear Ogorodnikov-Milne model. All of the distant (r>0.2 kpc) stars of the Hipparcos catalog have been found to rotate around the Galactic y axis with an angular velocity of M ₁₃ ^? =?0.36±0.09 mas yr^?1. One of the causes of this rotation may be an uncertainty in the lunisolar precession constant adopted when constructing the ICRS. In this case, the correction to the IAU (1976) lunisolar precession constant in longitude is shown to be Δp₁=?3.26±0.10 mas yr^?1. Based on the TRC catalog, we have determined the mean Oort constants: A=14.9±1.0 and B=?10.8±0.3 km s^?1 kpc^?1. The component of the model that describes the rotation of all TRC stars around the Galactic y axis is nonzero for all magnitudes, M ₁₃ ^? =?0.86±0.11 mas yr^?1.     

Stars with Discrepant <Emphasis Type="Italic">v</Emphasis> sin <Emphasis Type="Italic">i</Emphasis> as Derived from the Ca II 3933 and Mg II 4481 Å Lines. VII. HD9531 (SB), HD31592 (SB2), HD129174 (SB?)

In this paper of the series we analyze three stars listed among stars with discrepant v sin i: HD9531 and HD31592, which also show radial velocity variations inherent to spectroscopic binaries, and HD129174 which is an Mn-type star with a possible magnetic field. In HD9531 we confirm the radial velocity derived fromthe hydrogen lines as well as fromthe Ca II line at 3933 Å as variable. The profile of the calcium line also appears variable, and with the estimated magnetic induction B_e = ?630 ± 1340 G, this suggests that the abundance of calcium possibly varies over the surface of the star. We identified the lines of the secondary component in the spectrum of HD31592 revealing thus it is an SB2 binary with B9.5V and A0V components. While the primary star rotates with v sin i = 50 km s^?1, the secondary star is faster with v sin i = 170 km s^?1. We find that only 60% of the Mn lines identified in the spectrum of HD129174 can be fitted with a unique abundance value, whereas the remaining lines are stronger or fainter. We also identified two Xe II lines at 5339.33 Å and 5419.15 Å and estimated their log g f.     

Close binary systems in multiple stars: II. Eclipsing variables in visual binary and multiple systems

Eclipsing variables in visual binary and multiple stars are searched using data from GCVS, WDS, and CCDM catalogs. The list of 421 eclipsing variables is obtained. The masses of components of multiple systems from the list are estimated using the mass-luminosity relation for the main sequence stars. It is shown that, for 85% multiple systems from the list, the mass of visual components is smaller by a factor of 2 than the total mass of close binary systems. The distributions of orbital elements of visual binary systems are constructed and used for calculation of orbit semi-major axes for star from the list. The distributions of orbit semi-major axes and periods obtained from observations are approximated by Gaussian curves. The maxima of the curves correspond to a = 800 a.u. and P = 7600 years, respectively. The distribution of orbit semi-major axes larger than 800 a.u. is better described by Opik’s law; it is expected that this law describes the real a distribution in the region of small values as well. The frequency of eclipsing variables in multiple stars makes 12% of the total number of stars of this type in GCVS.     

Radial pulsations of helium stars with masses from 10 to 50<Emphasis Type="Italic">M</Emphasis><Subscript>⊙</Subscript>

We have performed hydrodynamic calculations of the radial pulsations of helium stars with masses 10M_⊙ ≤ M ≤ 50M_⊙, luminosity-to-mass ratios 5 × 10³L_⊙/M_⊙ ≤ L/M ≤ 2.5 × 10⁴L_⊙/M_⊙, and effective temperatures 2 × 10⁴ K ≤ T_eff ≤ 10⁵ K for helium and heavy-element mass fractions of Y=0.98 and Z=0.02, respectively. We show that the high-temperature boundary of the instability region for radial pulsations at L/M ? 10⁴L_⊙/M_⊙ extends to T_eff≈10⁵ K. The amplitude of the velocity variations for outer layers is several hundred km s^?1, while the brightness variations in the B band of the UBV photometric system are within the range from several hundredths to half a magnitude. At constant luminosity-to-mass ratio, the radial pulsation period is determined only by the effective temperature of the star. In the ranges of luminosity-to-mass ratios 10⁴L_⊙/M_⊙ ≤ L/M ≤ 2 × 10⁴L_⊙/M_⊙ and effective temperatures 5 × 10⁴ K ≤ T_eff ≤ 9 × 10⁴ K, the periods of the radial modes are within 6 min ?Π?10³ min.     

OB stars in the Tycho-2 and 2MASS catalogues

The Tycho-2 proper motions and five-band Tycho-2 and 2MASS photometry for approximately 2.5 million common stars have been used to select OB stars and to determine the extinction and photometric distance for each of them. We have selected 37 485 stars and calculated their reddenings based on their positions in the two-color V _T -H, J-Ks diagrams relative to the zero-age main sequence and the theoretical reddening line for B5 stars. Tests confirm that the selected stars belong to the spectral types O-B with a small admixture of later types. We calculate the extinction coefficient R and its variations with Galactic longitude based on the positions of the selected stars in the two-color B _T -V _T , V _T -Ks diagram. The interstellar extinction for each star is calculated as the product of the reddening found and the coefficient R. The extinction and its variations with Galactic longitude agree well with the extinction based on the model by Arenou et al. (1992). Calibration of the relation between the absolute magnitude and reduced proper motion V _T − + 5 + 5 log μ for Hipparcos stars has allowed us to calculate the absolute magnitudes and photometric distances for the selected stars. The distances found agree with those derived from the Hipparcos parallaxes within 500 pc. The distribution of the stars and the extinction variations with distance found show that the selected stars form an almost complete sample of stars with spectral types earlier than B5 within about 750 pc of the Sun. The sample includes many noticeably reddened stars in the first and second Galactic quadrants that are absent from the Hipparcos and Tycho Spectral Types Catalogues. This slightly changes the pattern of the distribution of OB stars compared to the classical pattern based on Hipparcos. Original Russian Text ? G.A. Goncharov, 2008, published in Pis’ma v Astronomicheskiĭ Zhurnal, 2008, Vol. 34, No. 1, pp. 10–20.     

Results of magnetic field measurements of CP-stars performed with the 6-m telescope. III. Observations in 2009

We present the results of measuring longitudinal magnetic fields (B_e), rotation velocities (v_e sin i), and radial velocities (V_r) of 44 stars observed with the Main Stellar Spectrograph (MSS) of the 6-m BTA telescope of the Special Astrophysical Observatory in 2009. For the first time, magnetic fields were detected for the stars HD5441, HD199180, HD225627, and BD+00° 4535. We show that for the same stars, the longitudinal fields B_e measured from the Hβ hydrogen line core and from metal lines can differ by 10% and up to a factor of 2–3. Except in rare cases, magnetic fields measured from the metal lines are stronger. We believe that this phenomenon is of a physical nature and depends on the magnetic field topology and the physical conditions inside a specific star. Observations of standard stars without magnetic fields confirm the absence of systematic errors capable of introducing distortions into the longitudinal-field measurement results. In this work we comment on the results for each of the stars.     

Planetary system formation through binary star merging

A scenario is considered for the formation of a planetary system through the merging of a binary star comprised of low-mass (0.5–1 M _⊙) stars in the stage of contracting towards the main sequence. According to our previous computations (Sirotkin and Karetnikov, 2006), under certain conditions, the destruction of the more massive component can result in the formation of a central star, an accretion disk, and an extended arm. The extended arm is fragmented to form clouds of planetary masses (<5M _J). The formed disk and clouds rotate in the same direction as the central star. The clouds are in elongated orbits (e > 0.3) lying in the orbital plane of the initial binary system. To test these earlier results, we repeated computations for the same system parameters but with higher accuracy. The new computations confirmed the earlier results and gave new information about the cloud and disk structure.