A novel look at the pulsar force-free magnetosphere

The stationary axisymmetric force-free magnetosphere of a pulsar is considered. We present an exact dipolar solution of the pulsar equation, construct the magnetospheric model on its basis and examine its observational support. The new model has toroidal rather than common cylindrical geometry, in line with that of the plasma outflow observed directly as the pulsar wind nebula at much larger spatial scale. In its new configuration, the axisymmetric magnetosphere consumes the neutron star rotational energy much more efficiently, implying re-estimation of the stellar magnetic field, \(B_{\mathrm{new}}^{0}=3.3\times10^{-4}B/P\), where \(P\) is the pulsar period. Then the 7-order scatter of the magnetic field derived from the rotational characteristics of the pulsars observed appears consistent with the \(\cot\chi\)-law, where \(\chi\) is a random quantity uniformly distributed in the interval \([0,\pi/2]\). Our result is suggestive of a unique actual magnetic field strength of the neutron stars along with a random angle between the magnetic and rotational axes and gives insight into the neutron star unification on the geometrical basis.     

Precision pulsar timing with the ORT and the GMRT and its applications in pulsar astrophysics

Radio pulsars show remarkable clock-like stability, which make them useful astronomy tools in experiments to test equation of state of neutron stars and detecting gravitational waves using pulsar timing techniques. A brief review of relevant astrophysical experiments is provided in this paper highlighting the current state-of-the-art of these experiments. A program to monitor frequently glitching pulsars with Indian radio telescopes using high cadence observations is presented, with illustrations of glitches detected in this program, including the largest ever glitch in PSR B0531+21. An Indian initiative to discover sub-\(\mu \)Hz gravitational waves, called Indian Pulsar Timing Array (InPTA), is also described briefly, where time-of-arrival uncertainties and post-fit residuals of the order of \(\mu \)s are already achievable, comparable to other international pulsar timing array experiments. While timing the glitches and their recoveries are likely to provide constraints on the structure of neutron stars, InPTA will provide upper limits on sub-\(\mu \)Hz gravitational waves apart from auxiliary pulsar science. Future directions for these experiments are outlined.     

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.     

Nature of Coherent Radio Emission from Pulsars

The pulsar radio emission originates from regions below 10% of the light cylinder radius. This requires a mechanism where coherent emission is excited in relativistic pair plasma with frequency \(\nu _{\mathrm{cr}}\) which is below the plasma frequency \(\nu _{\circ }\) i.e. \(\nu _{\mathrm{cr}} < \nu _{\circ }\). A possible model for the emission mechanism is charged bunches (charged solitons) moving relativistically along the curved open dipolar magnetic field lines capable of exciting coherent curvature radio emission. In this article, we review the results from high quality observations in conjunction with theoretical models to unravel the nature of coherent curvature radio emission in pulsars.     

Timing irregularities of PSR J1705–1906

Timing analysis of PSR J1705–1906 using data from Nanshan 25-m and Parkes 64-m radio telescopes, which span over fourteen years, shows that the pulsar exhibits significant proper motion, and rotation instability. We updated the astrometry parameters and the spin parameters of the pulsar. In order to minimize the effect of timing irregularities on measuring its position, we employ the Cholesky method to analyse the timing noise. We obtain the proper motion of \(-77(3)\) mas?yr^?1 in right ascension and \(-38(29)\) mas?yr^?1 in declination. The power spectrum of timing noise is analyzed for the first time, which gives the spectral exponent \(\alpha =-5.2\) for the power-law model indicating that the fluctuations in spin frequency and spin-down rate dominate the red noise. We detect two small glitches from this pulsar with fractional jump in spin frequency of \(\Delta \nu /\nu \sim 2.9 \times 10^{-10}\) around MJD 55199 and \(\Delta \nu /\nu \sim 2.7\times 10^{-10}\) around MJD 55953. Investigations of pulse profile at different time segments suggest no significant changes in the pulse profiles around the two glitches.     

Millisecond pulsar radiation properties

Two investigations of millisecond pulsar radiation are discussed: average total intensity pulse morphology and individual pulse to pulse fluctuations. The average emission profiles of millisecond pulsars are compared with those of slower pulsars in the context of polar cap models. In general the full widths of pulsar emission regions continue to widen inversely with periodP as P^-(0.30-0.5) as expected for dipole polar cap models. Many pulse components are very narrow. The period scaling of pulsar profiles -separations and widths -can tell us about the angular distribution of radiating currents. An investigation of individual pulses from two millisecond pulsars at 430 MHz shows erratic pulse to pulse variations similar to that seen in slow pulsars. PSR B1937+21 displays occasional strong pulses that are located in the trailing edge of the average profile with relative flux densities in the range of 100 to 400. These are similar to the giant pulses seen in the Crab pulsar.     

Revisiting Field Burial by Accretion onto Neutron Stars

The surface magnetic field strength of millisecond pulsars (MSPs) is found to be about 4 orders of magnitude lower than that of garden variety radio pulsars (with a spin of \({\sim }\)0.5–5 s and \(B{\sim }10^{12}\hbox { G}\)). The exact mechanism of the apparent reduction of field strength in MSPs is still a subject of debate. One of the proposed mechanisms is burial of the surface magnetic field under matter accreted from a companion. In this article we review the recent work on magnetic confinement of accreted matter on neutron stars poles. We present the solutions of the magneto-static equations with a more accurate equation of state of the magnetically confined plasma and discuss its implications for the field burial mechanism.     

Spectroscopic Variability of Supergiant Star HD14134, B3Ia

Profile variations in the \(\hbox {H}\alpha \) and \(\hbox {H}\beta \) lines in the spectra of the star HD14134 are investigated using observations carried out in 2013–2014 and 2016 with the 2-m telescope at the Shamakhy Astrophysical Observatory. The absorption and emission components of the \(\hbox {H}\alpha \) line are found to disappear on some observational days, and two of the spectrograms exhibit inverse P-Cyg profile of \(\hbox {H}\alpha \). It was revealed that when the \(\hbox {H}\alpha \) line disappeared or an inversion of the P-Cyg-type profile is observed in the spectra, the \(\hbox {H}\beta \) line is displaced to the longer wavelengths, but no synchronous variabilities were observed in other spectral lines (CII \( \lambda \) 6578.05 Å, \( \lambda \) 6582.88 Å  and HeI \( \lambda \) 5875.72 Å) formed in deeper layers of the stellar atmosphere. In addition, the profiles of the \(\hbox {H}\alpha \) and \(\hbox {H}\beta \) lines have been analysed, as well as their relations with possible expansion, contraction and mixed conditions of the atmosphere of HD14134. We suggest that the observational evidence for the non-stationary atmosphere of HD14134 can be associated in part with the non-spherical stellar wind.     

Interstellar communication: Short pulse duration limits of optical SETI

Previous and ongoing searches for extraterrestrial optical and infrared nanosecond laser pulses and narrow line-width continuous emissions have so far returned null results. At the commonly used observation cadence of \(\sim 10^{-9}\,\hbox {s}\), sky-integrated starlight is a relevant noise source for large field-of-view surveys. This can be reduced with narrow bandwidth filters, multipixel detectors, or a shorter observation cadence. We examine the limits of short pulses set by the uncertainty principle, interstellar scattering, atmospheric scintillation, refraction, dispersion and receiver technology. We find that optimal laser pulses are time-bandwidth limited Gaussians with a duration of \(\Delta t \approx \,10^{-12}\,\hbox {s}\) at a wavelength \(\lambda _{0}\approx 1\,\upmu \hbox {m}\), and a spectral width of \(\Delta \lambda \approx 1.5\,\hbox {nm}\). Shorter pulses are too strongly affected through Earth’s atmosphere. Given certain technological advances, survey speed can be increased by three orders of magnitude when moving from ns to ps pulses. Faster (and/or parallel) signal processing would allow for an all-sky-at-once survey of lasers targeted at Earth.     

Supernova neutrinos: Flavor conversion independent of their mass

In extremely dense neutrino environments like in supernova core, the neutrino-neutrino refraction may give rise to self-induced flavor conversion. These neutrino flavor oscillations are well understood from the idea of the exponentially growing modes of the interacting oscillators in the flavor space. Until recently, the growth rates of these modes were found to be of the order of the vacuum oscillation frequency \(\Delta m^2/2E\) [\(\mathcal {O}(1~\mathrm{km}^{-1})\)] and were considered slow growing. However, in the last couple of years it was found that if the system was allowed to have different zenith-angle distributions for the emitted \(\nu _e\) and \(\bar{\nu }_e\) beams then the fastest growing modes of the interacting oscillators grew at the order of \(\mu =\sqrt{2} G_\mathrm{F}n_{\nu }\), a typical \(\nu \)–\(\nu \) interaction energy [\(\mathcal {O}(10^5~\mathrm{km}^{-1})\)]. Thus the growth rates are very large in comparison to the so-called ‘slow oscillations’ and can result in neutrino flavor conversion on a much faster scale. In fact, the point that the growth rates are no longer dependent on the vacuum oscillation frequency \(\Delta m^2/2E\), makes these ‘fast flavor conversions’ independent of \(\Delta m^2\) (thus mass) and energy. This is a surprising result as neutrino flavor conversions are considered to be the ultimate proof of massive neutrinos. However, the importance of this effect in the realistic astrophysical scenarios still remains to be understood.     

A satellite relative motion model including $$J_2$$ and $$J_3$$J3 via Vinti’s intermediary

Vinti’s potential is revisited for analytical propagation of the main satellite problem, this time in the context of relative motion. A particular version of Vinti’s spheroidal method is chosen that is valid for arbitrary elliptical orbits, encapsulating \(J_2\), \(J_3\), and generally a partial \(J_4\) in an orbit propagation theory without recourse to perturbation methods. As a child of Vinti’s solution, the proposed relative motion model inherits these properties. Furthermore, the problem is solved in oblate spheroidal elements, leading to large regions of validity for the linearization approximation. After offering several enhancements to Vinti’s solution, including boosts in accuracy and removal of some singularities, the proposed model is derived and subsequently reformulated so that Vinti’s solution is piecewise differentiable. While the model is valid for the critical inclination and nonsingular in the element space, singularities remain in the linear transformation from Earth-centered inertial coordinates to spheroidal elements when the eccentricity is zero or for nearly equatorial orbits. The new state transition matrix is evaluated against numerical solutions including the \(J_2\) through \(J_5\) terms for a wide range of chief orbits and separation distances. The solution is also compared with side-by-side simulations of the original Gim–Alfriend state transition matrix, which considers the \(J_2\) perturbation. Code for computing the resulting state transition matrix and associated reference frame and coordinate transformations is provided online as supplementary material.     

An analysis of the convergence of Newton iterations for solving elliptic Kepler’s equation

In this note a study of the convergence properties of some starters \( E_0 = E_0(e,M)\) in the eccentricity–mean anomaly variables for solving the elliptic Kepler’s equation (KE) by Newton’s method is presented. By using a Wang Xinghua’s theorem (Xinghua in Math Comput 68(225):169–186, 1999) on best possible error bounds in the solution of nonlinear equations by Newton’s method, we obtain for each starter \( E_0(e,M)\) a set of values \( (e,M) \in [0, 1) \times [0, \pi ]\) that lead to the q-convergence in the sense that Newton’s sequence \( (E_n)_{n \ge 0}\) generated from \( E_0 = E_0(e,M)\) is well defined, converges to the exact solution \(E^* = E^*(e,M)\) of KE and further \( \vert E_n - E^* \vert \le q^{2^n -1}\; \vert E_0 - E^* \vert \) holds for all \( n \ge 0\). This study completes in some sense the results derived by Avendaño et al. (Celest Mech Dyn Astron 119:27–44, 2014) by using Smale’s \(\alpha \)-test with \(q=1/2\). Also since in KE the convergence rate of Newton’s method tends to zero as \( e \rightarrow 0\), we show that the error estimates given in the Wang Xinghua’s theorem for KE can also be used to determine sets of q-convergence with \( q = e^k \; \widetilde{q} \) for all \( e \in [0,1)\) and a fixed \( \widetilde{q} \le 1\). Some remarks on the use of this theorem to derive a priori estimates of the error \( \vert E_n - E^* \vert \) after n Kepler’s iterations are given. Finally, a posteriori bounds of this error that can be used to a dynamical estimation of the error are also obtained.     

Pulse phase-resolved analysis of SMC X-3 during its 2016–2017 super-Eddington outburst

The Be X-ray pulsar SMC X-3 underwent an extra long and ultraluminous giant outburst from 2016 August to 2017 March. The peak X-ray luminosity is up to \(\sim10^{39}~\mbox{erg/s}\), suggesting a mildly super-Eddington accretion onto the strongly magnetized neutron star. It therefore bridges the gap between the Galactic Be/X-ray binaries (\(L_{\mathrm{X}}^{\mathrm{peak}} \leq10^{38}~\mbox{erg/s}\)) and the ultraluminous X-ray pulsars (\(L_{\mathrm{X}}^{\mathrm{peak}} \geq10^{40}~\mbox{erg/s}\)) found in nearby galaxies. A number of observations were carried out to observe the outburst. In this paper, we perform a comprehensive phase-resolved analysis on the high quality data obtained with the Nustar and XMM-Newton, which were observed at a high and intermediate luminosity levels. In order to get a better understanding on the evolution of the whole extreme burst, we take the Swift results at the low luminosity state into account as well. At the early stage of outburst, the source shows a double-peak pulse profile, the second main peak approaches the first one and merges into the single peak at the low luminosity. The second main peak vanishes beyond 20 keV, and its radiation becomes much softer than that of the first main peak. The line widths of fluorescent iron line vary dramatically with phases, indicating a complicated geometry of accretion flows. In contrast to the case at low luminosity, the pulse fraction increases with the photon energy. The significant small pulse fraction detected below 1 keV can be interpreted as the existence of an additional thermal component located at far away from the central neutron star.     

GONG p-Mode Parameters Through Two Solar Cycles

We investigate the parameters of global solar p-mode oscillations, namely damping width \(\Gamma\), amplitude \(A\), mean squared velocity \(\langle v^{2}\rangle\), energy \(E\), and energy supply rate \(\mathrm{d}E/\mathrm{d}t\), derived from two solar cycles’ worth (1996?–?2018) of Global Oscillation Network Group (GONG) time series for harmonic degrees \(l=0\,\mbox{--}\,150\). We correct for the effect of fill factor, apparent solar radius, and spurious jumps in the mode amplitudes. We find that the amplitude of the activity-related changes of \(\Gamma\) and \(A\) depends on both frequency and harmonic degree of the modes, with the largest variations of \(\Gamma\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le3300~\upmu\mbox{Hz}\) and \(31\le l \le60\) with a minimum-to-maximum variation of \(26.6\pm0.3\%\) and of \(A\) for modes with \(2400~\upmu\mbox{Hz}\le\nu\le 3300~\upmu\mbox{Hz}\) and \(61\le l \le100\) with a minimum-to-maximum variation of \(27.4\pm0.4\%\). The level of correlation between the solar radio flux \(F_{10.7}\) and mode parameters also depends on mode frequency and harmonic degree. As a function of mode frequency, the mode amplitudes are found to follow an asymmetric Voigt profile with \(\nu_{\text{max}}=3073.59\pm0.18~\upmu\mbox{Hz}\). From the mode parameters, we calculate physical mode quantities and average them over specific mode frequency ranges. In this way, we find that the mean squared velocities \(\langle v^{2}\rangle\) and energies \(E\) of p modes are anticorrelated with the level of activity, varying by \(14.7\pm0.3\%\) and \(18.4\pm0.3\%\), respectively, and that the mode energy supply rates show no significant correlation with activity. With this study we expand previously published results on the temporal variation of solar p-mode parameters. Our results will be helpful to future studies of the excitation and damping of p modes, i.e., the interplay between convection, magnetic field, and resonant acoustic oscillations.     

Viscoelastic tides: models for use in Celestial Mechanics

In this note, by using Smale’s \(\alpha \)-theorem on the convergence of Newton’s method, the \(\alpha \)-sets of convergence of some starters of solving the elliptic Kepler’s equation are derived. For each starter we compute the exact \(\alpha \)-set in the eccentricity-main anomaly \((e,M)\in [0,1)\times [0,\pi ]\), showing that these sets are larger than those derived by Avendaño et al. (Celest Mech Dyn Astron 119:27–44, 2014). Further, new convergence tests based on the Newton–Kantorowitch theorem are given comparing with the derived from Smale’s \(\alpha \)-test.     

$$M_{\bullet } - \sigma $$ relation in spherical systems

To investigate the \(M_\bullet -\sigma \) relation, we consider realistic elliptical galaxy profiles that are taken to follow a single power-law density profile given by \(\rho (r) = \rho _{0}(r/ r_{0})^{-\gamma }\) or the Nuker intensity profile. We calculate the density using Abel’s formula in the latter case by employing the derived stellar potential; in both cases. We derive the distribution function f(E) of the stars in the presence of the supermassive black hole (SMBH) at the center and hence compute the line-of-sight (LoS) velocity dispersion as a function of radius. For the typical range of values for masses of SMBH, we obtain \(M_{\bullet } \propto \sigma ^{p}\) for different profiles. An analytical relation \(p = (2\gamma + 6)/(2 + \gamma )\) is found which is in reasonable agreement with observations (for \(\gamma = 0.75{-}1.4\), \(p = 3.6{-}5.3\)). Assuming that a proportionality relation holds between the black hole mass and bulge mass, \(M_{\bullet } =f M_\mathrm{b}\), and applying this to several galaxies, we find the individual best fit values of p as a function of f; also by minimizing \(\chi ^{2}\), we find the best fit global p and f. For Nuker profiles, we find that \(p = 3.81 \pm 0.004\) and \(f = (1.23 \pm 0.09)\times 10^{-3}\) which are consistent with the observed ranges.     

Formation of Double Neutron Stars,Millisecond Pulsars and Double Black Holes

The 1982 model for the formation of Hulse–Taylor binary radio pulsar PSR B1913+16 is described, which since has become the ‘standard model’ for the formation of the double neutron stars, confirmed by the 2003 discovery of the double pulsar system PSR J0737-3039AB. A brief overview is given of the present status of our knowledge of the double neutron stars, of which 15 systems are presently known. The binary-recycling model for the formation of millisecond pulsars is described, as put forward independently by Alpar et al. (1982), Radhakrishnan & Srinivasan (1982) and Fabian et al. (1983). This now is the ‘standard model’ for the formation of these objects, confirmed by the discovery in 1998 of the accreting millisecond X-ray pulsars. It is noticed that the formation process of close double black holes has analogies to that of close double neutron stars, extended to binaries with larger initial component masses, although there are also considerable differences in the physics of the binary evolution at these larger masses.     

Analysis of the exoplanet containing system Kepler-13

We have applied the close binary system analysis program WinFitter, with its physically detailed fitting function, to an intensive study of the complex multiple system Kepler-13 using photometry data from all 13 short cadence quarters downloaded from the NASA Exoplanet Archive (NEA) (http://exoplanetarchive.ipac.caltech.edu). The data-point error of our normalized, phase-sequenced and binned (380 points per bin: 0.00025 phase interval) flux values, at 14 ppm, allows the model’s specification for the mean reference flux level of the system to a precision better than 1 ppm. Our photometrically derived values for the mass and radius of KOI13.01 are \(6.8\pm0.6~\mbox{M}_{\mathrm{J}}\) and \(1.44\pm0.04~\mbox{R}_{\mathrm{J}}\). The star has a radius of \(1.67\pm0.05~\mbox{R}_{\odot}\). Our modelling sets the mean of the orbital inclination \(i\) at \(94.35\pm0.14^{\circ}\), with the star’s mean precession angle \(\phi_{p}\)—\(49.1\pm5.0^{\circ}\) and obliquity \(\theta_{o}\)\(67.9 \pm 3.0^{\circ}\), though there are known ambiguities about the sense in which such angles are measured.Our findings did not confirm secular variation in the transit modelling parameters greater than their full correlated errors, as argued by previous authors, when each quarter’s data was best-fitted with a determinable parameter set without prejudice. However, if we accept that most of the parameters remain the same for each transit, then we could confirm a small but steady diminution in the cosine of the orbital inclination over the 17 quarter timespan. This is accompanied by a slight increase of the star’s precession angle (less negative), but with no significant change in the obliquity of its spin axis. There are suggestions of a history of strong dynamical interaction with a highly distorted planet rotating in a 3:2 resonance with its revolution, together with a tidal lag of \(\sim30~\mbox{deg}\). The mean precessional period is derived to be about 1000 y, but at the present time the motion of the star’s rotation axis appears to be supporting the gravitational torque, rather than providing the balance against it that would be expected over long periods of time.The planet has a small but detectable backwarming effect on the star, which helps to explain the difference in brightness just after transit and just before occultation eclipses. In assessing these findings it is recognized that sources of uncertainty remain, notably with possible inherent micropulsational effects, variations from other components of the multiple star, stellar activity, differential rotation and the neglect of higher order terms (than \(r_{1}^{5}\)) in the fitting function, where \(r_{1}\) is the ratio of the radius of the star to the mean orbital separation of planet and host star.     

The Cadmium Zinc Telluride Imager on AstroSat

The Cadmium Zinc Telluride Imager (CZTI) is a high energy, wide-field imaging instrument on AstroSat. CZTI’s namesake Cadmium Zinc Telluride detectors cover an energy range from 20 keV to \(>200\) keV, with 11% energy resolution at 60 keV. The coded aperture mask attains an angular resolution of 17\(^\prime \) over a 4.6\(^\circ \) \(\times \) 4.6\(^\circ \)  (FWHM) field-of-view. CZTI functions as an open detector above 100 keV, continuously sensitive to GRBs and other transients in about 30% of the sky. The pixellated detectors are sensitive to polarization above \(\sim \)100 keV, with exciting possibilities for polarization studies of transients and bright persistent sources. In this paper, we provide details of the complete CZTI instrument, detectors, coded aperture mask, mechanical and electronic configuration, as well as data and products.