Radial pulsations of helium stars with masses from 10 to 50<Emphasis Type="Italic">M</Emphasis><Subscript>⊙</Subscript>

We have performed hydrodynamic calculations of the radial pulsations of helium stars with masses 10M_⊙ ≤ M ≤ 50M_⊙, luminosity-to-mass ratios 5 × 10³L_⊙/M_⊙ ≤ L/M ≤ 2.5 × 10⁴L_⊙/M_⊙, and effective temperatures 2 × 10⁴ K ≤ T_eff ≤ 10⁵ K for helium and heavy-element mass fractions of Y=0.98 and Z=0.02, respectively. We show that the high-temperature boundary of the instability region for radial pulsations at L/M ? 10⁴L_⊙/M_⊙ extends to T_eff≈10⁵ K. The amplitude of the velocity variations for outer layers is several hundred km s^?1, while the brightness variations in the B band of the UBV photometric system are within the range from several hundredths to half a magnitude. At constant luminosity-to-mass ratio, the radial pulsation period is determined only by the effective temperature of the star. In the ranges of luminosity-to-mass ratios 10⁴L_⊙/M_⊙ ≤ L/M ≤ 2 × 10⁴L_⊙/M_⊙ and effective temperatures 5 × 10⁴ K ≤ T_eff ≤ 9 × 10⁴ K, the periods of the radial modes are within 6 min ?Π?10³ min.     

Radial pulsations of helium stars with masses from 1 to 10 <Emphasis Type="Italic">M</Emphasis><Subscript>⊙</Subscript>

We present the results of our hydrodynamic calculations of radial pulsations in helium stars with masses 1 M_⊙ ≤ M ≤ 10 M_⊙, luminosity-to-mass ratios 1 × 10³L_⊙/M_⊙ ≤ L/M ≤ 2 × 10⁴L_⊙/M_⊙, and effective temperatures 2 × 10⁴ K ≤ T_eff ≤ 10⁵ K for mass fractions of helium Y=0.98 and heavy elements Z=0.02. We show that the lower boundary of the pulsation-instability region corresponds to L/M ～ 10³L_⊙/M_⊙ and that the instability region for L/M ? 5 × 10³L_⊙/M_⊙ is bounded by effective temperatures T_eff ? 3 × 10⁴ K. As the luminosity rises, the instability boundary moves into the left part of the Hertzsprung-Russell diagram and radial pulsations can arise in stars with effective temperatures T_eff ? 10⁵ K at L/M ? 7 × 10³L_⊙/M_⊙. The velocity amplitude for the outer boundary of the hydrodynamic model increases with L/M and lies within the range 200 ? ΔU ? 700 km s^?1 for the models under consideration. The periodic shock waves that accompany radial pulsations cause a significant change of the gas-density distribution in the stellar atmosphere, which is described by a dynamic scale height comparable to the stellar radius. The dynamic instability boundary that corresponds to the separation of the outer stellar atmospheric layers at a superparabolic velocity is roughly determined by a luminosity-to-mass ratio L/M ～ 3 × 10⁴L_⊙/M_⊙.     

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.     

Hydrodynamical models of type II supernovae

A series of hydrodynamical models of type-II supernova outbursts (SNII) has been calculated. Approximate relations connecting the total outburst energy ε, the mass of envelope ejectedM, the presupernova radiusR, and the amount of ionizing quanta radiated by the supernovaeN _H with such values as the duration of the light curve plateau Δt, and absolute magnitude in the wavelength bandV and photospheric velocityU _PH observed near the middle of the plateau have been established. Advantage has been taken of the relations to obtain a preliminary evaluation for the characteristics of the average SN II: ε=7×10⁵⁰ erg,M=6M _⊙,R=500R _⊙,N _H=2×10⁵⁸. The SNIIs with plateau-like light curves seem to be accounted for by thermonuclear explosions of degenerate cores of red giant stars and result in a total disruption of the star without any stellar remnant. To the contrary, SNIIs with linear light curves have substantially different properties (in particular, they throw considerably less massive envelopes off). These SNII must signify the birth of collapsed objects—neutron stars (pulsars) or black holes.     

Pulsations of intermediate–mass stars on the asymptotic giant branch

Evolutionary tracks from the zero age main sequence to the asymptotic giant branch were computed for stars with initial masses 2 M _⊙ ≤ M _ZAMS ≤ 5 M _⊙ and metallicity Z = 0.02. Some models of evolutionary sequences were used as initial conditions for equations of radiation hydrodynamics and turbulent convection describing radial stellar pulsations. The early asymptotic giant branch stars are shown to pulsate in the fundamental mode with periods 30 day ? Π ? 400day. The rate of period change gradually increases as the star evolves but is too small to be detected (Π?/Π < 10^?5 yr^?1). Pulsation properties of thermally pulsing AGB stars are investigated on time intervals comprising 17 thermal pulses for evolutionary sequences with initial masses M _ZAMS = 2 M _⊙ and 3 M _⊙ and 6 thermal pulses for M _ZAMS = 4 M _⊙ and 5 M _⊙. Stars with initial masses M _ZAMS ≤ 3 M _⊙ pulsate either in the fundamental mode or in the first overtone, whereas more massive red giants (M _ZAMS ≥ 4 M _⊙) pulsate in the fundamental mode with periods Π ? 10³ day. Most rapid pulsation period change with rate ?0.02 yr^?1 ? Π?/Π ? ?0.01 yr^?1 occurs during decrease of the surface luminosity after the maximum of the luminosity in the helium shell source. The rate of subsequent increase of the period is Π?/Π ? 5 × 10^?3 yr^?1.     

Evolution of low-mass helium stars in semidetached binaries

We systematically investigate the evolution of low-mass (0.35, 0.40, and 0.65M _⊙) helium donors in semidetached binaries with white-dwarf accretors. The initial periods of the binaries are chosen in such a way that the helium abundance in the center of the models at the time of Roche lobe overflow varies between Y _c = 0.98 and Y _c ? 0.1. The results of our calculations can be used to analyze the formation scenarios and evolutionary status of AM CVn stars. We show that the minimum orbital periods of the semidetached binaries depend weakly on the total mass of the components and the evolutionary phase of the donor at the time of Roche lobe overflow and are 9–10 min. The differences in the mass transfer rates after P _orb reaches its minimum in the range P _orb ≈ 10–40 min do not exceed a factor of ～2.5. For P _orb ? 20 min, the mass-losing stars are weakly degenerate homogeneous cooling objects; the He, C, N, O, and Ne abundances depend on the evolutionary phase at which Roche lobe overflow occurred. For the binaries that are currently believed to be the most probable candidates for AM CVn stars with helium donors, Y ? 0.4, X _C ? 0.3, X _O ? 0.25, and X _N ? 0.5 × 10^?2. In the binaries under consideration, once P _orb ≈ 40 min has been reached, the mass loss time scale begins to exceed the thermal time scale of the donors, the latter begin to contract, their matter becomes degenerate, and the populations of AMCVn stars with white-dwarf and helium-star progenitors of their donors probably merge together.     

Radial pulsations of red giant branch stars

We performed hydrodynamic computations of nonlinear stellar pulsations of population I stars at the evolutionary stages of the ascending red giant branch and the following luminosity drop due to the core helium flash. Red giants populating this region of the Hertzsprung–Russel diagram were found to be the fundamental mode pulsators. The pulsation period is the largest at the tip of the red giant branch and for stars with initial masses from 1.1 M _⊙ to 1.9 M _⊙ ranges from ∏ ≈ 254 day to ∏ ≈ 33 day , respectively. The rate of period change during the core helium flash is comparable with rates of secular period change in Mira type variables during the thermal pulse in the helium shell source. The period change rate is largest (∏?/∏ ≈ ?10^?2 yr^?1) in stars with initial mass M ^ZAMS = 1.1 M ^⊙ and decreases to ∏?/∏ ～ ?10^?3 yr^?1 for stars of the evolutionary sequence M ^ZAMS = 1.9 M ^⊙. Theoretical light curves of red giants pulsating with periods ∏ > 200 day show the presence of the secondary maximum similar to that observed in many Miras.     

Refinement of the parameters of three selected model Galactic potentials based on the velocities of objects at distances up to 200 kpc

This paper is a continuation of our recent paper devoted to refining the parameters of threecomponent (bulge, disk, halo) axisymmetric model Galactic gravitational potentials differing by the expression for the dark matter halo using the velocities of distant objects. In all models the bulge and disk potentials are described by the Miyamoto–Nagai expressions. In our previous paper we used the Allen–Santillán (I), Wilkinson–Evans (II), and Navarro–Frenk–White (III) models to describe the halo. In this paper we use a spherical logarithmic Binney potential (model IV), a Plummer sphere (model V), and a Hernquist potential (model VI) to describe the halo. A set of present-day observational data in the range of Galactocentric distances R from 0 to 200 kpc is used to refine the parameters of the listed models, which are employed most commonly at present. The model rotation curves are fitted to the observed velocities by taking into account the constraints on the local matter density ρ_⊙= 0.1 M _⊙pc^?3 and the force K _z=1.1/2πG = 77M _⊙pc^?2 acting perpendicularly to the Galactic plane. The Galactic mass within spheres of radius 50 and 200 kpc are shown to be, respectively, M ₅₀ = (0.409 ± 0.020) × 10¹² M _⊙ and M ₂₀₀ = (1.395 ± 0.082) × 10¹² M _⊙ in model IV, M ₅₀ = (0.417 ± 0.034) × 10¹² M _⊙ and M ₂₀₀ = (0.469 ± 0.038) × 10¹² M _⊙in model V, and M ₅₀ = (0.417 ± 0.032) × 10¹² M _⊙ and M ₂₀₀ = (0.641 ± 0.049)× 10¹² M _⊙ in model VI. Model VI looks best among the three models considered here from the viewpoint of the achieved accuracy of fitting the model rotation curves to the measurements. This model is close to the Navarro–Frenk–White model III refined and considered best in our previous paper, which is shown using the integration of the orbits of two globular clusters, Lynga 7 and NGC 5053, as an example.     

Rotation curve and mass distribution in the Galaxy from the velocities of objects at distances up to 200 kpc

Three three-component (bulge, disk, halo) model Galactic gravitational potentials differing by the expression for the dark matter halo are considered. The central (bulge) and disk components are described by the Miyamoto–Nagai expressions. The Allen–Santillán (I), Wilkinson–Evans (II), and Navarro–Frenk–White (III) models are used to describe the halo. A set of present-day observational data in the range of Galactocentric distances R from 0 to 200 kpc is used to refine the parameters of thesemodels. For the Allen–Santillán model, a dimensionless coefficient γ has been included as a sought-for parameter for the first time. In the traditional and modified versions, γ = 2.0 and 6.3, respectively. Both versions are considered in this paper. The model rotation curves have been fitted to the observed velocities by taking into account the constraints on the local matter density ρ _⊙ = 0.1 M _⊙ pc^?3 and the force K _z =1.1/2πG = 77 M _⊙ pc^?2 acting perpendicularly to the Galactic plane. The Galactic mass within a sphere of radius 50 kpc, M _G (R ≤ 50 kpc) ≈ (0.41 ± 0.12) × 10¹² M _⊙, is shown to satisfy all three models. The differences between the models become increasingly significant with increasing radius R. In model I, the Galactic mass within a sphere of radius 200 kpc at γ = 2.0 turns out to be greatest among the models considered, M _G (R ≤ 200 kpc) = (1.45 ±0.30)× 10¹² M _⊙, M _G (R ≤ 200 kpc) = (1.29± 0.14)× 10¹² M _⊙ at γ = 6.3, and the smallest value has been found in model II, M _G (R ≤ 200 kpc) = (0.61 ± 0.12) × 10¹² M _⊙. In our view, model III is the best one among those considered, because it ensures the smallest residual between the data and the constructed model rotation curve provided that the constraints on the local parameters hold with a high accuracy. Here, the Galactic mass is M _G (R ≤ 200 kpc) = (0.75 ± 0.19) × 10¹² M _⊙. A comparative analysis with the models by Irrgang et al. (2013), including those using the integration of orbits for the two globular clusters NGC 104 and NGC 1851 as an example, has been performed. The third model is shown to have subjected to a significant improvement.     

A Model of the Mira–Type Star T UMi

Stellar evolution calculations were carried out from the main sequence to the final stage of the asymptotic giant branch for stars with initial masses 1 M_⊙ ≤ M_ZAMS ≤ 2 M_⊙ and metallicity Z = 0.01. Selected models of evolutionary sequences were used as initial conditions for solution of the equations of radiation hydrodynamics and time–dependent convection describing radial stellar pulsations. The study was aimed to construct the hydrodynamic models of Mira–type stars that show the secular decrease in the pulsation period Π commenced in 1970th at Π = 315 day. We show that such a condition for the period change is satisfied with evolutionary sequences 1 M_⊙ ≤ M_ZAMS ≤ 1.2 M_⊙ and the best agreement with observations is obtained for M_ZAMS = 1.2 M_⊙. The pulsation period reduction is due to both the stellar radius decrease during the thermal pulse of the helium burning shell and mode switch from the fundamental mode to the first overtone. Theoretical estimates of the fundament parameters of the star at the onset of pulsation period reduction are as follows: the mass is M = 0.93 M_⊙, the luminosity is L = 4080 L_⊙, and the radius is R = 220 R_⊙. The mode switch occurs 35 years after the onset of period reduction.     

Ultraluminous quasars as a gravitational mesolensing effect

Two quasars SDSS J010013.02+280225.8 and J030642.51+185315.8 with redshifts z = 6.30 and z = 5.363 were recently discovered. Their apparent magnitudes in the standard cosmological model give the luminosities of L_bol ～ 4.3 × 10¹⁴L_⊙ and L_bol ～ 3.4 × 10¹⁴L_⊙. In the framework of modern concepts it is accepted that the energy release of quasars is provided by the accretion onto black holes with masses of 1.24 ± 0.19 × 10¹⁰M_⊙ and 1.07 ± 0.27 × 10¹⁰M_⊙. As within the standard cosmological model the ages of these objects are about one billion years, this creates serious difficulties for the scenario of formation of such objects. Here we interpret the ultra-high luminosities of quasars as the effect of lensing of their radiation by the foreground globular clusters or dwarf galaxies.     

Unsteady mass outflow from Wolf-Rayet stars

We show that hydrostatically equilibrium models for the thin photospheres of helium stars based on new opacities κ_R (OPAL and OP) can be constructed only for masses M<5M _⊙. The parameter Г=κL/4πGMc, defined as the ratio of light pressure to gravity, exceeds a critical value of 1.0 for larger masses, which must result in mass outflow under light pressure. This mass limit matches the observed lower limit for the masses of Wolf-Rayet stars (M _WR>5M _⊙)), which is an additional argument that the Wolf-Rayet stellar cores are actually helium stars. By solving the equation of radiative transfer in extended atmospheres, we construct a semiempirical model for a WN5 star (M _WN5=10M _⊙)) with a helium core and an expanding envelope, whose physical and geometric parameters are known mainly from light-curve solution for the eclipsing binary V444 Cyg (WN5+06): outflow rate $\dot M \approx 1.0 \times 10^{ - 5} M_ \odot yr^{ - 1} $ , terminal velocity V _∞≈2000 km s^?1, and expanding-envelope optical depth τ_env≈25. The temperature at the outer boundary of the photosphere of a helium star surrounded by such an envelope is approximately 130 kK higher than that in the absence of an envelope, being T _ph≈240 kK. Because of the high temperatures, the absorption coefficients at the corresponding photospheric levels are smaller than those in models with no envelope; therefore, the photosphere turns out to be in hydrostatic equilibrium and stable against light pressure (Г_max≈0.9). As a way out of this conflicting situation (an expanding envelope together with a hydrostatically equilibrium photosphere), we propose a model of discrete mass outflow, which is also supported by the observed cloudy structure of the envelopes in this type of stars. To quantitatively estimate parameters of the nonuniform outflow model requires detailed gas-dynamical calculations.     

Evolution and Period Change in RR Lyr Variables of the Globular Cluster M 3

The grid of evolutionary tracks of population II stars with initial masses 0.81 M_⊙ ≤ M_ZAMS ≤ 0.85 M_⊙ and chemical composition of the globular cluster M3 is computed. Selected models of horizontal branch stars were used as initial conditions for solution of the equations of radiation hydrodynamics and time–dependent convection describing radial stellar oscillations. The boundaries of the instability strip on the Herztsprung–Russel diagram were determined using ≈100 hydrodynamic models of RR Lyr pulsating variables. For each evolutionary track crossing the instability strip the pulsation period was determined as a function of evolutinary time. The rate of period change of most variables is shown to range within ?0.02 ≤ \(\dot{\Pi}\) ≤ 0.05 day/10⁶ yr. Theoretical estimate of the mean period change rate obtained by the population synthesis method is 〈\(\dot{\Pi}\)〉 = 6.0 × 10^?3 day/10⁶ yr and agrees well with observations of RR Lyr variables of the globular cluster M3.     

Graphite grain formation in the atmospheres of R coronae borealis stars

The results of calculations of graphite grain formation in the atmospheres of R CrB stars are given. The parameters for the models wereM=1M _⊙,M _bol=?6 mag. The effective temperature ranged from 5300K to 8300K. The chemical composition corresponded to the hydrogen-deficient carbon rich mixture:X=0,Y=0.9,Z _c=0.1. The results obtained show the existence of a critical mass loss rate which is ranged fromM _*≈10^?6 M _⊙yr^?1 forT _eff=5300 K toM _*≈10^?5 M _⊙ yr^?1 forT _eff=8300 K. As soon as the rate of mass loss exceedsM _* by 3–5 times the degree of condensation of carbon changes from 0 to 0.7. The finite radii of grains are about from 0.01 μm to 0.6 μm depending on the density near the condensation point, the velocity of matter outflow, and the stellar effective temperature. The duration of grain growth should amount to some dozens of days. It is supposed that the most probable explanation of dust-shell formation around R CrB stars is graphite condensation behind a shock wave arising from nonlinear stellar pulsation.     

Supersonic steady accretion of a rotating cold “Iron Gas” onto a neutron star after gravitational collapse

We analytically generalize the well-known solution of steady supersonic spherically symmetric gas accretion onto a star (Bondi 1952) for an iron atmosphere with completely degenerate electrons with an arbitrary degree of relativity. This solution is used for typical physical conditions in the vicinity of protoneutron stars produced by gravitational collapse with masses M ₀=(1.4?1.8)M _⊙ and over a wide range of nonzero “iron gas” densities at infinity, ρ_∞=(10⁴?5×10⁶)g cm^?3. Under these conditions, we determine all accretion parameters, including the accretion rate, whose value is ～(10?50)M _⊙s^?1 at M ₀=1.8M _⊙ (it is a factor of 1.7 lower for M ₀=1.4M _⊙, because the accretion rate is exactly ∝M ₀ ² ). We take into account the effect of accreting-gas rotation in a quasi-one-dimensional approximation, which has generally proved to be marginal with respect to the accretion rate.     

Mass—luminosity relation and pulsational properties of Wolf—Rayet stars

The evolution of Population I stars with initial masses 70M _⊙ ≤ M _ZAMS ≤ 130M _⊙ is considered. The computations were performed under various assumptions about the mass loss rate and were terminated at the phase of gravitational contraction after core helium exhaustion. The mass loss rate at the helium burning phase, ?_3α , is shown to be the main parameter that determines the coefficients of the mass—luminosity relation for Wolf—Rayet stars. Several more accurate mass—luminosity relations for mass loss rates ? = f _3α ?_3α , where 0.5 ≤ f _3α ≤ 3, are suggested, along with the mass—luminosity relation that combines all of the evolutionary sequences considered. The results of the stellar evolution computations were used as initial conditions in solving the hydrodynamic equations describing the spherically symmetric motions of a self-gravitating gas. The outer layers of massive Population I stars are unstable against radial oscillations throughout the helium burning phase. The oscillation amplitude is largest at enhanced carbon and oxygen abundances in the outer stellar layers, i.e., at a lower initial stellar mass M _ZAMS or a lower mass loss rate during the entire preceding evolution. In the course of evolution, the radial oscillation amplitude decreases and the small nonlinearity of the oscillations at M < 10M _⊙ allow the integral of mechanical work W done by an elementary spherical layer of gas in a closed thermodynamic cycle to be calculated with the necessary accuracy. The maximum of the radial dependence of W is shown to be located in layers with a gas temperature T ～ 2 × 10⁵ K, where the oscillations are excited by the iron Z-bump κ-mechanism. Comparison of the radial dependences of the integral of mechanical work W and the amplitude of the radiative flux variations suggests that the nonlinear radial oscillations of more massive Wolf—Rayet stars are also excited by the κ-mechanism.     

Spectroscopic study of the envelope of Nova Mon 2012, a γ-ray source

Based on our spectrophotometric observations, we have studied the envelope of the HeN Nova Mon 2012. The abundances of some chemical elements in the envelope and its mass have been estimated. Our results show that the helium, nitrogen, oxygen, and neon abundances in the Nova envelope exceed the solar ones by a factor of 1.5, 33, 9, and 95, respectively. The envelope mass has been found to be 2.3 × 10^?4 M _⊙.     

A family of exact solutions of Einstein-Maxwell field equations in isotropic coordinates: an application to optimization of quark star mass

In the present paper, we have obtained a class of charged super dense star models, starting with a static spherically symmetric metric in isotropic coordinates for perfect fluid by considering Hajj-Boutros (in J. Math. Phys. 27:1363, 1986) type metric potential and a specific choice of electrical intensity which involves a parameter K. The resulting solutions represent charged fluid spheres joining smoothly with the Reissner-Nordstrom metric at the pressure free interface. The solutions so obtained are utilized to construct the models for super-dense star like neutron stars (ρ _b =2 and 2.7×10¹⁴ g/cm³) and Quark stars (ρ _b =4.6888×10¹⁴ g/cm³). Our solution is well behaved for all values of n satisfying the inequalities \(4 < n \le4(4 + \sqrt{2} )\) and K satisfying the inequalities 0≤K≤0.24988, depending upon the value of n. Corresponding to n=4.001 and K=0.24988, we observe that the maximum mass of quark star M=2.335M _⊙ and radius R=10.04 km. Further, this maximum mass limit of quark star is in the order of maximum mass of stable Strange Quark Star established by Dong et al. (in arXiv:1207.0429v3, 2013). The robustness of our results is that the models are alike with the recent discoveries.     

A family of well behaved charge analogues of Durgapal’s perfect fluid exact solution in general relativity

This paper presents a new family of interior solutions of Einstein–Maxwell field equations in general relativity for a static spherically symmetric distribution of a charged perfect fluid with a particular form of charge distribution. This solution gives us wide range of parameter, K, for which the solution is well behaved hence, suitable for modeling of superdense star. For this solution the gravitational mass of a star is maximized with all degree of suitability by assuming the surface density equal to normal nuclear density, ρ _nm=2.5×10¹⁷ kg?m^?3. By this model we obtain the mass of the Crab pulsar, M _Crab, 1.36M _⊙ and radius 13.21 km, constraining the moment of inertia >?1.61×10³⁸ kg?m² for the conservative estimate of Crab nebula mass 2M _⊙. And M _Crab=1.96M _⊙ with radius R _Crab=14.38 km constraining the moment of inertia >?3.04×10³⁸ kg?m² for the newest estimate of Crab nebula mass, 4.6M _⊙. These results are quite well in agreement with the possible values of mass and radius of Crab pulsar. Besides this, our model yields moments of inertia for PSR J0737-3039A and PSR J0737-3039B, I _A =1.4285×10³⁸ kg?m² and I _B =1.3647×10³⁸ kg?m² respectively. It has been observed that under well behaved conditions this class of solutions gives us the overall maximum gravitational mass of super dense object, M _G(max)=4.7487M _⊙ with radius $R_{M_{\max}}=15.24~\mathrm{km}$ , surface redshift 0.9878, charge 7.47×10²⁰ C, and central density 4.31ρ _nm.     

Photometric elements and apsidal motion of the eclipsing binary NSV 18773

The photometric elements of the eclipsing binary NSV 18773 (HD 99898) have been determined for the first time by analyzing its V-and I-band light curves from the ASAS-2 and ASAS-3 catalogs. Based on these elements and using other published spectroscopic and photometric data, we constructed a consistent system of geometrical and physical parameters for the system that consists of two stars (M ₁ = 20M _⊙, Sp₁=B0V, R ₁ = 5.0R _⊙ and M ₂ = 14M _⊙, Sp₂ = B1V, R ₂ = 6.5R _⊙) in elliptical orbits (P = 5 _. ^d 049, e = 0.365, a = 40.1R _⊙). The distance to the system is d = 3.3 kpc, the interstellar extinction is A _V = 2 _. ^m 0, and the age is t = 2.8 × 10⁶ yr. NSV 18773 is a visual binary with components V _A = 9 _. ^m 9 and V _B = 10 _. ^m 3 separated by 0 _. ^" 8. The third light (L ₃ = 0.61) that we found by analyzing the light curves shows that the eclipsing binary is the system’s fainter component B. We confirmed the rapid apsidal motion of the star detected by Otero and Wils (2006) and refined its observed period: U _obs = 150 ± 6 yr. Our photometric elements and physical parameters allowed the apsidal parameter $\bar k_2^{obs} = 0.0135(14)$ , which reflects the density distribution along the radii of the component stars, to be determined. Within the error limits, the derived parameter agrees with its theoretically expected value, $\bar k_2^{th} = 0.0119(8)$ , from current evolutionary models of stars of the corresponding masses and ages.