Close binary systems among very low‐mass stars and brown dwarfs

Using Monte Carlo simulations and published radial velocity surveys we have constrained the frequency and separation (a ) distribution of very low‐mass star (VLM) and brown dwarf (BD) binary systems.We find that simple Gaussian extensions of the observed wide binary distribution, with a peak at 4AU and 0.6 < σ _{log(a /AU)} < 1.0, correctly reproduce the observed number of close binary systems, implying a close (a < 2.6 AU) binary frequency of 17–30% and overall frequency of 32–45%. N‐body models of the dynamical decay of unstable protostellar multiple systems are excluded with high confidence because they do not produce enough close binary VLMs/BDs. The large number of close binaries and high overall binary frequency are also completely inconsistent with published smoothed particle hydrodynamical modelling and argue against a dynamical origin for VLMs/BDs. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Absolute Parameters of Early-Type Close Binaries in the LMC

A project is presented which aims at high-precision determination of absolute parameters of close early-type eclipsing binaries in the LMC. We will use multi-object spectrographs (MOS) to measure RV curves of a large number of program stars selected from the MACHO archive. Spectroscopic mass ratios will be used as input for our light curve analysis code MORO. Application of the Simplex-based algorithm FITSB2 will achieve spectrum disentangling and fitting of NLTE model (TLUSTY) atmospheres, yielding orbital and atmospheric parameters. The method was extensively tested by application to time series of synthetic binary spectra for the expected range of S/N and MOS instrumental resolution. log g and T _eff were reestablished with a precision of better than 5%, radial velocity amplitudes with errors of <3%. An important by-product of our project will be the improvement of the distance modulus of the LMC, a topic which is still being intensely discussed.     

Radial Velocity and Light Curve Analysis of the Eclipsing Binary Nn Vir

The eclipsing binary NN Vir is a short period system showing an EW-type light curve. Photometric observations of NN Vir were done by Gomez–Ferrellad and Garcia–Melendo (1997) at Esteve Duran Observatory. The first spectroscopic observations of this system were obtained by Rucinski and Lu (1999). The radial velocity and light curves analysis was made with the latest version of the Wilson program (1998), and the geometric and physical elements of the system are derived. From the simultaneous solutions of the system, we determined the masses and radii of the components: 1.89 M _⊙ and 1.65 R _⊙ for the primary component; 0.93 M _⊙ and 1.23 R _⊙ for the secondary component. We estimated effective temperatures of 7030 K for the primary and 6977 K for the secondary component.     

Existence and stability of symmetric periodic simultaneous binary collision orbits in the planar pairwise symmetric four-body problem

We analytically prove the existence of a symmetric periodic simultaneous binary collision orbit in a regularized planar pairwise symmetric equal mass four-body problem. This is an extension of our previous proof of the analytic existence of a symmetric periodic simultaneous binary collision orbit in a regularized planar fully symmetric equal mass four-body problem. We then use a continuation method to numerically find symmetric periodic simultaneous binary collision orbits in a regularized planar pairwise symmetric 1, m, 1, m four-body problem for m between 0 and 1. Numerical estimates of the the characteristic multipliers show that these periodic orbits are linearly stability when 0.54 ≤ m ≤ 1, and are linearly unstable when 0 < m ≤ 0.53.     

A general catalogue of ephemerides and apparent orbits for visual binaries

Summary Informations on 736 pairs of visual binaries are given in the form of the Tables. The General Catalogue gives ephemerides (t,θ°,ρ″) in 20 years (1984–2003) for each pair with the figures of its apparent elliptical orbit where the positions of secondary component relative to the primary one at different epochs are indicated.The General Catalogue contains four parts: Part one — Source and Grading of Orbit; Part two — Ephemérides and Atlas of Apparent Orbits; Part three — Classifications of Visual Binary Stars; Part four — Indexes of Visual Binary Stars. The principle of calculation and the statistical data are presented in this paper. There are seven statistical tables, giving the elements distribution of true and apparent orbits, grading distribution of orbits, number distribution with different physical property of component of binary star. The number of binary stars in anyone constellation, the number of binary stars brighter than 6^m.5. The number of binary stars nearer than 25 parsec.     

Orbital elements of the RS CVn eclipsing binary,SV Camelopardalis

The extensiveUBV observations of SV Camelopardalis by Patkos (1982) have been analysed to derive the orbital elements of the system. The data were corrected for the effect of third body (Sarma, Sarma & Abhyankar 1985) and for the ‘RS CVn’ distortion wave (Sarma, Vivekanandarao & Sarma 1988). The cleaned data were used to obtain a preliminary solution by a modified version of Wellmann method (Sarma & Abhyankar 1979) from which we concluded that the primary eclipse is a transit. The final orbital elements of SV Cam were obtained by the modified version (Sarma 1988; Sarmaet al. 1987) of WINK program by Wood (1972). The colour and median brightness variation are discussed. From the spectroscopic mass functionf(m) = 0.118 M_⊙ (Hiltner 1953), the absolute dimensions of the components are found to be 0.826 M_bd & 0.592 M_⊙ and 1.236 R_⊙ & 0.778 R_⊙ for the primary and secondary components, respectively. The age of the binary system is estimated to be 6.0 ± 1.0 × 10⁸ years     

The Missing Intermediate Luminosity Galaxies in the Leo I Group

We present results from a large scale R-band survey of the nearby Leo I group to measure the luminosity function down to dwarf spheroidal luminosities and surface brightnesses. This program has utilized>7 square degrees of imaging with the Mosaic camera at the KPNO 0.9mand extensive, on-going spectroscopic follow-up. We estimate our completeness from simulations to be 90% at μ_R(0) ≃ 24.5, M _R ≃ –11.Leo I is a nearby, low-density group for which a steep faint-end slope has been suggested (Ferguson and Sandage, 1991).We find the luminosity function to be more similar to that of the Local Group at the faint-end (flat).We also find an unusual gap in the luminosity function at intermediate luminosities, –19.5 < M _R < –16. This gap cannot be explained by pure Poisson fluctuations around a typical Schechter function, nor is it likely to result from photometric incompleteness. This revised version was published online in July 2006 with corrections to the Cover Date.     

Stability criterion for a light binary attracted by a heavy body

Dynamical behaviour of a small binary with equal components, each of mass m, is considered under attraction of a heavy body of mass M. Differential equations of the general three-body problem are integrated numerically using the code by S. J. Aarseth (Aarseth, Zare 1974) for mass ratios m/M within 10⁻¹¹–10⁻⁴ range. The direct and retrograde orbits of light bodies about each other are considered which lie either in the plane of moving their center of mass or in the plane perpendicular to it. It is shown numerically that the critical separation between the binary components which leads to disruption of binary is proportional to (m/M)^1/3. The criterion can be used for studying (in the first approximation) the motion of double stars and binary asteroids or computing the parameters of magnetic monopol and antimonopol pairs.     

Generalized Chaplygin gas model with or without viscosity in the w–w′ plane

In this paper, we investigate the dynamics of generalized Chaplygin gas (GCG) model with or without viscosity in the w–w′ plane, which is defined by the equation of state parameter and its time derivative with respect to the logarithm of the scale factor. We show that GCG model without viscosity approaches to a late time de Sitter attractor (w _g =−1) and behaves like a “freezing” scalar field for the parameter α constrained by the latest observational data. However, introducing viscosity exerts an influence on the evolution of w and affects the location of the late time attractor (w _g >−1) in viscous GCG model. We also find numerically such a transition from w′>0 to w′<0 as the universe expands in viscous GCG model different from GCG model without viscosity (w′<0) in the w–w′ plane.     

H.E.S.S. observations of LS 5039

Recent observations of the binary system LS 5039 with the High Energy Stereoscopic System (H.E.S.S.) revealed that its Very High Energy (VHE) γ-ray emission is modulated at the 3.9 days orbital period of the system. The bulk of the emission is largely confined to half of the orbit, peaking around the inferior conjunction epoch of the compact object. The flux modulation provides the first indication of γ-ray absorption by pair production on the intense stellar photon field. This implies that the production region size must be not significantly greater than the gamma-gamma photosphere size (∼1 AU), thus excluding the large scale collimated outflows or jets (extending out to ∼1000 AU). A hardening of the spectrum is also observed at the same epoch between 0.2 and a few TeV which is unexpected under a pure absorption scenario and could rather arise from variation with phase in the maximum electron energy and/or the dominant VHE γ-ray production mechanism. This first-time observation of modulated γ-ray emission allows precise tests of the acceleration and emission models in binary systems. Mathieu de Naurois for the H.E.S.S. Collaboration.     

Interactions between brown-dwarf binaries and Sun-like stars

Several mechanisms have been proposed for the formation of brown dwarfs, but there is as yet no consensus as to which—if any—are operative in nature. Any theory of brown dwarf formation must explain the observed statistics of brown dwarfs. These statistics are limited by selection effects, but they are becoming increasingly discriminating. In particular, it appears (a) that brown dwarfs that are secondaries to Sun-like stars tend to be on wide orbits, a≳100 AU (the Brown Dwarf Desert), and (b) that these brown dwarfs have a significantly higher chance of being in a close (a≲10 AU) binary system with another brown dwarf than do brown dwarfs in the field. This then raises the issue of whether these brown dwarfs have formed in situ, i.e. by fragmentation of a circumstellar disc; or have formed elsewhere and subsequently been captured. We present numerical simulations of the purely gravitational interaction between a close brown-dwarf binary and a Sun-like star. These simulations demonstrate that such interactions have a negligible chance (<0.001) of leading to the close brown-dwarf binary being captured by the Sun-like star. Making the interactions dissipative by invoking the hydrodynamic effects of attendant discs might alter this conclusion. However, in order to explain the above statistics, this dissipation would have to favour the capture of brown-dwarf binaries over single brown-dwarfs, and we present arguments why this is unlikely. The simplest inference is that most brown-dwarf binaries—and therefore possibly also most single brown dwarfs—form by fragmentation of circumstellar discs around Sun-like protostars, with some of them subsequently being ejected into the field.     

Trigonometric parallaxes of 29 stars with large proper motions

While implementing the first stage of the Pulkovo program of research on stars with large proper motions, we determined the trigonometric parallaxes of 29 stars (12 ^m < V < 16 ^m ) based on CCD observations with a 26-inch refractor. The mean standard error was 3.7 mas. Comparison of the Pulkovo parallaxes with those obtained at the Observatory of the Yale University and the US Naval Observatory (USNO) has shown that the parallax differences (Pulkovo-Yale/USNO) lie within the limits of their measurement errors in an overwhelming majority of cases. On average, they are −0.6 ± 1.0 mas. No systematic dependences on stellar distance, magnitude, and color in this set of differences have been found. Our comparisons show that the observing and data reduction techniques used in the Pulkovo program of research on fast stars allow highly accurate trigonometric parallaxes of these objects to be obtained. All program stars are within 50 pc of the Sun; most of them belong to the immediate solar neighborhood (D < 25 pc). For two stars (J0522+3814 and J1202+3636), the trigonometric parallaxes have been determined for the first time.     

Models for the formation of binary and millisecond radio pulsars

The peculiar combination of a relatively short pulse period and a relatively weak surface dipole magnetic field strength of binary radio pulsars finds a consistent explanation in terms of (i) decay of the surface dipole component of neutron-star magnetic fields on a timescale of (2–5) × 10⁶ yr, in combination with (ii) spin-up of the rotation of the neutron star during a subsequent mass-transfer phase. The four known binary radio pulsars appear to fall into two different categories. Two of them, PSR 0655 + 64 and PSR 1913 + 16, have short orbital periods (<25 h) and high mass functions, indicating companion masses 0.7M⊙ (∼1 (± 0.3) M⊙ and 1.4 M⊙, respectively). The other two, PSR 0820 + 02 and PSR 1953 + 29, have long orbital periods (117d), nearly circular orbits, and low, almost identical mass functions of about 3×10^-3 M⊙, suggesting companion masses of about 0.3M⊙. It is pointed out that these two classes of systems are expected to be formed by the later evolution of binaries consisting of a neutron star and a normal companion star, in which the companion was (considerably) more massive than the neutron star, or less massive than the neutron star, respectively. In the first case the companion of the neutron star in the final system will be a massive white dwarf, in a circular orbit, or a neutron star in an eccentric orbit. In the second case the final companion to the neutron star will be a low-mass (∼ 0.3 M⊙) helium white dwarf in a wide and nearly circular orbit. In systems of the second type the neutron star was most probably formed by the accretion-induced collapse of a white dwarf. This explains in a natural way why PSR 1953 + 29 has a millisecond rotation period and PSR 0820 + 02 has not. Among the binary models proposed for the formation of the 1.5-millisecond pulsar, the only ones that appear to be viable are those in which the companion disappeared by coalescence with the neutron star. In such models the companion may have been a red dwarf of mass 0.03M⊙, a neutron star, or a massive (>0.7M⊙) white dwarf. Only in the last-mentioned case is a position of the pulsar close to the galactic plane a natural consequence. In the first-mentioned case the progenitor system most probably was a cataclysmic-variable binary in which the white dwarf collapsed by accretion.     

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.     

Supernovae-optical precursors of short gamma-ray bursts

The probability of observing “supernova-gamma-ray burst” (GRB) pair events and recurrent GRBs from one galaxy in a time interval of several years has been estimated. Supernova explosions in binary systems accompanied by the formation of a short-lived pair of compact objects can be the sources of such events. If a short GRB is generated during the collision of a pair, then approximately each of ∼300 short GRBs with redshift z must have an optical precursor—a supernova in the observer’s time interval <2(1 + z) yr. If the supernova explosion has the pattern of a hypernova, then a successive observation of long and short GRBs is possible. The scenario for the generation of multiple GRBs in collapsing galactic nuclei is also discussed.     

A new class of charged analogues of Vaidya–Tikekar type super-dense star

First ever closed form solution for charged fluid sphere expressed by a space time with its hypersurfaces t= constant as spheroid is obtained for the case 0<K<1. The same is utilized to construct a superdense star with surface density 2×10¹⁴ gm/cm³. The star is seen to satisfy the reality and causality conditions for 0<K≤0.045 and possesses maximum mass and radius to be 0.065216M _Θ and 1.137496 km respectively. Moreover the interior of the star satisfy strong energy condition. However in the absence of the causality condition, the reality conditions are valid for a wider range 0<K≤0.13. The maximum mass and radius for the later case are 1.296798M _Θ and 2.6107 km respectively for the strong energy condition, while the said parameters for the weak energy condition read as 1.546269M _Θ and 2.590062 km respectively.     

First CCD photometry for the contact binary VZ Trianguli

First CCD photometry is presented for the eclipsing binary VZ Trianguli, observed at the Sheshan Station of Shanghai Astronomical Observatory in 2008. Using the Wilson–Devinney Code, the photometric solution of VZ Tri was first deduced from the R-band observations. The results show that VZ Tri is an A-subtype late-type contact binary, with a mass ratio of q=0.350(±0.004) and a low contact degree of f=27.9%(±1.0%). Based on all available light minimum times covering over 40 years, it is found that the orbital period shows a long-term decrease at a rate of dP/dt=−1.52(±0.03)×10⁻⁷ d yr⁻¹, suggesting that VZ Tri is undergoing mass transfer from the more massive component to the less massive component, accompanied with angular momentum loss. With period decreasing, the inner and outer critical Roche lobes will shrink, and then cause the contact degree to increase. Therefore, the weak-contact binary VZ Tri with decreasing period may evolve into a deep-contact configuration.     

Mv Lyr: Transition from Low to High Brightness State

We present results of multicolour photometry of the nova-like binary MV Lyr in the years 2002–2003, corresponding to the transition of the star from its low (V = 17.8) to high (V = 12.3) brightness state. The transition lasted at most 200 days. During the first ∼50 days MV Lyr brightened at the rate 0.06 mag/day, and during the next ∼150 days at the rate 0.01 mag/day. The brightening was accompanied by a blue shift from 0.45 to 0.05 in V−R, interpreted as an increase of the accretion disk contribution to the total light. During the transition MV Lyr displayed well-known ‘quasi-orbital’ light variations and fast quasi-periodic oscillations with a typical time scale of tens of minutes.     

Examination of the Association between Quasars and Abell Clusters on All Sky

Using the Hewitt-Burbidge QSO Catalogue (1993) and all-sky catalogue of Abell clusters (ACO, 1989) at the region |b| > 40^° we analyze the cross correlation function and find anti-correlation between them at angular separations 3^° < θ < 10^° , which is mainly caused by optical-selected QSOs, rather than radio-selected QSOs. There is no such anti-correlation between QSOs and Abell clusters at smaller separations θ < 3^°. Considering that this phenomenon may be caused by different characters of the objects, we further estimate the correlation function with various subsamples. We find that the correlation is independent of the redshift of QSOs, but depends upon the type of Abell clusters: for the D ≤ 4 clusters there is an obvious tendency of overdensity of quasars at 0^° < θ < 5^°; around the R ≥ 2 Abell clusters there is about an 18.7% deficit of quasars in the region 3^° < θ < 7^°. K-S Test shows the overdensity or deficit of quasars around different types of clusters cannot be explained by the projection effect of background quasars. We get the enhancement factor of quasar overdensity (for D ≤ 4 clusters) q =1.13, and the extinction magnitude factor of QSO deficiency (for R ≥ 2 clusters) Av= 0.14. This revised version was published online in July 2006 with corrections to the Cover Date.     

Fast computation of complete elliptic integrals and Jacobian elliptic functions

As a preparation step to compute Jacobian elliptic functions efficiently, we created a fast method to calculate the complete elliptic integral of the first and second kinds, K(m) and E(m), for the standard domain of the elliptic parameter, 0 < m < 1. For the case 0 < m < 0.9, the method utilizes 10 pairs of approximate polynomials of the order of 9–19 obtained by truncating Taylor series expansions of the integrals. Otherwise, the associate integrals, K(1 − m) and E(1 − m), are first computed by a pair of the approximate polynomials and then transformed to K(m) and E(m) by means of Jacobi’s nome, q, and Legendre’s identity relation. In average, the new method runs more-than-twice faster than the existing methods including Cody’s Chebyshev polynomial approximation of Hastings type and Innes’ formulation based on q-series expansions. Next, we invented a fast procedure to compute simultaneously three Jacobian elliptic functions, sn(u|m), cn(u|m), and dn(u|m), by repeated usage of the double argument formulae starting from the Maclaurin series expansions with respect to the elliptic argument, u, after its domain is reduced to the standard range, 0 ≤ u < K(m)/4, with the help of the new method to compute K(m). The new procedure is 25–70% faster than the methods based on the Gauss transformation such as Bulirsch’s algorithm, sncndn, quoted in the Numerical Recipes even if the acceleration of computation of K(m) is not taken into account.