Devasthal Fast Optical Telescope Observations of Wolf–Rayet Dwarf Galaxy Mrk 996

The Devasthal Fast Optical Telescope (DFOT) is a 1.3 meter aperture optical telescope, recently installed at Devasthal, Nainital. We present here the first results using an Hα filter with this telescope on a Wolf–Rayet dwarf galaxy Mrk 996. The instrumental response and the Hα sensitivity obtained with the telescope are (3.3 ± 0.3) × 10^???15 erg s^?1 cm^?2/counts s^?1 and 7.5 × 10^???17 erg s^?1 cm^?2 arcsec^?2 respectively. The Hα flux and the equivalent width for Mrk 996 are estimated as (132 ± 37) × 10^?14 erg s^?1 cm^?2 and ～96 Å respectively. The star formation rate is estimated as 0.4 ± 0.1M _⊙ yr^?1. Mrk 996 deviates from the radio-FIR correlation known for normal star forming galaxies with a deficiency in its radio continuum. The ionized gas as traced by Hα emission is found in a disk shape which is misaligned with respect to the old stellar disk. This misalignment is indicative of a recent tidal interaction in the galaxy. We believe that galaxy–galaxy tidal interaction is the main cause of the WR phase in Mrk 996.     

A hydrodynamic model for asymmetric explosions of rapidly rotating collapsing supernovae with a toroidal atmosphere

We numerically solved the two-dimensional axisymmetric hydrodynamic problem of the explosion of a low-mass neutron star in a circular orbit. In the initial conditions, we assumed a nonuniform density distribution in the space surrounding the collapsed iron core in the form of a stationary toroidal atmosphere that was previously predicted analytically and computed numerically. The configuration of the exploded neutron star itself was modeled by a torus with a circular cross section whose central line almost coincided with its circular orbit. Using an equation of state for the stellar matter and the toroidal atmosphere in which the nuclear statistical equilibrium conditions were satisfied, we performed a series of numerical calculations that showed the propagation of a strong divergent shock wave with a total energy of ～0.2×10⁵¹ erg at initial explosion energy release of ～1.0×10⁵¹ erg. In our calculations, we rigorously took into account the gravitational interaction, including the attraction from a higher-mass (1.9M_⊙) neutron star located at the coordinate origin, in accordance with the rotational explosion mechanism for collapsing supernovae. We compared in detail our results with previous similar results of asymmetric supernova explosion simulations and concluded that we found a lower limit for the total explosion energy.     

MHD simulations of the long-term evolution of a dipolar magnetosphere surrounded by an accretion disk

The evolution of a stellar, initially dipole type magnetosphere interacting with an accretion disk is investigated using numerical ideal MHD simulations. The simulations follow several 1000 Keplerian periods of the inner disk (for animated movies see http://www.aip.de～cfendt).Our model prescribes a Keplerian disk around a rotating star as a fixed boundary condition. The initial magnetic field distribution remains frozen into the star and the disk. The mass flow rate into the corona is fixed for both components. The initial dipole type magnetic field develops into a spherically radial outflow pattern with two main components – a disk wind and a stellar wind – both evolving into a quasi-stationary final state. A neutral field line divides both components, along which small plasmoids are ejected in irregular time intervals. The half opening angle of the stellar wind cone varies from 30° to55° depending on the ratio of the mass flow rates of disk wind and stellar wind. The maximum speed of the outflow is about the Keplerian speed at the inner disk radius. An axial jet forms during the first decades of rotations. However, this feature does not survive on the very long time scale and a pressure driven low velocity flow along the axis evolves. Within a cone of 15° along the axis the formation of knots may be observed if the stellar wind is weak. With the chosen mass flow rates and field strength we see almost no indication for a flow self-collimation. This is due to the weak net poloidal electric current in the magnetosphere which is in difference to typical jet models.     

Photoionization of cool MHD disk winds

Recent ultraviolet observations point out that there is hot, dense plasma associated with the optical jet in some T Tauri stars. In this contribution, cool MHD disk wind physics is reviewed by means of a self-similar analytical model to analyze whether hot (Te ? 80,000 K) and dense (n_e ? 10⁹ cm^-3) plasma can be produced in disk winds. It is shown that these high densities can only be achieved at the base of the wind where the stellar X-rays radiation field is strong. The propagation of the X-rays radiation through the disk wind is analyzed: a cocoon of photoionized gas is generated around the star. However, it is difficult to foresee how temperatures as high as ～ 5 × 10⁴ can be reached unless a significant fraction of the X-rays radiation is produced by magnetic reconnection at the boundary between the stellar magnetosphere and the accretion disk.     

Neutral hydrogen around the oxygen-sequence Wolf-Rayet star WR 102 and the nebula G2.4+1.4

We carried out the first 21-cm line observations of an extended region around the Wolf-Rayet star WR 102 and the associated nebula G2.4+1.4 with the RATAN-600 radio telescope. An irregular H I shell was identified. Its maximum expansion velocity reaches ～50 km s^?1, and its outer diameter (at a distance of 3 kpc) is 56 pc. The mechanical luminosity of the stellar wind required to produce the observed shell is estimated to be ～0.8×10³⁸ erg s^?1; the age of the shell is ～3.4×10⁵ yr. We compare the inferred parameters of the H I shell with the structure and kinematics of the ionized nebula and with the dust distribution in the region.     

Acceleration of charged particles in the magnetosphere of a collapsing star and their nonthermal electromagnetic radiation

We investigate a transformation of a magnetic field and plasma in nonhomogeneous magnetospheres of collapsing stars with a dipole initial magnetic field and certain initial energy distributions of particles in the magnetosphere as the power low, relativistic Maxwell and Boltzmann. The betatron mechanism of the charged particles acceleration in a collapsing star’s magnetosphere is considered. When a magnetized star is compressed in the stage of the gravitational collapse, the magnetic field increases strongly. This variable magnetic field generates a vortical electric field. Our calculations show that this electric field will accelerate charged particles up to relativistic velocities. Thus, collapsing stars may be sources of high energy cosmic rays in our galaxy as in others. The acceleration of particles during the collapse happens mostly in polar regions of the magnetosphere that leads to polar relativistic streams (jets) formation. When moving in a magnetic field, these particles will generate nonthermal electromagnetic radiation in a broad electromagnetic wavelength band from radioto gamma rays. Thus, in the stage of the gravitational collapse, relativistic jets are formed in stellar magnetospheres. These jets are powerful sources of the nonthermal electromagnetic radiation.     

Thermal emission areas of heated neutron star polar caps

Observed hot spots on neutron stars are often associated with polar caps heated by the backflow of energetic electrons or positrons from accelerators on bundles of open magnetic field lines. Three effects are discussed that may be relevant to formation of hot spots and their areas. (1) The area of a polar cap is proportional to the ratio of the star’s surface dipole field to the local field at the polar cap. Because the field is coupled to the evolving spin of the superfluid core of the star, this ratio can depend on the stellar spin and its history. (2) The hot emission area may appear smaller to a distant observer when emitted X-rays propagate through electron-positron plasma created in the magnetosphere. The X-rays then change their energy spectrum because of cyclotron resonant scattering by pairs. (3) Hot spots may form on the star’s surface as a result of crust motions that are driven by the pull of core flux tubes pinned to the crust. Such motions twist the footprints of closed magnetic loops of the magnetosphere and induce an electric current in the loop, which will heat those footprints.     

A young cooling neutron star in the remnant of Supernova 1987A

We describe the cooling theory for isolated neutron stars that are several tens of years old. Their cooling differs greatly from the cooling of older stars that has been well studied in the literature. It is sensitive to the physics of the inner stellar crust and even to the thermal conductivity of the stellar core, which is never important at later cooling stages. The absence of observational evidence for the formation of a neutron star during the explosion of Supernova 1987A is consistent with the fact that the star was actually born there. It may still be hidden in the dense center of the supernova remnant. If, however, the star is not hidden, then it should have a low thermal luminosity (below ～10³⁴ erg s^?1) and a short internal thermal relaxation time (shorter than 13 yr). This requires that the star undergo intense neutrino cooling (e.g., via the direct Urca process) and have a thin crust with strong superfluidity of free neutrons and/or an anomalously high thermal conductivity.     

4U 0115+63 from RXTE and INTEGRAL data: Pulse profile and cyclotron line energy

We analyze the observations of the transient X-ray pulsar 4U 0115+63 with the RXTE and INTEGRAL observatories in a wide X-ray (3–100 keV) energy band during its intense outbursts in 1999 and 2004. The energy of the fundamental harmonic of the cyclotron resonance absorption line near the maximum of the X-ray flux from the source (luminosity range 5 × 10³⁷–2 × 10³⁸ erg s^?1) is ～11 keV. When the pulsar luminosity falls below ～5 × 10³⁷ erg s^?1, the energy of the fundamental harmonic is displaced sharply toward the high energies, up to ～16 keV. Under the assumption of a dipole magnetic field configuration, this change in cyclotron harmonic energy corresponds to a decrease in the height of the emitting region by ～2 km, while other spectral parameters, in particular, the cutoff energy, remain essentially constant. At a luminosity ～7 × 10³⁷ erg s^?1, four almost equidistant cyclotron line harmonics are recorded in the spectrum. This suggests that either the region where the emission originates is compact or the emergent spectrum from different (in height) segments of the accretion column is uniform. We have found significant pulse profile variations with energy, luminosity, and time. In particular, we show that the profile variations from pulse to pulse are not reduced to a simple modulation of the accretion rate specified by external conditions.     

Ultrahigh energy cosmic ray acceleration in Seyfert galactic nuclei

We propose a model for the particle acceleration to energy E≈10²¹ eV in Seyfert galactic nuclei. The model is based on the theory of active galactic nuclei by Vilkoviskij et al. (1999). The acceleration takes place in hot spots of relativistic jets, which decay in a dense stellar kernel at a distance of 1–3 pc from the center. The maximum energy and chemical composition of the accelerated particles depend on the jet magnetic-field strength. Fe nuclei acquire the largest energy, E≈8×10²⁰ eV, if the jet field strength is B≈16 G. At a field strength B～5–40 G, the nuclei with Z≥10 acquire energy E≥2×10²⁰ eV; the lighter nuclei are accelerated to E≤10²⁰ eV. In a field B～1000 G, only the particles with Z≥23 gain energy E≤10²⁰ eV. The protons are accelerated to E<4×10¹⁹ eV, and they do not fall within the energy range concerned at any field strength B. Interactions with infrared photons do not affect the accelerated-particle escape from the sources if the galactic luminosity L≤10⁴⁶ erg s^?1 and if the angle between the normal to the galactic plane and the line of sight is sufficiently small, i.e., if the galactic-disk axial ratio is comparatively large. The particles do not lose their energy through magnetodrift radiation if their deflection from the jet axis does not exceed 0.03–0.04 pc at a distance R≈40–50 pc from the center. The synchrotron losses are small, because the magnetic field frozen in the galactic wind at R≤40–50 pc is directed (as in the jet) predominantly along the motion. If this model is correct, then the detected protons are nuclear fragments or are accelerated in other sources. The jet magnetic fields can be estimated by using the cosmic-ray energy spectrum and chemical composition.     

WD + He star systems as the progenitors of Type Ia supernovae and their surviving companion stars

Employing Eggleton’s stellar evolution code with the optically thick wind assumption, we have systematically studied the WD + He star channel of Type Ia supernovae (SNe Ia), in which a carbon–oxygen WD accretes material from a He main-sequence star or a He subgiant to increase its mass to the Chandrasekhar mass. We mapped out the parameter spaces for producing SNe Ia. According to a detailed binary population synthesis approach, we find that the Galactic SN Ia birthrate from this channel is ～0.3×10^?3 yr^?1, and that this channel can produce SNe Ia with short delay times (～45–140 Myr). We also find that the surviving companion stars in this channel have a high spatial velocity (>400 km/s) after the SN explosion, which could be an alternative origin for hypervelocity stars (HVSs), especially for HVSs such as US 708.     

Electron-positron pairs production in the pulsar magnetosphere

The case of an aligned rotator magnetosphere is considered. Provided the ions ejection from the neutron stellar surface is absent, the pulsar magnetosphere consists of two polar electron caps. The upper parts of the caps are unstable. Electrons precipitate from these parts, fall onto the star and are accelerated to Lorentz-factor 10⁶–10⁷. Electrons radiate -quanta in the direction of the star. These -quanta are converted into electron-positron pairs. The region of size, about 10 stellar radii, around the star appears to be filled with electron-positron plasma. The inflow of electron-positron plasma interacts with the electron gas of the polar cap. For this reason longitudinal plasma vibrations arise, and bunched outflows of electron-positron plasma appear.     

High-mass X-ray binaries and recent star formation history of the Small Magellanic Cloud

We study the relation between high-mass X-ray binary (HMXB) population and recent star formation history (SFH) for the Small Magellanic Cloud (SMC). Using archival optical SMC observations, we have approximated the color-magnitude diagrams of the stellar population by model stellar populations and, in this way, reconstructed the spatially resolved SFH of the galaxy over the past 100 Myr. We analyze the errors and stability of this method for determining the recent SFH and show that uncertainties in the models of massive stars at late evolutionary stages are the main factor that limits its accuracy. By combining the SFH with the spatial distribution of HMXBs obtained from XMM-Newton observations, we have derived the dependence of the HMXB number on the time elapsed since the star formation event. The number of young systems with ages ? 10 Myr is shown to be smaller than the prediction based on the type-II supernova rate. The HMXB number reaches its maximum ～20–50 Myr after the star formation event. This may be attributable, at least partly, to a low luminosity threshold in the population of X-ray sources studied, L _min ～ 10³⁴ erg s^?1. Be/X systems make a dominant contribution to this population, while the contribution from HMXBs with black holes is relatively small.     

Relationship between the nonstationary radio-luminosity variations of the Crab Nebula and the activity of its pulsar

Based on the half-century-long history of radio observations of the Crab Nebula, we investigate the evolution of its radio luminosity. We found a secular decrease in the radio luminosity; it has decreased by 9% since the discovery of the radio source in 1948. Apart from the secular decrease in the luminosity of the Crab Nebula, we identified two time intervals, 1981–1987 and 1992–1998, when radio bursts with energy release ～10⁴¹ erg took place. In these years, the spectral indices of the instantaneous spectra decreased significantly due to the increase in the flux densities at short (centimeter and millimeter) wavelengths. These events were preceded by sudden increases in the pulsar’s rotation rates, the largest of which, with an amplitude of ΔΩ/Ω = 3 × 10^?8, occurred in 1975 and 1989. We show that the magnetospheric instability mechanism that accompanies strong glitches can provide the energetics of the excess luminosity of the Nebula through the ejection of relativistic electrons with a total energy higher than 6 × 10⁴² erg from the pulsar’s magnetosphere.     

The formation and evolution of massive stellar clusters in IC4662

We present a multiwavelength study of the formation of massive stellar clusters, their emergence from cocoons of gas and dust, and their feedback on surrounding matter. Using data that span from radio to optical wavelengths, including Spitzer and Hubble Space Telescope ACS observations, we examine the population of young star clusters in the central starburst region of the irregular Wolf–Rayet galaxy IC4662. We model the radio-to-infrared (IR) spectral energy distributions of embedded clusters to determine the properties of their Hii regions and dust cocoons (sizes, masses, densities, temperatures), and use near-IR and optical data with mid-IR spectroscopy to constrain the properties of the embedded clusters themselves (mass, age, extinction, excitation, abundance). The two massive star-formation regions in IC4662 are excited by stellar populations with ages of ～4 Myr and masses of ～3×10⁵ M_⊙ (assuming a Kroupa initial mass function). They have high excitation and subsolar abundances, and they may actually be comprised of several massive clusters rather than the single monolithic massive compact objects known as ‘super star clusters’ (SSCs). Mid-IR spectra reveal that these clusters have very high extinction values, A _V ～20–25 mag, and that the dust in IC4662 is well mixed with the emitting gas, not in a foreground screen.     

The mechanism of magnetospheric substorm and the MHD waves of the solar wind

A mechanism of the Earth's magnetospheric substorm is proposed. It is suggested that the MHD waves may propagate across the magnetopause from the magnetosheath into the magnetotail and will be dissipated in the plasma sheet, heating the plasma and accelerating the particles. When the solar wind parameters change, the Poynting flux of the waves transferred from the magnetosheath into the tail, may be greater than 10¹⁸ erg s^?1. The heated plasma and accelerated particles in the plasma sheet will be injected into the inner magnetosphere, and this may explain the process of the ring current formation and auroral substorm.The Alfvén wave can only propagate along the magnetic force line into the magnetosphere in the open magnetosphere, but the magnetosonic wave can propagate in both the open and closed magnetosphere. When the IMF turns southward, the configuration of the magnetosphere will change from a nearly closed model into some kind of open one. The energy flux of Alfvén waves is generally larger than that of the magnetosonic wave. This implies that it is easy to produce substorms when the interplanetary magnetic field (IMF) has a large southward component, but the substorm can also be produced even if the IMF is directed northward.     

QED can explain the non-thermal emission from SGRs and AXPs: variability

Owing to effects arising from quantum electrodynamics (QED), magnetohydrodynamical fast modes of sufficient strength will break down to form electron-positron pairs while traversing the magnetospheres of strongly magnetised neutron stars. The bulk of the energy of the fast mode fuels the development of an electron-positron fireball. However, a small, but potentially observable, fraction of the energy (～10³³ erg) can generate a non-thermal distribution of electrons and positrons far from the star. This paper examines the cooling and radiative output of these particles. Small-scale waves may produce only the non-thermal emission. The properties of this non-thermal emission in the absence of a fireball match those of the quiescent, non-thermal radiation recently observed non-thermal emission from several anomalous X-ray pulsars and soft-gamma repeaters. Initial estimates of the emission as a function of angle indicate that the non-thermal emission should be beamed and therefore one would expect this emission to be pulsed as well. According to this model the pulsation of the non-thermal emission should be between 90 and 180 degrees out of phase from the thermal emission from the stellar surface.     

Electron-positron and photon wind from the surface of a strange star

We consider the problem of strange-star (SS) radiation. The bare quark SS surface and electrons on the stellar surface generate an electric field that is strong enough for electron-positron pairs to be produced from a vacuum at a nonzero temperature. The luminosity in pairs is assumed to be within ?10⁴⁹ erg s^?1 from a surface with a characteristic radius of 10 km. We consider the energy transfer from pairs to photons by taking into account the well-studied reactions between e, e ⁺, γ and obtain a change in the photon spectrum with luminosity. Our analysis is restricted to the spherically symmetric case. The magnetic field is disregarded. To solve the problem, we developed a new numerical method of integrating the Boltzmann kinetic equations for pairs and photons. This method is used to calculate the problem up to a luminosity of 10⁴² erg s^?1 This region is difficult to investigate when the optical path for pairs or photons is considerably larger than unity but the two optical depths are not simultaneously much larger than unity (when hydrodynamics with heat conduction is applicable). It turns out that the mean photon energy is approximately equal to $\bar \in _\gamma \approx m_e c^2$ (the annihilation line for pairs) at a modest luminosity, L?1×10³⁷ erg s^?1, and decreases to ≈210 keV at L?10³⁸ erg s^?1. Hydrodynamic estimates point to an increase in the mean energy $\bar \in _\gamma$ to 1 MeV as the luminosity further increases to L?10⁴⁹ erg s^?1. Our calculations may prove to be useful in interpreting soft gamma repeaters (SGRs) and are of methodological interest.     

Energetics of galactic nuclei

A model of compact galactic nuclei in statistical equilibrium was developed in [L. Sh. Grigorian and G. S. Sahakian, Astrofizika (in press)]. It was shown that they should consist predominantly of neutron stars (pulsars) and white dwarfs. The problem of the energy reserves of galactic nuclei is discussed in terms of this concept. The mechanism of conversion of a white dwarf into a neutron star due to the accretion of interstellar matter is considered. This means that a galactic nucleus has an energy reserve of some 5·10⁶⁰ N₈ erg (N is the number of stars in the nucleus). It is shown that galactic nuclei are powerful sources of hard γ radiation [power L » 2·10⁴⁴µ₃₀N₈(Ω/50)^17/7 erg/sec, where µ is the magnetic moment and Ω is the angular rotation rate of a neutron star ] due to curvature radiation from relativistic electron fluxes flowing along channels of open magnetic field lines of pulsars. The x-ray and ultraviolet emission are due to synchrotron emission from the same electron fluxes in the magnetic field of the galactic nucleus (L » 10⁴²-10⁴⁴ erg/sec). The optical (visible and infrared) and radio emission are due to bremsstrahlung from electrons in the interstellar medium [L » 6·10⁴⁶N ₈ ² (5/R_pc)³ erg/sec, where R is the radius of the galactic nucleus]. An equation is obtained for the magnetic moment of a pulsar: µ ≈ 3.4·10^-5L_γP^17/7, where P is the pulsar’s period and L_03B3; is the luminosity of the pulsar’s y radiation.     

Study of individual geomagnetic storms in terms of the solar wind

This paper summarizes a study of the development of a large number of geomagnetic storms in terms of the solar wind—magnetosphere energy coupling function ε, the AE and Dst indices. It is shown that the maximum magnitude of the main phase decrease (¦Dst¦) is determined primarily by the peak value of ε; for ε < 10¹⁹ erg s^?1, ~10¹⁹ erg s^?1, 10¹⁹–10²⁰ erg s^?1, ? 10²⁰ erg s^?1, the maximum values of ¦Dst¦ are < 50γ, ~50γ, ~100γ and ? 200γ, respectively. A few examples for different peak values of ε (and thus of ¦Dst¦) are presented and examined in detail. Substorm activity during storms is well controlled by ε.