Cosmic rays from binary neutron stars

Within the more than 30 yr of cosmic ray astrophysics, neither their origin nor their precise mode of propagation have found undisputable explanations. Among the favoured boosters have been point sources, like supernovae and pulsars, as well as extended sources, like cosmic clouds and supernova remnants. Extended sources have been proposed by Fermi (1949), and pushed more recently by a number of investigators because of the huge available reservoirs, and because repetitive shock acceleration can generate power law spectra which are similar to the ones observed (Axfordet al., 1977; Bell, 1978; Blandford and Ostriker, 1978; Krymsky, 1977). Yet the shock acceleration model cannot easily be adjusted to achieve particle energies in excess of some critical energy, of order 10^4±1 GeV (Völket al., 1981). For this and several other reasons, the suggestion is revived that neutron stars are the dominant source of high-energy cosmic rays. To be more precise: the (relativistic) ionic component of the cosmic rays is argued to be injected by young binary neutron stars (?10⁵ yr) whose rotating magnetospheres act like grindstones in the wind of their companion (Kundt, 1976). The high-energy (?30 GeV) electron-positron component may be generated by young pulsars (?10⁵ yr) and by collision processes, and the electron component below 30 GeV predominantly by supernova remnants.     

Magnetic field evolution in neutron stars

Neutron stars contain persistent, ordered magnetic fields that are the strongest known in the Universe. However, their magnetic fluxes are similar to those in magnetic A and B stars and white dwarfs, suggesting that flux conservation during gravitational collapse may play an important role in establishing the field, although it might also be modified substantially by early convection, differential rotation, and magnetic instabilities. The equilibrium field configuration, established within hours (at most) of the formation of the star, is likely to be roughly axisymmetric, involving both poloidal and toroidal components. The stable stratification of the neutron star matter (due to its radial composition gradient) probably plays a crucial role in holding this magnetic structure inside the star. The field can evolve on long time scales by processes that overcome the stable stratification, such as weak interactions changing the relative abundances and ambipolar diffusion of charged particles with respect to neutrons. These processes become more effective for stronger magnetic fields, thus naturally explaining the magnetic energy dissipation expected in magnetars, at the same time as the longer-lived, weaker fields in classical and millisecond pulsars. (© 2007 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The spin evolution of nascent neutron stars

The loss of angular momentum owing to unstable r-modes in hot young neutron stars has been proposed as a mechanism for achieving the spin rates inferred for young pulsars. One factor that could have a significant effect on the action of the r-mode instability is fallback of supernova remnant material. The associated accretion torque could potentially counteract any gravitational-wave-induced spin-down, and accretion heating could affect the viscous damping rates and hence the instability. We discuss the effects of various external agents on the r-mode instability scenario within a simple model of supernova fallback on to a hot young magnetized neutron star. We find that the outcome depends strongly on the strength of the magnetic field of the star. Our model is capable of generating spin rates for young neutron stars that accord well with initial spin rates inferred from pulsar observations. The combined action of r-mode instability and fallback appears to cause the spin rates of neutron stars born with very different spin rates to converge, on a time-scale of approximately 1 year. The results suggest that stars with magnetic fields ≤10¹³ G could emit a detectable gravitational wave signal for perhaps several years after the supernova event. Stars with higher fields (magnetars) are unlikely to emit a detectable gravitational wave signal via the r-mode instability. The model also suggests that the r-mode instability could be extremely effective in preventing young neutron stars from going dynamically unstable to the bar-mode.     

Broadband gravitational-wave pulses from binary neutron stars in eccentric orbits

The gravitational-wave radiation from binary stars in elliptical orbits peaks at times close to the periastron passage. For a stationary distribution of binary neutron stars in the Galaxy, there are several systems with large orbital eccentricities and periods in the range from several tens of minutes to several days from which gravitational-wave radiation at periastron will be observed as a broad pulse in the frequency range 1–100 mHz. The LISA space interferometer will be able to record pulsed signals from these systems at a signal-to-noise ratio $S/N > 5\sqrt 5$ in the frequency range ～10^?3–10^?1 Hz. Algorithms for detecting such signals are discussed.     

Properties of Wolf-Rayet stars in eclipsing binary systems

The results of investigations of a number of eclipsing Wolf-Rayet binaries are presented. The ‘core’ radiuses, the ‘core’ temperatures and masses of WR stars in the eclipsing WR+OB binary systems V 444 Cyg, CX Cep, CQ Cep, and CV Ser are obtained (see Table I). The results obtained from the light curves analysis of the V 444 Cyg in the range λλ2460 Å-3.5μ give strong evidence for the Beals (1944) model of WR phenomenon. The chromospheric-coronal effects in the WN5 extended atmosphere are not observed up to a distance ofr?20R _⊙. In the Hertzsprung—Russell diagram all the WR stars lie on the left side from the main sequence between the main sequence and the sequence of uniform helium stars (see Figure 9). Their locations are close to those of the helium remnants formed as a result of mass exchange in massive close binary systems. The period variations in the systems V 444 Cyg and CQ Cep have been discovered and a reliable value of the mass loss rateM=10^?5 M _⊙yr^?1 is obtained, for the two WR stars. The results of the photometric and spectroscopic investigations of the WR stars with low mass companions (post X-ray binary stage?) are presented too (see Table II). The masses of the companions are (1–2)M _⊙, their optical luminosity is ～10³⁶, erg s^?1 which implies that these companions cannot be the normal stars. It is possible that these companions are neutron stars accreting from the stellar wind of the WR stars. Low values of the X-ray luminosities of such WR stars with low mass companions imply that the accretion of matter in such systems is distinct from the accretion process in classical X-ray binary systems. It is noted also that the parameters of low massive companions coupled with WR stars are close to those of helium stars.     

Magnetic and spin evolution of neutron stars in close binaries

The evolution of neutron stars in close binary systems with a low-mass companion is considered, assuming the magnetic field to be confined within the solid crust. We adopt the standard scenario for the evolution in a close binary system, in which the neutron star passes through four evolutionary phases ('isolated pulsar'–'propeller'– accretion from the wind of a companion – accretion resulting from Roche-lobe overflow). Calculations have been performed for a great variety of parameters characterizing the properties of both the neutron star and the low-mass companion. We find that neutron stars with more or less standard magnetic field and spin period that are processed in low-mass binaries can evolve to low-field rapidly rotating pulsars. Even if the main-sequence life of a companion is as long as 10¹⁰ yr, the neutron star can maintain a relatively strong magnetic field to the end of the accretion phase. The model that is considered can account well for the origin of millisecond pulsars.     

The evolution of the magnetic fields of neutron stars

Observational evidence, and theoretical models of the magnetic field evolution of neutron stars is discussed. Observational data indicates that the magnetic field of a neutron star decays significantly only if it has been a member of a close interacting binary. Theoretically, the magnetic field evolution has been related to the processing of a neutron star in a binary system through the spin evolution of the neutron star, and also through the accretion of matter on the neutron star surface. I describe two specific models, one in which magnetic flux is expelled from the superconducting core during spin-down, via a copuling between Abrikosov fluxoids and Onsager-Feynman vortices; and another in which the compression and heating of the stellar crust by the accreted mass drastically reduces the ohmic decay time scale of a magnetic field configuration confined entirely to the crust. General remarks about the behaviour of the crustal field under ohmic diffusion are also made.     

Thermal evolution of neutron stars with decaying magnetic fields

Rotochemical heating originates in the deviation from beta equilibrium due to spin-down compression, which is closely related to the dipole magnetic field. We numerically calculate the deviation from chemical equilibrium and thermal evolution of neutron stars with decaying magnetic fields. We find that the power-law long term decay of the magnetic field slightly affects the deviation from chemical equilibrium and surface temperature. However, the magnetic decay leads to older neutron stars that could have a different surface temperature with the same magnetic field strength. That is, older neutron stars with a low magnetic field(108G) could have a lower temperature even with rotochemical heating in operation, which probably explains the lack of other observations on older millisecond pulsars with higher surface temperature,except millisecond pulsar J0437–4715.     

Statistical study of Am stars in spectroscopic binary systems

Using the Am components in spectroscopic binaries, the region of metallicism in the (logm, logR) diagram is statistically discussed. The metallic-line characteristics appear in a slightly evolved stage near the Main Sequence within a belt characterized by logg=3.8–4.2, for which the radiiR/R =1.5–2.7. The distribution in the diagram indicates that appreciable differences should exist in the abundances as well as in the mixing lengths of the atmospheres of the Am components.Synchronism between rotation and revolution for the Am components is discussed, and it is clearly found that the synchronism holds accurately for the Am binaries with periods less than about six days, but for those with periods between six and ten days some Am components corotate and some do not, perhaps depending upon their ages.Using the Am spectroscopic binaries with periods less than twelve days, we find that the excess of metallicity is statistically correlated with the rotational velocity; and it is confirmed that the metallicity grossly decreases with increasing rotational velocity. If the effect of a non-Am secondary component is eliminated from a [m ₁]-value for the combined light, the correlation between metallicity excess and rotational velocity for the Am components only should approach [m ₁]/V (km s^–1)=–0.00077.     

Magnetic field evolution of accreting neutron stars – III

The evolutionary scenario of a neutron star magnetic field is examined assuming a spin-down induced expulsion of magnetic flux originally confined to the core, in a case in which the expelled flux undergoes ohmic decay. The nature of field evolution, for accreting neutron stars, is investigated incorporating the crustal microphysics and material movement resulting from accretion. This scenario may explain the observed field strengths of neutron stars but only if the crustal lattice contains a large amount of impurity, which is in direct contrast to the models that assume an original crustal field.     

The gravity-darkening of highly distorted stars in close binary systems

From a quantitative analysis of the photometric ellipticity effect, values of the exponent of gravity-darkening have been empirically deduced for the components of sixteen W UMa systems with accurate cross-correlation spectroscopic mass ratios. In the analysis, the exponents for both components of a W UMa system are assumed to be the same. In view of the generally unstable nature of their light curves and the consequent difficulty to find an unique set of the fractional radii of both components, -values have been derived for various equipotential surfaces with boundary cases set at the inner and outer Roche lobes.The results indicate that if the fill-out factors are less than 0.5 (as is so for most of the actual cases), the deduced -values are all greater than Lucy's conventional value of 0.32. The empirical -values as a whole are found to increase with the mean effective temperature of the components and, furthermore, they show a tendency to decrease with increasing mass ratio, particularly for A-type systems. Thus, a W UMa system consisting of two components of quite different masses has a large empirical -value. These observational evidences indicate that the structure of the photosphere of contact components of W UMa systems is significantly different from that of normal Main-Sequence stars of similar spectral type and may provide a clue to solve puzzles of W UMa stars.     

Tidal evolution in close binary systems

The aim of the present paper will be to give a mathematical outline of the theory of tidal evolution in close binary systems of secularly constant total momentum — an evolution activated by viscous friction of dynamical tides raised by the two components on each other. The first section contains a general outline of the problem; and in Section 2 we shall establish the basic expressions for the energy and momenta of close binaries consisting of components of arbitrary internal structure. In Section 3 we shall investigate the maximum and minimum values of the energy (kinetic and potential) which such systems can attain for given amount of total momentum; while in Section 4 we shall compare these results with the actual facts encountered in binaries with components whose internal structure (and, therefore, rotational momenta) are known to us from evidence furnished by the observed rates of apsidal advance.The results show that all such systems — be these of detached or semi-detached type — disclose that more than 99% of their total momenta are stored in the orbital momentum. The sum of the rotational momenta of the constituent components amounts to less than a percent of the total — a situation characteristic of a state close to the minimum energy for given total momentum. This appears, moreover, to be true not only of the systems with both components on the Main Sequence, but also of those possessing evolved components in contact with their Roche limits.Under such conditions, a synchronism between rotation and revolution (characteristic of both extreme states of maximum and minimum energy) is not only possible, but appears to have been actually approached — if not attained — in the majority of cases. In other words, it would appear that — in at least a large majority of known cases — the existing close binaries have already attained orbits of maximum distension consistent with their momenta; and tidal evolution alone can no longer increase the present separations of the components to any appreciable extent.The virtual absence, in the sky, of binary systems intermediate between the stages of maximum and minimum energy for given momentum leads us to conjecture that the process of dynamical evolution activated by viscous tides may enroll on a time-scale which is relatively short in comparison with their total age — even for systems like Y Cygni or AG Persei, whose total age can scarcely exceed 10⁷ yr. A secular increase of the semi-major axes of relative orbits is dynamically coupled with a corresponding variation in the velocity of axial rotation of both components through the tidal lag arising from the viscosity of stellar material. The differential equations of so coupled a system are given in Section 5; but their solution still constitutes a task for the future.The Lunar Science Institute Contribution No. 90. The Lunar Science Institute is operated by the Universities Space Research Association under Contract No. NSR 09-051-001 with the National Aeronautics and Space Administration.     

Analytic formulas for the mass-transfer rate and the evolution of a close binary system of neutron (degenerate) stars

We derive approximate analytic relations between the mass-transfer rate in a close binary system described in terms of the Roche potential and its basic parameters, such as the total mass of the binary, the radius of its circular orbit, the mass of the mass-losing component, and the degree of its Roche lobe overfilling. Using simplifying assumptions (conservative mass transfer, a short relaxation time of matter on the mass-gaining component compared to the mass-transfer time scale, adiabaticity and quasi-stationarity of the mass flow through the Lagrangian point L ₁) allows the evolution of a binary system of neutron (degenerate) stars to be described in terms of two ordinary differential equations. This makes it possible to qualitatively analyze the evolution process, which is useful in those cases where the evolution of a close binary system must be investigated in general terms, for example, in terms of the scenario for the transformation of the collapse of a rotating presupernova core into a supernova explosion proposed by Imshennik and Nadyozhin (1992) and Imshennik (1992).     

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.     

Signal of quark deconfinement in thermal evolution neutron stars with deconfinement heating

As neutron stars spin-down and contract, the deconfinement phase transition can continue to occur, resulting in energy release (so-called deconfinement heating) in case of the first-order phase transition. The thermal evolution of neutron stars is investigated to combine phase transition and the related energy release self-consistently. We find that the appearance of deconfinement heating during spin-down result in not only the cooling delay but also the increase of surface temperature of stars. For stars characterized by intermediate and weak magnetic field strength, a period of increasing surface temperature could exist. Especially, a sharp jump in surface temperature can be produced as soon as quark matter appears in the core of stars with a weak magnetic field. We think that this may serve as evidence for the existence of deconfinement quark matter. The results show that deconfinement heating facilitates the emergence of such characteristic signature during the thermal evolution process of neutron stars.     

Neutrinos in neutron stars

It is shown that in neutron stars some of the neutrinos cannot escape from the star. For incompressible fluid sphere models, this part, in the optimum case, is 25%, and for the modelV (B) integrated byTsuruta andCameron (1966b) 1%. As a result a neutrino sea is formed and some radioactive isotopes may become stable.     

Masses and ages of delta scuti stars in eclipsing binary systems

By using data mainly from Frolovet al. (1982) for four Delta Scuti stars in eclipsing binary systems, AB Cas, Y Cam, RS Cha, and AI Hya, their physical parameters, distances, and radial pulsation modes are determined. The evolutionary track systems of Iben (1967), Paczyski (1970), and Maeder and Meynet (1988) are interpolated, in order to estimate evolutionary massesM _eand agest of these variables. Their pulsation massesM _Qare estimated from the fitting formulae of Faulkner (1977) and Fitch (1981). Our estimates of evolutionary massesM _eand pulsation massesM _Qare close to the massesM determined by Frolovet al. (1982) from the star binarity. The only exception is AB Cas, for which there is no agreement between certain star parameters.Another, independent approach is also applied to the stars RS Cha and AI Hya: by using their photometric indicesb — y andc ₁ from the catalogue of López de Cocaet al. (1990) and appropriate photometric calibrations, other sets of physical parameters, distances, modes, ages, evolutionary and pulsation masses of both variables are obtained.