Diagnostics of a flare on EQ Peg B from optical pulsations

Based on the methods of coronal seismology, we have investigated the ten-second quasi-periodic pulsations of the optical flare emission from the active red dwarf EQ Peg B detected with the William Herschel Telescope on La Palma. We propose and analyze a model in which they could be produced by sausage oscillations of a coronal flare loop. The amplitude and phase relations between the displacement components of the radial oscillations and the conditions for their excitation in loops with footpoints frozen into the photosphere are considered. The temperature (≈6 × 10⁷ K), plasma density (≈2.7 × 10¹¹ cm⁻³), and magnetic field strength (≈540 G) in the region of energy release have been determined. Our estimate of the flare loop length (≈0.4R _⋆) provides evidence for the existence of extended coronae on red dwarf stars.     

MMIV: de SGR 1806–20 Anno Mirabili

On 27th December 2004 SGR 1806–20, one of the most active Soft γ-ray Repeaters (SGRs), displayed an extremely rare event, also known as giant flare, during which up to 10⁴⁷ ergs were released in the ∼1–1000 keV range in less than 1 s. Before and after the giant flare we carried out IR observations by using adaptive optics (NAOS-CONICA) mounted on VLT which provided images of unprecedented quality (FWHM better than 0.1″). We discovered the likely IR counterpart to SGR 1806–20 based on positional coincidence with the VLA uncertainty region and flux variability of a factor of about 2 correlated with that at higher energies. Moreover, by analysing the Rossi-XTE/PCA data we have discovered rapid Quasi-Periodic Oscillations (QPOs) in the pulsating tail of the 27th December 2004 giant flare of SGR 1806–20. QPOs at ∼92.5 Hz are detected in a 50 s interval starting 170 s after the onset of the giant flare. These QPOs appear to be associated with increased emission by a relatively hard unpulsed component and are seen only over phases of the 7.56 s spin period pulsations away from the main peak. QPOs at ∼18 and ∼30 Hz are also detected ∼200–300 s after the onset of the giant flare. This is the first time that QPOs are unambiguously detected in the flux of a Soft Gamma-ray Repeater, or any other isolated neutron star. We interpret the highest QPOs in terms of the coupling of toroidal seismic modes with Alfvén waves propagating along magnetospheric field lines. The lowest frequency QPO might instead provide indirect evidence on the strength of the internal magnetic field of the neutron star.      

WR 140 (=V1687 Cyg): Infrared photometry, 2001–2010

We present the results of our infrared observations of WR 140 (=V1687 Cyg) in 2001–2010. Analysis of the observations has shown that the J brightness at maximum increased near the periastron by about 0 ^m .3; the M brightness increased by ∼2 ^m in less than 50 days. The minimum J brightness and the minimum L and M brightnesses were observed 550–600 and 1300–1400 days after the maximum, respectively. The JHKLM brightness minimum was observed in the range of orbital phases 0.7–0.9. The parameters of the primary O5 component of the binary have been estimated to be the following: R(O5) ≈ 24.7R _⊙, L(O5) ≈ 8 × 10⁵ L _⊙, and M _bol(O5) ≈ −10 ^m . At the infrared brightness minimum, T _g ∼ 820–880 K, R _g ≈ 2.6 × 10⁵ R _⊙, the optical depth of the shell at 3.5 μm is ∼5.3 × 10⁻⁶, and its mass is ≈1.4 × 10⁻⁸ M _⊙. At the maximum, the corresponding parameters are ∼1300 K, 8.6 × 10⁴ R _⊙, ∼2 × 10⁻⁴, and ∼6 × 10⁻⁸ M _⊙; the mean rate of dust inflow (condensation) into the dust structure is ∼3.3 × 10⁻⁸ M _⊙ yr⁻¹. The mean escape velocity of the shell from the heating source is ∼10³ km s⁻¹ and the mean dispersal rate of the shell is ∼1.1 × 10⁻⁸ M _⊙ yr⁻¹.     

Chromospheric Height and Density Measurements in a Solar Flare Observed with RHESSI – II. Data Analysis

We present an analysis of hard X-ray imaging observations from one of the first solar flares observed with the Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) spacecraft, launched on 5 February 2002. The data were obtained from the 22 February 2002, 11:06 UT flare, which occurred close to the northwest limb. Thanks to the high energy resolution of the germanium-cooled hard X-ray detectors on RHESSI we can measure the flare source positions with a high accuracy as a function of energy. Using a forward-fitting algorithm for image reconstruction, we find a systematic decrease in the altitudes of the source centroids z(ε) as a function of increasing hard X-ray energy ε, as expected in the thick-target bremsstrahlung model of Brown. The altitude of hard X-ray emission as a function of photon energy ε can be characterized by a power-law function in the ε=15–50 keV energy range, viz., z(ε)≈2.3(ε/20 keV)^−1.3 Mm. Based on a purely collisional 1-D thick-target model, this height dependence can be inverted into a chromospheric density model n(z), as derived in Paper I, which follows the power-law function n _e(z)=1.25×10¹³(z/1 Mm)^−2.5 cm⁻³. This density is comparable with models based on optical/UV spectrometry in the chromospheric height range of h≲1000 km, suggesting that the collisional thick-target model is a reasonable first approximation to hard X-ray footpoint sources. At h≈1000–2500 km, the hard X-ray based density model, however, is more consistent with the `spicular extended-chromosphere model' inferred from radio sub-mm observations, than with standard models based on hydrostatic equilibrium. At coronal heights, h≈2.5–12.4 Mm, the average flare loop density inferred from RHESSI is comparable with values from hydrodynamic simulations of flare chromospheric evaporation, soft X-ray, and radio-based measurements, but below the upper limits set by filling-factor insensitive iron line pairs.     

On origin of ultra high energy cosmic rays

Propagation of UHE protons through CMB radiation leaves the imprint on energy spectrum in the form of Greisen–Zatsepin–Kuzmin (GZK) cutoff, bump (pile-up protons) and dip. The dip is a feature in energy range 1×10¹⁸–4×10¹⁹ eV, caused by electron-positron pair production on CMB photons. Calculated for power-law generation spectrum with index γ _g =2.7, the shape of the dip is confirmed with high accuracy by data of Akeno—AGASA, HiRes, Yakutsk and Fly’s Eye detectors. The predicted shape of the dip is robust: it is valid for the rectilinear and diffusive propagation, for different discretenesses in the source distribution, for local source overdensity and deficit etc. This property of the dip allows us to use it for energy calibration of the detectors. The energy shift λ for each detector is determined by minimum χ ² in comparison of observed and calculated dip. After this energy calibration the absolute fluxes, measured by AGASA, HiRes and Yakutsk detectors remarkably coincide in energy region 1×10¹⁸–1×10²⁰ eV. Below the characteristic energy E _c ≈1×10¹⁸ eV the spectrum of the dip flattens for both diffusive and rectilinear propagation, and more steep galactic spectrum becomes dominant at E<E _c . The energy of transition E _tr<E _c approximately coincides with the position of the second knee E _2kn , observed in the cosmic ray spectrum. The dip-induced transition from galactic to extragalactic cosmic rays at the second knee is compared with traditional model of transition at ankle, the feature observed at energy ∼1×10¹⁹ eV.      

The Evaporation Regime in a Confined Flare

We studied the evolution of a small eruptive flare (GOES class C1) from its onset phase using multi-wavelength observations that sample the flare atmosphere from the chromosphere to the corona. The main instruments involved were the Coronal Diagnostic Spectrometer (CDS) aboard SOHO and facilities at the Dunn Solar Tower of the National Solar Observatory/Sacramento Peak. Transition Region and Coronal Explorer (TRACE) together with Ramaty High-Energy Spectroscopic Imager (RHESSI) also provided images and spectra for this flare. Hα and TRACE images display two loop systems that outline the pre-reconnection and post-reconnection magnetic field lines and their topological changes revealing that we are dealing with an eruptive confined flare. RHESSI data do not record any detectable emission at energies ≥25 keV, and the observed count spectrum can be well fitted with a thermal plus a non-thermal model of the photon spectrum. A non-thermal electron flux F ≈ 5 × 10¹⁰ erg cm⁻² s⁻¹ is determined. The reconstructed images show a very compact source whose peak emission moves along the photospheric magnetic inversion line during the flare. This is probably related to the motion of the reconnection site, hinting at an arcade of small loops that brightens successively. The analysis of the chromospheric spectra (Ca II K, He I D₃ and Hγ, acquired with a four-second temporal cadence) shows the presence of a downward velocity (between 10 and 20 km s⁻¹) in a small region intersected by the spectrograph slit. The region is included in an area that, at the time of the maximum X-ray emission, shows upward motions at transition region (TR) and coronal levels. For the He I 58.4 and O v 62.97 lines, we determine a velocity of ≈−40 km s⁻¹ while for the Fe XIX 59.22 line a velocity of ≈−80 km s⁻¹ is determined with a two-component fitting. The observations are discussed in the framework of available hydrodynamic simulations and they are consistent with the scenario outlined by Fisher (1989). No explosive evaporation is expected for a non-thermal electron beam of the observed characteristics, and no gentle evaporation is allowed without upward chromospheric motion. It is suggested that the energy of non-thermal electrons can be dissipated to heat the high-density plasma, where possibly the reconnection occurs. The consequent conductive flux drives the evaporation process in a regime that we can call sub-explosive.     

Solar Flares and Coronal Mass Ejections: A Statistically Determined Flare Flux – CME Mass Correlation

In an effort to examine the relationship between flare flux and corresponding CME mass, we temporally and spatially correlate all X-ray flares and CMEs in the LASCO and GOES archives from 1996 to 2006. We cross-reference 6733 CMEs having well-measured masses against 12 050 X-ray flares having position information as determined from their optical counterparts. For a given flare, we search in time for CMEs which occur 10 – 80 minutes afterward, and we further require the flare and CME to occur within ± 45° in position angle on the solar disk. There are 826 CME/flare pairs which fit these criteria. Comparing the flare fluxes with CME masses of these paired events, we find CME mass increases with flare flux, following an approximately log-linear, broken relationship: in the limit of lower flare fluxes, log (CME mass)∝0.68×log (flare flux), and in the limit of higher flare fluxes, log (CME mass)∝0.33×log (flare flux). We show that this broken power-law, and in particular the flatter slope at higher flare fluxes, may be due to an observational bias against CMEs associated with the most energetic flares: halo CMEs. Correcting for this bias yields a single power-law relationship of the form log (CME mass)∝0.70×log (flare flux). This function describes the relationship between CME mass and flare flux over at least 3 dex in flare flux, from ≈ 10⁻⁷ – 10⁻⁴ W m⁻².     

Energy-Dependent Timing of Thermal Emission in Solar Flares

We report solar flare plasma to be multi-thermal in nature based on the theoretical model and study of the energy-dependent timing of thermal emission in ten M-class flares. We employ high-resolution X-ray spectra observed by the Si detector of the “Solar X-ray Spectrometer” (SOXS). The SOXS onboard the Indian GSAT-2 spacecraft was launched by the GSLV-D2 rocket on 8 May 2003. Firstly we model the spectral evolution of the X-ray line and continuum emission flux F(ε) from the flare by integrating a series of isothermal plasma flux. We find that the multi-temperature integrated flux F(ε) is a power-law function of ε with a spectral index (γ)≈−4.65. Next, based on spectral-temporal evolution of the flares we find that the emission in the energy range E=4 – 15 keV is dominated by temperatures of T=12 – 50 MK, while the multi-thermal power-law DEM index (δ) varies in the range of −4.4 and −5.7. The temporal evolution of the X-ray flux F(ε,t) assuming a multi-temperature plasma governed by thermal conduction cooling reveals that the temperature-dependent cooling time varies between 296 and 4640 s and the electron density (n _e) varies in the range of n _e=(1.77 – 29.3)×10¹⁰ cm⁻³. Employing temporal evolution technique in the current study as an alternative method for separating thermal from nonthermal components in the energy spectra, we measure the break-energy point, ranging between 14 and 21±1.0 keV.     

Flare Plasma Cooling from 30 MK down to 1 MK modeled from Yohkoh,GOES, and TRACE observations during the Bastille Day Event (14 July 2000)

We present an analysis of the evolution of the thermal flare plasma during the 14 July 2000, 10 UT, Bastille Day flare event, using spacecraft data from Yohkoh/HXT, Yohkoh/SXT, GOES, and TRACE. The spatial structure of this double-ribbon flare consists of a curved arcade with some 100 post-flare loops which brighten up in a sequential manner from highly-sheared low-lying to less-sheared higher-lying bipolar loops. We reconstruct an instrument-combined, average differential emission measure distribution dEM(T)/dT that ranges from T=1 MK to 40 MK and peaks at T ₀=10.9 MK. We find that the time profiles of the different instrument fluxes peak sequentially over 7 minutes with decreasing temperatures from T≈30 MK to 1 MK, indicating the systematic cooling of the flare plasma. From these temperature-dependent relative peak times t _peak(T) we reconstruct the average plasma cooling function T(t) for loops observed near the flare peak time, and find that their temperature decrease is initially controlled by conductive cooling during the first 188 s, T(t)∼[1+(t/τ_cond)]^−2/7, and then by radiative cooling during the next 592 s, T(t)∼[1−(t/τ_rad)]^3/5. From the radiative cooling phase we infer an average electron density of n _e=4.2×10¹¹ cm⁻³, which implies a filling factor near 100% for the brightest observed 23 loops with diameters of ∼1.8 Mm that appear simultaneously over the flare peak time and are fully resolved with TRACE. We reproduce the time delays and fluxes of the observed time profiles near the flare peak self-consistently with a forward-fitting method of a fully analytical model. The total integrated thermal energy of this flare amounts to E _thermal=2.6×10³¹ erg. Supplementary material to this paper is available in electronic form at http://dx.doi.org/10.1023/A:1014257826116     

Well-behaved relativistic charged super-dense star models

A new class of charged super-dense star models is obtained by using an electric intensity, which involves a parameter, K. The metric describing the model shares its metric potential g ₄₄ with that of Durgapal’s fourth solution (J. Phys. A, Math. Gen. 15:2637, 1982). The pressure-free surface is kept at the density ρ _b =2×10¹⁴ g/cm³ and joins smoothly with the Reissner-Nordstrom solution. The charge analogues are well-behaved for a wide range, 0≤K≤59, with the optimum value of X=0.264 i.e. the pressure, density, pressure–density ratio and velocity of sound are monotonically decreasing and the electric intensity is monotonically increasing in nature for the given range of the parameter K. The maximum mass and the corresponding radius occupied by the neutral solution are 4.22M _Θ and 20 km, respectively for X=0.264. For the charged solution, the maximum mass and radius are defined by the expressions M≈(0.0059K+4.22)M _Θ and r _b ≈−0.021464K+20 km respectively.     

Shrinking and Cooling of Flare Loops in a Two-Ribbon Flare

We analyze the evolution of the flare/postflare-loop system in the two-ribbon flare of November 3, 2003, utilizing multi-wavelength observations that cover the temperature range from several tens of MK down to 10⁴ K. A non-uniform growth of the loop system enables us to identify analogous patterns in the height–time, h(t), curves measured at different temperatures. The “knees,” “plateaus,” and “bends” in a higher-temperature curve appear after a certain time delay at lower heights in a lower-temperature curve. We interpret such a shifted replication as a track of a given set of loops (reconnected field lines) while shrinking and cooling after being released from the reconnection site. Measurements of the height/time shifts between h(t) curves of different temperatures provide a simultaneous estimate of the shrinkage speed and cooling rate in a given temperature domain, for a period of almost ten hours after the flare impulsive phase. From the analysis we find the following: (a) Loop shrinkage is faster at higher temperatures – in the first hour of the loop-system growth, the shrinkage velocity at 5 MK is 20 – 30 km s⁻¹, whereas at 1 MK it amounts to 5 km s⁻¹; (b) Shrinking becomes slower as the flare decays – ten hours after the impulsive phase, the shrinkage velocity at 5 MK becomes 5 km s⁻¹; (c) The cooling rate decreases as the flare decays – in the 5 MK range it is 1 MK min⁻¹ in the first hour of the loop-system growth, whereas ten hours later it decreases to 0.2 MK min⁻¹; (d) During the initial phase of the loop-system growth, the cooling rate is larger at higher temperatures, whereas in the late phases the cooling rate apparently does not depend on the temperature; (e) A more detailed analysis of shrinking/cooling around one hour after the impulsive phase reveals a deceleration of the loop shrinkage, amounting to ā ≈ 10 m s⁻² in the T < 5 MK range; (f) In the same interval, conductive cooling dominates down to T ≈ 3 MK, whereas radiation becomes dominant below T ≈ 2 MK; (g) A few hours after the impulsive phase, radiation becomes dominant across the whole T < 5 MK range. These findings are compared with results of previous studies and discussed in the framework of relevant models.     

Jupiter's Synchrotron Emission Induced by the Collision of Comet Shoemaker-Levy 9

From July 13 to August 21, 1994, we observed Jupiter at 1420 MHz using one of the 30-m single dishes of the Instituto Argentino de Radioastronomía. After the impact of fragment G, we detected a rapid increase of the 21cm-continuum flux, which reached the maximum (≈ 20% of Jupiter's flux) at the end of the impact period. The nature of this radiation is clearly synchrotron. We interpret it in terms of a new population of relativistic electrons (≈ 2 × 10²⁹) injected into the Jovian magnetosphere as a consequence of the impact explosions. The proposed mechanism is that the relativistic plasma was blown as magnetic clouds that flowed along the magnetic lines of force towards the jovimagnetic equator. We constructed a model in which the energies of the fresh electrons, generated within the magnetized clouds with a power law energy spectrum, were highly degraded by the comet dust grains attached to the magnetized plasma. The model can account for the spectral shape based on observations at several frequencies (de Pater et al., 1995, Science 268, 1879; Venturi et al., 1996, Astron. Astrophys. 316, 243). The energy released by the explosions under the form of relativistic electrons is of ≈ 2 × 10²⁵ erg, which represents a fraction of about 1–3 per cent of the explosion energy. The efficiency in converting the explosion energy into the relativistic electron energy is, therefore, of the same order of magnitude as that of supernova explosions. An alternative model is considered. This gives figures for the total energy and number of relativistic electrons that are similar to the corresponding ones of the favoured model. Finally, we suggest that the behavior of the flux decay in the various observed frequencies is the result of the diffusion of electrons into the loss-cone due to the resonant scattering of the electrons by Alfven waves. This revised version was published online in July 2006 with corrections to the Cover Date.     

Evaluation of a Selected Case of the Minimum Dissipative Rate Method for Non-Force-Free Solar Magnetic Field Extrapolations

The minimum dissipative rate (MDR) method for deriving a coronal non-force-free magnetic field solution is partially evaluated. These magnetic field solutions employ a combination of three linear (constant-α) force-free-field solutions with one being a potential field (i.e., α=0). The particular case of the solutions where the other two α’s are of equal magnitude but of opposite sign is examined. This is motivated by studying the SOLIS (Synoptic Optical Long-term Investigation of the Sun (SOLIS), a National Solar Observatory facility) vector magnetograms of AR 10987, which show a global α value consistent with an α=0 value as evaluated by (∇×B) _z /B _z over the region. Typical of the current state of the observing technology, there is no definitive twist for input into the general MDR method. This suggests that the special α case, of two α’s with equal magnitudes and opposite signs, is appropriate given the data. Only for an extensively twisted active region does a dominant, nonzero α normally emerge from a distribution of local values. For a special set of conditions, is it found that (i) the resulting magnetic field is a vertically inflated magnetic field resulting from the electric currents being parallel to the photosphere, similar to the results of Gary and Alexander (Solar Phys. 186:123, 1999), and (ii) for α≈(α _max/2), the Lorentz force per unit volume normalized by the square of the magnetic field is on the order of 1.4×10⁻¹⁰ cm⁻¹. The Lorentz force (F _L) is a factor of ten higher than that of the magnetic force d(B ²/8π)/dz, a component of F _L. The calculated photospheric electric current densities are an order of magnitude smaller than the maximum observed in all active regions. Hence both the Lorentz force density and the generated electric current density seem to be physically consistent with possible solar dynamics. The results imply that the field could be inflated with an overpressure along the neutral line. However, the implementation of this or any other extrapolation method using the electric current density as a lower boundary condition must be done cautiously, with the current magnetography.     

Vertical structure of the outer accretion disk in persistent low-mass X-ray binaries

We have investigated the influence of X-ray irradiation on the vertical structure of the outer accretion disk in low-mass X-ray binaries by performing a self-consistent calculation of the vertical structure and X-ray radiation transfer in the disk. Penetrating deep into the disk, the field of scattered X-ray photons with energy E ≳ 10 keV exerts a significant influence on the vertical structure of the accretion disk at a distance R ≳ 10¹⁰ cm from the neutron star. At a distance R ∼ 10¹¹ cm, where the total surface density in the disk reaches Σ₀ ∼ 20 g cm⁻², X-ray heating affects all layers of an optically thick disk. The X-ray heating effect is enhanced significantly in the presence of an extended atmospheric layer with a temperature T _atm ≈ (2–3) × 10⁶ K above the accretion disk. We have derived simple analytic formulas for the disk heating by scattered X-ray photons using an approximate solution of the transfer equation by the Sobolev method. This approximation has a ≲10% accuracy in the range of X-ray photon energies E < 20 keV.     

Multiwavelength Observations of a Failed Flux Rope in the Eruption and Associated M-Class Flare from NOAA AR 11045

We present the multiwavelength observations of a flux rope that was trying to erupt from NOAA AR 11045 and the associated M-class solar flare on 12 February 2010 using space-based and ground-based observations from TRACE, STEREO, SOHO/MDI, Hinode/XRT, and BBSO. While the flux rope was rising from the active region, an M1.1/2F class flare was triggered near one of its footpoints. We suggest that the flare triggering was due to the reconnection of a rising flux rope with the surrounding low-lying magnetic loops. The flux rope reached a projected height of ≈0.15R _⊙ with a speed of ≈90 km s⁻¹ while the soft X-ray flux enhanced gradually during its rise. The flux rope was suppressed by an overlying field, and the filled plasma moved towards the negative polarity field to the west of its activation site. We found the first observational evidence of the initial suppression of a flux rope due to a remnant filament visible both at chromospheric and coronal temperatures that evolved a couple of days earlier at the same location in the active region. SOHO/MDI magnetograms show the emergence of a bipole ≈12 h prior to the flare initiation. The emerged negative polarity moved towards the flux rope activation site, and flare triggering near the photospheric polarity inversion line (PIL) took place. The motion of the negative polarity region towards the PIL helped in the build-up of magnetic energy at the flare and flux rope activation site. This study provides unique observational evidence of a rising flux rope that failed to erupt due to a remnant filament and overlying magnetic field, as well as associated triggering of an M-class flare.     

Do EUV Nanoflares Account for Coronal Heating?

Recent observations with EUV imaging instruments such as SOHO/EIT and TRACE have shown evidence for flare-like processes at the bottom end of the energy scale, in the range of E _th≈10²⁴–10²⁷ erg. Here we compare these EUV nanoflares with soft X-ray microflares and hard X-ray flares across the entire energy range. From the observations we establish empirical scaling laws for the flare loop length, L(T)∼T, the electron density, n _e(T)∼T ², from which we derive scaling laws for the loop pressure, p(T)∼T ³, and the thermal energy, E _th∼T ⁶. Extrapolating these scaling laws into the picoflare regime we find that the pressure conditions in the chromosphere constrain a height level for flare loop footpoints, which scales with h _eq(T)∼T ^−0.5. Based on this chromospheric pressure limit we predict a lower cutoff of flare loop sizes at L _∖min≲5 Mm and flare energies E _∖min≲10²⁴ erg. We show evidence for such a rollover in the flare energy size distribution from recent TRACE EUV data. Based on this energy cutoff imposed by the chromospheric boundary condition we find that the energy content of the heated plasma observed in EUV, SXR, and HXR flares is insufficient (by 2–3 orders of magnitude) to account for coronal heating.     

Investigation of Quasi-periodic Variations in Hard X-Rays of Solar Flares. II. Further Investigation of Oscillating Magnetic Traps

In our recent paper (Jakimiec and Tomczak, Solar Physics 261, 233, 2010) we investigated quasi-periodic oscillations of hard X-rays during the impulsive phase of solar flares. We have come to the conclusion that they are caused by magnetosonic oscillations of magnetic traps within the volume of hard-X-ray (HXR) loop-top sources. In the present paper we investigate four flares that show clear quasi-periodic sequences of the HXR pulses. We also describe our phenomenological model of oscillating magnetic traps to show that it can explain the observed properties of the HXR oscillations. The main results are the following: i) Low-amplitude quasi-periodic oscillations occur before the impulsive phase of some flares. ii) The quasi-periodicity of the oscillations can change in some flares. We interpret this as being due to changes of the length of oscillating magnetic traps. iii) During the impulsive phase a significant part of the energy of accelerated (non-thermal) electrons is deposited within a HXR loop-top source. iv) The quick development of the impulsive phase is due to feedback between the pressure pulses by accelerated electrons and the amplitude of the magnetic-trap oscillation. v) The electron number density and magnetic field strength values obtained for the HXR loop-top sources in several flares fall within the limits of N≈(2 – 15)×10¹⁰ cm⁻³, B≈(45 – 130) gauss. These results show that the HXR quasi-periodic oscillations contain important information about the energy release in solar flares.     

The Magnetic Field, Horizontal Motion and Helicity in a Fast Emerging Flux Region which Eventually forms a Delta Spot

We studied the behavior of magnetic field, horizontal motion and helicity in a fast emerging flux region NOAA 10488 which eventually forms a δ spot. It is found that the rotation of photospheric footpoints forms in the earlier stage of magnetic flux emergence and the relative shear motion of different magnetic flux systems appears later in this active region (AR). Therefore the emerging process of the AR can be separated into two phases: rotation and shear. We have computed the magnetic helicity injected into the corona using the local correlation tracking (LCT) technique. Furthermore we determined the vertical component of current helicity density and the vertical component of induction electric fields E_z = (V× B)_z in the photosphere. Particularly we have presented the comparison of the injection rate of magnetic helicity and the variation of the current helicity density. The main results are as follows: (1) The strong shear motion (SSM) between the new emerging flux system and the old one brings more magnetic helicity into the corona than the twisting motions. (2) After the maturity of the main bipolar spots, their twist decreases and the SSM becomes dominant and the major contributor of magnetic non-potentiality in the solar atmosphere in this AR. (3) The positions of the maxima of E_z (about 0.1 ∼ 0.2 V cm⁻¹) shift from the twisting areas to the areas showing SSMs as the AR evolved from the rotation phase to the shear one, but no obvious correlation is found between the kernels of Hα flare and E_z for the M1.6 flare in this AR. (4) The coronal helicity inferred from the horizontal motion of this AR amounts to −6 × 10⁴³ Mx². It is comparable with the coronal helicity of ARs producing flares with coronal mass ejections (CMEs) or helicity carried away by magnetic clouds (MCs) reported in previous studies (Nindos, Zhang, and Zhang, 2003; Nindos and Andrews, 2004). In addition, the formation of the δ configuration in this AR belongs to the third formation type indicated by Zirin and Liggett (1987), i.e., collision of opposite polarities from different dipoles, and can be naturally explained by the SSM.     

Variable synchrotron emission from BL Lacertae objects. II. Optical and X-ray flares in HBL source 1ES 1959+650

This paper presents the results of the optical R band and 1.5–12 keV band X-ray monitoring of the high-energy peaked BL Lacertae source 1ES 1959+650 performed during 2002–2007 with the 70 cm Meniscus Telescope of Abastumani Astrophysical Observatory (Georgia) and the All-Sky Monitor on board the Rossi X-ray Time Explorer, respectively. The observed long- and short-term outbursts are fitted with the lightcurves obtained by means of the modeling of synchrotron flares that are assumed to be the result of a propagation of the relativistic shock waves through the jet of 1ES 1959+650, pointed to the observer. Different values of the input parameters (shock velocity, particles’ spectral index, sizes of emission region, minimum and maximum Lorentz factors of the particles etc.) are used in order to fit the simulated lightcurves whose constructed by means of observational data. This investigation shows that both shock velocity and physical conditions in the jet of 1ES 1959+650 should be variable from flare to flare. The shocks are found to be mildly relativistic with the apparent speeds β=0.46–0.85, expressed in the units of c. Spectral index of the particle energy distribution varied from 2.10 to 2.17 for the long-term flares while it is higher in the case of short-term outbursts: s=2.32–2.45 that is suggested to be a result of the deceleration of shock front during its passage through the shell situated downstream the Mach disc. The average strength of a turbulent magnetic field ranged from 0.025 gauss to 0.04 gauss for different long-term flares while the values of 0.07–0.14 gauss were adopted for the different short-term outbursts. The lengths of variable jet area found to be of 0.13–0.47 pc with the transverse extents of (0.5–1.0)×10¹⁷ cm in the case of long-term flares. The same characteristics for short-term outbursts were (2.74–5.5)×10¹⁶ cm and (0.2–04)×10¹⁷ cm, respectively. We conclude that both shock velocity and properties of pre-shocked plasma were not the same in 1ES 1959+650 for the different flaring epochs.     

Variety of well behaved exact solutions of Einstein–Maxwell field equations: an application to Strange Quark stars, Neutron stars and Pulsars

We present a variety of well behaved classes of Charge Analogues of Tolman’s iv (1939). These solutions describe charged fluid balls with positively finite central pressure, positively finite central density; their ratio is less than one and causality condition is obeyed at the centre. The outmarch of pressure, density, pressure-density ratio and the adiabatic speed of sound is monotonically decreasing, however, the electric intensity is monotonically increasing in nature. These solutions give us wide range of parameter for every positive value of n for which the solution is well behaved hence, suitable for modeling of super dense stars. keeping in view of well behaved nature of these solutions, one new class of solutions is being studied extensively. Moreover, this class of solutions gives us wide range of constant K (0.3≤K≤0.91) for which the solution is well behaved hence, suitable for modeling of super dense stars like Strange Quark stars, Neutron stars and Pulsars. For this class of solutions the mass of a star is maximized with all degree of suitability, compatible with Quark stars, Neutron stars and Pulsars. By assuming the surface density ρ _b =2×10¹⁴ g/cm³ (like, Brecher and Caporaso in Nature 259:377, 1976), corresponding to K=0.30 with X=0.39, the resulting well behaved model has the mass M=2.12M _Θ, radius r _b ≈15.27 km and moment of inertia I=4.482×10⁴⁵ g cm²; for K=0.4 with X=0.31, the resulting well behaved model has the mass M=1.80M _Θ, radius r _b ≈14.65 km and moment of inertia I=3.454×10⁴⁵ g cm²; and corresponding to K=0.91 with X=0.135, the resulting well behaved model has the mass M=0.83M _Θ, radius r _b ≈11.84 km and moment of inertia I=0.991×10⁴⁵ g cm². For n=0 we rediscovered Pant et al. (in Astrophys. Space Sci. 333:161, 2011b) well behaved solution. These values of masses and moment of inertia are found to be consistent with other models of Neutron stars and Pulsars available in the literature and are applicable for the Crab and the Vela Pulsars.