Pulsar simulations for the Fermi Large Area Telescope

Pulsars are among the prime targets for the Large Area Telescope (LAT) aboard the recently launched Fermi observatory. The LAT will study the gamma-ray Universe between 20 MeV and 300 GeV with unprecedented detail. Increasing numbers of gamma-ray pulsars are being firmly identified, yet their emission mechanisms are far from being understood. To better investigate and exploit the LAT capabilities for pulsar science, a set of new detailed pulsar simulation tools have been developed within the LAT collaboration. The structure of the pulsar simulator package (PulsarSpectrum) is presented here. Starting from photon distributions in energy and phase obtained from theoretical calculations or phenomenological considerations, gamma-rays are generated and their arrival times at the spacecraft are determined by taking into account effects such as barycentric effects and timing noise. Pulsars in binary systems also can be simulated given orbital parameters. We present how simulations can be used for generating a realistic set of gamma-rays as observed by the LAT, focusing on some case studies that show the performance of the LAT for pulsar observations.     

Study of X-ray and γ-ray flare from Mkn-421 during February 2010 using Swift and Fermi observations

In this work, we report on the intense flaring activity from Mkn-421 in X-ray and γ-ray regimes simultaneously observed by Swift-XRT/BAT and Fermi-LAT satellite telescopes in February 2010. With the aim of understanding the underlying physics of the flaring state in Mkn-421, we have performed a detailed spectral analysis of Swift/XRT and Fermi/LAT observations of Mkn-421 during February 12–25, 2010 (MJD 55239–55252). Over this period, we study the daily light curves and spectral variability of the source in 1–10 keV, 0.1–1 GeV and 1–100 GeV energy bands. We have performed the spectral analysis of Swift-XRT and Fermi/LAT observations to study the spectral evolution in the X-ray and gamma-ray energy domains respectively. We also compute the fractional variability amplitude in both the energy bands during the above period. We study trends between spectral parameters and physical insights provided by the parameter responsible for X-ray and γ-ray emission from the source. We search for energetic features phenomenologically linked to the single zone SSC model for blazar emission. We also produce the broad band SED with a leptonic single zone SSC model for the source.     

IACT observations of gamma-ray bursts: prospects for the Cherenkov Telescope Array

Gamma rays at rest frame energies as high as 90 GeV have been reported from gamma-ray bursts (GRBs) by the Fermi Large Area Telescope (LAT). There is considerable hope that a confirmed GRB detection will be possible with the upcoming Cherenkov Telescope Array (CTA), which will have a larger effective area and better low-energy sensitivity than current-generation imaging atmospheric Cherenkov telescopes (IACTs). To estimate the likelihood of such a detection, we have developed a phenomenological model for GRB emission between 1 GeV and 1 TeV that is motivated by the high-energy GRB detections of Fermi-LAT, and allows us to extrapolate the statistics of GRBs seen by lower energy instruments such as the Swift-BAT and BATSE on the Compton Gamma-ray Observatory. We show a number of statistics for detected GRBs, and describe how the detectability of GRBs with CTA could vary based on a number of parameters, such as the typical observation delay between the burst onset and the start of ground observations. We also consider the possibility of using GBM on Fermi as a finder of GRBs for rapid ground follow-up. While the uncertainty of GBM localization is problematic, the small field-of-view for IACTs can potentially be overcome by scanning over the GBM error region. Overall, our results indicate that CTA should be able to detect one GRB every 20–30 months with our baseline instrument model, assuming consistently rapid pursuit of GRB alerts, and provided that spectral breaks below ~100 GeV are not a common feature of the bright GRB population. With a more optimistic instrument model, the detection rate can be as high as 1 to 2 GRBs per year.     

Minimal spanning tree algorithm for γ-ray source detection in sparse photon images: cluster parameters and selection strategies

The minimal spanning tree (MST) algorithm is a graph-theoretical cluster-finding method. We previously applied it to γ-ray bidimensional images, showing that it is quite sensitive in finding faint sources. Possible sources are associated with the regions where the photon arrival directions clusterize. MST selects clusters starting from a particular “tree” connecting all the point of the image and performing a cut based on the angular distance between photons, with a number of events higher than a given threshold. In this paper, we show how a further filtering, based on some parameters linked to the cluster properties, can be applied to reduce spurious detections. We find that the most efficient parameter for this secondary selection is the magnitude M of a cluster, defined as the product of its number of events by its clustering degree. We test the sensitivity of the method by means of simulated and real Fermi-Large Area Telescope (LAT) fields. Our results show that $\sqrt{M}$ is strongly correlated with other statistical significance parameters, derived from a wavelet based algorithm and maximum likelihood (ML) analysis, and that it can be used as a good estimator of statistical significance of MST detections. We apply the method to a 2-year LAT image at energies higher than 3 GeV, and we show the presence of new clusters, likely associated with BL Lac objects.     

A Heliospheric Imager for Deep Space: Lessons Learned from Helios, SMEI, and STEREO

The zodiacal-light photometers on the twin Helios spacecraft, the Solar Mass Ejection Imager (SMEI) on the Coriolis spacecraft, and the Heliospheric Imagers (HIs) on the Solar-TErrestrial RElations Observatory (STEREO) twin spacecraft all point the way to optimizing future remote-sensing Thomson-scattering observations from deep space. Such data could be provided by wide-angle viewing instruments on Solar Orbiter, Solar Probe, or other deep-space probes. Here, we present instrument specifications required for a successful heliospheric imager, and the measurements and data-processing steps that make the best use of this remote-sensing system. When this type of instrument is properly designed and calibrated, its data are capable of determining zodiacal-dust properties, and of three-dimensional reconstructions of heliospheric electron density over large volumes of the inner heliosphere. Such systems can measure fundamental properties of the inner heliospheric plasma, provide context for the in-situ monitors on board spacecraft, and enable physics-based analyses of this important segment of the Sun-spacecraft connection.     

Exploring the Capabilities of the Anti-Coincidence Shield of the International Gamma-Ray Astrophysics Laboratory (INTEGRAL) Spectrometer to Study Solar Flares

The International Gamma-Ray Astrophysics Laboratory (INTEGRAL) is a European Space Agency hard X-ray/γ-ray observatory for astrophysics, covering photon energies from 15 keV to 10 MeV. It was launched in 2002, and since then the Bismuth Germanate (BGO) detectors of the Anti-Coincidence Shield (ACS) of the Spectrometer on INTEGRAL (SPI) have detected many hard X-ray (HXR) bursts from the Sun, producing light curves at photon energies above ≈?100 keV. The spacecraft has a highly elliptical orbit, providing long uninterrupted observing (about 90 % of the orbital period) with nearly constant background due to the shorter time needed to cross Earth’s radiation belts. However, because of technical constraints, INTEGRAL cannot be pointed at the Sun, and high-energy solar photons are always detected in nonstandard observation conditions. To make the data useable for solar studies, we have undertaken a major effort to specify the observing conditions through Monte Carlo simulations of the response of ACS for several selected flares. We checked the performance of the model employed for the Monte Carlo simulations using the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI) observations for the same sample of solar flares. We conclude that although INTEGRAL was not designed to perform solar observations, ACS is a useful instrument for solar-flare research. In particular, its relatively large effective area allows determining good-quality HXR/γ-ray light curves for X- and M-class solar flares and, in some cases, probably also for C-class flares.     

Uncovering the Bias in Low-Degree p-Mode Linewidth Fitting

Obtaining reliable estimates of linewidths in the power spectra of low-degree p modes is problematic at low frequency. In this regime, the mode coherence time increases with decreasing frequency, often causing the modes to be unresolved in relatively long duration spectra. The signal-to-noise ratio is also less favourable at low frequency, resulting in fits to power spectra underestimating the true linewidth of the p modes owing to the tails of the Lorentzian peaks becoming dominated by the background noise. We use a numerical simulation approach to assess the effect of this bias on the fitted widths of p-mode peaks and calculate observational duration limits required to obtain an unbiased estimate of the p-mode linewidth as a function of frequency. This is done in four different cases, where the precision of the artificial data is set at 0.25, 0.50, 0.75, and 1.00 m?s^?1 by adding random?scatter to increase the sample standard deviation per 40-second measurement. In all cases, the observational duration required to accurately obtain width estimates increases beyond that required for sufficient spectral resolution below a certain threshold frequency. For modes at ≈?1500 μHz, with an amplitude of approximately ten times the background, observations of up to 972 days are required to obtain an unbiased estimate of the linewidth. This is equivalent to ≈?18 times the coherence time of the corresponding p modes.     

Gamma ray bursts

Gamma-ray bursts have been detected at photon energies up to tens of GeV. We review some recent developments in the X-ray to GeV photon phenomenology in the light of Swift and Fermi observations, and some of the theoretical models developed to explain them, with a view towards implications for C.T.A.     

Multi-spacecraft Study of the 8 November 2000 SEP Event: Electron Injection Histories 100° Apart

We present the analysis of a large solar near-relativistic electron event observed by the Ulysses and the ACE spacecraft on 8 November 2000, when Ulysses was located at a heliocentric distance of 2.4 AU and at a heliographic latitude of ??80° S. We use a particle propagation model to infer the local interplanetary transport conditions and the injection histories of the near-relativistic electrons observed by both spacecraft. We find different local transport conditions for each set of observations. The inferred injection profiles for both spacecraft extend for several hours; but the injection at Ulysses was smaller and started later. The association with type II radio emission suggests that the heliospheric electrons were provided by coronal shock acceleration. An analysis of the in situ magnetic field and plasma measurements indicates that the global configuration of the heliosphere (disturbed by transient structures) could play a role in shaping the characteristics of solar energetic particle events observed from different locations.     

Long-term observations of the evolution of the southern seasonal cap of Mars: Neutron measurements by the HEND instrument onboard the 2001 <Emphasis Type="Italic">Mars Odyssey</Emphasis> spacecraft

We present the results of five-year observations of the southern seasonal cap of Mars based on neutron spectroscopy of the surface fulfilled by the Russian HEND instrument onboard the NASA 2001 Mars Odyssey spacecraft. The numerical modeling of the observational data allowed us to reconstruct the curves of the variations of the total mass of the southern seasonal cap of Mars for different years (three Martian years) and to find the year-to-year variations of the seasonal cycle.     

Inclusion of In-Situ Velocity Measurements into the UCSD Time-Dependent Tomography to Constrain and Better-Forecast Remote-Sensing Observations

The University of California, San Diego (UCSD) three-dimensional (3-D) time-dependent tomography program has been used successfully for a decade to reconstruct and forecast coronal mass ejections from interplanetary scintillation observations. More recently, we have extended this tomography technique to use remote-sensing data from the Solar Mass Ejection Imager (SMEI) on board the Coriolis spacecraft; from the Ootacamund (Ooty) radio telescope in India; and from the European Incoherent SCATter (EISCAT) radar telescopes in northern Scandinavia. Finally, we intend these analyses to be used with observations from the Murchison Widefield Array (MWA), or the LOw Frequency ARray (LOFAR) now being developed respectively in Australia and Europe. In this article we demonstrate how in-situ velocity measurements from the Advanced Composition Explorer (ACE) space-borne instrumentation can be used in addition to remote-sensing data to constrain the time-dependent tomographic solution. Supplementing the remote-sensing observations with in-situ measurements provides additional information to construct an iterated solar-wind parameter that is propagated outward from near the solar surface past the measurement location, and throughout the volume. While the largest changes within the volume are close to the radial directions that incorporate the in-situ measurements, their inclusion significantly reduces the uncertainty in extending these measurements to global 3-D reconstructions that are distant in time and space from the spacecraft. At Earth, this can provide a finely-tuned real-time measurement up to the latest time for which in-situ measurements are available, and enables more-accurate forecasting beyond this than remote-sensing observations alone allow.     

GLAST large area telescope multiwavelength planning

Gamma-ray astrophysics depends in many ways on multiwavelength studies. The Gamma-ray Large Area Space Telescope (GLAST) Large Area Telescope (LAT) Collaboration has started multiwavelength planning well before the scheduled 2007 launch of the observatory. Some of the high-priority multiwavelength needs include: (1) availability of contemporaneous radio and X-ray timing of pulsars; (2) expansion of blazar catalogs, including redshift measurements; (3) improved observations of molecular clouds, especially at high galactic latitudes; (4) simultaneous broad-band blazar monitoring; (5) characterization of gamma-ray transients, including gamma ray bursts; (6) radio, optical, X-ray and TeV counterpart searches for reliable and effective sources identification and characterization. Several of these activities are needed to be in place before launch.      

Masses of the Main Asteroid Belt and the Kuiper Belt from the Motions of Planets and Spacecraft

Dynamicalmass estimates for the main asteroid belt and the trans-Neptunian Kuiper belt have been found from their gravitational influence on the motion of planets. Discrete rotating models consisting ofmovingmaterial points have been used tomodel the total attraction fromsmall or as yet undetected bodies of the belts. The masses of the model belts have been included in the set of parameters being refined and determined and have been obtained by processing more than 800 thousand modern positional observations of planets and spacecraft. We have processed the observations and determined the parameters based on the new EPM2017 version of the IAA RAS planetary ephemerides. The large observed radial extent of the belts (more than 1.2 AU for the main belt and more than 8 AU for the Kuiper belt) and the concentration of bodies in the Kuiper belt at a distance of about 44 AU found from observations have been taken into account in the discrete models. We have also used individual mass estimates for large bodies of the belts as well as for objects that spacecraft have approached and for bodies with satellites. Our mass estimate for the main asteroid belt is (4.008 ± 0.029) × 10^?4/m_⊕ (3σ). The bulk of the Kuiper belt objects are in the ring zone from 39.4 to 47.8 AU. The estimate of its total mass together with the mass of the 31 largest trans-Neptunian Kuiper belt objects is (1.97 ± 0.30) × 10^?2m_⊕ (3σ), which exceeds the mass of the main asteroid belt almost by a factor of 50. The mass of the 31 largest trans-Neptunian objects (TNOs) is only about 40% of the total one.     

Assessing the Effect of Spacecraft Motion on Single-Spacecraft Solar Wind Tracking Techniques

Recent advances in wide-angle imaging by the Solar Mass Ejection Imager (SMEI) on board the Coriolis spacecraft and more recently by the Heliospheric Imagers (HI) aboard NASA’s Solar TErrestrial RElations Observatory (STEREO), have enabled solar wind transients to be imaged and tracked from the Sun to 1 AU and beyond. In this paper we consider two of the techniques that have been used to determine the propagation characteristics of solar wind transients based on single-spacecraft observations, in particular propagation direction and radial speed. These techniques usually assume that the observing spacecraft remains stationary for the duration of observation of the solar wind transient. We determine the inaccuracy introduced by this assumption for the two STEREO spacecraft and find that it can be significant, and it can lead to an overestimation of the transient velocity as seen from STEREO-A and an underestimation as seen by STEREO-B. This has implications for the prediction or solar wind transients at 1 AU and hence is important for the study of space weather.     

Astrometric observations of planetary satellites at the Abastumani Astrophysical Observatory

We present and discuss the results of the astrometry project during which we observed the satellites of Mars, Jupiter, Saturn, Uranus, and Neptune at the Abastumani Astrophysical Observatory (Georgia) between 1983 and 1994. Observations at the Abastumani Observatory were performed with the double Zeiss astrograph (DZA: D/F = 400/3024 mm) and AZT-11 telescope (F = 16 m). We processed a large array of observations and determined exact coordinates of the planets and their satellites in a system of reference stars of modern catalogues as well as relative coordinates of the satellites. The results were compared with modern ephemerides using the MULTI-SAT software. The comparison enabled us to estimate the accuracy of observations (their random and systematic uncertainties) and the accuracy of modern theories of the motion of planets and their satellites. Random uncertainties of observations are estimated to be 0.10??C0.40?? for various objects and observational conditions. Observational results obtained for Uranus, Neptune and the satellites Titania and Oberon were shown to deviate appreciably and systematically from theories of their motion. The results of observations are presented in the Pulkovo database for Solar System bodies that is available at the website http://www.puldb.ru.     

The equatorial shape and gravity field of Mercury from MESSENGER flybys 1 and 2

On 14 January and 6 October 2008 the MESSENGER spacecraft passed within 200 km of the surface of Mercury. These flybys by MESSENGER provided the first observations of Mercury from a spacecraft since the Mariner 10 flybys in 1974 and 1975. Data from the Mercury Laser Altimeter (MLA) provided new information on the equatorial shape of Mercury, and Doppler tracking of the spacecraft through the flybys provided new data on the planet’s gravity field. The MLA passes were on opposite hemispheres of the planet and span collectively ∼40% of the equatorial circumference. The mean elevation of topography observed during flyby 1, in the longitude range 0-90°E, is greater than that seen during flyby 2 in the longitude range 180-270°E, indicating an offset between centers of mass and figure having a magnitude and phase in general agreement with topography determined by Earth-based radar. Both MLA profiles are characterized by slopes of ∼0.015° downward to the east, which is consistent with a long-wavelength equatorial shape defined by a best-fitting ellipse. The Doppler tracking data show sensitivity to the gravitational structure of Mercury. The equatorial ellipticity of the gravitational field, C_2,2, is well determined and correlates with the equatorial shape. The S_2,2 coefficient is ∼0, as would be expected if Mercury’s coordinate system, defined by its rotational state, is aligned along its principal axes of inertia. The recovered value of the polar flattening of the gravitational potential, J₂, is considerably lower in magnitude than the value obtained from Mariner 10 tracking, a result that is problematic for internal structure models. This parameter is not as well constrained as the equatorial ellipticity because the flyby trajectories were nearly in the planet’s equatorial plane. The residuals from the Doppler tracking data suggest the possibility of mascons on Mercury, but flyby observations are of insufficient resolution for confident recovery. For a range of assumptions on degree of compensation and crustal and mantle densities, the allowable crustal thickness is consistent with the upper limit of about 100 km estimated from the inferred depth of faulting beneath a prominent lobate scarp, an assumed ductile flow law for crustal material, and the condition that temperature at the base of the crust does not exceed the solidus temperature. The MESSENGER value of C_2,2 has allowed an improved estimate of the ratio of the polar moment of inertia of the mantle and crust to the full polar moment (C_m/C), a refinement that strengthens the conclusion that Mercury has at present a fluid outer core.     

Monitoring of energy spectra of particle bursts in the plasma sheet and magnetosheath

Energetic ion (E ? 290 keV) and electron (E_e ? 220 keV) burst intensities were simultaneously monitored at various regions of the plasma sheet and magnetosheath by the CPME JHU/APL instruments on board the IMP-7 and 8 s/c during an extended period from day 250, 1975 to day 250, 1976 when the two spacecraft were closely trailing each other in crossing the geomagnetotail. The energy spectra of the energetic particle populations of different regions in the magnetotail were also computed and monitored simultaneously at the positions of the two spacecraft. The results indicate that the energetic particle intensities are higher and the energy spectra in general considerably softer inside the plasma sheet than the adjacent magnetosheath. The spectral index γ of a power law fit in the computed energy spectrum inside the plasma sheet occasionally exceeds γ > 10 for the ions and γ > 6 for the electrons. Furthermore simultaneous monitoring of particle intensities in the vicinity of the neutral sheet and the high latitude plasma sheet shows higher intensities in the former region. The observations suggest that the energetic particles escape to the magnetosheath from their source inside the plasma sheet by a rigidity dependent process. A dawn-dusk asymmetry in the particle acceleration and escape processes is implied in the observations and discussed in detail.     

Characterization of the particle radiation environment at three potential landing sites on Mars using ESA’s MEREM models

The ‘Mars Energetic Radiation Environment Models’ (dMEREM and eMEREM) recently developed for the European Space Agency are herein used to estimate, for the first time, background Galactic Cosmic Ray (GCR) radiation and flare related solar energetic particle (SEP) events at three candidate martian landing sites under conditions where particle arrival occurred at solar minimum (December, 2006) and solar maximum (April, 2002) during Solar Cycle 23. The three landing sites were selected on the basis that they are characterized by significantly different hydrological conditions and soil compositions. Energetic particle data sets recorded on orbit at Mars at the relevant times were incomplete because of gaps in the measurements due to operational constraints. Thus, in the present study, comprehensive near-Earth particle measurements made aboard the GOES spacecraft were used as proxies to estimate the overall particle doses at each perspective landing site, assuming in each case that the fluxes fell off as 1/r² (where r is the helio-radial distance) and that good magnetic connectivity always prevailed. The results indicate that the particle radiation environment on Mars can vary according to the epoch concerned and the landing site selected. Particle estimations obtained using MEREM are in reasonable agreement, given the inherent differences between the models, with the related NASA Heavy Ion–Nucleon Transport Code for Space Radiation/HZETRN. Both sets of results indicated that, for short (30 days) stays, the atmosphere of Mars, in the cases of the SEPs studied and the then prevailing background galactic cosmic radiation, provided sufficient shielding at the planetary surface to maintain annual skin and blood forming organ/BFO dose levels below currently accepted ionizing radiation exposure limits. The threat of occurrence of a hard spectrum SEP during Cruise-Phase transfers to/from Mars over 400 days, combined with the associated cumulative effect of prolonged GCR exposure, poses an as yet unsolved hazard to prospective onboard personnel.     

Three-Dimensional (3-D) Reconstructions of EISCAT IPS Velocity Data in the Declining Phase of Solar Cycle 23

The European Incoherent SCATter (EISCAT) radar has been used for remote-sensing observations of interplanetary scintillation (IPS) for a quarter of a century. During the April/May 2007 observing campaign, a large number of observations of IPS using EISCAT took place to give a reasonable spatial and temporal coverage of solar wind velocity structure throughout this time during the declining phase of Solar Cycle 23. Many co-rotating and transient features were observed during this period. Using the University of California, San Diego three-dimensional (3-D) time-dependent computer assisted tomography (C.A.T.) solar-wind reconstruction analysis, we show the velocity structure of the inner heliosphere in three dimensions throughout the time interval of 20 April through 20 May 2007. We also compare to white-light remote-sensing observations of an interplanetary coronal mass ejection (ICME) seen by the STEREO Ahead spacecraft inner Heliospheric Imager on 16 May 2007, as well as to in-situ solar-wind measurements taken with near-Earth spacebourne instrumentation throughout this interval. The reconstructions show clear co-rotating regions during this period, and the time-series extraction at spacecraft locations compares well with measurements made by the STEREO, Wind, and ACE spacecraft. This is the first time such clear structures have been revealed using this 3-D technique with EISCAT IPS data as input.