The formation and evolution of accretion disks in binary systems

In this paper, we consider a close binary system consisting of a compact star and a optical Main-Sequence star which fills its critical Roche lobe and transfers matter through the inner Lagrangian pointL ₁ toward the compact object.We use the Hill's problem as the dynamical model of the binaries. The following binaries are calculated by Duncanet al.'s mapping approximation: RY Per, RZ Sct, RS Vul, and Tau. Figures 3–6 show that the trajectories of accretion disk particles in the binaries. The relation between the dimensionless semi-major axis and times of conjunctionN are presented in Figures 7–10.This work was supported by National Natural Science Foundation of China.     

Measurement of the orbital period of the X-ray burster GS 1826-238 based on observations of its optical brightness variations

The variability of the optical and X-ray fluxes from the binary GS 1826-238 is investigated. An epoch-folding analysis of the optical data obtained with the RTT-150 telescope in 2003–2004 has revealed periodic brightness variations in the source with a period P _orb = 2.24940 ± 0.00015 h with a high statistical significance. When estimating the detection significance of the periodic signal, we have specially taken into account the presence of a powerful aperiodic component (“red noise”) in the source’s brightness variability. The source’s power density spectra in the frequency range ∼10⁻⁵–0.01 Hz have been obtained. We have detected a statistically significant break in the power density spectrum of GS 1826-238 at a frequency ν _br ≈ (8.48 ± 0.14) × 10⁻⁵ Hz in both optical and X-ray energy bands. We have estimated the orbital period of the binary GS 1826-238 using the correlation between the break frequency in the power density spectrum and the orbital period of binaries, P _orb ∝ 1/ν _br, found by Gilfanov and Arefiev (2005): P _orb = 3.7 ± 0.8 h and P _orb = 11.3 ± 5.9 h when using Sco X-1 and 1H 16267-273, respectively, as reference sources. It seems to us that the method for estimating the orbital periods of low-mass X-ray binaries using the correlation P _orb ∝ 1/ν _br may turn out to be very promising, especially for persistent low-luminosity X-ray binaries.     

The origin of intergalactic thermonuclear supernovae

The population synthesis method is used to study the possibility of explaining the appreciable fraction of the intergalactic type-Ia supernovae (SN Ia), 20 ₋₁₅ ⁺¹² %, observed in galaxy clusters (Gal-Yam et al. 2003) when close white dwarf binaries merge in the cores of globular clusters. In a typical globular cluster, the number of merging double white dwarfs does not exceed ∼10⁻¹³ per year per average cluster star in the entire evolution time of the cluster, which is a factor of ∼3 higher than that in a Milky-Way-type spiral galaxy. From 5 to 30% of the merging white dwarfs are dynamically expelled from the cluster with barycenter velocities up to 150 km s⁻¹. SN Ia explosions during the mergers of double white dwarfs in dense star clusters may account for ∼1% of the total rate of thermonuclear supernovae in the central parts of galaxy clusters if the baryon mass fraction in such star clusters is ∼0.3%.     

Light-Time Effect in the Eclipsing Binaries GO Cyg, GW Cep, AR Aur and V505 Sgr

Orbital period changes of the eclipsing binaries GO Cyg and GW Cep are explained by the light-time effect for the first time. New minima of the eclipsing binary AR Aur improve the predicted light-time orbit. The light-time orbit with the quadratic ephemeris of the binary matches the new observations of V505 Sgr better than the linear one. As the light-time effect fits in corresponding O − C diagrams of all four systems have been reaching extreme values, the observations of minima times in forthcoming years are highly desirable.     

Chirikov diffusion in the asteroidal three-body resonance (5, −2, −2)

The theory of diffusion in many-dimensional Hamiltonian system is applied to asteroidal dynamics. The general formulation developed by Chirikov is applied to the Nesvorny-Morbidelli analytic model of three-body (three-orbit) mean-motion resonances (Jupiter-Saturn-asteroid). In particular, we investigate the diffusion along and across the separatrices of the (5, −2, −2) resonance of the (490) Veritas asteroidal family and their relationship to diffusion in semi-major axis and eccentricity. The estimations of diffusion were obtained using the Melnikov integral, a Hadjidemetriou-type sympletic map and numerical integrations for times up to 10⁸ years.     

Nature and distribution of dark matter: 2. binaries,groups and clusters

We study the mass-radius relationship for aggregates of galaxies, viz. binaries, small groups and clusters. The data are subjected to a simple best-fit analysis similar to the one carried out earlier for individual field galaxies. The analysis shows that: (i) The data on binary galaxies are consistent with the assumption that binaries are just two galaxies, each with an individual isothermal (M ∫R) dark matter halo, moving under the mutual gravitational attraction, (ii) The data on the groups of galaxies are too scattered to obey a single power-law relation of the formM = kR ⁿ with any degree of reliability, (iii) The data on groups and clusters fit better with a law of the formM = AR ³ +BR. This form suggests the existence of two components in dark matter—one which is clustered around the galaxies (M ∫R) and another which is distributed smoothly (M ∫R ³ ). The smooth distributions becomes significant only at scales ≥ 1 Mpc and hence does not affect binaries significantly. We briefly discuss the theoretical implications of this analysis     

Constraining Coronal Heating: Employing Bayesian Analysis Techniques to Improve the Determination of Solar Atmospheric Plasma Parameters

One way of revealing the nature of the coronal heating mechanism is by comparing simple theoretical one-dimensional hydrostatic loop models with observations at the temperature and/or density structure along these features. The most well-known method for dealing with comparisons like that is the χ ² approach. In this paper we consider the restrictions imposed by this approach and present an alternative way for making model comparisons using Bayesian statistics. In order to quantify our beliefs we use Bayes factors and information criteria such as AIC and BIC. Two datasets (Ugarte-Urra et al. 2005; Priest et al. 2000) are reanalyzed using the method described above. For the dataset of Ugarte-Urra et al. (2005), we conclude to apex dominant heating as the likely heating candidate, whereas the dataset of Priest et al. (2000) implies basal heating. Note that these new results are different from those obtained using the chi-squared statistic. For this we suggest that proper usage of Classical and Bayesian statistics should be applied in order to make safe assumptions about the nature of the coronal heating mechanisms.     

Magnetic field of the young star RW Aur

Results of ourmeasurements of the longitudinal magnetic field B _z for the young star RWAur A are presented. B _z measured from the so-called narrow component of the He I 5876 line varies in the range from −1.47 ± 0.15 to +1.10 ± 0.15 kG. Our data are consistent with a stellar rotation period of }~5.6 days and the model of two hot spots with opposite magnetic field polarities spaced about 180° apart in longitude. Relative to the Earth, the spot with B _z < 0 lies in the hemisphere above the midplane of the accretion disk, while the spot with B _z > 0 is below the midplane. The upper limit for B _z (at the 3σ level) obtained by averaging all observations is 180 G for the photosphere and 220 and 230 G for the Hα and [OI] 6300 line formation regions, respectively. We have also failed to detect a field in the formation region of broad emission line components: the upper limit for B _z is 600 G. In two of 11 cases, we have detected a magnetic field in the formation region of the blue absorption wing of the Na I D doublet lines, i.e., in the wind from RW Aur A: B _z = −180 ± 50 and −810 ± 80 G. The radial velocity of the photospheric lines in RW Aur A averaged over all our observations is }~+10.5 km s⁻¹, i.e., a value lower than that obtained by Petrov et al. (2001) ten years earlier by 5.5 km s⁻¹. Therefore, we discuss the possibility that RW Aur is not a binary but a triple system.     

Period Change of CU Pegasi

On the basis of the published times of minima and our own observations, we analysed the period change of the Algol-type eclipsing binary CU Pegasi. Over almost seventy years of observations, the parabolic period change has been clearly seen as dP/dt = 1.38 × 10⁻⁶ d/year. The estimated mass transfer in the system is about 1 × 10⁻⁷ M_M⊙/year.     

Galactic rotation curve from the space velocities of selected masers

Based on currently available observations of 28 maser sources in 25 star-forming regions with measured trigonometric parallaxes, proper motions, and radial velocities, we have constructed the rotation curve of the Galaxy. Taking different distances to the Galactic center R ₀, we have estimated the peculiar velocity of the Sun, the angular velocity of Galactic rotation, and its three derivatives. For R ₀ = 8 kpc, we have found the circular velocity of the Sun to be V ₀ = 243 ± 16 km s⁻¹, which corresponds to a revolution period of 202 ± 10 Myr. We have obtained the Oort constants A = 16.9 ± 1.2 km s⁻¹ kpc⁻¹ and B = −13.5 ± 1.4 km s⁻¹ kpc⁻¹. Our simulation of the influence of a spiral density wave has shown that the peculiar velocity of the Sun with respect to the local standard of rest and the component (V _⊙)_LSR depend significantly on the Sun’s phase in the spiral wave.     

Parameters of the local warp of the stellar-gaseous galactic disk from the kinematics of Tycho-2 nearby red giant clump stars

We analyze the three-dimensional kinematics of about 82 000 Tycho-2 stars belonging to the red giant clump (RGC). First, based on all of the currently available data, we have determined new, most probable components of the residual rotation vector of the optical realization of the ICRS/HIPPARCOS system relative to an inertial frame of reference, (ω _x , ω _y , ω _z ) = (−0.11, 0.24, −0.52) ± (0.14, 0.10, 0.16) mas yr⁻¹. The stellar proper motions in the form μ_α cos δ have then be corrected by applying the correction ω _z = −0.52 mas yr⁻¹. We show that, apart from their involvement in the general Galactic rotation described by the Oort constants A = 15.82 ± 0.21 km s⁻¹ kpc⁻¹ and B = −10.87 ± 0.15 km s⁻¹ kpc⁻¹, the RGC stars have kinematic peculiarities in the Galactic yz plane related to the kinematics of the warped stellar-gaseous Galactic disk. We show that the parameters of the linear Ogorodnikov-Milne model that describe the kinematics of RGC stars in the zx plane do not differ significantly from zero. The situation in the yz plane is different. For example, the component of the solid-body rotation vector of the local solar neighborhood around the Galactic x axis is M ₃₂⁻ = −2.6 ± 0.2 km s⁻¹ kpc⁻¹. Two parameters of the deformation tensor in this plane, namely M ₂₃⁺ = 1.0 ± 0.2 km s⁻¹ kpc⁻¹ and M ₃₃ − M ₂₂ = −1.3 ± 0.4 km s⁻¹ kpc⁻¹, also differ significantly from zero. On the whole, the kinematics of the warped stellar-gaseous Galactic disk in the local solar neighborhood can be described as a rotation around the Galactic x axis (close to the line of nodes of this structure) with an angular velocity −3.1 ± 0.5 km s⁻¹ kpc⁻¹ ≤ Ω_W ≤ −4.4 ± 0.5 km s⁻¹ kpc⁻¹.     

INTEGRAL serendipitous detection of the gamma-ray microquasar LS 5039

LS 5039 is the only X-ray binary persistently detected at TeV energies by the Cherenkov HESS telescope. It is moreover a γ-ray emitter in the GeV and possibly MeV energy ranges. To understand important aspects of jet physics, like the magnetic field content or particle acceleration, and emission processes, such as synchrotron and inverse Compton (IC), a complete modeling of the multiwavelength data is necessary. LS 5039 has been detected along almost all the electromagnetic spectrum thanks to several radio, infrared, optical and soft X-ray detections. However, hard X-ray detections above 20 keV have been so far elusive and/or doubtful, partly due to source confusion for the poor spatial resolution of hard X-ray instruments. We report here on deep (∼300 ks) serendipitous INTEGRAL hard X-ray observations of LS 5039, coupled with simultaneous VLA radio observations. We obtain a 20–40 keV flux of 1.1±0.3 mCrab (5.9 (±1.6) ×10⁻¹² erg cm⁻² s⁻¹), a 40–100 keV upper limit of 1.5 mCrab (9.5×10⁻¹² erg cm⁻² s⁻¹), and typical radio flux densities of ∼25 mJy at 5 GHz. These hard X-ray fluxes are significantly lower than previous estimates obtained with BATSE in the same energy range but, in the lower interval, agree with extrapolation of previous RXTE measurements. The INTEGRAL observations also hint to a break in the spectral behavior at hard X-rays. A more sensitive characterization of the hard X-ray spectrum of LS 5039 from 20 to 100 keV could therefore constrain key aspects of the jet physics, like the relativistic particle spectrum and the magnetic field strength. Future multiwavelength observations would allow to establish whether such hard X-ray synchrotron emission is produced by the same population of relativistic electrons as those presumably producing TeV emission through IC.     

Present Evidence for Intermediate Mass Black Holes in ULXs and Future Prospects

In a number of the most luminous ULXs (those with L_X ∼ 10⁴⁰ erg s⁻¹) in nearby galaxies, observations with XMM-Newton and Chandra are revealing evidence which suggests that these ULXs may harbor intermediate-mass black holes (IMBHs). The detection of accretion disk spectral components with temperatures 5–10 times lower than the temperatures observed in stellar-mass black hole binaries near to their Eddington limit may be particularly compelling evidence for IMBH primaries, since T ∝ M^−1/4 for disks around black holes. In some sources, X-ray timing diagnostics also hint at IMBHs. Evidence for IMBHs in a subset of the most luminous ULXs, a discussion of the robustness of this evidence and alternatives to the IMBH interpretation, and prospects for better determining the nature of these sources in the future, are presented in this work.     

Statistical properties of spectroscopic binary stars

As part of a study of the mass-ratio distribution of spectroscopic binary stars, the statistical properties of the systems in theEighth Catalogue of the Orbital Elements of Spectroscopic Binary Stars, compiled by Battenet al. (1989), are investigated.Histograms are presented of the distributions of various parameters of the systems in the catalogue as a whole, and compared to those of the previous edition.Histograms of binaries of various spectral types are presented. It is noted that the early-type binaries in the catalogue have on average higher radial-velocity amplitudes, shorter periods, and smaller semi-major axes than late-type binaries. Late-type binaries have relatively more eccentric orbits. Whether the differences noticed between the early- and late-type binaries have any significance with respect to the population of binary stars in the Galaxy is not clear, because it is very hard to assess the extent to which the catalogue is a statistically representative and complete sample.The distribution of semi-major axesa ₁ sini varies considerably among different subsamples.The mass-ratio distribution of single-lined spectroscopic binaries is, for all (sub-)samples, characterized by a decrease in the number of systems according to a power law asq1, forq=M _sec/M _prim>0.25.The mass-ratio distribution of double-lined spectroscopic binaries (SBII) is, for all (sub-)samples but one, characterized by an increase in numbers according to a power law, asq1. The exception to this general behaviour is the sample of SBII systems with magnitudesm5 ^m , which has its maximum atq0.65.The distributions are presented as they are, without corrections for selection effects.     

Theoretical models for wide binary evolution

Summary Wide binary evolution under the influence of stars, giant molecular clouds, and the galactic tidal field is reviewed. Recent results show that binary lifetimes ata≳0.1 pc are determined by an interplay between relatively gentle evolution due to stellar encounters and catastrophic collisions with cloud subclumps. It is emphasized that binary semimajor axis distributions may not be understood from lifetimes alone. To determine the theoretical distribution binary ‘birthrate functions’ must be convolved with probability distributions for the evolution of an ensemble of binaries from given initial separations. Simple models show no sharp breaks below the tidal cut-off imposed by the Galaxy.     

Photometric Observations and Apsidal Motion Study of V1143 Cyg

Photometric observations of the eccentric eclipsing binary V1143 Cyg were performed during Aug.–Sep. 2000 and July 2002, in Johnson B and V bands. The analysis of both light curves was made separately using the 1998 version of Wilson’s LC code. In order to find a new observed rate of apsidal motion, we followed the procedure described by Guinan and Maloney (1985). A new observed rate of apsidal motion of 3^∘.72/100 yr was computed, which is close to the one reported earlier by Khaliullin (1983), Gimenez and Margrave (1985), and Burns et al. (1996).     

New CCD photometry for the extreme lower mass-ratio binary AW Ursae Majoris

This paper presents charge-couple device (CCD) photometric observations for the eclipsing binary AW UMa. The V-band light curve in 2007 was analyzed using the 2003 version of the Wilson–Devinney code. It is confirmed that AW UMa is a total eclipsing binary with a higher degree of contact f=80.2% and a lower mass ratio of q=0.076. From the (O−C) curve, the orbital period shows a continuous period decrease at a rate of dP/dt=−2.05×10⁻⁷ d yr⁻¹. The long-term period decrease suggested that AW UMa is undergoing the mass transfer from the primary component to the secondary one, accompanied by angular momentum loss due to mass outflow L ₂. Weak evidence indicates that there exists a cyclic variation with a period of 17.6 yr and a small amplitude of A=0. ^d 0019, which may be attributed to the light-time effect via the third body. If the existence of an additional body is true, it may remove a great amount of angular momentum from the central system. For this kind of contact binary, as the orbital period decreases, the shrinking of the inner and outer critical Roche lobes will cause the contact degree f to increase. Finally, this kind of binary will merge into a single rapid-rotation star.     

Observations of double stars at Pulkovo with 65 cm Zeiss refractor

Some results of the photographic observations of double stars with 65 cm refractor of Pulkovo observatory are presented. We use the apparent motion parameters (AMP) method which allows to determine the orbits and to carry out the dynamical investigation of wide binaries on the basis of a short arc of their orbital motion. We have determined more than 40 orbits for wide pairs and also the sum of masses and in some cases—the mass-ratio of components. The references to our works and the basic results of observations are contained in Kisselev et al. [2004. Catalogue of relative positions of visual double stars made on the observations with 26^″ refractor of Pulkovo observatory. Strassbourg, I/297]. We apply two ways of revealing the hidden mass of our stars, namely: revealing of possible perturbations from comparison of observational and calculated positions using differences O-C (for instance, perturbations in the orbital motion of ADS 15571) and also by means of comparison of the sum of the masses obtained by us and the sum of the masses obtained by means of the mass-luminosity relation. An excess of masses of about 1-3 solar masses is detected for binaries: ADS 497, ADS 8450 and ADS 10329 by means of last method.The estimations of the masses for some binaries are discussed. Also we justify the necessity of precise parallaxes and relative radial velocities of stars, which could be measured by space telescopes such as the GAIA as the additional parameters for determination of orbits of binaries.     

The Large-Scale, Decelerating X-ray Jets from the Microquasar Xte J1550−564: Evidence for External Shocks Caused by the Jet-Ism Interaction?

Large-scale, decelerating, relativistic X-ray jets from microquasar XTE J1550−564 has been recently discovered with Chandra by Corbel et al. (2002). We find that the dynamical evolution of the approaching jet at the late time is consistent with the well-known Sedov evolutionary phase R∝ t ^2/5. A trans-relativistic external shock dynamic model by analogy with the evolution of gamma-ray burst remnants, is shown to be able to fit the proper-motion data of the approaching jet reasonably well. The inferred interstellar medium density around the source is well below the canonical value n _ISM∼1 cm⁻³. The rapidly fading X-ray emission can be interpreted as synchrotron radiation from the non-thermal electrons in the adiabatically expanding ejecta. These electrons were accelerated by the reverse shock (moving back into the ejecta) which becomes important when the inertia of the swept external matter leads to an appreciable slowing down of the original ejecta.     

WR 140 (=V1687 Cyg): Infrared photometry, 2001–2010

We present the results of our infrared observations of WR 140 (=V1687 Cyg) in 2001–2010. Analysis of the observations has shown that the J brightness at maximum increased near the periastron by about 0 ^m .3; the M brightness increased by ∼2 ^m in less than 50 days. The minimum J brightness and the minimum L and M brightnesses were observed 550–600 and 1300–1400 days after the maximum, respectively. The JHKLM brightness minimum was observed in the range of orbital phases 0.7–0.9. The parameters of the primary O5 component of the binary have been estimated to be the following: R(O5) ≈ 24.7R _⊙, L(O5) ≈ 8 × 10⁵ L _⊙, and M _bol(O5) ≈ −10 ^m . At the infrared brightness minimum, T _g ∼ 820–880 K, R _g ≈ 2.6 × 10⁵ R _⊙, the optical depth of the shell at 3.5 μm is ∼5.3 × 10⁻⁶, and its mass is ≈1.4 × 10⁻⁸ M _⊙. At the maximum, the corresponding parameters are ∼1300 K, 8.6 × 10⁴ R _⊙, ∼2 × 10⁻⁴, and ∼6 × 10⁻⁸ M _⊙; the mean rate of dust inflow (condensation) into the dust structure is ∼3.3 × 10⁻⁸ M _⊙ yr⁻¹. The mean escape velocity of the shell from the heating source is ∼10³ km s⁻¹ and the mean dispersal rate of the shell is ∼1.1 × 10⁻⁸ M _⊙ yr⁻¹.