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.     

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     

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.     

Gamma-ray emission expected from Kepler’s SNR

Nonlinear kinetic theory of cosmic ray (CR) acceleration in supernova remnants (SNRs) is used to investigate the properties of Kepler’s SNR and, in particular, to predict the γ-eay spectrum expected from this SNR. Observations of the nonthermal radio and X-ray emission spectra as well as theoretical constraints for the total supernova (SN) explosion energy E _sn are used to constrain the astronomical and particle acceleration parameters of the system. Under the assumption that Kepler’s SN is a type Ia SN we determine for any given explosion energy E _sn and source distance d the mass density of the ambient interstellar medium (ISM) from a fit to the observed SNR size and expansion speed. This makes it possible to make predictions for the expected γ-eay flux. Exploring the expected distance range we find that for a typical explosion energy E _sn=10⁵¹ erg the expected energy flux of TeV γ-rays varies from 2×10⁻¹¹ to 10⁻¹³ erg/(cm² s) when the distance changes from d=3.4 kpc to 7 kpc. In all cases the γ-eay emission is dominated by π ⁰-decay γ-rays due to nuclear CRs. Therefore Kepler’s SNR represents a very promising target for instruments like H.E.S.S., CANGAROO and GLAST. A non-detection of γ-rays would mean that the actual source distance is larger than 7 kpc.     

Upper limits on neutrino fluxes from point-like sources with AMANDA-II

The AMANDA-II telescope, operated by the IceCube collaboration, is currently the world’s most sensitive telescope to fluxes of neutrinos from individual sources. A data sample of 4282 neutrino induced events collected in 1001 days of detector livetime during the years 2000–2004 have now been analyzed looking for a neutrino signal from point-like sources. A sensitivity to fluxes of of d Φ/dE=1.0×10⁻¹⁰(E/TeV)⁻² TeV⁻¹ cm⁻²s⁻¹ was reached in the energy range between 1.7 TeV and 2.4 PeV. So far no statistically significant localized excess of events over the background of atmospheric neutrinos has been found, which would be ascribed to a neutrino source. However, the flux upper limits derived from the non-observation of a signal are comparable to observed fluxes of high energy gamma rays from blazars and within the range of current models for neutrino emission from selected sources. Possible constraints on these models are discussed.      

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.     

Optical imaging and spectroscopic observation of the galactic supernova remnant G85.9-0.6

Optical CCD imaging with Hα and [SII] filters and spectroscopic observations of the galactic supernova remnant G85.9-0.6 have been performed for the first time. The CCD image data are taken with the 1.5 m Russian-Turkish Telescope (RTT150) at TüBİTAK National Observatory (TUG) and spectral data are taken with the Bok 2.3 m telescope on Kitt Peak, AZ. The images are taken with narrow-band interference filters Hα, [SII] and their continuum. [SII]/Hα ratio image is performed. The ratio obtained from [SII]/Hα is found to be ∼0.42, indicating that the remnant interacts with HII regions. G85.9-0.6 shows diffuse-shell morphology. [SII]λ λ6716/6731 average flux ratio is calculated from the spectra, and the electron density N _e is obtained to be 395 cm⁻³. From [OIII]/Hβ ratio, shock velocity has been estimated, pre-shock density of n _c =14 cm⁻³, explosion energy of E=9.2×10⁵⁰ ergs, interstellar extinction of E(B−V)=0.28, and neutral hydrogen column density of N(HI)=1.53×10²¹ cm⁻² are reported.     

Peak-Size Distributions of Proton Fluxes and Associated Soft X-Ray Flares

A database combining information about solar proton enhancements (SPEs) near the Earth and soft X-ray flares (GOES measurements) has been used for the study of different correlations through the period from 1975 to May 2006. The emphasis of this work is on the treatment of peak-size distributions of SXR flares and SPEs. The frequency of SXR flares and solar proton events (>10 and >100 MeV, respectively) for the past three solar cycles has been found to follow mainly a power-law distribution over three to five orders of magnitude of fluxes, which is physically correct beyond the “sensitivity” problem with the smallest peak values. The absence of significant spectral steepening in the domain of the highest peak values demonstrates that during the period considered, lasting 30 years, the limit of the highest flare’s energy release has not yet been achieved. The power-law exponents were found to be −2.19±0.04, −1.34±0.02, and −1.46±0.04, for the total SXR flare distribution and the total SPE distributions (for both E _P>10 MeV and E _P>100 MeV), respectively. For SPEs associated with flares located to the West of 20° W, the exponents are −1.22±0.05 (E _P>10 MeV) and −1.26±0.03 (E _P>100 MeV). The size distribution for corresponding flares follows a power law with a slope of −1.29±0.12. Thus, X-ray and proton fluxes produced in the same solar events have very similar distribution shapes. Moreover, the derived slopes are not incompatible with a linear dependence between X-ray flare power and proton fluxes near the Earth. A similar statistical relation is obtained independently from the direct comparison of the X-ray and proton fluxes. These all argue for a statistically significant relationship between X-ray and proton emissions.     

Stability of strange dwarfs II. Computational results

The stability of strange dwarfs for quark cores with M _0core /M _⨀ = 10⁻⁴, has been studied by calculating, in each individual case, a series of strange dwarfs with configurations in which 5 ⋅ 10⁻⁴, 10⁻³, 5 ⋅ 10⁻³, 10⁻², 1.31 ⋅ 10⁻², 1.6 ⋅ 10⁻², 1.7 ⋅ 10⁻², 2 ⋅ 10⁻², ranges from the values in white dwarfs to ρ _drip = 4.3 ⋅ 10¹¹ g/cm³, at which free neutrons are produced in the crust. For the series with M _0core /M _⨀ < 0.0131, stability is lost when ρ _tr < ρ _drip . For the series with M _0core /M _⨀ > 0.0131, the equality ρ _tr = ρ _drip is reached before the strange dwarf attains its maximum mass. Although the frequency of the radial pulsations in the fundamental mode obeys ω₀² > 0 for these configurations, they are unstable with respect to transitions into a strange star state with the same total number of baryons and a radius on the order of that of neutron stars. An energy on the order of the energy in a supernova explosion is released during these transitions. It is shown that the gravitational red shift of white and strange dwarfs are substantially different for low and limiting (high) masses.     

The physics of E × B-drifting jets

E × B-drifting jets have been generally ignored for the past 25 years even though they may well describe all the astrophysical jet sources, both on galactic and stellar scales. Here we present closed-form solutions for their joint field-and-particle distribution, argue that the observed jets are near equipartition, with extremely relativistic, monoenergetic e^±-pairs of bulk Lorentz factor γ ≲ 10⁴, and are first-order stable. We describe plausible mechanisms for the jets’ (i) formation, (ii) propagation, and (iii) termination. Wherever a beam meets with resistance, its frozen-in Poynting flux transforms the delta-shaped energy distribution of the pairs into an almost white power law,E ² N _E ∼E ^−∫ with ∫ ≳ 0, via single-step falls through the huge convected potential.     

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⁻¹.     

Tritium chain of the hydrogen cycle of thermonuclear reactions in the solar interior

The tritium chain of the hydrogen cycle on the Sun including the reactions ³He(e⁻, ν _e) ³H(p, γ)⁴He is considered. The flux of tritium neutrinos at a distance of 1 AU is 8.1 × 10⁴ cm⁻² s⁻¹. It exceeds the neutrino flux from the (hep)-reaction by one order of magnitude. The radial distribution of the yield of ³H neutrinos inside the Sun and their energy spectrum, which has the form of a line at an energy of 2.5–3.0 keV, have also been calculated. The flux of thermal tritium neutrinos is accompanied by a very weak flux of antineutrinos (∼10³ cm⁻² yr⁻¹) with an energy below 18.6 keV. These antineutrinos are produced in the URCA processes ³He ⇆ ³H.     

New evidence for strong nonthermal effects in Tycho’s supernova remnant

For the case of Tycho’s supernova remnant (SNR) we present the relation between the blast wave and contact discontinuity radii calculated within the nonlinear kinetic theory of cosmic ray (CR) acceleration in SNRs. It is demonstrated that these radii are confirmed by recently published Chandra measurements which show that the observed contact discontinuity radius is so close to the shock radius that it can only be explained by efficient CR acceleration which in turn makes the medium more compressible. Together with the recently determined new value E _sn=1.2×10⁵¹ erg of the SN explosion energy this also confirms our previous conclusion that a TeV γ-ray flux of (2–5)×10⁻¹³ erg/(cm² s) is to be expected from Tycho’s SNR. Chandra measurements and the HEGRA upper limit of the TeV γ-ray flux together limit the source distance d to 3.3≤d≤4 kpc.     

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.     

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.      

Structure of Radiative Shock Waves in Stellar Atmospheres

A global iteration method to determine the self-consistent structure of steady plane-parallel radiative shock waves is shown to converge to the stable solution with upstream front velocities of 15 km/s ≤ U ₁≤ 60 km/s and for hydrogen gas of unperturbed temperature T= 3000 K and density ρ = 10⁻¹⁰gcm⁻³. This revised version was published online in July 2006 with corrections to the Cover Date.     

Acoustic Events in the Solar Atmosphere from <Emphasis Type="Italic">Hinode</Emphasis>/SOT NFI Observations

We investigate the properties of acoustic events (AEs), defined as spatially concentrated and short duration energy flux, in the quiet Sun, using observations of a 2D field of view (FOV) with high spatial and temporal resolution provided by the Solar Optical Telescope (SOT) onboard Hinode. Line profiles of Fe i 557.6 nm were recorded by the Narrow-band Filter Imager (NFI) on a 82″×82″ FOV during 75 min with a time step of 28.75 s and 0.08″ pixel size. Vertical velocities were computed at three atmospheric levels (80, 130, and 180 km) using the bisector technique, allowing the determination of energy flux to be made in the range 3 – 10 mHz using two complementary methods (Hilbert transform and Fourier power spectrum). Horizontal velocities were computed using local correlation tracking (LCT) of continuum intensities providing divergences. We found that the net energy flux is upward. In the range 3 – 10 mHz, a full FOV space and time averaged flux of 2700 W m⁻² (lower layer 80 – 130 km) and 2000 W m⁻² (upper layer 130 – 180 km) is concentrated in less than 1 % of the solar surface in the form of narrow (0.3″) AE. Their total duration (including rise and decay) is of the order of 10³ s. Inside each AE, the mean flux is 1.6×10⁵ W m⁻² (lower layer) and 1.2×10⁵ W m⁻² (upper). Each event carries an average energy (flux integrated over space and time) of 2.5×10¹⁹ J (lower layer) to 1.9×10¹⁹ J (upper). More than 10⁶ events could exist permanently on the Sun, with a birth and decay rate of 3500 s⁻¹. Most events occur in intergranular lanes, downward velocity regions, and areas of converging motions.     

The Radio Recombination Line Intensities Of Atoms And Ions: Helium,Carbon And Oxygen

The radio recombination line intensities of heavy elements of helium, carbon and oxygen are calculated with accounting for dielectronic recombination. Dielectronic recombination rates are determined accurate to the second order of a perturbation theory and the rates are described as function of principal quantum number for helium-like atom or ion. Balance equations are solved for the departure coefficients from LTE b_n. The collision and spontaneous transition rates are accounted for the balance equations, in which non-equilibrium distribution source is dielectronic recombination. Non-equilibrium amplification coefficients are found as functions of a medium temperature, density and ion charge z = 1–3 for radio recombination lines. Optical depths are calculated for the heavy element low-frequency lines with the numbers 300 > n > 1200. For the chosen electronic temperatures and densities T_e = 0.8× 10⁴–10× 10⁴ K, N_e = 0.05–0.1 cm⁻³ the line optical depth is determined by the values τ_L∼ 0.1× 10⁻⁴–100× 10⁻⁴. Calculated for free-free transition rates, the optical depth is given by using the value τ_ff∼ 10⁻²τ_L.     

Relativistic and non-relativistic dust grains

Charged dust grains of radiia≃3×10⁻⁶∼3×10⁻⁵ cm could be a help in understanding high energy particles in extensive air showers (EAS). It is suggested that the dust grains in the intergalactic medium may attain relatistic energy (≥10²⁰ eV), and may be responsible for the apparent ‘bump’ in the energy spectrum. The relativistic and non-relativistic dust grains may help to understand the anomalies in the energy spectrum as regards the slope and intensity.     

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.