Propagation and 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 discreteness in the source distribution, for local source overdensity, 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.     

A linear stability analysis of large-Prandtl-number thermocapillary liquid bridges

A linear stability analysis is applied to determine the onset of oscillatory thermocapillary convection in cylindrical liquid bridges of large Prandtl numbers (4  Pr  50). We focus on the relationships between the critical Reynolds number Re_c, the azimuthal wave number m, the aspect ratio Γ and the Prandtl number Pr. A detailed Re_c–Pr stability diagram is given for liquid bridges with various Γ. In the region of Pr > 1, which has been less studied previously and where Re_c has been usually believed to decrease with the increase of Pr, we found Re_c exhibits an early increase for liquid bridges with Γ around one. From the computed surface temperature gradient, it is concluded that the boundary layers developed at both solid ends of liquid bridges strengthen the stability of basic axisymmetric thermocapillary convection at large Prandtl number, and that the stability property of the basic flow is determined by the “effective” part of liquid bridge.     

Equatorial F-region vertical ion drifts during quiet solar maximum

Median values of ionosonde h′F data acquired at Ibadan (Geographic:7.4°N, 3.9°E, Magnetic: dip 6°S, and magnetic declination, 3°W), Nigeria, West Africa, have been used to determine vertical ion drift (electric field) characteristics in the postsunset ionosphere in the African region during a time of high solar activity (average F_10.7 −208). The database spans from January and December 1958 during the era of International Geophysical Year (IGY) for geomagnetic quiet conditions. Bimonthly averaged diurnal variations patterns are very similar, but differ significantly in magnitude and in the evening reversal times. Also, monthly variations of F-region vertical ion drift reversal times inferred from the time of h′F maximum indicates early reversal during equinoxes and December solstice months except for the month of April. Late reversal is observed during the June solstice months. The equatorial evening prereversal enhancement in vertical ion drift (V_zp) occurs largely near 1900 LT with typical values 20–45 m/s. Comparison of Ibadan ionosonde V_zp with the values of prereversal peak velocity reported for Jicamarca (South America), Kodaikanal (India), and Scherliess and Fejer global model show considerable disparity. The changes of postsunset peak in virtual height of F-layer (h′F_P) with prereversal velocity peak V_zp are anti-correlated. Investigation of solar effects on monthly values of V_zp and h′F_P revealed that these parameters are independent of monthly averaged solar flux intensity during quiet-time sunspot maximum conditions.     

Geomagnetic effects on mid-latitude railways: A statistical study of anomalies in the operation of signaling and train control equipment on the East-Siberian Railway

Intense geomagnetically induced currents (GIC) can hamper rail traffic by disturbing signaling and train control systems. GIC threats have been a concern for technological systems at high-latitude locations due to geomagnetic disturbances driven by substorm expansion electrojet or convection electrojet intensifications. However, other geomagnetic storm processes such as storm sudden commencement (SSC) and geomagnetic pulsations can also cause GIC concerns for technological systems. We present in this paper the first evidence based on statistical data for links between geomagnetic disturbances and faulty operations (anomalies) in the functioning of railway automatics and telemetry. We analyze anomalies of automatic signaling and train control equipment which occurred in 2004 on the East-Siberian Railway (corrected geomagnetic latitude _m = 46–51°N and longitude λ_m = 168–187°E). Our results reveal a seasonal effect in the number of anomalies per train similar to the one observed in geomagnetic activity (Kp, Ap, Dst indices). We also found an increase by a factor of 3 in the total duration of daily anomalies during intense geomagnetic storms (local geomagnetic index specific to Siberian Observatory A_max > 30), with a significant correlation between the daily sum of durations of anomalies with geomagnetic activity. Special attention was paid to failures not related to recognized technical malfunctions. We found that the probability of these failures occurring in geomagnetically disturbed periods was 5–7 times higher than the average anomaly occurrence.     

X-ray spectra and time variability of active galactic nuclei

The X-ray spectra of broad line active galactic nuclei of all types (Seyfert I's, NELG's, Broadline radio galaxies) is well fit by a power law in the .5–100 keV band of mean energy slope α = .68±.15. There is, as yet, no strong evidence for time variability of this slope in a given object. The constraints that this places on simple models of the central energy source are discussed. BL Lac objects have quite different X-ray spectral properties and show pronounced X-ray spectral variability.On time scales longer than 12 hours most radio quiet AGN do not show strong, ΔI/I > .5, variability. The probability of variability of these AGN seems to be inversely related to their luminosity. However characteristic timescales for variability have not been measured for many objects. This general lack of variability may imply that most AGN are well below the Eddington limit. Radio bright AGN tend to be more variable on long, τ 6 month, timescales.     

ETALON spin period determination from kHz SLR data

Using kHz Satellite LASER Ranging (SLR) data of the SLR station Graz only, we determined the spin periods of the two ETALON satellites – launched into high orbits of about 20,000 km – and their spin period increase during 3 years. The determined spin period values and spin period increase rates at 2004-01-01 are: T_ET1 = 63 s + 0.484 s/year, and T_ET2 = 65.5 s + 0.401 s/year.     

The iron Kα-line diagnostics of the rotating black hole metric in Seyfert galaxies

Observations of Seyfert galaxies in X-ray region reveal the wide emissive lines in their spectra, which can arise in inner parts of accretion disks, where the effects of General Relativity (GR) must be taken into account. A spectrum of a solitary emission line of a hot spot in Kerr accretion disk is simulated depending on the radial coordinate r and the angular momentum a = J/M of a black hole, under the assumption of equatorial circular motion of a hot spot. It is shown that the characteristic two-peak line profile with the sharp edges arises at a large distance (about r ≈ (3–10) r_g). The inner regions emit the line, which is observed with one maximum and extremely wide red wing. We present results of simulations for the isothermal and Shakura–Sunayev disks.     

Analytical approach to cosmic ray ionization by nuclei with charge Z in the middle atmosphere – Distribution of galactic CR effects

The effects of galactic and solar cosmic rays (CR) in the middle atmosphere are considered in this work. A new analytical approach for CR ionization by protons and nuclei with charge Z in the lower ionosphere and middle atmosphere is developed in this paper. For this purpose the ionization losses (dE/dh) according to the Bohr–Bethe–Bloch formula for the energetic charged particles are approximated in three different energy intervals. More accurate expressions for energy decrease E(h) and electron production rate profiles q(h) are derived. The obtained formulas allow comparatively easy computer programming. The integrand in q(h) gives the possibility for application of adequate numerical methods – such as Romberg method or Gauss quadrature, for the solution of the mathematical problem. On this way the process of interaction of cosmic ray particles with the upper, middle and lower atmosphere will be described much more realistically. Computations for cosmic ray ionization in the middle atmosphere are made. The full CR composition is taken into account: protons, Helium (α-particles), light L, medium M, heavy H and very heavy VH group of nuclei.     

High-energy neutrinos from powerful radio galaxies

In order to explore mechanisms for the production of radio lobes from radio galaxies, we propose observational tests involving neutrinos at E ≥ 4 TeV. Among the mechanisms that have been suggested are: an explosive burst of energy that has been stored in or near the galactic nucleus; diffusive escape of particles from the vicinity of the galactic nucleus into plasmons; a beam from the galactic core that interacts with the circumgalactic medium; and black holes or spinars ejected from the galactic nucleus by a gravitational slingshot mechanism. As an example, we estimate neutrino fluxes from Cen A. The annual neutrino event rate at energies above 4 TeV is ≈ 10²–10³ for a DUMAND-type dectector if relativistic particles are temporarily stored near the galactic nucleus, and if these have a differential energy spectrum with exponent –2.0. With a similar exponent, but allowing free escape from the galaxy, the fluxes are about 10 times lower.     

Asymmetric precipitation in a coronal loop as explanation of a singular observed spectrum

Almost 10 years of solar submillimeter observations have shown new aspects of solar activity, such as the presence of rapid solar spikes associated with the launch of coronal mass ejections and an increasing submillimeter spectral component in flares. We analyse the singular microwave–submillimeter spectrum of an M class solar flare on 20 December, 2002. Flux density observations measured by Sun patrol telescopes and the Solar Submillimeter Telescope are used to build the radio spectrum, which is fitted using Ramaty’s code. At submillimeter frequencies the spectrum shows a component different from the microwave classical burst. The fitting is achieved proposing two homogeneous sources of emission. This theoretical fitting is in agreement with differential precipitation through a magnetically asymmetric loop or set of loops. From a coronal magnetic field model we infer an asymmetric magnetic structure at the flare location. The model proposed to quantify the differential precipitation rates due to the asymmetry results in a total precipitation ratio Q₂/Q₁≈10⁴–10⁵, where Q₁(Q₂) represents the total precipitation in the loop foot with the high (low) magnetic field intensity. This ratio agrees with the electron total number ratio of the two sources proposed to fit the radio spectrum.     

Study of correlations between waves and particle fluxes measured on board the DEMETER satellite

The topic of relativistic electron dynamics in the outer radiation belt has received considerable attention for many years. Nevertheless, the problem of understanding the physical phenomenon involved is far from being resolved. In this paper, we use DEMETER observations to examine the variations of the energetic electron fluxes and ELF/VLF wave intensities in the inner magnetosphere during the intense 8 November 2004 magnetic storm. Electron flux spectra and associated wave intensity spectra are analysed throughout the magnetic storm and common characteristics or differences to other storm events are retained. The overall objective of this study is to identify and derive parameters that are relevant for particle flux modelling; the time constant characterizing the persistent decay after particle enhancement was found to be one of these important model parameters.The analysis of the 8 November 2004 event reveals that for L-shell parameter higher than 4, an electron flux dropout is observed during the storm’s main phase for electrons in the energy range 0.1–1 MeV, as has been reported from other measurements. Characteristic wave spectra accompanying this phase are analysed. They show a typical enhancement in the frequency range 0.3–10 kHz at onset for all L-shell values under consideration (2 < L < 5). During the first stage of the recovery phase, the electron fluxes are increased to a level higher than the pre-storm level, whereas the level of wave intensity in the frequency range observed below 300 Hz is at its highest. In the second stage, the particle flux decrease goes hand in hand with a global wave activity decline, the relaxation time of the latter being smaller than the former’s one. In some other cases, long-lasting electron enhancement associated with constant wave activity has been observed during this latter stage. For the above mentioned storm, while at low L values the decay time constants are higher for low energy electrons than for high energy electrons, this order is reversed at high L values. At about L = 3.6 the time constant is independent of electron energy.     

H2 and H∞ control options for the combined attitude and thermal control system (CATCS)

The combined attitude and thermal control system (CATCS) combines the conventional attitude control and thermal control subsystems. Its principle is based on circulating a heat conducting fluid inside a closed duct wielding the excess onboard heat in order to produce the attitude control torques. Previously only the proportional-integral (PI) controller has been tested for CATCS. In this paper two other control options for CATCS were designed based on the H₂ and H_∞ control methods to improve the attitude control performance of a small satellite. The control gain matrix with the minimum cost function is obtained by solving the Riccati equation and fed back to the system in order to achieve the system’s performance. The designed controllers can efficiently control the roll, pitch and yaw satellite attitudes. Simulations for the two techniques were carried out using Matlab and Simulink for ideal and non-ideal system models. Results show that the H₂ controller has a better attitude control performance over the H_∞ controller and PI controller itself.     

Spectral characteristics of electron fluxes at L < 2 under the Radiation Belts

The paper presents the analysis of experimental data on electron fluxes with energies 10 keV–10 MeV. Data were obtained during 1978–2005 years in different space experiments (COSMOS-900, MIR Space Station, ACTIVE, SAMPEX, CORONAS-I, CORONAS-F, NOAA POES-17, TATYANA and others). Two areas of electron flux enhancements are studied in the paper: the near-equatorial (L < 1.2) zone and the middle-latitude (1.2 < L < 1.9) zone. It is shown that electron flux enhancements are regularly registered at L < 2 and the observed formations have some typical features. Electron peaks at L < 1.2 appear sporadically while peaks at 1.2 < L < 1.9 are observed regularly. The approximations of spectra by several functions including kappa-function are presented.     

Low energy observations of Cygnus X-2 by Ariel VI

Ariel VI observations of Cygnus X-2 have revealed a rather flat spectrum between 0.1 and 1.5 keV with variable emission at low energy. Of the two conflicting interpretations of this object in terms of i) a distant high-luminosity (L_x 10³⁸ ergs s⁻¹) binary and ii) a nearby low-luminosity (L_x 10³⁵ ergs s⁻¹) degenerate dwarf system, our measurements support the latter.     

Navigating BepiColombo during the weak-stability capture at Mercury

BepiColombo is scheduled for launch in August 2013 and to arrive after a nearly six-year long transfer at Mercury in June 2019. The trajectory has a number of challenging elements: a launch with Soyuz/Fregat into a geostationary transfer orbit, followed by a lunar flyby, long low-thrust arcs and five more planetary flybys (one at the Earth, two at Venus and two at Mercury). At arrival the low thrust arcs reduce the approach velocity so much that BepiColombo passes by the Sun–Mercury Lagrange points L₁ and L₂ and gets weakly captured in a highly eccentric orbit around Mercury in case the orbit insertion manoeuvre would fail.This paper describes the navigation strategy during the final phase. Five trajectory correction manouevres during the last 65 days requiring up to 20 m/s (3σ) are proposed. With this strategy it is possible to navigate BepiColombo safely through the weak-stability boundary of Mercury and to reach the target periherm with a precision of 11 km.     

Evaluation of gravity-dependent membrane potential shift in Paramecium

It is still debated whether or not gravity can stimulate unicellular organisms. This question may be settled by revealing changes in the membrane potential in a manner depending on the gravitational forces imposed on the cell. We estimated the gravity-dependent membrane potential shift to be about 1 mV G⁻¹ for Paramecium showing gravikinesis at 1–5 G, on the basis of measurements of gravity-induced changes in active propulsion and those of propulsive velocity in solutions, in which the membrane potential has been measured electrophysiologically. The shift in membrane potential to this extent may occur from mechanoreceptive changes in K⁺ or Ca²⁺ conductance by about 1% and might be at the limit of electrophysiological measurement using membrane potential-sensitive dyes. Our measurements of propulsive velocity vs membrane potential also suggested that the reported propulsive force of Paramecium measured in a solution of graded densities with the aid of a video centrifuge microscope at 350 G was 11 times as large as that for −29 mV, i.e., the resting membrane potential at [K⁺]_o = 1 mM and [Ca²⁺]_o = 1 mM, and, by extrapolation, that Paramecium was hyperpolarized to −60 mV by gravity stimulation of 100- G equivalent, the value corrected by considering the reduction of density difference between the interior and exterior of the cell in the graded density solution. The estimated shift of the membrane potential from −29 mV to −60 mV by 100- G equivalent stimulation, i.e., 0.3 mV G⁻¹, could reach the magnitude entirely feasible to be measured more directly.     

On the possibility to check the magnetosphere’s model by CR: The strong geomagnetic storm in November 2003

We report on a study of cosmic ray cutoff rigidity variations during the strong geomagnetic storm of 18–24 November 2003. We employed the most recent Tsyganenko magnetospheric model to represent the very strong disturbed Magnetosphere. We used this magnetic field for the cosmic ray trajectory calculations to determine the geomagnetic cutoff rigidity throughout this period of severe geomagnetic disturbance. We determine the cutoff rigidity changes during this period by two methods, by trajectory calculations and by the spectrographic global survey method. The values of the change of cutoff rigidities obtained by two different methods are correlated with the D_st and interplanetary magnetic field and plasma parameters and result in correlation coefficients in the range 0.63–0.84 for the various cosmic ray stations. The result of this study indicates that the most significant contributions to the cutoff rigidity changes are due to D_st variation although the influence of solar wind density and B_z and B_y components of IMF variations is significant.     

Investigations of artificial HF plasma turbulence features using stimulated electromagnetic emission

New results on growth and decay times of the stimulated electromagnetic emission (SEE), obtained using a special quasi-continuous schedule of HF pump wave radiation, are presented. The employed technique allows to study the evolution of HF plasma turbulence with a high time resolution ( 0.25 ms). It is established that for the broad continuum SEE component the emission intensity starts to decrease within the delay time τ_D 0.5−3.2 ms after pump wave turn-off, where longer times correspond to lower SEE frequencies. The effect is explained in terms of transfer of HF-excited plasma waves throughout their spectrum due to induced scattering by thermal ions.     

Expanding shells of young pulsars as sources of high-energy neutrinos

The initial power outputs P_o of pulsars are not yet well known, but these seem to follow approximately a distribution law N(> P_o) ∝ P_o⁻ⁿ where 0.5 ≤ n ≤ 1.0. It seems likely that P_o ≥ 10³⁸ ergs/sec. With these assumptions, we estimate that the DUMAND detector can record ≥ 10³ high-energy (> 4 TeV) neutrino events in a four-month period per Galactic supernova; (in our Galaxy, these are estimated to occur at the rate of about 8 per century.) Neutrinos from supernova shells in the Virgo supercluster would be marginally detectable (one very bright supernova per decade at about 20 Mpc) if N(> P_o) ∝ P_o^−0.5, but undetectable if N(> P_o) ∝ P_o⁻¹. The diffuse flux summed over distant extragalactic supernovae is likely to be well below the detection threshold.     

Gamma rays from accreting matter onto collapsed objects

In the spherical accretion onto massive objects, the matter may be heated up to temperatures as high as 10¹² °K. In such a hot plasma, the thermal bremsstrahlung (e-e and e-p) and π° decay from inelastic collisions of protons are the main γ-ray sources. We determined the γ -ray production spectra from the π° decay and from bremsstrahlung for different temperatures. The expected γ-ray spectra were evaluated too in order to fit experimental data. We have fitted COS B data from 3C 273 using a two temperatures plasma model. The best fit is for

Full-size image (1K)

(M₈ is the black hole mass in 10⁸ M) which gives . The hard X-ray measurements do not contradict the bremsstrahlung spectrum.