Magnetic field decay of magnetars in supernova remnants

In this paper, we modify our previous research carefully, and derive a new expression of electron energy density in superhigh magnetic fields. Based on our improved model, we re-compute the electron capture rates and the magnetic fields’ evolutionary timescales t of magnetars. According to the calculated results, the superhigh magnetic fields may evolve on timescales ～(10⁶?10⁷) yrs for common magnetars, and the maximum timescale of the field decay, t≈2.9507×10⁶ yrs, corresponding to an initial internal magnetic field B ₀=3.0×10¹⁵ G and an initial inner temperature T ₀=2.6×10⁸ K. Motivated by the results of the neutron star-supernova remnant (SNR) association of Zhang and Xie (2011), we calculate the maximum B ₀ of magnetar progenitors, B _max～(2.0×10¹⁴?2.93×10¹⁵) G when T ₀=2.6×10⁸ K. When T ₀～2.75×10⁸?1.75×10⁸ K, the maximum B ₀ will also be in the range of ～10¹⁴?10¹⁵ G, not exceeding the upper limit of magnetic field of a magnetar under our magnetar model. We also investigate the relationship between the spin-down ages of magnetars and the ages of their SNRs, and explain why all AXPs associated with SNRs look older than their real ages, whereas all SGRs associated with SNRs appear younger than they are.     

On the origin of the high-Q high-frequency oscillations of magnetars

A model for the high-frequency (20–2400 Hz) quasi-periodic oscillations (QPOs) of magnetars based on the representation of coronal magnetic loops as equivalent electric RLC circuits is proposed. The observed periods of the QPOs and their high Q-factor (Q ≈ 10⁴–10⁵) are explained. It follows from the model that the QPOs can be excited not only in the tail of a flare but also before the main pulse. The parameters of the QPO source at the “ringing tail” stage of the flare from SGR 1806–20 on December 27, 2004, have been estimated: electric current I ≈ 3 × 10¹⁹ A, minimum magnetic field strength B _min ≈ 10¹³ G, and electron density n ≈ 2 × 10¹⁶ cm⁻³.     

Stokes Diagnostics of 2D MHD-Simulated Solar Magnetogranulation

The properties of solar magnetic fields on scales less than the spatial resolution of solar telescopes are studied. A synthetic infrared spectropolarimetric diagnostic based on a 2D MHD simulation of magnetoconvection is used for this. Analyzed are two time sequences of snapshots that likely represent two regions of the network fields with their immediate surroundings on the solar surface with unsigned magnetic flux densities of 300 and 140 G. In the first region from the probability density functions of the magnetic field strength it is found that the most probable field strength at log τ ₅=0 is equal to 250 G. Weak fields (B<500 G) occupy about 70% of the surface, whereas stronger fields (B>1000 G) occupy only 9.7% of the surface. The magnetic flux is −28 G and its imbalance is −0.04. In the second region, these parameters are correspondingly equal to 150 G, 93.3%, 0.3%, −40 G, and −0.10. The distribution of line-of-sight velocities on the surface of log τ ₅=−1 is estimated. The mean velocity is equal to 0.4 km s⁻¹ in the first simulated region. The average velocity in the granules is −1.2 km s⁻¹ and in the intergranules it is 2.5 km s⁻¹. In the second region, the corresponding values of the mean velocities are equal to 0, −1.8, and 1.5 km s⁻¹. In addition the asymmetry of synthetic Stokes V profiles of the Fe i 1564.8 nm line is analyzed. The mean values of the amplitude and area asymmetry do not exceed 1%. The spatially smoothed amplitude asymmetry is increased to 10% whereas the area asymmetry is only slightly varied.     

Magnetic fields in the nearby Universe,as observed in solar and planetary realms,stars, and interstellar starforming nurseries

Here I present a survey of planetary, stellar and clumpy interstellar magnetic fields, with an emphasis on discoveries made in the last decade. Ground-based and Earth–satellite observations of magnetic fields in astronomical objects are cataloged, as obtained at all multiple electromagnetic wavelengths.For each type of object, the basic observational properties are summarized, and the best theoretical scenario which accounts for the large body of observations is discussed. There is no single unifying magnetic model for all objects (although dynamo theory comes close for many compact objects, and helical theory for many extended objects).The strength of the observed magnetic fields reaches exagauss values (10¹⁸ G) in compact objects such as magnetars, and decreases to microgauss values (10^?6 G) in the clumpy interstellar medium.     

Equation of state and anisotropy of pressure of magnetars

Magnetars are the neutron stars with the highest magnetic fields up to 10¹⁵–10¹⁶ G. It has been proposed that they are also responsible for a variety of extra-galactic phenomena, ranging from giant flares in nearby galaxies to fast radio bursts. Utilizing a relativistic mean field model and a variable magnetic field configuration, we investigate the effects of strong magnetic fields on the equation of state and anisotropy of pressure of magnetars. It is found that the mass and radius of low-mass magnetars are weakly enhanced under the action of the strong magnetic field, and the anisotropy of pressure can be ignored. Unlike other previous investigations, the magnetic field is unable to violate the mass limit of the neutron stars.     

Magnetic fields of active galactic nuclei and quasars from the SDSS catalog

We analyze the spectropolarimetric observations of 12 candidates for quasars from the spectroscopic database of the SDSS Catalog. The magnetic fields of these objects are estimated in the context of a theory that includes the Faraday rotation of the polarization plane on the mean free path of a photon in the outflow from an accretion disk. As a result, we have determined the column density in the outflow, N _H ∼ 6 × 10²³ cm⁻², and the radial, B ∼ 1 G, and toroidal, B ∼ 600 G, magnetic fields.     

Fundamentals of collisionless shocks for astrophysical application, 2. Relativistic shocks

In this concise review of the recent developments in relativistic shock theory in the Universe we restrict ourselves to shocks that do not exhibit quantum effects. On the other hand, emphasis is given to the formation of shocks under both non-magnetised and magnetised conditions. We only briefly discuss particle acceleration in relativistic shocks where much of the results are still preliminary. Analytical theory is rather limited in predicting the real shock structure. Kinetic instability theory is briefed including its predictions and limitations. A recent self-similar relativistic shock theory is described which predicts the average long-term shock behaviour to be magnetised and to cause reasonable power-law distributions for energetic particles. The main focus in this review is on numerical experiments on highly relativistic shocks in (i) pair and (ii) electron-nucleon plasmas and their limitations. These simulations do not validate all predictions of analytic and self-similar theory and so far they do not solve the injection problem and the self-modification by self-generated cosmic rays. The main results of the numerical experiments discussed in this review are: (i) a confirmation of shock evolution in non-magnetised relativistic plasma in 3D due to either the lepton-Weibel instability (in pair plasmas) or to the ion-Weibel instability; (ii) the sensitive dependence of shock formation on upstream magnetisation which causes suppression of Weibel modes for large upstream magnetisation ratios σ>10⁻³; (iii) the sensitive dependence of particle dynamics on the upstream magnetic inclination angle θ _Bn , where particles of θ _Bn >34° cannot escape upstream, leading to the distinction between ‘subluminal’ and ‘superluminal’ shocks; (iv) particles in ultra-relativistic shocks can hardly overturn the shock and escape to upstream; they may oscillate around the shock ramp for a long time, so to speak ‘surfing it’ and thereby becoming accelerated by a kind of SDA; (v) these particles form a power-law tail on the downstream distribution; their limitations are pointed out; (vi) recently developed methods permit the calculation of the radiation spectra emitted by the downstream high-energy particles; (vii) the Weibel-generated downstream magnetic fields form large-amplitude vortices which could be advected by the downstream flow to large distances from the shock and possibly contribute to an extended strong field region; (viii) if cosmic rays are included, Bell-like modes can generate upstream magnetic turbulence at short and, by diffusive re-coupling, also long wavelengths in nearly parallel magnetic field shocks; (ix) advection of such large-amplitude waves should cause periodic reformation of the quasi-parallel shock and eject large-amplitude magnetic field vortices downstream where they contribute to turbulence and to maintaining an extended region of large magnetic fields.     

The Correlation between the Magnetic and Velocity Fields on the Full Solar Disk

A possible correlation between the magnetic and velocity fields has been analyzed based on the SOHO/MDI magnetograms and Dopplergrams. It is found that the observed large-scale weak magnetic field (weaker than 50 G (gauss)) is correlated with the velocity statistically. The curves of u⋅b with latitude, where u and b are the velocity and magnetic fields in a rectangular region (±15^○ in longitude, ±45^○ in latitude) on the Sun, show the same patterns in the years 2000, 2004, and 2007. The patterns indicate that u and b are positively correlated near the equator but are anti-correlated at the middle latitudes. For a strong magnetic field between 50 G and 3000 G, the curves of u⋅b with latitude show the same tendencies at the middle latitudes. Near the equator, however, the slope of the curve is positive in 2000 and is negative in 2004 and 2007. In addition, we give an estimation for the amplitude of the cross helicity h _χ (h_c=[`(u·b)]h_{\chi}=\overline{\mathbf{u}\cdot\mathbf{b}}) inferred from the MDI data, which is of the order of 10³ G m s⁻¹ near the center of the solar disk.     

Cooling of magnetars with internal layer heating

We model thermal evolution of magnetars with a phenomenological heat source in a spherical internal layer and compare the results with observations of persistent thermal radiation from magnetars. We show that the heat source should be located in the outer magnetar’s crust, at densities ρ≲5×10¹¹ g cm⁻³, and the heating rate should be ∼10²⁰ erg cm⁻³ s⁻¹. Heating deeper layers is extremely inefficient because the thermal energy is mainly radiated away by neutrinos and does not warm up the surface to the magnetar’s level. This deep heating requires too much energy; it is inconsistent with the energy budget of neutron stars.      

Dirac neutrino magnetic moment and possible time evolution of the neutrino signal from a supernova

We analyze the influence of neutrino helicity conversion, ν _L → ν _R, on the neutrino flux from a supernova attributable to the interaction of the Dirac neutrino magnetic moment with a magnetic field.We show that if the neutrino has a magnetic moment in the interval 10⁻¹³μ_B < μ_ν < 10⁻¹²μ_B and provided that a magnetic field of ∼10¹³–10¹⁴ G exists in the supernova envelope, a peculiar kind of time evolution of the neutrino signal from the supernova attributable to the resonance transition ν _L → ν _R in the magnetic field of the envelope can appear.     

Orbital period analysis of eclipsing dwarf novae: U Geminorum

A new orbital period analysis for U Geminorum is made by means of the standard O–C technique based on 187 times of light minima including the three newest CCD data from our observation. Although there are large scatter near 70,000 cycles in its O–C diagram, there is strong evidence (>99.9% confidence level) to show the secular increase of orbital period with a rate  s⁻¹. Using the physical parameters recently derived by Echevarría et al. (Astron. J. 134:262, 2007), the range of mass transfer rate for U Geminorum is estimated as from −3.5(5)×10⁻⁹ M _⊙ yr⁻¹ to −1.30(6)×10⁻⁸ M _⊙ yr⁻¹. Moreover, the data before 60,000 cycles shows the obvious quasi-period variations. The least square estimation of a ∼17.4 yr quasi-periodic variation superimposed on secular orbital period increase is derived. Considering the possibility that solar-type magnetic activity cycles in the secondary star of U Geminorum may produce the quasi-period variations of the orbital period, Applegate’s mechanism is discussed and the results indicate such mechanism has difficulty explaining the quasi-period variation for U Geminorum. Hence, we attempted to apply the light-travel time effect to interpret the quasi-period variation and found the perturbation of ∼17.4 yr quasi-period may result from a brown dwarf. If the orbital inclination is assumed as i∼15°, corresponding to the upper limit of mass of a brown dwarf, then its orbital radii is ∼7.7 AU.     

The ordinary URCA process in the presence of positrons and large magnetic fields

The effect of positron capture on the ordinary URCA neutrino luminosity in a zero magnetic field is investigated for several values of the degeneracy parameter and the range of temperatures 5×10⁸K–5×10¹⁰K. The rate for this process is then compared with those in large magnetic fields (on the order ofH _c =m ² c ³/eh=4.414×10¹⁰ G). The results indicate that positron capture reduces the effect of large magnetic fields on this process at high temperatures.     

Two-photon positron–electron annihilation in a strong magnetic field

We consider the two-photon positron and electron annihilation in flight, it means the annihilating particles exhibit the middly relativistic momenta in a super strong magnetic field. Such particles are present in the corona of pulsars and magnetars. The paper presents how the total emission rate for the two-photon process is affected not only by magnetic field but also by the relativistic momentum of the annihilating particles. We found that the momenta influence significantly the total emission rate and the directions of the emitted photons. Additionally, the total emission for the two-photon process is comparable to that for the one-photon process at the momentum of annihilating particles of about m₀, where m₀ is the electron mass, and the magnetic field close to the critical Schwinger value of 4.41 × 10¹³ G. The latter is reported as a main annihilation channel in a super strong magnetic field. We calculated also the energetic spectra of annihilating photons emitted, which are also affected by the magnetic field and the momenta of the annihilating particles.     

Newborn magnetars as sources of gravitational radiation: constraints from high energy observations of magnetar candidates

Two classes of high-energy sources, the Soft Gamma Repeaters and the Anomalous X-ray Pulsars are believed to contain slowly spinning “magnetars,” i.e. neutron stars the emission of which derives from the release of energy from their extremely strong magnetic fields (>10¹⁵ G). The enormous energy liberated in the 2004 December 27 giant flare from SGR 1806-20 (～5×10⁴⁶ erg), together with the likely recurrence time of such events, points to an internal magnetic field strength of ≥10¹⁶ G. Such strong fields are expected to be generated by a coherent α?Ω dynamo in the early seconds after the Neutron Star (NS) formation, if its spin period is of a few milliseconds at most. A substantial deformation of the NS is caused by such fields and, provided the deformation axis is offset from the spin axis, a newborn millisecond-spinning magnetar would thus radiate for a few days a strong gravitational wave signal the frequency of which (～0.5–2 kHz range) decreases in time. This signal could be detected with Advanced LIGO-class detectors up to the distance of the Virgo cluster, where ≥1 yr^?1 magnetars are expected to form. Recent X-ray observations revealed that SNRs around magnetar candidates do not appear to have received a larger energy input than in standard SNRs (see Vink and Kuiper, Mon. Not. Roy. Astron. Soc. 319, L14 (2006)). This is at variance with what would be expected if the spin energy of the young, millisecond NS were radiated away as electromagnetic radiation and/or relativistic particle winds. In fact, such energy would be transferred quickly and efficiently to the expanding gas shell. This may thus suggest that magnetars did not form with the expected very fast initial spin. We show here that these findings can be reconciled with the idea of magnetars being formed with fast spins, if most of their initial spin energy is radiated through GWs. In particular, we find that this occurs for essentially the same parameter range that would make such objects detectable by Advanced LIGO-class detectors up to the Virgo Cluster. If our argument holds for at least a fraction of newly formed magnetars, then these objects constitute a promising new class of gravitational wave emitters.     

Search for nearby host galaxies of short gamma-ray bursts detected and well localized by BATSE/IPN

We have collected the observational data accumulated before the Swift experiment to check the possible connection of short gamma-ray bursts (GRBs) with low-redshift galaxies. The BATSE/IPN experimental data on well-localized short GRBs and the SDSS DR5 and PSCz catalogs of galaxies are used. The PSCz sky coverage has allowed us to search for host galaxies for a sample of 34 short GRBs. One or more galaxies have been found in the error boxes of six bursts, but the probability of a chance coincidence for each of them is high. No excess of nearby galaxies in the total sample has been detected. The 90% confidence limit corresponds to the fact that no more than 7%of the short GRBs could originate in nearby galaxies of the PSCz sample. The estimated upper limit of several percent may be considered to be valid in the volume z = 0.015–0.025. Based on the results of our search, we have estimated the lower limits for the isotropic energies E _γiso of 31 short bursts from our sample. Their values lie within the range 1.0 × 10⁴⁷–2.7 × 10⁴⁹ erg. The possible fraction of the flares from magnetars in our sample of short GRBs is discussed. The SDSS sky coverage is currently insufficient to perform a similar analysis.     

Evidence for a coronal magnetic bottle at 10 solar radii

The Faraday rotation of a radio source (Pioneer 6) occulted by the solar corona has been measured by Levy et al. (1969). During the course of these measurements, three large-scale transient phenomena were observed. These events were preceded by subflares and class 1 flares. These transient events are interpreted as evidence for a coronal magnetic bottle at 10 R _⊙. The velocity of propagation for the disturbance is set at 200 km/sec; the dimension of the region, 10 R _⊙; field strength at 10 R _⊙, 0.02 G; particle density, 2.0 × 10⁴/cm³; Alfvén speed, 320 km/sec. From the nature of the observations and the lack of related effects from similar flares on the interplanetary sector pattern observed at 1 AU, it is suggested that such coronal magnetic bottles expand to perhaps 10–30 R _⊙ and then contract to a few solar radii. Such a phenomena is evidence for an expansion of the corona with a sub-Alfvénic velocity. It is further suggested that such magnetic bottles may be important in the storage and diffusion of solar generated cosmic ray particles. NAS-NRC Postdoctoral Resident Research Associate.     

Relativistic model of radiating massive fluid sphere

In this paper we present a new class of nonsingular solutions representing time dependent balls of perfect fluid with matter-radiation in general relativity. The solution of the class is suitable for interior modeling of a quasar i.e. a massive radiating star. The interior solution is matched with a zero pressure Vaidya metric. From this solution we constructed a quasar model by assuming the life time of the quasar of ≈10⁷ year. We obtained a mass of the quasar of ≈10⁹ M _θ , linear dimension ≈10¹⁷ km and a rate of emission L _∞≈10⁴⁷ erg/s.     

Observations of magnetic flux compression in jet impact experiments

The astrophysical jet experiment at Caltech generates a T=2–5 eV, n=10²¹–10²² m⁻³ plasma jet using coplanar disk electrodes linked by a poloidal magnetic field. A 100 kA current generates a toroidal magnetic field; the toroidal field pressure inflates the poloidal flux surface, magnetically driving the jet. The jet travels at up to 50 km/s for ∼20–25 cm before colliding with a cloud of initially neutral gas. We study the interaction of the jet and the cloud in analogy to an astrophysical jet impacting a molecular cloud. Diagnostics include magnetic probe arrays, a 12-channel spectroscopic system and a fast camera with optical filters. When a hydrogen plasma jet collides with an argon target cloud, magnetic measurements show the magnetic flux compressing as the plasma jet deforms. As the plasma jet front slows and the plasma piles up, the density of the frozen-in magnetic flux increases.     

The Variation of Solar Wind Correlation Lengths Over Three Solar Cycles

We present the results of a study of solar wind velocity and magnetic field correlation lengths over the last 35 years. The correlation length of the magnetic field magnitude λ _|B| increases on average by a factor of two at solar maxima compared to solar minima. The correlation lengths of the components of the magnetic field l_BXYZ\lambda_{B_{XYZ}} and of the velocity l_VYZ\lambda_{V_{YZ}} do not show this change and have similar values, indicating a continual turbulent correlation length of around 1.4×10⁶ km. We conclude that a linear relation between λ _|B|, VB ², and Kp suggests that the former is related to the total magnetic energy in the solar wind and an estimate of the average size of geoeffective structures, which is, in turn, proportional to VB ². By looking at the distribution of daily correlation lengths we show that the solar minimum values of λ _|B| correspond to the turbulent outer scale. A tail of larger λ _|B| values is present at solar maximum causing the increase in mean value.     

How to Reach Superequipartition Field Strengths in Solar Magnetic Flux Tubes

A number of independent arguments indicate that the toroidal flux system responsible for the sunspot cycle is stored at the base of the convection zone in the form of flux tubes with field strength close to 10⁵ G. Although the evidence for such strong fields is quite compelling, how such field strength can be reached is still a topic of debate. Flux expulsion by convection should lead to about the equipartition field strength, but the magnetic energy density of a 10⁵-G field is two orders of magnitude larger than the mean kinetic energy density of convective motions. Line stretching by differential rotation (i.e., the “Ω effect” in the classical mean-field dynamo approach) probably plays an important role, but arguments based on energy considerations show that it does not seem feasible that a 10⁵-G field can be produced in this way. An alternative scenario for the intensification of the toroidal flux system in the overshoot layer is related to the explosion of rising, buoyantly unstable magnetic flux tubes, which opens a complementary mechanism for magnetic-field intensification. A parallelism is pointed out with the mechanism of “convective collapse” for the intensification of photospheric magnetic flux tubes up to field strengths well above equipartition; both mechanisms, which are fundamentally thermal processes, are reviewed.