Determination of a temporally and spatially resolved supernova rate from OB stars within 5 kpc

We spatially and temporally resolve the future Supernova (SN) rate in the Solar vicinity and the whole Galaxy by comparing observational parameters of massive stars with theoretical models for estimating age and mass and, hence, the remaining lifetime until the SN explosion. Our SN rate derived in time and space for the future (few Myr) should be the same as in the last few Myr by assuming a constant rate. From BVRIJHK photometry, parallax, spectral type, and luminosity class we compile a Hertzsprung‐Russell diagram (HRD) for 25027 massive stars and derive extinction, and luminosity, then mass, age, and remaining lifetime from evolutionary models. Within 600 pc our sample of SN progenitors and, hence, SN prediction, is complete, and all future SN events of our sample stars take place in 8 % of the area of the sky, whereas 90 % of the events take place in 7 % of the area of the sky. The current SN rate within 600 pc is increased by a factor of 5–6 compared with the Galactic rate. For a distance of 5 kpc our sample is incomplete, nevertheless 90 % of those SN events take place in only 12 % of the area of the projected sky. If the SN rate in the near future is the same as the recent past, there should be unknown young neutron stars concentrated in those areas. Our distribution can be used as input for constraints of gravitational waves detection and for neutron star searches. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Gamma-ray bursts from nearby neutron stars

We interpret the puzzling-ray bursts as emitted by cooling sparks from the surface of spasmodically accreting, old neutron stars. Their spiky, anisotropic radiation is oriented w.r.t. the galactic disk via interstellar accretion, whose orbital angular momentum tends to counteralign with the galactic spin; in this way, larger source numbers in directions of the galactic disk are compensated by smaller beaming probabilities, resulting in a near-isotropic arrival distribution, as observed by BATSE. The source distances range between 10 pc and 500 pc. Their radiated energies are of order 10³⁵ erg, corresponding to accreted clumps (blades) of typical mass 10¹⁵ g per burst. Magnetic surface field strengths range between 10¹⁰ and 10¹² G, somewhat weaker than those of newborn neutron stars.     

The Galactic Center: a laboratory for AGN?

Summary. This paper reviews the physical state of stars and Interstellar Matter in the Galactic Bulge (radius kpc from the dynamical center of the Galaxy), in the Nuclear Bulge (kpc) and in the Sgr A Radio and GMC Complex, i.e. the central \,pc of our Galaxy. The Galactic Bulge is devoid of cold Interstellar Matter and consists mainly of old stars, while the Nuclear Bulge accounts for of the mass of all of the Interstellar Matter in the Galaxy. A similar ratio holds for the formation rate of medium and high mass stars in Bulge and Disk. The metal abundance of the Interstellar Matter in the Galactic Bulge is found to be . The H-to-CO conversion factors to be applied to molecular gas in the Central Region are by factors 3 (Arimoto et al. 1996) to 10 (Sodroski et al. 1995) lower than in the solar vicinity. Hence, most H masses derived for the Central Region appear to be considerably overestimated. The Nuclear Bulge is pervaded by a thermal plasma (K) which is responsible for the diffuse free-free emission. Lyman continuum photon and dust IR luminosity of the Nuclear Bulge again account for of the respective total luminosities of the Galaxy. Magnetic fields in the Nuclear Bulge are strong (up to mG) as compared with the Galactic Disk (a few tens of G). The field lines are oriented parallel to the galactic plane inside giant molecular clouds and perpendicular to the plane in the intercloud medium. The compact source Sgr A* is close to or at the dynamical center of the Galaxy. Its radio spectrum with a high frequency cut-off at GHz, a low frequency turnover at GHz and a flux density dependence in between can be explained by synchrotron emission from quasi-monoenergetic relativistic electrons. Due to an extinction between Sun and Galactic Center corresponding to , an intrinsic weakness of this source in the near infrared, and a strong background emission from warm dust there are only upper limits available for the flux density of Sgr A* in the far, mid and near infrared and X-ray regime. The size of Sgr A* in the radio regime is cm, its dereddened K-band flux density is mJy, its luminosity has upper limits of (if radiation comes from an Accretion Disk) and (if black-body radiation from an object with a single temperature of K is assumed). If anyone of the soft X-ray sources detected by ROSAT actually coincides with Sgr A*, its X-ray luminosity would be less than a few . With a dark mass of Sgr A* is the best candidate for a starving black hole, although there are no observational indications for the presence of a (Standard) Accretion Disk. While the radio/IR spectrum of Sgr A* is purely nonthermal, the spectrum integrated over the central parsec resembles that of a Seyfert galaxy. Sgr A* is embedded in the Hii region Sgr A West with part of the ionized gas forming a minispiral. Sgr A West is surrounded by the Circum Nuclear Disk, an irregular shaped assembly of molecular gas which extends from pc and rotates around the Galactic Center with an estimated dynamical time scale of \,yr. The total luminosity of of the central parsec is due to the radiation of early-type stars of which have now been directly identified as luminous blue supergiants. It is still debated, however, if these stars can also account for all of the ionization of Sgr A West. In addition, the central parsec contains red giants, AGB stars, and a few super giants of which the brightest are now identified by direct imaging. These stars – together with a few million low mass main sequence stars – account for the bulk of the 2.2\,m emission. The spatial distributions of the three stellar populations in the central pc are remarkably different. Sgr A* is – along the line-of-sight – presumably located close to the center of the Hii region Sgr A West, which in turn is located in front of the extended (pc) synchrotron source Sgr A East, which appears to be the remnant of a gigantic explosion (of the order of the energy of a single supernova explosion) which took place yr ago inside the GMC Sgr A East Core. X-ray observations show within pc a pervasive hot (keV) plasma of expansion age of yr. Both phenomena – as well as the formation of the Circum Nuclear Disk – may have the same origin. Influx of material is observed within the Nuclear Bulge on all distance scales. In the Nuclear Bulge (pc) as well as in the Circum Nuclear Disk (pc) inflow towards the Galactic Center occurs primarily in the galactic plane and amounts to a few . The accretion rate into the central Black Hole, deduced from the luminosity of Sgr A*, however, appears to be lower by at least five orders of magnitude (assuming standard disk accretion). But in an equilibrium state only part of the infalling mass which is not accreted by the Black Hole can be consumed by star formation. A mass inflow rate varying with time is a more natural explanation. Comparing the physical state of the Center of our Galaxy with that of Active Galactic Nuclei derived from observations and modelling, we find that most of the basic characteristics of an AGN are also present in the Galactic Center. Lacking are, however, both the evidence for a standard Accretion Disk and a hard UV spectrum with accompanying high excitation emission lines in the Galactic Center which are characteristic for AGN. The luminosity of the central parsec, , amounts to only of the total luminosity of the Galaxy of . Seen from a distance of M31 (kpc) with an angular resolution of (corresponding to a linear size of pc) the Center of our Galaxy would appear as a mildly active nucleus with some starburst activity and would probably be classified as a weak Seyfert galaxy. The synchrotron spectrum of Sgr A*, however, would be completely masked by reprocessed stellar light (i.e. free-free and dust emission). Received: October 21, 1996     

A General Relativistic Calculation of Boundary Layer and Disc Luminosity for Accreting Non-magnetic Neutron Stars in Rapid Rotation

We calculate the disc and boundary layer luminosities for accreting rapidly rotating neutron stars with low magnetic fields in a fully general relativistic manner. Rotation increases the disc luminosity and decreases the boundary layer luminosity. A rapid rotation of the neutron star substantially modifies these quantities as compared with the static limit. For a neutron star rotating close to the centrifugal mass shed limit, the total luminosity has contribution only from the extended disc. For such maximal rotation rates, we find that well before the maximum stable gravitational mass configuration is reached, there exists a limiting central density, for which particles in the innermost stable orbit will be more tightly bound than those at the surface of the neutron star. We also calculate the angular velocity profiles of particles in Keplerian orbits around the rapidly rotating neutron star. The results are illustrated for a representative set of equation of state models of neutron star matter.     

Stellar associations and complexes in M33

About 460 OB associations were selected by a comparison of theUBV plates. TheUBV photographic photometry of 1944 blue stars in the associations was made. The new associations appear like cores within Humphreys and Sandage's associations. Their star content, size distribution, and mean size 80 pc confirm their identity with the OB associations in the Galaxy and in the Magellanic clouds. The boundaries of the associations are delineated mainly by the density of the resolved stars on theU plates. It is impossible to divide them into smaller areas. The genuine OB associations form groups of two or more members with a length scale of 250 pc. Their boundaries were delineated independently, but they coincide with the OB associations of Humphreys and Sandage (1980). These groups represent real concentration of blue massive stars with a large age dispersion. The star complexes unify a group of associations,Hii regions, andHi peak distribution. Their mean size is 570 pc. The extensiveHi clouds with a mean size of 1.2 kpc contain two or more star complexes. The questions related to star formation are briefly discussed.     

Analysis of the Z distribution of young objects in the Galactic thin disk

We have obtained new estimates of the Sun’s distance from the symmetry plane Z _⊙ and the vertical disk scale height h using currently available data on stellar OB associations, Wolf–Rayet stars, HII regions, and Cepheids. Based on individual determinations, we have calculated the mean Z _⊙ = ?16 ± 2 pc. Based on the model of a self-gravitating isothermal disk for the density distribution, we have found the following vertical disk scale heights: h = 40.2 ± 2.1 pc from OB associations, h = 47.8 ± 3.9 pc from Wolf–Rayet stars, h = 48.4 ± 2.5 pc from HII regions, and h = 66.2 ± 1.6 pc from Cepheids. We have estimated the surface, Σ = 6 kpc^?2, and volume, D(Z _⊙) = 50.6 kpc^?3, densities from a sample of OB associations. We have found that there could be ～5000 OB associations in the Galaxy.     

Population of the galactic disc in the solar neighbourhood

Spatial density parameters ofz-distribution and luminosity dispersion in the HR diagram for nine stellar spectral groups (A-K III, A-G IV, A-F V) were determined. It was shown that all stars considered can be divided into two subsets according to their scale-height . For the younger stars is of the order of 100–150 pc, and for older ones of 300–400 pc. The accuracy of the parameters derived was evaluated with use of numerical simulations. We found that for subgiants and dwarfs the relative error in the local density makes 10 to 15%, and for giants it is a little higher (20 to 25%). For other parameters the typical error is of the order of 20–30%.     

Low mass stars

Low mass stars contribute an important fraction to the mass of our Galaxy. Due to the faintness of these stars a direct investigation of their space distribution and kinematics can be carried out only in the immediate solar neighbourhood. This fact emphasizes the importance of the Third Catalogue of Nearby Stars (CNS3) as a probe of the stellar content of our galaxy.A preliminary version (Gliese and Jahreiss, 1991) of the CNS3 was recently released. Based on this version the spatial distribution of the nearby red dwarf stars is discussed. An infrared and a bolometric luminosity function is presented and compared with independent determinations from photometric surveys. An outlook is given on the expectation for the next decade due to the various surveys presently carried out or planned for the near future.     

Reconstruction of the Three-Dimensional Density Distribution of Flare Stars in the Pleiades

This paper deals with the determination of the three-dimensional distribution of flare stars in the Pleiades cluster. For this purpose a one-dimensional distribution is first constructed from the observed two-dimensional distribution of the stars. It is shown that reliable construction of one dimensioned distribution requires solution of the Abel equation. The last one used to determine the dependence of the true three-dimensional distribution of the density of flare stars on distance from the center of the cluster. This reveals a spherical layer of width 0.5 pc (2.8 R 3.5 pc) with a deficit in the number of flare stars. A profile of the three-dimensional density distribution of flare stars is constructed in the region of deficit. The characteristics of this region are described.     

Constraining the number of compact remnants near Sgr A*

Due to dynamical friction stellar mass black holes and neutron stars are expected to form high-density cusps in the inner parsec of our Galaxy. These compact remnants, expected to number around 20 000, may be accreting cold dense gas present there, and give rise to potentially observable X-ray emission. Here we build a simple but detailed time-dependent model of such emission. The possibility that these accretion flows are radiatively inefficient is taken into account and brings in some uncertainty in the conclusions. Despite this uncertainty, we find that at least several X-ray sources of this nature should be detectable with Chandra at any one time. Turning this issue around, we also ask a question of what current observational constraints might be telling us about the total number of compact remnants. In our 'best guess' model, a cusp of ∼40 000 remnants overpredicts the number of discrete sources and the total X-ray luminosity of the inner parsec, and is hence ruled out. In the most radiatively inefficient scenario that we consider, the radiative efficiency is set to be as small as  ɛ= 10⁻⁵  . In this rather unlikely scenario, a cusp of ∼40 000 black holes would be allowed by the data, but several individual sources should still be visible. Future observations of the distribution and orbits of the cold ionized gas in the inner parsec of our Galaxy will put tighter constraints on the cusp of compact remnants.     

Something about the structure of the Galaxy

We analyse a sample of 507 evolved (OH/IR) stars in the region (10°>ℓ>−45°), (| b |<3°). We derive average ages for subsets of this sample, and use those sets as beacons for the evolution of the Galaxy. In the bulge, the oldest OH/IR stars in the plane are 7.5 Gyr (1.3 M_⊙), and in the disc 2.7 Gyr (2.3 M_⊙). The vertical distribution of almost all AGB stars in the disc is found to be nearly exponential, with scaleheight increasing from 100 pc for stars ≲1 Gyr old to 500 pc for stars ≳5 Gyr old. There may be a small, disjunct population of OH/IR stars. The radial distribution of AGB stars is dictated by the metallicity gradient. Unequivocal morphological evidence is presented for the existence of a central bar, but parameters can be constrained only for a given spatial-density model. Using a variety of indicators, we identify the radii of the inner ultraharmonic (2.5 kpc) and corotation resonance (3.5 kpc). We show that the 3-kpc arm is likely to be an inner ring, as observed in other barred galaxies, by identifying a group of evolved stars that is connected to the 3-kpc H  i filament. Also, using several observed features, we argue that an inner-Lindblad resonance exists, at ∼1–1.5 kpc. The compositions of OH/IR populations within 1 kpc of the Galactic Centre give insight into the bar-driven evolution of the inner regions. We suggest that the bar is ∼8 Gyr old, relatively weak (SAB), and may be in a final stage of its existence.     

Accretion by isolated neutron stars

Accretion of interstellar material by an isolated neutron star is discussed. The point I address here is the interaction between the accretion flow and the stellar magnetosphere. I show that the interchange instabilities of the magnetospheric boundary under the conditions of interest are basically suppressed. The entry of the material into the magnetosphere is governed by diffusion. Due to this reason the persistent accretion luminosity of isolated neutron stars is limited to <4×10²⁶ erg s⁻¹. These objects can also appear as X-ray bursters with the burst durations of ∼30 min and repetition time of ∼10⁵ yr. This indicates that the number of the accreting isolated neutron stars which could be observed with recent and modern X-ray missions is a few orders of magnitude smaller than that previously estimated.      

Sequential explosions of supernovae in an ob association and formation of a superbubble

From the standpoint of view that the early type stars are formed sequentially at an OB association, it is expected that the supernova explosions will also occur sequentially. We study the expansion law of a supernova remnant, which is formed by sequential explosions of supernovae. The superbubbles and supershells with the radii 2001000 pc are naturally explained by this model. Assuming that the sequential explosion of supernovae occurs at every OB association, we deduce the star formation rate in our Galaxy.     

Relation between the X-ray and optical luminosities in binary systems with accreting nonmagnetic white dwarfs

We investigate the relation between the optical (g-band) and X-ray (0.5–10 keV) luminosities of accreting nonmagnetic white dwarfs. According to the present-day counts of the populations of star systems in our Galaxy, these systems have the highest space density among the close binary systems with white dwarfs. We show that the dependence of the optical luminosity of accreting white dwarfs on their X-ray luminosity forms a fairly narrow one-parameter curve. The typical half-width of this curve does not exceed 0.2–0.3 dex in optical and X-ray luminosities, which is essentially consistent with the amplitude of the aperiodic flux variability for these objects. At X-ray luminosities L _x ～ 10³² erg s^?1 or lower, the optical g-band luminosity of the accretion flow is shown to be related to its X-ray luminosity by a factor ～2–3. At even lower X-ray luminosities (L _x ? 10³⁰ erg s^?1), the contribution from the photosphere of the white dwarf begins to dominate in the optical spectrum of the binary system and its optical brightness does not drop below M _g ～ 13–14. Using the latter fact, we show that in current and planned X-ray sky surveys, the family of accreting nonmagnetic white dwarfs can be completely identified to the distance determined by the sensitivity of an optical sky survey in this region. For the Sloan Digital Sky Survey (SDSS) with a limiting sensitivity m _g ～ 22.5, this distance is ～400–600 pc.     

Neutron stars in globular clusters: Formation and observational manifestations

Population synthesis is used to model the number of neutron stars in globular clusters that are observed as low-mass X-ray sources and millisecond radio pulsars. The dynamical interactions between binary and single stars in a cluster are assumed to take place only with a continuously replenished “background” of single stars whose properties keep track of the variations in parameters of the cluster as a whole and the evolution of single stars. We use the hypothesis that the neutron stars forming in binary systems from components with initial masses of ～8–12 M _⊙ during the collapse of degenerate O-Ne-Mg cores through electron captures do not acquire a high space velocity. The remaining neutron stars (from single stars with masses >8 M _⊙ or from binary components with masses >12 M _⊙) are assumed to be born with high space velocities. According to this hypothesis, a sizeable fraction of the forming neutron stars remain in globular clusters (about 1000 stars in a cluster with a mass of 5 × 10⁵ M _⊙). The number of millisecond radio pulsars forming in such a cluster in the case of accretion-driven spinup in binary systems is found to be ～10, in agreement with observations. Our modeling also reproduces the observed shape of the X-ray luminosity function for accreting neutron stars in binary systems with normal and degenerate components and the distribution of spin periods for millisecond pulsars.     

The spiral structure in the inner parts of the Galaxy

The spiral structure of the inner parts of the Galaxy is studied using 21 cm line data and stellar data. To study the neutral hydrogen distribution in the galactic layer a parameter =(dV/dr) proportional to the mean densities is calculated using a first approximation for the velocity gradients due to differential rotation.The obtained distribution (R, Z) shows spiral features completely consistent with the early star distribution and with the Hii regions. The corrugation effect of the galactic layer is observed in all the studied zones in neutral hydrogen and in the distribution of the OB stars in the Carina zone.The pattern obtained indicates four spiral arms for the inner parts of the Galaxy, three of which are identified also in the stellar data (arms -I, -II, and -III) and the more distant -IV in Hii regions.The local arm according to the stellar data of Kilkennyet al. forms a feature completely similar to the arms -I and -II and there are no indications that this arm is a special material branch between two main spiral arms as has been supposed in order to conciliate the neutral hydrogen pattern with the stellar distribution.The pitch angles for the spiral arms are approximately 13°–17°.The observed wave form distribution of the hydrogen cloud layer is completely consistent with the theoretical predictions of Nelson (1976) but there are no indications of such an effect in the intercloud hydrogen. The corrugated cloud layer has a width of 100 pc, a wave amplitude of 70 pc, and a wavelength which grows with the galactic center distance (approx. 2 kpc in the zones next to the galactic nucleus and 2.6–3.0 kpc in the zones next to the Sun). To each wavelength correspond two spiral arms. The spiral features in our Galaxy show characteristics quite similar to the features in the Andromeda nebula, not only in the component materials (neutral hydrogen, Hii regions and possibly also dust and stars) but also in their kinematics.     

Masses and luminosities of O‐ and B‐type stars and red supergiants

Massive stars are of interest as progenitors of supernovae, i.e. neutron stars and black holes, which can be sources of gravitational waves. Recent population synthesis models can predict neutron star and gravitational wave observations but deal with a fixed supernova rate or an assumed initial mass function for the population of massive stars. Here we investigate those massive stars, which are supernova progenitors, i.e. with O‐ and early B‐type stars, and also all supergiants within 3 kpc. We restrict our sample to those massive stars detected both in 2MASS and observed by Hipparcos, i.e. only those stars with parallax and precise photometry. To determine the luminosities we calculated the extinctions from published multi‐colour photometry, spectral types, luminosity class, all corrected for multiplicity and recently revised Hipparcos distances. We use luminosities and temperatures to estimate the masses and ages of these stars using different models from different authors. Having estimated the luminosities of all our stars within 3 kpc, in particular for all O‐ and early B‐type stars, we have determined the median and mean luminosities for all spectral types for luminosity classes I, III, and V. Our luminosity values for supergiants deviate from earlier results: Previous work generally overestimates distances and luminosities compared to our data, this is likely due to Hipparcos parallaxes (generally more accurate and larger than previous ground‐based data) and the fact that many massive stars have recently been resolved into multiples of lower masses and luminosities. From luminosities and effective temperatures we derived masses and ages using mass tracks and isochrones from different authors. From masses and ages we estimated lifetimes and derived a lower limit for the supernova rate of ≈20 events/Myr averaged over the next 10 Myr within 600 pc from the sun. These data are then used to search for areas in the sky with higher likelihood for a supernova or gravitational wave event (like OB associations) (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Search for relativistic companions in ‘non-X-ray’ binary systems

Consideration is given to a search for relativistic objects in massive close binary systems without strong X-ray emission (L _x <10³⁴ erg s^–1). It is pointed out that, according to the present-day theory on the evolution of massive close binaries, the number of neutron stars and black holes in non-X-ray binary systems must be 100 times the number of the known X-ray binaries comprising OB supergiant stars; that is why, in studying non-X-ray binary systems, the chances are to detect about a hundred of black holes in the Galaxy.Criteria are formulated for the relativistic nature of companions in the binary systems, such as high spatial velocity values and height Z over the galactic plane for OB stars (runaway stars) and for Wolf-Rayet stars. As reported by Tutukov and Yungelson (1973), as well as by van den Heuvel (1976), the presence of ring-type nebulae can serve as another indication of a relativistic nature of companions in the case of Wolf-Rayet stars.Data are collected on Wolf-Rayet stars with low-mass companions (Table I), which can be relativistic objects accreting within a strong stellar wind from Wolf-Rayet stars. Presented are new findings in respect of spectral examination of the runaway OB-stars (Table II), bringing together data on eight OB stars which can represent binary systems with relativistic companions (Table III).A list of 28 OB-stars (Table IV) which offer a good chance for finding relativistic companions is given.     

Symmetry of the neutron and proton superfluidity effects in cooling neutron stars

We investigate the combined effect of neutron and proton superfluidities on the cooling of neutron stars whose cores consist of nucleons and electrons. We consider the singlet state paring of protons and the triplet pairing of neutrons in the cores of neutron stars. The critical superfluid temperatures T _c are assumed to depend on the matter density. We study two types of neutron pairing with different components of the total angular momentum of a Cooper pair along the quantization axis (|m _J |=0 or 2). Our calculations are compared with the observations of thermal emission from isolated neutron stars. We show that the observations can be interpreted by using two classes of superfluidity models: (1) strong proton superfluidity with a maximum critical temperature in the stellar core T _c ^max ?4×10⁹ K and weak neutron superfluidity of any type (T _c ^max ?2×10⁸ K); (2) strong neutron superfluidity (pairing with m _J =0) and weak proton superfluidity. The two types of models reflect an approximate symmetry with respect to an interchange of the critical neutron and proton pairing temperatures.     

Small-scale Gradients and Large-scale Diffusion of Charged Particles in the Heliospheric Magnetic Field

We present a simple spin-evolution model that predicts that rapidly rotating accreting neutron stars will be confined mainly to a narrow range of spin frequencies: P=1.5-5 ms. This is in agreement with current observations of neutron stars in both the low-mass X-ray binaries and the millisecond radio pulsars. The main ingredients in the model are (1) the instability of r-modes above a critical spin rate, (2) the thermal runaway that is due to the heat released as viscous damping mechanisms counteract the r-mode growth, and (3) a revised estimate of the strength of the dissipation that is due to the presence of a viscous boundary layer at the base of the crust in an old and relatively cold neutron star. We discuss the gravitational waves that are radiated during the brief r-mode-driven spin-down phase. We also briefly touch on how the new estimates affect the predicted initial spin periods of hot young neutron stars.