Hard X-ray footpoint motions in solar flares: Comparing magnetic reconnection models with observations

Hard X-ray observations from the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) of the October 29, 2003 GOES X10 two-ribbon flare are used together with magnetic field observations from the Michelson Doppler Imager (MDI) onboard SoHO to compare footpoint motions with predictions from magnetic reconnection models. The temporal variations of the velocity v of the hard X-ray footpoint motions and the photospheric magnetic field strength B in footpoints are investigated. The underlying photospheric magnetic field strength is generally higher (B ∼ 700–1200 G) in the slower moving (v ∼ 20–50 km s⁻¹) western footpoint than in the faster (v ∼ 20–100 km s⁻¹) moving eastern source (∼100–600 G). Furthermore, a rough temporal correlation between the HXR flux and the product vB² is observed.     

Creating an artificial Geminid meteor shower: Correlation between ejecta velocity and observability

One of the interesting arguments for a space impact mission to asteroid 3200 Phaethon is to create an artificial Geminid meteor shower. In this work we investigate the artificial shower’s dates of observability and dependence on ejecta velocity using dust trail theory. We find that when the dust ejecta velocities are 200 m/s the artificial meteor showers start to be visible in 2204 and continue for about 30 years. If the dust ejecta velocity is 20 m/s they only last 10 years from 2215 to 2225. Thus, the onset of artificial shower activity begins sooner and lasts longer with higher ejecta velocities. To produce an artificial meteor shower with 3200 Phaethon as the parent will require higher impact energy than the Deep Impact spacecraft delivered to 9P/Tempel 1.     

Gamma-ray observations of explosive nucleosynthesis products

In this review I discuss the various γ-ray emission lines that can be expected and, in some cases have been observed, from radioactive explosive nucleosynthesis products. The most important γ-ray lines result from the decay chains of ⁵⁶Ni, ⁵⁷Ni, and ⁴⁴Ti. ⁵⁶Ni is the prime explosive nucleosynthesis product of Type Ia supernovae, and its decay determines to a large extent the Type Ia light curves. ⁵⁶Ni is also a product of core-collapse supernovae, and in fact, γ-ray line emission from its daughter product, ⁵⁶Co, has been detected from SN1987A by several instruments. The early occurrence of this emission was surprising and indicates that some fraction of ⁵⁶Ni, which is synthesized in the innermost supernova layers, must have mixed with the outermost supernova ejecta.Special attention is given to the γ-ray line emission of the decay chain of ⁴⁴Ti (⁴⁴Ti → ⁴⁴Sc → ⁴⁴Ca), which is accompanied by line emission at 68, 78, and 1157 keV. As the decay time of ⁴⁴Ti is ∼86 yr, one expects this line emission from young supernova remnants. Although the ⁴⁴Ti yield (typically 10⁻⁵–10⁻⁴M_⊙) is not very high, its production is very sensitive to the energetics and asymmetries of the supernova explosion, and to the mass cut, which defines the mass of the stellar remnant. This makes ⁴⁴Ti an ideal tool to study the inner layers of the supernova explosion. This is of particular interest in light of observational evidence for asymmetric supernova explosions.The γ-ray line emission from ⁴⁴Ti has so far only been detected from the supernova remnant Cas A. I discuss these detections, which were made by COMPTEL (the 1157 keV line) and BeppoSAX (the 68 and 78 keV lines), which, combined, give a flux of (2.6 ± 0.4 ± 0.5) × 10⁻⁵ ph cm⁻² s⁻¹ per line, suggesting a ⁴⁴Ti yield of (1.5 ± 1.0) × 10⁻⁴M_⊙. Moreover, I present some preliminary results of Cas A observations by INTEGRAL, which so far has yielded a 3σ detection of the 68 keV line with the ISGRI instrument with a flux that is consistent with the BeppoSAX detections. Future observations by INTEGRAL-ISGRI will be able to constrain the continuum flux above 90 keV, as the uncertainty about the continuum shape, is the main source of systematic error for the 68 and 78 keV line flux measurements. Moreover, with the INTEGRAL-SPI instrument it will be possible to measure or constrain the line broadening of the 1157 keV line. A preliminary analysis of the available data indicates that narrow line emission (i.e., Δv < 1000 km s⁻¹) can be almost excluded at the 2σ level, for an assumed line flux of 1.9 × 10⁻⁵ ph cm⁻² s⁻¹.     

OPTIS – A Satellite test of Special and General Relativity

OPTIS has been proposed as a small satellite platform in a high elliptical orbit (apogee 40,000 km, perigee 10,000 km) and is designed for high precision tests of foundations of Special and General Relativity. The experimental set-up consists of two ultrastable Nd:YAG lasers, three crossed optical resonators (monolithic cavities), an ensemble of atomic clocks, an optical comb generator, laser tracking devices and a drag-free control system. OPTIS enables improved tests of (1) the isotropy and (2) constancy of the speed of light, (3) special relativistic time dilation, (4) the universality of the gravitational redshift by comparison of various clocks, can measure (5) the absolute value of the gravitational redshift, (6) the Lense–Thirring effect and (7) the perigee advance and (8) can make a test of a hypothetical Yukawa part in the gravitational potential. To avoid any influence from atmospheric drag, solar radiation, or Earth albedo, the satellite needs drag-free control to depress the residual acceleration down to 10⁻¹⁴ m/s² in the frequency range between 10⁻² and 10⁻³ Hz. Precise thermal control must be used to stabilize the cavity temperature to within one part in 10⁷ at time scales of 100 s and to one part in 10⁵ on the orbit time scale.     

Supernova remnant 1987A: The latest report from the Chandra X-ray Observatory

We continue monitoring supernova remnant (SNR) 1987A with the Chandra X-ray Observatory. As of 2004 January, bright X-ray spots in the northwest and the southwest are now evident in addition to the bright eastern ring. The overall X-ray spectrum, since 2002 December, can be described by a planar shock with an electron temperature of ∼2.1 keV. The soft X-ray flux is now 8 × 10⁻¹³ ergs cm⁻² s⁻¹, which is about five times higher than four years ago. This flux increase rate is consistent with our prediction based on an exponential density distribution along the radius of the SNR between the HII region and the inner ring. We still have no direct evidence of a central point source, and place an upper limit of L_X = 1.3 × 10³⁴ ergs s⁻¹ on the 3–10 keV band X-ray luminosity.     

Study of >100 GeV electrons with BETS detector using a long duration balloon flight in Antarctica

We have observed cosmic-ray electrons from 10 to 1000 GeV by a long duration balloon flight using Polar Patrol Balloon (PPB) in Antarctica. The observation was carried out for 13 days at an altitude of 35 km in January 2004. The detector is an imaging calorimeter composed of scintillating-fiber belts and plastic scintillators inserted between lead plates. The geometrical factor of detector is about 600 cm²sr and the total thickness of lead absorber is 9 radiation lengths. The performance of the detector has been confirmed by the CERN-SPS beam test and also investigated by Monte-Carlo simulations. New telemetry system using a commercial satellite of iridium, power supply by solar batteries, and automatic level control using CPU have successfully been developed and operated during the flight. We have collected 5.7 × 10³ events over 100 GeV including nearly 100 candidates of primary electrons.     

The superbubble origin for galactic cosmic rays

The source of energy for cosmic-ray acceleration appears to be shock waves driven by supernova (SNe) ejecta. The great majority (80–90%) of SNe (SNII and SNIb) are formed by the core collapse of young, massive O and B stars. However, it has been known for more than forty years that the births of such massive stars in stellar clumps, termed OB associations, are correlated in space and time. The combined ejecta of core-collapse SNe, occurring at the deaths of these massive stars, create low-density (∼5.0 × 10⁻³ cm⁻³) superbubbles that reach dimensions of several hundred pc. The occurrence of correlated SNe in superbubbles affects not just the source of cosmic-ray energy, SNe shock waves, it impinges as well on the elemental and isotopic source abundances of cosmic-ray nuclei. We argue that the well-known anomalous cosmic-ray ²²Ne/²⁰Ne ratio, a factor of five times the Solar System ratio, results from a mixing of freshly synthesized nucleosynthetic material in supernova active cores of superbubbles. Although diluted by mixing with older, lower metallicity interstellar gas, the mean metallicities in the superbubble, SNe-dominated cores are high ∼3 times the Solar System value.     

XMM-Newton observations of the Mouse,SLX 1744–299 and SLX 1744–300

We observed the radio and X-ray source G359.23–0.82, also known as “the Mouse”, with XMM-Newton. The X-ray image of this object shows a point-like source at the Mouse’s “head”, accompanied by a “tail” that extends for about 40″ westward. The morphology is consistent with that observed recently with Chandra [Gaensler, B.M., van der Swaluw, E., Camilo, F., et al. The Mouse that soared: high resolution X-ray imaging of the pulsar-powered bow shock G359.23–0.82, ApJ 616, 383–402, 2004]. The spectrum of the head can be described by a power-law model with a photon index Γ ≃ 1.9. These results confirm that the Mouse is a bow-shock pulsar wind nebula (PWN) powered by PSR J1747–2958. We found that the hydrogen column density toward the Mouse, N_H = (2.60 ± 0.09) × 10²² cm⁻², is 20%–40% lower than those toward two serendipitously detected X-ray bursters, SLX 1744–299 and SLX 1744–300. At a plausible distance of 5 kpc, the X-ray luminosity of the Mouse, L(0.5–10 keV) = 3.7 × 10³⁴ erg s⁻¹, is 1.5% of the pulsar’s spin-down luminosity. We detected a Type I X-ray burst from SLX 1744–300 and found a possible decrease of N_H and persistent luminosity for this source, in comparison with those observed with ROSAT in 1992.     

Numerical modeling of the pulsar wind interaction with ISM

Time dependent numerical simulation of relativistic wind interaction with interstellar medium was performed. The winds are ejected from magnetosphere of rotation powered pulsars. The particle flux in the winds is assumed to be isotropic. The energy flux is taken as strongly anisotropic in accordance with prediction of the MHD theory of the relativistic winds. The modeling has been performed for the wind magnetization in the range 3 × 10⁻³–10⁻¹. The numerical solutions reproduce the most spectacular features observed in the central part of plerions: toroidal structure and jet-like features. Increase of the wind’s magnetization results in decrease of the size of the synchrotron nebula.     

The temperature map of the loop and coronal sources of an occulted solar flare

An occulted solar flare occurred at about 06:07 UT on 2002, November 2. The RHESSI X-ray images show two separate parts. The lower part consists of a complete loop and the upper part a coronal source which well extends above the solar limb. The loop source shrank for about 3 min with a speed of ∼24 km s⁻¹ during the early impulsive phase and then expanded at ∼7 km s⁻¹, while the coronal source presented an upward motion at about 6 km s⁻¹. We obtained the temperature map of the loop source from RHESSI image spectrum. The temperature of the loop increases with altitude, indicating that the reconnection X-point of this flare is located above the loop source. However, the apparent coronal source is the top of another independent large-scale loop.     

Features of the 11-year variation of galactic cosmic rays in different periods of solar magnetic cycles

Variations of galactic cosmic ray intensity have been studied based on the neutron monitors and interplanetary magnetic field experimental data for different ascending and descending epochs of solar activity. The dependence of the diffusion coefficient on the cosmic ray particles rigidity R is stronger in the maxima epoch than in the minima epoch of solar activity. For the period of 1977–1981 (qA > 0) the diffusion coefficient for the minimum epoch is, χ_min ∝ R^{0.7 ± 0.04} and for the maximum χ_max ∝ R^{1.3 ± 0.05}; for the period of 1987–1990 (qA < 0), χ_min ∝ R^{0.8 ± 0.05} and χ_max ∝ R^{1.1 ± 0.04}. The exponents ν_y and ν_z of the power spectral density of the B_y and B_z components of the IMF in the region of the frequencies (10⁻⁶– 4 × 10⁻⁶) Hz are larger for the minimum epoch of 1987 (ν_y ≈ 2.0 and ν_z ≈ 1.93) than for the maximum epoch of 1990 (ν_y ≈ 1.43 and ν_z ≈ 1.27).     

The nature of nonthermal X-ray filaments near the galactic center

Recent Chandra and XMM-Newton observations reported evidence of two X-ray filaments G359.88−0.08 (SgrA-E) and G359.54+0.18 (the ripple filament) near the Galactic center. The X-ray emission from these filaments has a nonthermal spectrum and coincides with synchrotron emitting radio sources. Here, we report the detection of a new X-ray feature coincident with a radio filament G359.90−0.06 (SgrA-F) and show more detailed VLA, Chandra and BIMA observations of the radio and X-ray filaments. In particular, we show that radio emission from the nonthermal filaments G359.90−0.06 (SgrA-F) and G359.54+0.18 (the ripple) has a steep spectrum whereas G359.88−0.08 (SgrA-E) has a flat spectrum. The X-ray emission from both these sources could be due to synchrotron radiation. However, given that the 20 km s⁻¹ molecular cloud, with its intense 1.2 mm dust emission, lies in the vicinity of SgrA-F, it is possible that the X-rays could be produced by inverse Compton scattering of far-infrared photons from dust by the relativistic electrons responsible for the radio synchrotron emission. The production of X-ray emission from ICS allows an estimate of the magnetic field strength of 0.08 mG within the nonthermal filament. This should be an important parameter for any models of the Galactic center nonthermal filaments.     

Giant radio pulses from the Crab pulsar

Individual giant radio pulses (GRPs) from the Crab pulsar last only a few microseconds. However, during that time they rank among the brightest objects in the radio sky reaching peak flux densities of up to 1500 Jy even at high radio frequencies. Our observations show that GRPs can be found in all phases of ordinary radio emission including the two high frequency components (HFCs) visible only between 5 and 9 GHz [Moffett, D.A., Hankins, T.H. Multifrequency radio observations of the Crab pulsar. Astrophys. J. 468, 779–783, 1996]. This leads us to believe that there is no difference in the emission mechanism of the main pulse (MP), inter pulse (IP) and HFCs. High resolution dynamic spectra from our recent observations of giant pulses with the Effelsberg telescope at a center frequency of 8.35 GHz show distinct spectral maxima within our observational bandwidth of 500 MHz for individual pulses. Their narrow band components appear to be brighter at higher frequencies (8.6 GHz) than at lower ones (8.1 GHz). Moreover, there is an evidence for spectral evolution within and between those structures. High frequency features occur earlier than low frequency ones. Strong plasma turbulence might be a feasible mechanism for the creation of the high energy densities of ∼6.7 × 10⁴ erg cm⁻³ and brightness temperatures of ∼10³¹ K.     

Impacts of the PAH size and the radiation intensity on the IR features of illuminated dust within the reflection nebulae

Interstellar dust grains are illuminated in the reflection nebulae. Under conditions of the PAH size and the intensity of the interstellar radiation field, we follow their impact on the PAH aromatic infrared bands using the numerical DustEM code. For a dust model consisting of PAH, amorphous C and amorphous silicate, the PAH size varies in a range from 0.31 to 4.9 nm while the radiation intensity varies by a scale factor from 0.1 to 10⁴. Various trends of the results show the effect of varying both the PAH size and the radiation intensity on the strength of the aromatic mid-IR bands. Through small PAH sizes less than 0.7 nm, the grain temperature distribution of PAHs shows a small variation within 2–3 K at low radiation intensity while it increases to 15 and 8 K for PAH⁰ and PAH⁺, respectively, at higher radiation intensity. In final the variability in these results reveals the evolution of the dust grains under the physical space conditions of the reflection nebulae. In the mid-IR region, the contributions of PAH⁰ and PAH⁺ in the total SED intensity agree with the proportions of these PAHs observed in some reflection nebulae having higher radiation intensities.     

Evaluation of solar proton spectra using balloon cosmic ray observations and Monte Carlo simulation results

We have developed a method to evaluate the spectrum of solar energetic protons at the top of the Earth’s atmosphere from the measurements of our balloon cosmic ray experiment. By using the Monte Carlo PLANETOCOSMICS code based on Geant4 we compute the interaction of solar protons [10 MeV–10 GeV] with the Earth’s atmosphere. We obtain the angular and energy distributions of secondary particles (p, e⁻, e⁺, photons, muons) at different atmospheric levels as a function of primary proton spectra. By comparing the calculated depth dependence of the particle flux with the data obtained by our balloon experiment we can deduce the parameters of the solar proton spectrum that best fit the observations. In this paper we discuss our solar proton spectrum estimation method, and present results of its application to selected solar proton events from 2001 to 2005.     

Global MGS TES data and Mars-GRAM validation

Mars Global Reference Atmospheric Model (Mars-GRAM 2001) is an engineering-level Mars atmosphere model widely used for many Mars mission applications. From 0 to 80 km, it is based on NASA Ames Mars General Circulation Model (MGCM), while above 80 km it is based on University of Michigan Mars Thermospheric General Circulation Model. Mars-GRAM 2001 and MGCM use surface topography from Mars Global Surveyor Mars Orbiter Laser Altimeter (MOLA). Validation studies are described comparing Mars-GRAM with a global summary data set of Mars Global Surveyor Thermal Emission Spectrometer (TES) data. TES averages and standard deviations were assembled from binned TES data which covered surface to ∼40 km, over more than a full Mars year (February 1999–June 2001, just before start of a Mars global dust storm). TES data were binned in 10° × 10° latitude–longitude bins (36 longitude bins, centered at 5°–355°, by 18 latitude bins, centered at −85° to +85°), and 12 seasonal bins (based on 30° increments of L_s angle). Bin averages and standard deviations were assembled at 23 data levels (temperature at 21 pressure levels, plus surface temperature and surface pressure). Two time-of-day bins were used: local time near 2 or 14 h. Two dust optical depth bins were used: infrared optical depth, either less than or greater than 0.25 (which corresponds to visible optical depth less than or greater than about 0.5). For interests in aerocapture and precision entry and landing, comparisons focused on atmospheric density. TES densities versus height were computed from TES temperature versus pressure, using assumptions of perfect gas law and hydrostatics. Mars-GRAM validation studies used density ratio (TES/Mars-GRAM) evaluated at data bin center points in space and time. Observed average TES/Mars-GRAM density ratios were generally 1 ± 0.05, except at high altitudes (15–30 km, depending on season) and high latitudes (>45°N), or at most altitudes in the southern hemisphere at L_s ∼ 90° and 180°. Compared to TES averages for a given latitude and season, TES data had average density standard deviation about the mean of ∼2.5% for all data, or ∼1–4%, depending on time of day and dust optical depth. Average standard deviation of TES/Mars-GRAM density ratio was 8.9% for local time 2 h and 7.1% for local time 14 h. Thus standard deviation of observed TES/Mars-GRAM density ratio, evaluated at matching positions and times, is about three times the standard deviation of TES data about the TES mean value at a given position and season.     

Propagation of secondary antiprotons and cosmic rays in the Galaxy

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratios (e.g., B/C) produce too few antiprotons. In the present paper, we discuss one possibility to overcome these difficulties. Using the measured antiproton flux and B/C ratio to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local “unprocessed” component at low energies perhaps associated with the Local Bubble, thus decreasing the measured secondary to primary nuclei ratio. The independent evidence for SN activity in the solar vicinity in the last few Myr supports this idea. The model reproduces antiprotons, B/C ratio, and elemental abundances up to Ni (Z ⩽ 28). Calculated isotopic distributions of Be and B are in perfect agreement with CR data. The abundances of three “radioactive clock” isotopes in CR, ¹⁰Be, ²⁶Al, ³⁶Cl, are all consistent and indicate a halo size z_h ∼ 4 kpc based on the most accurate data taken by the ACE spacecraft.     

Scaling and universality of the thermal conductivity of liquid 4He near the superfluid transition and in restricted geometries

We present measurements of the thermal conductivity λ(t, P, L) = l/ρ(t, P, L) near the superfluid transition of ⁴He at saturated vapor pressure and confined in cylindrical geometries with radii L = 0.5 and 1.0 μm (t ≡ T/T_λ(P) − 1). For L = 1.0 μm measurements at six pressures P are presented. At and above T_λ the data are consistent with a universal scaling function F(X) = (L/ξ_o)^x/ν(ρ/ρ₀), X = (L/ξ_o)^1/νt valid for all P (ρ₀ and x are the pressure-dependent amplitude and effective exponent of the bulk resistivity ρ(t, P, ∞) = ρ₀t^x and ξ = ξ₀t^−ν is the correlation length). Indications of breakdown of scaling and universality are observed below T_λ.     

Observations of “zebra” pattern in cm-range with spatial resolution

We present the results of the first observations of the solar microwave burst with fine spectral structure of zebra type at the frequency about 5.7 GHz. The burst has been detected simultaneously by the Siberian Solar Radio Telescope and by the spectropolarimeter of the National Astronomical Observatory of China. Zebra pattern consisted of three parallel stripes with complex frequency drift. The degree of circular polarization of emission reached 100%, the polarization sense corresponded to the extraordinary wave (X-mode). We have determined the plasma parameters in the emission source: plasma density about 10¹¹ cm⁻³, magnetic field strength 60–80 G. We argue that in the given event the most probable mechanism of the zebra pattern generation is non-linear coupling of harmonics of Bernstein modes.     

Nonlinear processes in cosmic-ray precursor of strong supernova shock: Maximum energy and average energy spectrum of accelerated particles

The instability in the cosmic-ray precursor of a supernova shock is studied. The level of turbulence in this region determines the maximum energy of accelerated particles. The consideration is not limited by the case of weak turbulence. It is assumed that the Kolmogorov type nonlinear wave interactions together with the ion-neutral collisions restrict the amplitude of random magnetic field. As a result, the maximum energy of accelerated particles strongly depends on the age of a SNR. The average spectrum of cosmic rays injected in the interstellar medium in the course of adiabatic SNR evolution takes the approximate form E⁻² at energies larger than 10–30 GeV/nucleon with the maximum energy that is close to the position of the knee in cosmic-ray spectrum at 4 × 10¹⁵ eV. At an earlier stage of SNR evolution – the ejecta-dominated stage, the particles are accelerated to higher energies and have a rather steep power-law distribution. These results suggest that the knee may mark the transition from the ejecta-dominated to the adiabatic evolution of SNR shocks which accelerate cosmic rays.