The solar modulation of electrons in the heliosphere

The propagation and modulation of electrons in the heliosphere play an important part in improving our understanding and assessment of the modulation processes. A full three-dimensional numerical model is used to study the modulation of galactic electrons, from Earth into the inner heliosheath, over an energy range from 10 MeV to 30 GeV. The modeling is compared with observations of 6–14 MeV electrons from Voyager 1 and observations at Earth from the PAMELA mission. Computed spectra are shown at different spatial positions. Based on comparison with Voyager 1 observations, a new local interstellar electron spectrum is calculated. We find that it consists of two power-laws: In terms of kinetic energy E, the results give E ^?1.5 below ～500 MeV and E ^?3.15 at higher energies. Radial intensity profiles are computed also for 12 MeV electrons, including a Jovian source, and compared to the 6–14 MeV observations from Voyager 1. Since the Jovian and galactic electrons can be separated in the model, we calculate the intensity of galactic electrons below 100 MeV at Earth. The highest possible differential flux of galactic electrons at Earth with E=12 MeV is found to have a value of 2.5×10^?1 electrons m^?2?s^?1?sr^?1?MeV^?1 which is significantly lower (a factor of 3) than the Jovian electron flux at Earth. The model can also reproduce the extraordinary increase of electrons by a factor of 60 at 12 MeV in the inner heliosheath. A lower limit for the local interstellar spectrum at 12 MeV is estimated to have a value of (90±10) electrons m^?2?s^?1?sr^?1?MeV^?1.     

Measurement of the radial density gradient of cosmic ray in the heliosphere by the GRAPES-3 experiment

A radial anisotropy in the flux of cosmic rays in heliosphere was theoretically predicted by Parker and others within the framework of the diffusion–convection mechanism. The solar wind is responsible for sweeping out the galactic cosmic rays, creating a radial density gradient within the heliosphere. This gradient coupled with the interplanetary magnetic field induces a flow of charged particles perpendicular to the ecliptic plane which was measured and correctly explained by Swinson, and is hereafter referred as ‘Swinson flow’. The large area GRAPES-3 tracking muon telescope offers a powerful probe to measure the Swinson flow and the underlying radial density gradient of the galactic cosmic rays at a relatively high rigidity of ∼100 GV. The GRAPES-3 data collected over a period of six years (2000–2005) were analyzed and the amplitude of the Swinson flow was estimated to be (0.0644 ± 0.0008)% of cosmic ray flux which was an ∼80σ effect. The phase of the maximum flow was at a sidereal time of (17.70 ± 0.05) h which was 18 min earlier than the expected value of 18 h. This small 18 min phase difference had a significance of ∼6σ indicating the inherent precision of the GRAPES-3 measurement. The radial density gradient of the galactic cosmic rays at a median rigidity of 77 GV was found to be 0.65% AU⁻¹.     

An investigation of the field-aligned currents associated with a large-scale ULF wave in the morning sector

Previous work by Scoffield, H.C., Yeoman, T.K., Wright, D.M., Milan, S.E., Wright, A.N., Strangeway, R.J. [2005. An investigation of the field aligned currents associated with a large scale ULF wave using data from CUTLASS and FAST. Ann. Geophys. 23, 487–498) investigated a large-scale ULF wave, occurring in the dusk sector (∼1900 MLT). The wave had a period of ∼800 s (corresponding to 1.2 mHz frequency), an azimuthal wave number of ∼7 and a full-width at half-maximum (FWHM) across the resonance of 350 km. IMAGE ground magnetometer and SuperDARN radar observations of the wave's spatial and temporal characteristics were used to parameterise a simple, two-dimensional field line resonance (FLR) model. The model-calculated field-aligned current (FAC) was compared with FACs derived from the FAST energetic particle spectra and magnetic field measurement. Here the authors use the same method to investigate the FAC structure of a second large-scale ULF wave, with a period of ∼450 s, occurring the dawn sector (∼0500 MLT) with an opposite sense background region 1–region 2 current system. This wave has a much larger longitudinal scale (m∼4.5) and a smaller latitude scale (FWHM=150 km). Unlike the dusk sector wave, which was dominated by upward FAC, FAST observations of the dawn sector wave show an interval of large-scale downward FAC of ∼1.5 μA m⁻². Downgoing magnetospheric electrons with energies of a few keV were observed, which are associated with upward FACs of ∼1 μA m⁻². For both wave studies, downward currents appear to be carried partially by upgoing electrons below the FAST energy detection threshold (5 eV), but also consist of a mixture of hotter downgoing magnetospheric electrons and upgoing ionospheric electrons of energies 30 eV–1 keV. Strong intervals of upward current show that small-scale structuring of scale ∼50 km has been imposed on the current carriers. In general, this study confirms the findings of Scoffield, H.C., Yeoman, T.K., Wright, D.M., Milan, S.E., Wright, A.N., Strangeway, R.J. [2005. An investigation of the FACs associated with a large-scale ULF wave using data from CUTLASS and FAST. Ann. Geophys. 23, 487–498).     

Sources of stellar energy,Einstein–Eddington timescale of gravitational contraction and eternally collapsing objects

We point out that although conventional stars are primarily fed by burning of nuclear fuel at their cores, in a strict sense, the process of release of stored gravitational energy, known as, Kelvin–Helmholtz (KH) process is either also operational albeit at an arbitrary slow rate, or lying in wait to take over at the disruption of the nuclear channel. In fact, the latter mode of energy release is the true feature of any self-gravity bound object including stars. We also highlight the almost forgotten fact that Eddington was the first physicist to introduce special relativity into the problem and correctly insist that, actually, total energy stored in a star is not the mere Newtonian energy but the total mass energy (E = Mc²). Accordingly, Eddington defined an “Einstein time scale” of Evolution where the maximum age of the Sun turned out to be t_E ≈ 1.4 × 10¹³ yr. This concept has a fundamental importance though we know now that Sun in its present form cannot survive for more than 10 billion years. We extend this concept by introducing general relativity and show that the minimum value of depletion of total mass–energy is t_E = ∞ not only for Sun but for and sufficiently massive or dense object. We propose that this time scale be known in the name of “Einstein–Eddington”. We also point out that, recently, it has been shown that as massive stars undergo continued collapse to become a Black Hole, first they become extremely relativistic radiation pressure supported stars. And the life time of such relativistic radiation pressure supported compact stars is indeed dictated by this Einstein–Eddington time scale whose concept is formally developed here. Since this observed time scale of this radiation pressure supported quasistatic state turns out to be infinite, such objects are called eternally collapsing objects (ECO). Further since ECOs are expected to have strong intrinsic magnetic field, they are also known as “Magnetospheric ECO” or MECO.     

Gamma–hadron separation methods for the VERITAS array of four imaging atmospheric Cherenkov telescopes

Ground-based arrays of imaging atmospheric Cherenkov telescopes have emerged as the most sensitive γ-ray detectors in the energy range of about 100 GeV and above. The strengths of these arrays are a very large effective collection area on the order of 10⁵ m², combined with excellent single photon angular and energy resolutions. The sensitivity of such detectors is limited by statistical fluctuations in the number of Cosmic-ray initiated air showers that resemble γ-ray air showers in many ways. In this paper, we study the performance of simple event reconstruction methods when applied to simulated data of the Very Energetic Radiation Imaging Telescope Array System (VERITAS) experiment. We review methods for reconstructing the arrival direction and the energy of the primary photons, and examine means to improve on their performance. For a software threshold energy of 300 GeV (100 GeV), the methods achieve point source angular and energy resolutions of σ_63% = 0.1° (0.2°) and σ_68% = 15% (22%), respectively. The main emphasis of the paper is the discussion of γ–hadron separation methods for the VERITAS experiment. We find that the information from several methods can be combined based on a likelihood ratio approach and the resulting algorithm achieves a γ–hadron suppression with a quality factor that is substantially higher than that achieved with the standard methods used so far.     

The major upgrade of the MAGIC telescopes,Part II: A performance study using observations of the Crab Nebula

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma, Spain. During summer 2011 and 2012 it underwent a series of upgrades, involving the exchange of the MAGIC-I camera and its trigger system, as well as the upgrade of the readout system of both telescopes. We use observations of the Crab Nebula taken at low and medium zenith angles to assess the key performance parameters of the MAGIC stereo system. For low zenith angle observations, the standard trigger threshold of the MAGIC telescopes is ∼ 50  GeV. The integral sensitivity for point-like sources with Crab Nebula-like spectrum above 220 GeV is (0.66 ± 0.03)% of Crab Nebula flux in 50 h of observations. The angular resolution, defined as the σ of a 2-dimensional Gaussian distribution, at those energies is ≲ 0.07°, while the energy resolution is 16%. We also re-evaluate the effect of the systematic uncertainty on the data taken with the MAGIC telescopes after the upgrade. We estimate that the systematic uncertainties can be divided in the following components: < 15% in energy scale, 11%–18% in flux normalization and ± 0.15 for the energy spectrum power-law slope.     

Comparison of Fermi-LAT and CTA in the region between 10–100 GeV

The past decade has seen a dramatic improvement in the quality of data available at both high (HE: 100 MeV to 100 GeV) and very high (VHE: 100 GeV to 100 TeV) gamma-ray energies. With three years of data from the Fermi Large Area Telescope (LAT) and deep pointed observations with arrays of Cherenkov telescope, continuous spectral coverage from 100 MeV to ∼10 TeV exists for the first time for the brightest gamma-ray sources. The Fermi-LAT is likely to continue for several years, resulting in significant improvements in high energy sensitivity. On the same timescale, the Cherenkov Telescope Array (CTA) will be constructed providing unprecedented VHE capabilities. The optimisation of CTA must take into account competition and complementarity with Fermi, in particularly in the overlapping energy range 10–100 GeV. Here we compare the performance of Fermi-LAT and the current baseline CTA design for steady and transient, point-like and extended sources.     

A Study of Potential Earth-like Planets in the GJ 876 Planetary System

We perform numerical simulations to investigate potential Earth-like planets in the GJ 876 planetary system. We show that the secular resonances ν₁ and ν₂ (resulting respectively from the inner and outer giant planets) can excite the eccentricities of Earth-like planets with orbits 0.21 AU ≤ a < 0.50 AU and cause them to be ejected out of the system in a short time. However, in the dynamical sense, Earth-like planets potentially exist in the region 0.50 AU ≤ a ≤ 1.00 AU, in stable low-eccentricity orbits which may last up to 10⁵ yr.     

The extragalactic background light and the gamma-ray opacity of the universe

The extragalactic background light (EBL) is one of the fundamental observational quantities in cosmology. All energy releases from resolved and unresolved extragalactic sources, and the light from any truly diffuse background, excluding the cosmic microwave background (CMB), contribute to its intensity and spectral energy distribution. It therefore plays a crucial role in cosmological tests for the formation and evolution of stellar objects and galaxies, and for setting limits on exotic energy releases in the universe. The EBL also plays an important role in the propagation of very high energy γ-rays which are attenuated en route to Earth by pair producing γ–γ interactions with the EBL and CMB. The EBL affects the spectrum of the sources, predominantly blazars, in the ∼10 GeV–10 TeV energy regime. Knowledge of the EBL intensity and spectrum will allow the determination of the intrinsic blazar spectrum in a crucial energy regime that can be used to test particle acceleration mechanisms and very high energy (VHE) γ-ray production models. Conversely, knowledge of the intrinsic γ-ray spectrum and the detection of blazars at increasingly higher redshifts will set strong limits on the EBL and its evolution. This paper reviews the latest developments in the determination of the EBL and its impact on the current understanding of the origin and production mechanisms of γ-rays in blazars, and on energy releases in the universe. The review concludes with a summary and future directions in Cherenkov Telescope Array techniques and in infrared ground-based and space observatories that will greatly improve our knowledge of the EBL and the origin and production of very high energy γ-rays.     

Robust limits on Lorentz violation from gamma-ray bursts

We constrain the possibility of a non-trivial refractive index in free space corresponding to an energy-dependent velocity of light: c(E) ≃ c₀(1 − E/M), where M is a mass scale that might represent effect of quantum-gravitational space-time foam, using the arrival times of sharp features observed in the intensities of radiation with different energies from a large sample of gamma-ray bursters (GRBs) with known redshifts. We use wavelet techniques to identify genuine features, which we confirm in simulations with artificial added noise. Using the weighted averages of the time-lags calculated using correlated features in all the GRB light curves, we find a systematic tendency for more energetic photons to arrive earlier. However, there is a very strong correlation between the parameters characterizing an intrinsic time-lag at the source and a distance-dependent propagation effect. Moreover, the significance of the earlier arrival times is less evident for a subsample of more robust spectral structures. Allowing for intrinsic stochastic time-lags in these features, we establish a statistically robust lower limit: M > 0.9 × 10¹⁶ GeV on the scale of violation of Lorentz invariance.     

The Contribution from Radio Quiet Quasars to γ-ray Background

The EGRET observations have confirmed and proposed the new isotropically distributed γ-ray background, but the known objects radiating γ-rays can not supply so much radiation. Meanwhile, EGRET also reveals a population of γ-ray sources with no radio counterparts which are isotropically distributed in the sky, indicating their possible cosmological origins. Wang et al.^[13] proposed a new γ-ray radiation process driven by the radiation feedback of AGNs. The energy of the radiation peaks around 1 GeV—0.1 TeV with the typical luminosity of 10⁴²—10⁴³ ergs · s^?1. This kind of radiation process in the radio quiet quasars make them the potential γ-ray radiation sources as well as the contributors to the γ-ray background. We consider two cases in which the seed photons in the inverse-Compton processes are from the accretion disks of quasars and cosmic microwave background (CMB), respectively. We find that the former contributes 78%—92% of the background radiation around 1 GeV, while the contribution from the latter is negligible. The radio quiet quasars are highly likely to become the objects which contribute the most energy to the γ-ray background around 1 GeV.     

All-particle cosmic ray energy spectrum measured with 26 IceTop stations

We report on a measurement of the cosmic ray energy spectrum with the IceTop air shower array, the surface component of the IceCube Neutrino Observatory at the South Pole. The data used in this analysis were taken between June and October, 2007, with 26 surface stations operational at that time, corresponding to about one third of the final array. The fiducial area used in this analysis was 0.122 km². The analysis investigated the energy spectrum from 1 to 100 PeV measured for three different zenith angle ranges between 0° and 46°. Because of the isotropy of cosmic rays in this energy range the spectra from all zenith angle intervals have to agree. The cosmic-ray energy spectrum was determined under different assumptions on the primary mass composition. Good agreement of spectra in the three zenith angle ranges was found for the assumption of pure proton and a simple two-component model. For zenith angles θ < 30°, where the mass dependence is smallest, the knee in the cosmic ray energy spectrum was observed at about 4 PeV, with a spectral index above the knee of about −3.1. Moreover, an indication of a flattening of the spectrum above 22 PeV was observed.     

Milagro limits and HAWC sensitivity for the rate-density of evaporating Primordial Black Holes

Primordial Black Holes (PBHs) are gravitationally collapsed objects that may have been created by density fluctuations in the early universe and could have arbitrarily small masses down to the Planck scale. Hawking showed that due to quantum effects, a black hole has a temperature inversely proportional to its mass and will emit all species of fundamental particles thermally. PBHs with initial masses of ∼5.0 × 10¹⁴ g should be expiring in the present epoch with bursts of high-energy particles, including gamma radiation in the GeV–TeV energy range. The Milagro high energy observatory, which operated from 2000 to 2008, is sensitive to the high end of the PBH evaporation gamma-ray spectrum. Due to its large field-of-view, more than 90% duty cycle and sensitivity up to 100 TeV gamma rays, the Milagro observatory is well suited to perform a search for PBH bursts. Based on a search on the Milagro data, we report new PBH burst rate density upper limits over a range of PBH observation times. In addition, we report the sensitivity of the Milagro successor, the High Altitude Water Cherenkov (HAWC) observatory, to PBH evaporation events.     

Outer Main Belt asteroids: Identification and distribution of four 3-μm spectral groups

This paper examines the distribution and the abundance of hydrated minerals (any mineral that contains H₂O or OH) on outer Main Belt asteroids spanning the 2.5 < a < 4.0 AU region. The hypothesis we are testing is whether planetesimals that accreted closer to the Sun experienced a higher degree of aqueous alteration. We would expect then to see a gradual decline of the abundance of hydrated minerals among the outer Main Belt asteroids with increasing heliocentric distance (2.5 < a < 4.0 AU). We measured spectra (0.8–2.5 μm and 1.9–4.1 μm) of 28 outer Main Belt asteroids using the SpeX spectrograph/imager at the NASA Infrared Telescope Facility (IRTF). We identified four groups on the basis of the shape and the band center of the 3-μm feature. The first group, which we call “sharp”, exhibits a sharp 3-μm feature, attributed to hydrated minerals (phyllosilicates). Most asteroids in this group are located in the 2.5 < a < 3.3 AU region. The second group, which we call “Ceres-like”, consists of 10 Hygiea and 324 Bamberga. Like Asteroid Ceres, these asteroids exhibit a 3-μm feature with a band center of 3.05 ± 0.01 μm that is superimposed on a broader absorption feature from ～2.8 to 3.7 μm. The third group, which we call “Europa-like”, includes 52 Europa, 31 Euphrosyne, and 451 Patientia. Objects in this group exhibit a 3-μm feature with a band center of 3.15 ± 0.01 μm. Both the Ceres-like and Europa-like groups are concentrated in the 2.5 < a < 3.3 AU region. The fourth group, which we call “rounded”, is concentrated in the 3.4 < a < 4.0 AU region. Asteroids in this group are characterized by a rounded 3-μm feature, attributed to H₂O ice. A similar rounded 3-μm feature was also identified in 24 Themis and 65 Cybele. Unlike the sharp group, the rounded group did not experience aqueous alteration. Of the asteroids observed in this study, 140 Siwa, a P-type, is the only one that does not exhibit a 3-μm feature. These results are important to constrain the nature and the degree of aqueous alteration in outer Main Belt asteroids.     

Search for relativistic magnetic monopoles with the Baikal neutrino telescope

We present results of a search for relativistic magnetic monopoles with the Baikal neutrino telescope NT200, using data taken between April 1998 and February 2003. No monopole candidates have been found. We set an upper limit 4.6 × 10⁻¹⁷ cm⁻² s⁻¹ sr⁻¹ for the flux of monopoles with β_m = 1. This is a factor of 20 below the Chudakov–Parker bound which is inferred from the very existence of large-scale galactic magnetic fields.     

New wind measurements in Venus’ lower mesosphere from visible spectroscopy

The dynamics of Venus’ mesosphere (70–110 km) is characterized by the superposition of two different wind regimes: (1) Venus’ retrograde superrotation; (2) a sub-solar to anti-solar (SS–AS) flow pattern, driven by solar EUV heating on the sunlit hemisphere. Here, we report on new ground-based velocity measurements in the lower part of the mesosphere. We took advantage of two essentially symmetric Venus elongations in 2001 and 2002 to perform high-resolution Doppler spectroscopy (R=120,000) in ¹²C¹⁶O₂ visible lines of the 5ν₃ band and in a few solar Fraunhofer lines near 8700 Å. These measurements, mapped over several points on Venus’ illuminated hemisphere, probe the region of cloud tops. More precisely, the solar Fraunhofer lines sample levels a few kilometers below the UV features (i.e. near ∼67 km), while the CO₂ lines probe an altitude higher by about 7 km. The wind field over Venus’ disk is retrieved with an rms uncertainty of 15–25 m s⁻¹ on individual measurements. Kinematical fit to a one- or two-component circulation model indicates the dominance of the zonal retrograde flow with a mean equatorial velocity of ∼75 m s⁻¹, exhibiting very strong day-to-day variations (±65 m s⁻¹). Results are very consistent for the two kinds of lines, suggesting a negligible vertical wind shear over 67–74 km. The SS–AS flow is not detected in single-day observations, but combining the results from all data suggests that this component may invade the lower mesosphere with a ∼40 m s⁻¹ velocity.     

The outcome of the protoplanetary disk of very massive stars

We suggest that planets, brown dwarfs, and even low mass stars can be formed by fragmentation of protoplanetary disks around very massive stars (M ? 100 M_⊙). We discuss how fragmentation conditions make the formation of very massive planetary systems around very massive stars favorable. Such planetary systems are likely to be composed of brown dwarfs and low mass stars of ～0.1–0.3 M_⊙, at orbital separations of ～ few × 100–10⁴ AU. In particular, scaling from solar-like stars suggests that hundreds of Mercury-like planets might orbit very massive stars at ～10³ AU where conditions might favor liquid water. Such fragmentation objects can be excellent targets for the James Webb Space Telescope and other large telescopes working in the IR bands. We predict that deep observations of very massive stars would reveal these fragmentation objects, orbiting in the same orbital plane in cases where there are more than one object.     

The orbital distribution of Near-Earth Objects inside Earth’s orbit

Canada’s Near-Earth Object Surveillance Satellite (NEOSSat), set to launch in early 2012, will search for and track Near-Earth Objects (NEOs), tuning its search to best detect objects with a < 1.0 AU. In order to construct an optimal pointing strategy for NEOSSat, we needed more detailed information in the a < 1.0 AU region than the best current model (Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J.M., Levison, H.F., Michel, P., Metcalfe, T.S. [2002]. Icarus 156, 399–433) provides. We present here the NEOSSat-1.0 NEO orbital distribution model with larger statistics that permit finer resolution and less uncertainty, especially in the a < 1.0 AU region. We find that Amors = 30.1 ± 0.8%, Apollos = 63.3 ± 0.4%, Atens = 5.0 ± 0.3%, Atiras (0.718 < Q < 0.983 AU) = 1.38 ± 0.04%, and Vatiras (0.307 < Q < 0.718 AU) = 0.22 ± 0.03% of the steady-state NEO population. Vatiras are a previously undiscussed NEO population clearly defined in our integrations, whose orbits lie completely interior to that of Venus. Our integrations also uncovered the unexpected production of retrograde orbits from main-belt asteroid sources; this retrograde NEA population makes up ?0.1% of the steady-state NEO population. The relative NEO impact rate onto Mercury, Venus, and Earth, as well as the normalized distribution of impact speeds, was calculated from the NEOSSat-1.0 orbital model under the assumption of a steady-state. The new model predicts a slightly higher Mercury impact flux.     

Spectral behavior of the symbiotic nova AG Pegasi observed with IUE and HST

Ultraviolet spectra from the International Ultraviolet Explorer (IUE) and from the Hubble Space Telescope (HST) of the symbiotic novae AG Peg during the period 1978–1996 are analyzed. Some spectra showing the modulations of spectral lines at different times are presented. We determined the reddening from the 2200 Å feature, finding that E(B−V) = 0.10 ± 0.02. We studied N IV] at 1486 Å, C IV 1550 Å, and O III] at 1660 Å, produced in the fast wind from the hot white dwarf. The mean wind velocity of the three emission lines is 1300 km s⁻¹ (FWHM). The mean wind mass loss rate is ∼6 × 10⁻⁷ M_⊙ yr⁻¹. The mean temperature is ∼6.5 × 10⁵ K. The mean ultraviolet luminosity is ∼5 × 10³³ erg s⁻¹. The modulations of line fluxes in the emitting region at different times are attributed to the variations of density and temperature of the ejected matter as a result of variations in the rate of mass loss.