Effects of the position of the solar wind termination shock and the heliopause on the heliospheric modulation of cosmic rays

The effects of changing the position of the solar wind termination shock and the position of the heliopause, and therefore the extent of the heliosheath, on the modulation of cosmic ray protons are illustrated. An improved numerical model with diffusive termination shock acceleration, a heliosheath and drifts is used. The modulation is computed in the equatorial plane and at 35 heliolatitude using recently derived diffusion coefficients applicable to a number of cosmic ray species during both magnetic polarity cycles of the Sun. It was found that qualitatively the modulation results for the different heliopause positions are similar although they differ quantitatively, e.g., clearly different radial gradients are predicted for the regions beyond the termination shock compared to inside the shock. The difference between the modulation for the two solar polarity cycles are less significant at a heliolatitude of 35° than in the equatorial plane. We found that moving the termination shock from 90 to 100 AU, with the heliopause fixed at 120 AU, caused only quantitative differences so that the exact position of the TS in the outer heliosphere seems not crucially important to global modulation. Moving the heliopause outwards, to represent the modulation in the tail region of the heliosphere, causes overall decreases in the cosmic ray intensities but not linearly as a function of energy, e.g., at 1 GeV the effect is insignificant. We conclude from this modelling that the modulation of protons in the heliospheric nose and tail regions are qualitatively similar although, clear quantitative and interesting differences occur.     

The heliospheric modulation of cosmic rays: Effects of a latitude dependent solar wind termination shock

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock and that it has begun its merger with the local interstellar medium. The compression ratio of this shock affects galactic cosmic rays when they enter the heliosphere. Hydrodynamic (HD) models show that the compression ratio can have a significant latitude dependence; with the largest value in the nose direction of the heliosphere, becoming significantly less towards the polar regions. The modulation effects of such large latitude dependence are studied, using a well-established numerical drift and shock modulation model. We focus on computing the modulated spectra for galactic protons with emphasis on the radial and polar gradients in the equatorial plane and at a polar angle of θ = 55°, corresponding to the heliolatitude of Voyager 1. Two sets of solutions are computed and compared each time; with and without a latitude dependence for the compression ratio. All computations are done for the two magnetic field polarity cycles assuming solar minimum conditions. Including the termination shock in the model allows the study of the re-acceleration of galactic protons in the outer heliosphere. We find that for the A < 0 polarity cycle the intensity between ∼200 MeV and ∼1 GeV in the vicinity of the shock in the heliospheric equatorial plane may exceed the local interstellar value specified at the heliopause. Unfortunately, at θ = 55°, the effect is reduced. This seems not possible during an A > 0 cycle because significant modulation is then predicted between the heliopause and the termination shock, depending on how strong global gradient and curvature drifts are in the heliosheath. The overall effect of the shock on galactic protons in the equatorial plane is to reduce the total modulation as a function of radial distance with respect to the interstellar spectrum. Making the compression ratio latitude dependent enhances these effects at energies E < 200 MeV in the equatorial plane. At larger heliolatitudes these effects are even more significant. The differences in the modulation between the two drift cycles are compelling when the compression ratio is made latitude dependent but at Earth this effect is insignificant. A general result is that the computed radial gradient changes for galactic protons at and close to the TS and that these changes are polarity dependent. In line with previous work, large polarity dependent effects are predicted for the inner heliosphere and also close to the shock’s position in the equatorial plane. In contrast, at θ = 55°, the largest polarity effect occurs in the middle heliosphere (50 AU), enhanced by the latitude dependence of the compression ratio. At this latitude, the amount of proton modulation between the heliopause and the termination shock is much reduced. If galactic cosmic rays were to experience some diffusive shock acceleration over the 100–1000 MeV range at the shock, the radial gradient should change its sign in the vicinity of the shock, how large, depends on the compression ratio and the amount of drifts taking place in the outer heliosphere. The effective polar gradient shows a strong polarity dependence at Earth but this dissipates at θ = 55°, especially with increasing radial distance. This tendency is enhanced by making the compression ratio latitude dependent.     

The interstellar heliopause probe technology reference study

The interstellar heliopause probe (IHP) is one of ESA’s technology reference studies (TRS). The TRS aim to focus the development of strategically important technologies of relevance to future science missions by studying technologically demanding and scientifically interesting missions that are currently not part of the science mission programme.

Equipped with a highly integrated payload suite (HIPS), the IHP will perform in situ exploration of the heliopause and the heliospheric interface. The HIPS, which is a standard element in all TRSs, miniaturize payloads through resource reduction by using miniaturized components and sensors, and by sharing common structures and payload functionality.

To achieve the scientific requirements of the mission, the spacecraft is to leave the heliosphere as close to the heliosphere nose as possible and reach a distance of 200 AU from the Sun within 25 years. This is possible by using a trajectory with two solar flybys and a solar sail with characteristic acceleration of 1.1 mm/s², which corresponds to a 245 × 245 m² solar sail and a sail thickness of 1–2 μm. The trajectory facilitates a modest sail design that could potentially be developed in a reasonable timeframe.

In this paper, an update to the results of studies being performed on this mission will be given and the current mission baseline and spacecraft design will be described. Furthermore, alternative solar sail systems and enabling technologies will be discussed.     

A first step towards proton flux forecasting

We present a preliminary version of a potential tool for real time proton flux prediction which provides proton flux profiles and cumulative fluence profiles at 0.5 and 2 MeV of solar energetic particle events, from their onset up to the arrival of the interplanetary shock at the spacecraft position (located at 1 or 0.4 AU). Based on the proton transportation model by Lario et al. [Lario, D., Sanahuja, B., Heras, A.M. Energetic particle events: efficiency of interplanetary shocks as 50 keV E < 100 MeV proton accelerators. Astrophys. J. 509, 415–434, 1998] and the magnetohydrodynamic shock propagation model of Wu et al. [Wu, S.T., Dryer, M., Han, S.M. Non-planar MHD model for solar flare-generated disturbances in the Heliospheric equatorial plane. Sol. Phys. 84, 395–418, 1983], we have generated a database containing “synthetic” profiles of the proton fluxes and cumulative fluences of 384 solar energetic particle events. We are currently validating the applicability of this code for space weather forecasting by comparing the resulting “synthetic” flux profiles with those of several real events.     

Modulation of galactic cosmic rays in a north–south asymmetrical heliosphere

Observations made with the two Voyager spacecraft confirmed that the solar wind decelerates to form the heliospheric termination shock. Voyager 1 crossed this termination shock at ∼94 AU in 2004, while Voyager 2 crossed it in 2007 at a different heliolatitude, about 10 AU closer to the Sun. These different positions of the termination shock confirm the dynamic and cyclic nature of the shock’s position. Observations from the two Voyager spacecraft inside the heliosheath indicate significant differences between them, suggesting that apart from the dynamic nature caused by changing solar activity there also may exist a global asymmetry in the north–south (polar) dimensions of the heliosphere, in addition to the expected nose–tail asymmetry. This relates to the direction in which the heliosphere is moving in interstellar space and its orientation with respect to the interstellar magnetic field. In this paper we focus on illustrating the effects of this north–south asymmetry on modulation of galactic cosmic ray Carbon, between polar angles of 55° and 125°, using a numerical model which includes all four major modulation processes, the termination shock and the heliosheath. This asymmetry is incorporated in the model by assuming a significant dependence on heliolatitude of the thickness of the heliosheath. When comparing the computed spectra between the two polar angles, we find that at energies E < ∼1.0 GeV the effects of the assumed asymmetry on the modulated spectra are insignificant up to 60 AU from the Sun but become increasingly more significant with larger radial distances to reach a maximum inside the heliosheath. In contrast, with E > ∼1.0 GeV, these effects remain insignificant throughout the heliosphere even very close to the heliopause. Furthermore, we find that a higher local interstellar spectrum for Carbon enhances the effects of asymmetric modulation between the two polar angles at lower energies (E < ∼300 MeV). In conclusion, it is found that north–south asymmetrical effects on the modulation of cosmic ray Carbon depend strongly on the extent of the geometrical asymmetry of the heliosheath together with the assumed value of the local interstellar spectrum.     

Acceleration and transport of energeticparticles at CME-driven shocks

Historically, solar energetic particle (SEP) events are classified in two classes as “impulsive” and “gradual”. Whether there is a clear distinction between the two classes is still a matter of debate, but it is now commonly accepted that in large “gradual” SEP events, Fermi acceleration, also known as diffusive shock acceleration, is the underlying acceleration mechanism. At shock waves driven by coronal mass ejections (CMEs), particles are accelerated diffusively at the shock and often reach > MeV energies (and perhaps up to GeV energies). As a CME-driven shock propagates, expands and weakens, the accelerated particles can escape ahead of the shock into the interplanetary medium. These escaping energized particles then propagate along the interplanetary magnetic field, experiencing only weak scattering from fluctuations in the interplanetary magnetic field (IMF). In this paper, we use a Monte-Carlo approach to study the transport of energetic particles escaping from a CME-driven shock. We present particle spectra observed at 1 AU. We also discuss the particle “crossing number” at 1AU and its implication to particle anisotropy. Based on previous models of particle acceleration at CME-driven shocks, our simulation allows us to investigate various characteristics of energetic particles arriving at various distances from the sun. This provides us an excellent basis for understanding the observations of high-energy particles made at 1 AU by ACE and WIND.     

Cosmic ray anisotropies in the outer heliosphere

The observation of the directional distribution of energetic and cosmic ray particles has been done with the Voyager spacecraft over a long period. Since 2002, when the first flux enhancements of charged particles associated with the approach of Voyager 1 to the solar wind termination shock were observed, these anisotropy measurements have become of special interest. They play an important role to understand the magnetic field and shock structure and the basics of the modulation of cosmic ray and anomalous particles at and beyond the termination shock. They also serve as motivation to study the spatial behavior of galactic and anomalous cosmic ray anisotropies with numerical modulation models in order to illustrate how the radial anisotropy, at different energies, change from upstream to downstream of the termination shock. Observations made by Voyager 1 indicate that the termination shock is a complicated region than previously thought, hence the effects of the latitude dependence of the termination shock’s compression ratio and injection efficiency on the radial anisotropies of galactic and anomalous protons will be illustrated. We find that the magnitude and direction of the radial anisotropy strongly depends on the position in the heliosphere and the energy of particles. The effect of the TS on the radial anisotropy is to abruptly increase its value in the heliosheath especially in the A > 0 cycle for galactic protons and in both polarity cycles for anomalous protons. Furthermore, the global effect of the latitude dependence of the shock’s compression ratio is to increase the radial anisotropy for galactic protons throughout the heliosphere, while when combined with the latitude dependence of the injection efficiency this increase depends on modulation factors for anomalous protons and can even alter the direction of the radial anisotropy.     

On the non-occurrence of Type I X-ray bursts from the black hole candidates

It has been justifiably questioned if the black hole candidates (BHCs) have “hard surface” why Type I X-ray bursts are not seen from them [Narayan, R., Black holes in astrophysics, New J. Phys, 7, 199–218, 2005]. It is pointed out that a “physical surface” need not always be “hard” and could be “gaseous” in case the compact object is sufficiently hot [Mitra, A., The day of the reckoning: the value of the integration constant in the vacuum Schwarzschild solution, physics/0504076, p1–p6, 2005; Mitra, A., BHs or ECOs: A review of 90 years of misconceptions, in: Focus on Black Holes Research, Nova Science Pub., NY, p1–p94, 2005]. Even if a “hard surface” would be there, presence of strong intrinsic magnetic field could inhibit Type I X-ray burst from a compact object as is the case for Her X-1. Thus, non-occurrence of Type I bursts actually rules out those alternatives of BHs which are either non-magnetized or cold and, hence, is no evidence for existence of Event Horizons (EHs). On the other hand, from the first principle, we again show that the BHCs being uncharged and having finite masses cannot be BHs, because uncharged BHs have a unique mass M = 0. Thus the previous results that the so-called BHCs are actually extremely hot, ultramagnetized, Magnetospheric Eternally Collapsing Objects (ECOs) [Robertson, S., Leiter, D., Evidence for intrinsic magnetic moment in black hole candidates, Astrophys. J., 565, 447–451, (astro-ph/0102381), 2002 ; Robertson, S., Leiter, D., MECO model of galactic black hole candidates and active galactic nuclei, in: New Developments in Black Hole Research, Nova Science Pub., NY, p1–p44, astro-ph/0602453, 2005] rather than anything else get reconfirmed by non-occurrence of Type I X-ray bursts in BHCs.     

Solar wind dependence of the Pc1 wave activity

There are a host of factors influencing the excitation of Pc1 geomagnetic pulsations, which are ULF waves in the frequency range between 0.2 and 5 Hz. We have studied carefully the dependence of the pearl-type Pc1 activity at Sodankylä, Finland (L = 5.1) on the plasma density N in front of the magnetosphere, the bulk velocity V of the solar wind, and the intensity B of the IMF. The result is as follows: high values of N and reduced values of V are favorable to appearance of Pc1, whereas the dependence of Pc1 activity on B is practically absent. We also show that the probability of Pc1 occurrence decreases with the interplanetary electric field, and increases with solar wind impact pressure and with the plasma to magnetic pressure ratio “beta”.     

Heliospheric asymmetries due to the action of the interstellar magnetic field

We discuss the asymmetry of the heliospheric discontinuities obtained from the analysis of 3D modeling of the solar wind (SW) interaction with local interstellar medium (LISM). The flow of charged particles is governed by the ideal MHD equations and the flow of neutral particles is described by the Boltzmann equation. The emphasis is made on the asymmetries of the termination shock (TS) and the heliopause under the combined action of the interstellar and interplanetary magnetic fields (ISMF and IMF) in the presence of neutral hydrogen atoms whose transport through the heliosphere is modeled kinetically, using a Monte Carlo approach. We show that the deflection of neutral hydrogen flow from its original direction in the unperturbed LISM is highly anisotropic and evaluate a possible angle between the hydrogen deflection plane measured in the SOHO SWAN experiment and the plane containing the ISMF and LISM velocity vectors for different ISMF strengths. It is shown that the ISMF of a strength greater than 4 μG can account for the 10 AU difference in the TS heliocentric difference observed during its crossing by the Voyager 1 and Voyager 2 spacecraft, which however results in a larger discrepancy between the calculated and observed velocity distributions. The effect of a strong ISMF on the distribution of plasma quantities in the inner heliosheath and on 2–3 kHz radio emission is discussed.     

Two-velocity structure observed in the inner solar wind

Measurements of the motion of plasma density inhomogeneities in the inner solar wind are presented. The speeds were estimated using a cross-correlation analysis of radio frequency fluctuations of the Galileo spacecraft measured simultaneously at widely spaced ground stations. The radial projections of the correlation baselines on the pattern plane were of the order of several thousand kilometers. For cross-correlation functions calculated with comparatively short averaging times, we find that a pronounced two-velocity configuration is occasionally observed over the range of heliocentric distances 20 R < R < 40 R. The typical mean speed for such observations is about 300–400 km/s and the difference between the two predominant speeds is about 150–200 km/s. These results may indicate that the density fluctuations are associated with slow magnetosonic waves propagating in opposite directions at the local speed of sound in the reference frame moving with the mean solar wind speed. Quite reasonable estimates of the solar wind speed and speed of sound are obtained from this model. Another possible explanation of the two-velocity structures is that two independent solar wind streams are present simultaneously along different segments of the radio ray path.     

X-ray spectroscopy and variability of AGN detected in the 2 Ms Chandra Deep Field-North Survey

We investigate the nature of the faint X-ray source population through X-ray spectroscopy and variability analyses of 136 active galactic nuclei (AGN) detected in the 2 Ms Chandra Deep Field-North Survey with >200 background-subtracted 0.5–8.0 keV counts [F_{0.5–8.0 keV} = (1.4−200) × 10⁻¹⁵ erg cm⁻² s⁻¹]. Our preliminary spectral analyses yield median spectral parameters of Γ = 1.61 and intrinsic N_H = 6.2 × 10²¹ cm⁻² (z = 1 assumed when no redshift available) when the AGN spectra are fitted with a simple absorbed power-law model. However, considerable spectral complexity is apparent (e.g., reflection, partial covering) and must be taken into account to model the data accurately. Moreover, the choice of spectral model (i.e., free vs. fixed photon index) has a pronounced effect on the derived JVH distribution and, to a lesser extent, the X-ray luminosity distribution. We also find that among the 136 AGN, 10 (≈7%) show significant Fe K emission-line features with equivalent widths in the range 0.1–1.3 keV. Two of these emission-line AGN could potentially be Compton thick (i.e., Γ < 1.0 and large Fe K equivalent width). Finally, we find that 81 (≈60%) of the 136 AGN show signs of variability, and that this fraction increases significantly (≈80–90%) when better photon statistics are available.     

Near infrared studies of the massive X-ray binary 2S 0114+65

We present near-infrared spectroscopy of the massive X-ray binary 2S 0114+650. These observations covering the spectral range 1.08–2.35 μm span the region where Paschen and Brackett series recombination lines of hydrogen are expected to be seen, namely, Paβ , Brγ and Br 10–17 lines. The absence of any of these lines in emission supports earlier inferences that the optical component in 2S 0114+650 is unlikely to be a Be star but rather a B type supergiant. Near-IR photometry gives J = 8.78, H = 8.53 and K = 7.96; these values show marginal variations from earlier reported measurements.     

Solar wind-driven electron radiation belt response functions at 100-min time scales

We present a simple yet numerically robust technique, using autoregressive linear filters, to remove unwanted “colored noise” from solar wind and radiation belt electron data at sub-daily resolution. The remaining signal is then studied using finite impulse response linear prediction filters to represent the driven portion of the linear dynamics that describe the coupling between solar wind speed and electron flux. Sub-daily resolution response profiles covering magnetic L-shells between 1.1 and 8.0 R_E are presented which are consistent with daily resolution response functions. Namely, while there is strong global coherence governing electron flux dynamics, there are at least two distinct responses. The first response is an immediate dropout of electrons between L = 4 and L = 7 that is at least a partly adiabatic effect associated with enhancements in the ring current. This is followed by a 1–2 day delayed enhancement across the same L-shells that is likely a result of increased radial diffusion. The second response is an immediate enhancement seen between L = 3 and L = 4 with a typical duration of less than one day. Plausible explanations for this second response are briefly discussed, but neither empirical nor theoretical evidence can establish conclusively a definite physical cause. Finally, the response profiles show significant solar cycle and seasonal dependencies, indicating that better model output might be achieved with: (1) additional simultaneous solar wind inputs; (2) more sophisticated dynamical model structures capable of incorporating non-linear feedback; and/or (3) time-adaptive linear filters that can track non-stationary dynamics in time.     

Investigation of magnetospheric processes with the use of a source of strong magnetic field in the ionosphere

More than 20 years ago V.P. Shabansky suggested that the magnetic system installed aboard the satellite, could be used as a physical instrument for studying the processes which occur in the near Earth space. The corresponding space scales of an artificial “magnetosphere”—“magnisphere”—are 10 m in the experiment with relatively small magnets in the ionosphere and 100 m in the solar wind. The corresponding similarity criteria are estimated. The possible scheme of the experiment with a superconducting magnet (magnetic moment 10⁵ A · m²) installed aboard the satellite is considered. The experimental complex includes a number of systems for measuring the fluxes of charged particles in a wide energy range, DC electric and magnetic fields, the electromagnetic fields in different frequency bands (from X-rays to radio). The scientific objectives are discussed in detail.     

X-ray measurements of the dark matter distribution in clusters of galaxies with Chandra

We investigate the dark matter distributions in the central region of two clusters of galaxies (A1835 and MKW3S) using Chandra data. N-body simulations in the standard cold dark matter (CDM) model predict the dark matter distribution shows a cuspy dark matter profile: ρ(r) ∝ r, with in the range 1–2, while observations of dwarf and low surface brightness galaxies seem to favor the presence of a relatively flat core: 0 <  < 1. To investigate the dark matter distributions in the central region of clusters of galaxies, we analyze the Chandra data of A1835 and MKW3S with a deprojection method. We derive the mass profiles without the assumption of analytical models. We examine the inner slope of derived mass profiles assuming the dark matter profile is described with a power-law expression. The values of the slope are 0.95 ± 0.10 for A1835 and 1.33 ± 0.12 for MKW3S within the radius of 200 kpc. These are consistent with the result of the CDM simulations. However, within the radius of 100 kpc, the value of is less than unity for A1835 (0.47 ± 0.31). Our result implies that the central dark matter profile of some clusters cannot be described by CDM halos.     

Cosmic rays in the dynamic heliosheath

Cosmic ray modulation in the outer heliosphere is discussed from a modeling perspective. Emphasis is on the transport and acceleration of these particles at and beyond the solar wind termination shock in the inner heliosheath region and how this changes over a solar cycle. We will show that by using numerical models, and by comparing results to spacecraft observations, much can be learned about the dependence of cosmic ray modulation on solar cycle changes in the solar wind and heliospheric magnetic field. While the first determines the heliospheric geometry and shock structure, the latter results in a time-dependence of the transport coefficients. Depending on energy, both these effects contribute to cosmic ray intensities in the inner heliosheath changing over a solar cycle.     

Cosmic ray modulation in the outer heliosphere: Predictions for cosmic ray intensities up to the heliopause along Voyager 1 and 2 trajectories

Time dependent cosmic ray modulation in the outer heliosphere is calculated and results are compared to Voyager 1 and 2 observations using a two-dimensional time-dependent cosmic ray transport model. We predict possible future 133–242 MeV proton observations along the Voyager 1 and 2 spacecraft trajectories. Recent theoretical advances in cosmic ray transport parameters are introduced in order to provide a time-dependence for the assumed transport parameters used in the model. This leads to results that are in general compatible with the spacecraft observations in the inner and outer heliosphere over multiple solar cycles. However, for the outer heliosphere, we find that the Voyager 1 and 2 spacecraft observations cannot be fitted with an identical set of parameters along both trajectories. This indicates a possible asymmetric heliosphere or a symmetric heliosphere but with different diffusion parameters in the northern and southern hemispheres, respectively. Furthermore, results indicate that Voyager 2 observations are still under the influence of solar cycle related changes because of the large modulation volume between the heliopause and spacecraft location in contrast to Voyager 1 which shows a steady increase in cosmic ray intensities.     

Acceleration of interplanetary pick-up ions and anomalous cosmic rays

It may not be doubted anymore that anomalous cosmic rays (ACRs) are produced in the heliosphere from interplanetary pick-up ions through their acceleration at the solar wind termination shock. However, there is no general agreement in the community of heliospheric researchers concerning the mechanism of injection of the pick-up ions into the shock acceleration. We discuss here three possible ways for pick-up ions to be involved into the acceleration process at the termination shock: (1) preacceleration of pick-up ions in the whole region from the Sun up to the termination shock by solar wind turbulences and interplanetary shock waves, (2) local preacceleration of pick-up ions in a vicinity of the termination shock by shock surfing, and (3) formation of high-velocity tails in pick-up ion spectra consisting of secondary pick-up ions which are produced in the supersonic solar wind due to ionization of energetic neutral atoms entering from the inner heliosheath.