Some static relativistic compact charged fluid spheres in general relativity

In this work a new family of relativistic models of electrically charged compact star has been obtained by solving Einstein–Maxwell field equations with preferred form of one of the metric potentials and a suitable form of electric charge distribution function. The resulting equation of state (EOS) has been calculated. The relativistic stellar structure for matter distribution obtained in this work may reasonably models an electrically charged compact star whose energy density associated with the electric fields is on the same order of magnitude as the energy density of fluid matter itself (e.g. electrically charged bare strange stars). Based on the analytic model developed in the present work, the values of the relevant physical quantities have been calculated by assuming the estimated masses and radii of some well known strange star candidates like X-ray pulsar Her X-1, millisecond X-ray pulsar SAX J 1808.4-3658, and 4U 1820-30.     

Modeling the luminosity function of galactic low-mass X-ray binaries

The evolution of the family of binaries with a low-mass star and a compact neutron star companion (low-mass X-ray binaries (LMXBs) with neutron stars) ismodeled by the method of population synthesis. Continuous Roche-lobe filling by the optical star in LMXBs is assumed to be maintained by the removal of orbital angular momentum from the binary by a magnetic stellar wind from the optical star and the radiation of gravitational waves by the binary. The developed model of LMXB evolution has the following significant distinctions: (1) allowance for the effect of the rotational evolution of a magnetized compact remnant on themass transfer scenario in the binary, (2) amore accurate allowance for the response of the donor star to mass loss at the Roche-lobe filling stage. The results of theoretical calculations are shown to be in good agreement with the observed orbital period-X-ray luminosity diagrams for persistent Galactic LMXBs and their X-ray luminosity function. This suggests that the main elements of binary evolution, on the whole, are correctly reflected in the developed code. It is shown that most of the Galactic bulge LMXBs at luminosities L _x > 10³⁷ erg s^?1 should have a post-main-sequence Roche-lobe-filling secondary component (low-mass giants). Almost all of the models considered predict a deficit of LMXBs at X-ray luminosities near ～10^36.5 erg s^?1 due to the transition of the binary from the regime of angular momentum removal by a magnetic stellar wind to the regime of gravitational waves (analogous to the widely known period gap in cataclysmic variables, accreting white dwarfs). At low luminosities, the shape of the model luminosity function for LMXBs is affected significantly by their transient behavior-the accretion rate onto the compact companion is not always equal to the mass transfer rate due to instabilities in the accretion disk around the compact object. The best agreement with observed binaries is achieved in the models suggesting that heavy neutron stars with masses 1.4–1.9M _⊙ can be born.     

X-ray Observations of Eight Young Open Star Clusters: I. Membership and X-ray Luminosity

We present a detailed investigation of X-ray source contents of eight young open clusters with ages between 4 to 46 Myr using archival X-ray data from XMM-Newton. The probable cluster memberships of the X-ray sources have been established on the basis of multi-wavelength archival data, and samples of 152 pre-main sequence (PMS) low mass (<2M _⊙), 36 intermediate mass (2–10M _⊙) and 16 massive (>10M _⊙) stars have been generated. X-ray spectral analyses of high mass stars reveal the presence of high temperature plasma with temperature <2 keV, and mean L _X/L _bol of 10^???6.9. In the case of PMS low mass stars, the plasma temperatures have been found to be in the range of 0.2 keV to 3 keV with a median value of ~1.3 keV, with no significant difference in plasma temperatures during their evolution from 4 to 46 Myr. The X-ray luminosity distributions of the PMS low mass stars have been found to be similar in the young star clusters under study. This may suggest a nearly uniform X-ray activity in the PMS low mass stars of ages ~4–14 Myr. These observed values of L _X/L _bol are found to have a mean value of 10^??3.6±0.4, which is below the X-ray saturation level. The L _X/L _bol values for the PMS low mass stars are well correlated with their bolometric luminosities, that implies its dependence on the internal structure of the low mass stars. The difference between the X-ray luminosity distributions of the intermediate mass stars and the PMS low mass stars has not been found to be statistically significant. Their L _X/L _bol values, however have been found to be significantly different from each other with a confidence level greater than 99.999% and the strength of X-ray activity in the intermediate mass stars is found to be lower compared to the low mass stars. However, the possibility of X-ray emission from the intermediate mass stars due to a low mass star in close proximity of the intermediate mass star can not be ruled out.     

Late stages of the evolution of close binaries

Close binaries can evolve through various ways of interaction into compact objects (white dwarfs, neutron stars, black holes). Massive binary systems (mass of the primaryM ₁ larger than 14 to 15M ₀) are expected to leave, after the first stage of mass transfer a compact component orbiting a massive star. These systems evolve during subsequent stages into massive X-ray binaries. Systems with initial large periode evolve into Be X-ray binaries.Low mass X-ray sources are probably descendants of lower mass stars, and various channels for their production are indicated. The evolution of massive close binaries is examined in detail and different X-ray stages are discussed. It is argued that a first X-ray stage is followed by a reverse extensive mass transfer, leading to systems like SS 433, Cir X1. During further evolution these systems would become Wolf-Rayet runaways. Due to spiral in these system would then further evolve into ultra short X-ray binaries like Cyg X-3.Finally the explosion of the secondary will in most cases disrupt the system. In an exceptional case the system remains bound, leading to binary pulsars like PSR 1913+16. In such systems the orbit will shrink due to gravitational radiation and finally the two neutron stars will coalesce. It is argued that the millisecond pulsar PSR 1937+214 could be formed in this way.A complete scheme starting from two massive ZAMS stars, ending with a millisecond pulsar is presented.Paper presented at the Lembang-Bamberg IAU Colloquium No. 80 on Double Stars: Physical Properties and Generic Relations, held at Bandung, Indonesia 3–7 June, 1983.     

Terrestrial planet formation surrounding close binary stars

Most stars reside in binary/multiple star systems; however, previous models of planet formation have studied growth of bodies orbiting an isolated single star. Disk material has been observed around both components of some young close binary star systems. Additionally, it has been shown that if planets form at the right places within such disks, they can remain dynamically stable for very long times. Herein, we numerically simulate the late stages of terrestrial planet growth in circumbinary disks around ‘close’ binary star systems with stellar separations 0.05 AU?aB?0.4 AU and binary eccentricities 0?eB?0.8. In each simulation, the sum of the masses of the two stars is 1 M_⊙, and giant planets are included. The initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet formation within our Solar System by Chambers [Chambers, J.E., 2001. Icarus 152, 205-224], and around each individual component of the α Centauri AB binary star system by Quintana et al. [Quintana, E.V., Lissauer, J.J., Chambers, J.E., Duncan, M.J., 2002. Astrophys. J. 576, 982-996]. Multiple simulations are performed for each binary star system under study, and our results are statistically compared to a set of planet formation simulations in the Sun-Jupiter-Saturn system that begin with essentially the same initial disk of protoplanets. The planetary systems formed around binaries with apastron distances QB≡aB(1+eB)?0.2 AU are very similar to those around single stars, whereas those with larger maximum separations tend to be sparcer, with fewer planets, especially interior to 1 AU. We also provide formulae that can be used to scale results of planetary accretion simulations to various systems with different total stellar mass, disk sizes, and planetesimal masses and densities.     

Millisecond radio pulsars with known masses: Parameter values and equation of state models

The recent fast growth of a population of millisecond pulsars with precisely measured mass provides an excellent opportunity to characterize these compact stars at an unprecedented level. This is because the stellar parameter values can be accurately computed for known mass and spin rate and an assumed equation of state (EoS) model. For each of the 16 such pulsars and for a set of EoS models from nucleonic, hyperonic, strange quark matter and hybrid classes, we numerically compute fast spinning stable stellar parameter values considering the full effect of general relativity. This first detailed catalogue of the computed parameter values of observed millisecond pulsars provides a testbed to probe the physics of compact stars, including their formation, evolution and EoS. We estimate uncertainties on these computed values from the uncertainty of the measured mass, which could be useful to quantitatively constrain EoS models. We note that the largest value of the central density ρ_c in our catalogue is ∼5.8 times the nuclear saturation density ρ_sat, which is much less than the expected maximum value 13ρ_sat. We argue that the ρ_c-values of at most a small fraction of compact stars could be much larger than 5.8ρ_sat. Besides, we find that the constraints on EoS models from accurate radius measurements could be significantly biased for some of our pulsars, if stellar spinning configurations are not used to compute the theoretical radius values.     

A Statistical Study of the Mass Distribution of Neutron Stars

By reviewing the methods of mass measurements of neutron stars in four different kinds of systems, i.e., the high-mass X-ray binaries (HMXBs), low-mass X-ray binaries (LMXBs), double neutron star systems (DNSs) and neutron star-white dwarf (NS-WD) binary systems, we have collected the orbital parameters of 40 systems. By using the boot-strap method and the Monte-Carlo method, we have rebuilt the likelihood probability curves of the measured masses of 46 neutron stars. The statistical analysis of the simulation results shows that the masses of neutron stars in the X-ray neutron star systems and those in the radio pulsar systems exhibit different distributions. Besides, the Bayes statistics of these four different kind systems yields the most-probable probability density distributions of these four kind systems to be (1.340 ± 0.230)M8, (1, 505 ± 0.125)M8,(1.335 ± 0.055)M8 and (1.495 ± 0.225)M8, respectively. It is noteworthy that the masses of neutron stars in the HMXB and DNS systems are smaller than those in the other two kind systems by approximately 0.16M8. This result is consistent with the theoretical model of the pulsar to be accelerated to the millisecond order of magnitude via accretion of approximately 0.2M8. If the HMXBs and LMXBs are respectively taken to be the precursors of the BNS and NS-WD systems, then the influence of the accretion effect on the masses of neutron stars in the HMXB systems should be exceedingly small. Their mass distributions should be very close to the initial one during the formation of neutron stars. As for the LMXB and NS-WD systems, they should have already under- gone the process of suffcient accretion, hence there arises rather large deviation from the initial mass distribution.     

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.     

Star formation rate in starburst galaxies

The statistical investigation of integrated characteristics of the blue compact galaxies from the Second Byurakan survey has been carried out. The tight correlation between the luminosities in the far-infrared range and radio range (radio-continuum at 21 cm) has been revealed. It implies the common origin of this radiation in star forming regions. The current star formation rates have been derived both from the luminosity in the far-infrared range (SFR_FIR) and from the luminosity in the radio range (SFR_{21 cm}). The Salpeter shape of the initial mass function (IMF) and low-mass cut-off of 0.1 M_⊙ in IMF has been assumed for the SFR_FIR calculation. We derive the SFR_{21 cm} assuming that radio free–free emission at λ=21 cm is thermal emission of gas ionised by massive stars. The SFR_{21 cm} is shown to be about 3.5 times as large as SFR_FIR for our sample of galaxies. These star formation rates are proposed to be considered as lower and upper limits of true galaxy star formation rates.     

Bimodality of anomalous pulsars?

At present, it is widely believed that anomalous X-ray pulsars (AXPs), soft gamma-ray repeaters (SGRs), rotational radio transients (RRATs), compact central objects (CCOs) in supernova remnants, and X-ray dim isolated neutron stars (XDINSs) belong to different classes of anomalous objects in which the central bodies are isolated neutron stars. Previously, we have shown that AXPs and SGRs can be described in terms of the drift model for parameters of the central neutron star typical of radio pulsars (rotation periods P ～ 0.1–1 s and surface magnetic fields B ～ 10¹¹–10¹³ G). Here, we show that some of the peculiarities of the sources under consideration can be explained by their geometry (in particular, by the angle β between the rotation axis and the magnetic moment). If β ? 10° (an aligned rotator), the drift waves in the outer layers of the neutron star magnetosphere can account for the observed periodicity in the radiation. For large β (a nearly orthogonal rotator), the observed modulation of the radiation and its short bursts can be explained by mass accretion from the ambient medium (e.g., a relic disk).     

A class of exact isotropic solutions of Einstein’s equations and relativistic stellar models in general relativity

In this paper we have studied a particular class of exact solutions of Einstein’s gravitational field equations for spherically symmetric and static perfect fluid distribution in isotropic coordinates. The Schwarzschild compactness parameter, GM/c ² R, can attain the maximum value 0.1956 up to which the solution satisfies the elementary tests of physical relevance. The solution also found to have monotonic decreasing adiabatic sound speed from the centre to the boundary of the fluid sphere. A wide range of fluid spheres of different mass and radius for a given compactness is possible. The maximum mass of the fluid distribution is calculated by using stellar surface density as parameter. The values of different physical variables obtained for some potential strange star candidates like Her X-1, 4U 1538–52, LMC X-4, SAX J1808.4?3658 given by our analytical model demonstrate the astrophysical significance of our class of relativistic stellar models in the study of internal structure of compact star such as self-bound strange quark star.     

Red giants in the blue compact galaxy IZw18

Archival Hubble Space Telescope (HST) data have been used for the photometry of stars in blue compact dwarf (BCD) galaxies of the IZw18 system. Applying the spatial selection of stars, we have detected red giants, stars older than 1–2 Gyr, in the galaxies. These red giants have allowed the distance to IZw18 to be reliably determined for the first time: D = 13.9 ± 1.2 Mpc. The presence of old stars in the galaxies of the IZw18 system refutes the hypothesis about the observed primary star formation in these galaxies.     

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

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

Radius of the neutron star magnetosphere during disk accretion

The dependence of the spin frequency derivative \(\dot \nu \) of accreting neutron stars with a strongmagnetic field (X-ray pulsars) on the mass accretion rate (bolometric luminosity, L_bol) has been investigated for eight transient pulsars in binary systems with Be stars. Using data from the Fermi/GBM and Swift/BAT telescopes, we have shown that for seven of the eight systems the dependence \(\dot \nu \) (L_bol) can be fitted by the model of angular momentum transfer through an accretion disk, which predicts the relation \(\dot \nu \) ～ L^6/7_bol. Hysteresis in the dependence \(\dot \nu \) (L_bol) has been confirmed in the system V 0332+53 and has been detected for the first time in the systems KS 1947+300, GRO J1008-57, and 1A 0535+26. Estimates for the radius of the neutron star magnetosphere in all of the investigated systems have been obtained. We show that this quantity varies from pulsar to pulsar and depends strongly on the analytical model and the estimates for the neutron star and binary system parameters.     

Evolution of massive close binaries

The further evolution of a massive X-ray binary consisting of a compact object and an OB supergiant is outlined. The supergiant exceeds its critical Roche lobe and a second stage of mass transfer starts. The remnant of the mass losing star — a pure helium star — develops a collapsing iron core and finally undergoes a supernova explosion. If the compact companion is a black hole the system remains bound; if the compact companion is a neutron star the system is disrupted unless an extra kick allowing an asymmetric explosion is given. Computations were performed for the massive binary 22.5M _⊙+2M _⊙. The possible final evolutionary products are: (1) a black hole and a compact object, in a binary system, (2) two run-away pulsars, (3) a binary pulsar. As final parameters for the described system the eccentricity and period for the recently discovered binary pulsar 1913+16 may be found. An orbital inclination ofi=40° may be derived. The probability for the generation of binary pulsars is very low; in most cases the system is disrupted during the supernova explosion.     

The correlations and anticorrelations in QPO data

Double peak kHz QPO frequencies in neutron star sources varies in time by a factor of hundreds Hz while in microquasar sources the frequencies are fixed and located at the line ν ₂ = 1.5ν ₁ in the frequency‐frequency plot. The crucial question in the theory of twin HFQPOs is whether or not those observed in neutron‐star systems are essentially different from those observed in black holes. In black hole systems the twin HFQPOs are known to be in a 3:2 ratio for each source. At first sight, this seems not to be the case for neutron stars. For each individual neutron star, the upper and lower kHz QPO frequencies, ν ₂ and ν ₁, are linearly correlated, ν ₂ = Aν ₁ + B , with the slope A < 1.5, i.e., the frequencies definitely are not in a 1.5 ratio. In this contribution we show that when considered jointly on a frequency‐frequency plot, the data for the twin kHz QPO frequencies in several (as opposed to one) neutron stars uniquely pick out a certain preferred frequency ratio that is equal to 1.5 for the six sources examined so far. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Effects of external magnetic field on the outflows of an accretion disc

The aim of this work is to study the effects of an external magnetic field generated by a magnetized compact star on the outflows of its accretion disc. For this purpose, we solve a set of magneto-hydrodynamic (MHD) equations for an accretion disc in spherical coordinates to consider the disc structure along the θ-direction. We also consider the magnetic field of a compact star beyond its surface as a dipolar field, producing a toroidal magnetic field inside the disc. We convert the equations to a set of ordinary differential equations (ODEs) as a function of the θ only by applying self-similar assumptions in the radial direction. Then, this set of equations is solved under symmetrical boundary conditions in the equatorial plane to obtain the velocity field. The results are considered in the gas-pressure-dominated (GPD) region and radiation-pressure-dominated (RPD) region as well. The dipolar field of the compact stars can significantly enhance the speed of outflows. It also can change the structure of the disc. The results of this work would be useful in the study of X-ray binaries, the origin of ultra-relativistic outflows, and jet formation around the compact stars.     

The location and X-ray emission of shocks around isolated be stars

All hot stars are observed to have X-ray emission: O stars haveL _X /L _bol 10^–7, whilst B stars' emission drops off with spectral subtype. Dynamical instability of OB star radiatively driven winds generates shocked regions which may be responsible for the bulk of the X-rays observed. The wind-compressed disc model of Bjorkman & Cassinelli (1993) presents another site for X-ray emission. The disc formed in the equatorial plane of a fast rotating Be star from equatorward drift of wind streamlines is confined on both sides by a shock which may also generate X-rays. As the X-ray emission originating from the wind shocking is ubiquitous amongst B and Be stars then the wind-compressed disc model näively predicts that Be stars should generate more X-rays than B stars of equivalent spectral subtype.The X-ray emission from the shocks confining compression discs has been calculated and compared to a limited set of observations. The excess X-ray emission from the Be star disc shocks is found to be undetectable over the inherent wind shocking emission.     

The HYPER-MUCHFUSS project—the constant high-velocity population

The HYPER-MUCHFUSS (HYPER-velocity stars or Massive Unseen Companions of Hot Faint Underluminous Stars from SDSS) project targets a population of high-velocity subluminous B stars to discover either close binaries with massive unseen companions or hyper-velocity stars. We re-observed high-velocity subdwarf selected candidates from the SDSS spectroscopic Data Release 6. Starting in 2007 we used several instruments and have now reached a completion level of 33% (from 265 targets), whereas we found at least 16 close binaries. Here we present results for two of our 39 hyper-velocity star candidates. From the available Digitized Sky Surveys photographic plates we measured a significant proper motion for 14 stars. Combining this information with a detailed spectroscopic analysis allows for the first time a complete determination of the 3D-trajectories for a high-velocity sample. We present our preliminary results for the two subdwarfs J1644+4523 and J1211+1437. Assuming the Standard Allen and Santillan (Rev. Mex. Astron. Astrofis. 22:255, 1991) potential the first one is bound and originates in the central region of the Galaxy. The subdwarf B star J1211+1437 is possibly unbound and seems to originate in the Galactic rim. We also performed numerical kinematical experiments with increased dark matter halo mass. and found that the origin of J1644+4523 in the central region is not changed but the time-of-flight is drastically shortened. J1211+1437 would be bound and probably belongs to population II.     

Stellar adiabatic mass loss model and applications

Roche-lobe overflow and common envelope evolution are very important in binary evolution, which is believed to be the main evolutionary channel to hot subdwarf stars. The details of these processes are difficult to model, but adiabatic expansion provides an excellent approximation to the structure of a donor star undergoing dynamical time scale mass transfer. We can use this model to study the responses of stars of various masses and evolutionary stages as potential donor stars, with the urgent goal of obtaining more accurate stability criteria for dynamical mass transfer in binary population synthesis studies. As examples, we describe here several models with the initial masses equal to 1 M _⊙ and 10 M _⊙, and identify potential limitations to the use of our results for giant-branch stars.