Goddard Robotic Telescope – Optical follow-up of GRBs and coordinated observations of AGNs

Since it is not possible to predict when a Gamma-Ray Burst (GRB) will occur or when Active Galactic Nucleus (AGN) flaring activity starts, follow-up/monitoring ground telescopes must be located as uniformly as possible all over the world in order to collect data simultaneously with Fermi and Swift detections. However, there is a distinct gap in follow-up coverage of telescopes in the eastern U.S. region based on the operations of Swift. Motivated by this fact, we have constructed a 14″ fully automated optical robotic telescope, Goddard Robotic Telescope (GRT), at the Goddard Geophysical and Astronomical Observatory. The aims of our robotic telescope are (1) to follow-up Swift/Fermi GRBs and (2) to perform the coordinated optical observations of Fermi Large Area Telescope (LAT) AGN. Our telescope system consists of off-the-shelf hardware. With the focal reducer, we are able to match the field of view of Swift narrow instruments (20′ × 20′). We started scientific observations in mid-November 2008 and GRT has been fully remotely operated since August 2009. The 3σ upper limit in a 30 s exposure in the R filter is ∼15.4 mag; however, we can reach to ∼18 mag in a 600 s exposures. Due to the weather condition at the telescope site, our observing efficiency is 30–40% on average.     

Using gamma-ray burst prompt emission to probe relativistic shock acceleration

It is widely accepted that the prompt transient signal in the 10 keV–10 GeV band from gamma-ray bursts (GRBs) arises from multiple shocks internal to the ultra-relativistic expansion. The detailed understanding of the dissipation and accompanying acceleration at these shocks is a currently topical subject. This paper explores the relationship between GRB prompt emission spectra and the electron (or ion) acceleration properties at the relativistic shocks that pertain to GRB models. The focus is on the array of possible high-energy power-law indices in accelerated populations, highlighting how spectra above 1 MeV can probe the field obliquity in GRB internal shocks, and the character of hydromagnetic turbulence in their environs. It is emphasized that diffusive shock acceleration theory generates no canonical spectrum at relativistic MHD discontinuities. This diversity is commensurate with the significant range of spectral indices discerned in prompt burst emission. Such system diagnostics are now being enhanced by the broad-band spectral coverage of bursts by the Fermi Gamma-Ray Space Telescope; while the Gamma-Ray Burst Monitor (GBM) provides key diagnostics on the lower energy portions of the particle population, the focus here is on constraints in the non-thermal, power-law regime of the particle distribution that are provided by the Large Area Telescope (LAT).     

GRB observed by IBIS/PICsIT in the MeV energy range

We present the preliminary results of a systematic search for GRB and other transients in the publicly available data for the IBIS/PICsIT (0.2–10 MeV) detector on board INTEGRAL. Lightcurves in 2–8 energy bands with time resolution from 1 to 62.5 ms have been collected and an analysis of spectral and temporal characteristics has been performed. This is the nucleus of a forthcoming first catalog of GRB observed by PICsIT.     

Can a double component outflow explain the X-ray and optical lightcurves of Swift Gamma-Ray Bursts?

An increasing sample of Gamma-Ray Bursts (GRBs) observed by Swift show evidence of ‘chromatic breaks’, i.e. breaks that are present in the X-ray but not in the optical. We find that in a significant fraction of these GRB afterglows the X-ray and the optical emission cannot be produced by the same component. We propose that these afterglow lightcurves are the result of a two-component jet, in which both components undergo energy injection for the whole observation and the X-ray break is due to a jet break in the narrow outflow. Bursts with chromatic breaks also explain another surprising finding, the paucity of late achromatic breaks. We propose a model that may explain the behaviour of GRB emission in both X-ray and optical bands. This model can be a radical and noteworthy alternative to the current interpretation for the ‘canonical’ XRT and UVOT lightcurves, and it bears fundamental implications for GRB physics.     

Testing the blast wave model with Swift GRBs

The complex structure of the light curves of Swift GRBs (e.g. superimposed flares and shallow decay) has made their interpretation and that of the blast wave caused by the burst, more difficult than in the pre-Swift era. We aim to constrain the blast wave parameters: electron energy distribution, p, density profile of the circumburst medium, k, and the continued energy injection index, q. We do so by comparing the observed multi-wavelength light curves and X-ray spectra of a Swift sample to the predictions of the blast wave model.     

Modeling and data analysis of a Forbush decrease

We study the Forbush decrease of the galactic cosmic ray intensity observed in 9–25 September 2005 using the experimental data and a newly developed time-dependent three dimensional modeling. We analyze neutron monitors and muon telescopes, and the interplanetary magnetic field data. We demonstrate a clear relationship between the rigidity (R) spectrum exponent (γ) of the Forbush decrease and the exponent (ν) of the power spectral density of the components of the interplanetary magnetic field in the frequency range of ∼ 10⁻⁶–10 ⁻⁵ Hz. We confirm that an inclusion of the time-dependent changes of the exponent ν makes the newly developed nonstationary three dimensional model of the Forbush decrease compatible with the experimental data. Also, we show that the changes of the rigidity spectrum exponent γ does not depend on the level of convection of the galactic cosmic rays stream by solar wind; depending on the changes of the exponent ν, i.e. on the state of the turbulence of the interplanetary magnetic field.     

X-ray study of accretion flow in narrow-line Seyfert 1 galaxies

We present the results of a systematic study of narrow-line Seyfert 1 galaxies (NLS1s) observed with XMM-Newton. The 2–12 keV X-ray spectra of NLS1s are well represented by a single power law with a photon index Γ ∼ 2. When this hard power law continuum is extrapolated into the low energy band, we found that all objects in our sample show prominent soft excess emission. This excess emission is well parameterized by the thermal emission expected from an optically thick accretion disk, and we found the following three peculiar features: (1) The derived disk temperatures are significantly higher than the expectation from a standard Shakura-Sunyaev accretion disk, if we assume a central mass of a black hole to be 10^6–8M_⊙. (2) The temperatures are distributed within narrow range (ΔkT ∼ 0.08 keV) with an average temperature of 0.18 keV in spite of the range of four orders of magnitude in luminosity (10^41–45 erg s⁻¹). (3) We found a peculiar temperature–luminosity relation, where the luminosity seems to be almost saturated in spite of the significant change in temperature, during the observations of the most luminous NLS1 PKS 0558-504. These results strongly suggest that the standard accretion disk picture is no longer appropriate in the nuclei of NLS1s. We discuss a possible origin for the soft excess component, and suggest that a slim disk may be able to explain the observational results, if the photon trapping effect is properly taken into account.     

Blazar nuclei in radio-loud narrow-line Seyfert 1?

It has been suggested that some radio-loud narrow-line Seyfert 1 contain relativistic jets, on the basis of their flat-spectrum radio nuclei and studies on variability. We present preliminary results of an ongoing investigation of the X-ray and multiwavelength properties of 5 radio-loud NLS1 based on archival data from Swift and XMM-Newton. Some sources present interesting characteristics, very uncharacteristic for a radio-quiet narrow-line Seyfert 1, such as very hard X-ray spectra, and correlated optical and ultraviolet variability. However, none of the studied sources show conclusive evidence for relativistic jets. γ-Ray observations with Fermi are strongly recommended to definitely decide on the presence or not of relativistic jets.     

Solar cycle and latitude dependence of high-beta suprathermal plasma conditions in interplanetary space between 1.3 and 5.4 AU

The analysis of energetic particles and magnetic field measurements from the Ulysses spacecraft has shown that in a series of events, the energy density contained in the suprathermal tail particle distribution is comparable to or larger than that of the magnetic field, creating conditions of high-beta plasma. In this work we analyze periods of high-beta suprathermal plasma occurrences (β_ep > 1) in interplanetary space, using the ratio (β_ep) of the energetic particle (20 keV to ∼5 MeV) and magnetic field energy densities from measurements covering the entire Ulysses mission lifetime (1990–2009) in order to reveal new or to reconfirm some recently defined interesting characteristics. The main key-results of the work are summarized as follows: (i) we verify that high-beta events are detected within well identified regions corresponding mainly to the vicinity of shock surfaces and magnetic structures, and associated with energetic particle intensity enhancements due to (a) reacceleration at shock-fronts and (b) unusually large magnetic field depressions. (ii) We define three considerable features for the high-beta events, concentrated on the next points: (a) there is an appreciable solar-activity influence on the high-beta events, during the maximum and middle solar-cycle phase, (b) the annual peak magnitude and the number of occurrences of high events are well correlated with the sunspot number, (c) the high-beta suprathermal plasma events present a spatial distribution in heliographic latitudes (HL) up to ∼±80°, and a specific important concentration on the low (−25° ? HL < −6°, 6° < HL ? 25°) and median (−45° ? HL < −25°, 25° < HL ? 45°) latitudes. We also reconfirm by a statistical analysis the results of Marhavilas and Sarris (2011), that the high-beta suprathermal plasma (β_ep > 1) events are characterized by a very large parameter β_ep (up to 1732.5), a great total duration (406 days) and a large percentage of the Ulysses-mission lifetime (which is equal to 6.34% of the total duration with usable measurements, and 11.3% of the duration in presence of suprathermal particles events).     

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.     

X-ray redshifts with the International X-ray Observatory (IXO)

We explore the capabilities of the future space science mission IXO (International X-ray Observatory) for obtaining cosmological redshifts of distant Active Galactic Nuclei (AGNs) using the X-ray data only. We first find in which regions of the X-ray luminosity (L_X) versus redshift (z) plane the weak but ubiquitous Fe Kα narrow emission line can deliver an accurate redshift (δz < 5%) as a function of exposure time, using a CCD-based Wide Field Imager (IXO/WFI) as the one baselined for IXO. Down to a 2–10 keV X-ray flux of 10⁻¹⁴ erg cm⁻² s⁻¹ IXO/WFI exposures of 100 ks, 300 ks and 1 Ms will deliver 20%, 40% and 60% of the redshifts. This means that in a typical 18′ × 18′ IXO/WFI field of view, 4, 10 and 25 redshifts will be obtained for free from the X-ray data alone, spanning a wide range up to z ∼ 2–3 and fairly sampling the real distribution. Measuring redshifts of fainter sources will indeed need spectroscopy at other wavebands.     

Spectral analysis of time-integrated Konus–Wind GRBs: Implication on radiative mechanisms

We have analysed a sample of 328 time-integrated GRB prompt emission spectra taken via the Konus instrument on board the US GGS-Wind spacecraft between 2002 and 2004 using a couple of two-components models, Cut-off Power Law (CPL) + Power Law (PL) and blackbody (BB) + PL. The spectra show clear deviation from the Band function. The PL term is interpreted as the low energy tail of a nonthermal emission mechanism. The distributions of corresponding index β give values β < −2/3 consistent with synchrotron and synchrotron self-Compton mechanisms. The distribution of low energy index α associated with the CPL term shows clear discordance with synchrotron models for 31.4% of the analysed GRBs with values exceeding that for the line of death, α = −2/3. Then, a set of nonthermal radiation mechanisms producing harder slopes, i.e., α > −2/3, are presented and discussed. For the remaining majority (68.6%) of GRBs with CPL index α < −2/3, we show that optically thin synchrotron produced by a power law electron distribution of type, N(γ) ∼ γ^−p, γ₁ < γ < γ₂, for finite energy range (γ₂ ≠ ∞) is a likely emission mechanism with α ∼−(p + 1)/2 in the frequency range ν₁ ? ν ? ν₂ (where ν₂ = η²ν₁ with η = γ₂/γ₁), such that for p > 1/3, one gets α < −2/3. We also show that corresponding spectra in terms of F_ν and νF_ν functions are peaked around frequency ν₂ instead of ν₁, respectively for p < 1 and p < 3. Besides, thermal emission is examined taking a single Planck function for fitting the low energy range. It can be interpreted as an early emission from the GRB fireball photosphere with observed mean temperature, kT′ ∼ 16.8 keV. Furthermore, we have performed a statistical comparison between the CPL + PL and BB + PL models finding comparable χ²-values for an important fraction of GRBs, which makes it difficult to distinguish which model and specific radiation mechanism (possible thermal or nonthermal γ-ray emissions) are best suitable for describing the reported data. Therefore, additional information for those bursts, such as γ-ray polarization, would be highly desirable in future determinations of GRBs observational data.     

Probing the method of correcting Faraday rotation in vector magnetograms

By analyzing the vector magnetograms of Huairou Solar Observing Station (HSOS) taken at the line center (0.0 Å) and the line wing (−0.12 Å) of FeI λ5324.19 Å, we make an estimate of the measured errors in transversal azimuths (δ?) caused by Faraday rotation. Since many factors, such as the magnetic saturation and scattered light, can affect the measurement accuracy of the longitudinal magnetic field in the umbrae of sunspots, we limit our study in the region ∣B_z∣ < 800 G. The main mean azimuth rotations are about 4°, 6°, 7° and 9°, while ∣B_z∣ are in the ranges of 400–500 G, 500–600 G, 600–700 G and 700–800 G, respectively. Moreover, we find there is also an azimuth rotation of about 8° at the wavelength offset −0.12 Å of the line compared against a previous numerical simulation.     

Quantifying the relationship between the measurement precision and specifications of a UV/visible sensor on a geostationary satellite

To investigate the feasibility of new satellite observations, including air quality (AQ) observations from geostationary (GEO) orbit, it is essential to link the measurement precision (ε) with sensor specifications in advance. The present study attempts to formulate the linkage between ε and specifications of a UV/visible sensor (signal-to-noise ratio (SNR), full width at half maximum (FWHM) of the slit function, and sampling ratio (SR)) on a GEO satellite. A sophisticated radiative transfer model (JACOSPAR) is used to calculate synthetic radiance spectra that would be measured by a UV/visible sensor observing the atmosphere over Tokyo (35.7°N, 139.7°E) from GEO orbit at 120°E longitude. The spectra, modified according to given sensor specifications, are analyzed by the differential optical absorption spectroscopy technique to estimate the ε for slant column densities of O₃ and NO₂. We find clear relationships: for example, the ε of the O₃ slant column density (molecules cm⁻²) and SNR at 330 nm are linked by the equation log(ε) = −1.06 · log(SNR) + 20.71 in the UV region, and the ε of the NO₂ slant column density and SNR at 450 nm are linked by log(ε) = −0.98 · log(SNR) + 18.00, at a FWHM = 0.6 nm (for the Gaussian slit function) and SR = 4. The relationships are mostly independent of other specifications (e.g., horizontal and temporal resolutions), as they affect ε primarily through SNR, providing constraints in determining the optimal SNR (and alternatively FWHM and SR) for similar UV/visible sensors dedicated for AQ studies.     

Cowpeas and pinto beans: Performance and yields of candidate space crops in the laboratory biosphere closed ecological system

An experiment utilizing cowpeas (Vigna unguiculata L.), pinto beans (Phaseolus vulgaris L.) and Apogee ultra-dwarf wheat (Triticum sativa L.) was conducted in the soil-based closed ecological facility, Laboratory Biosphere, from February to May 2005. The lighting regime was 13 h light/11 h dark at a light intensity of 960 μmol m⁻² s⁻¹, 45 mol m⁻² day⁻¹ supplied by high-pressure sodium lamps. The pinto beans and cowpeas were grown at two different planting densities. Pinto bean production was 341.5 g dry seed m⁻² (5.42 g m⁻² day⁻¹) and 579.5 dry seed m⁻² (9.20 g m⁻² day⁻¹) at planted densities of 32.5 plants m⁻² and 37.5 plants m⁻², respectively. Cowpea yielded 187.9 g dry seed m⁻² (2.21 g m⁻² day⁻¹) and 348.8 dry seed m⁻² (4.10 g m⁻² day⁻¹) at planted densities of 20.8 plants m⁻² and 27.7 plants m⁻², respectively. The crop was grown at elevated atmospheric carbon dioxide levels, with levels ranging from 300–3000 ppm daily during the majority of the crop cycle. During early stages (first 10 days) of the crop, CO₂ was allowed to rise to 7860 ppm while soil respiration dominated, and then was brought down by plant photosynthesis. CO₂ was injected 27 times during days 29–71 to replenish CO₂ used by the crop during photosynthesis. Temperature regime was 24–28 °C day/deg 20–24 °C night. Pinto bean matured and was harvested 20 days earlier than is typical for this variety, while the cowpea, which had trouble establishing, took 25 days more for harvest than typical for this variety. Productivity and atmospheric dynamic results of these studies contribute toward the design of an envisioned ground-based test bed prototype Mars base.     

Study of the 27-day variations of the galactic cosmic ray intensity and anisotropy

We demonstrate that the general features of the radial and azimuthal components of the anisotropy of galactic cosmic rays can be studied by the harmonic analysis method using data from an individual neutron monitor with cut off rigidity <5 GV. In particular, we study the characteristics of the 27-day (solar rotation period) variations of the galactic cosmic ray intensity and anisotropy, solar wind velocity, interplanetary magnetic field strength and sunspot number. The amplitudes of the 27-day variations of the galactic cosmic ray anisotropy are greater, and the phases more clearly established, in A > 0 polarity periods than in A < 0 polarity periods at times of minimum solar activity. The phases of the 27-day variations of the galactic cosmic rays intensity and anisotropy are opposite with respect to the similar changes of the solar wind velocity in A > 0 polarity periods. No significant dependence of the amplitude of the 27-day variation of the galactic cosmic ray anisotropy on the tilt angle of the heliospheric neutral sheet is found. Daily epicyclegrams obtained by Chree’s method show that the 27-day variations of the galactic cosmic ray anisotropy during A > 0 polarity periods follow elliptical paths with the major axes oriented approximately along the interplanetary magnetic field. The paths are more irregular during A < 0 polarity periods.     

Rigidity spectrum of the 27-day variation of the galactic cosmic ray intensity in different epochs of solar activity

We study the temporal evolution of the power rigidity spectrum of the first (27 days) and the second (14 days) harmonics of the 27-day variation of the galactic cosmic ray intensity measured by neutron monitors in the period of 1965–2002. The rigidity spectrum of these variations can be approximated by a power law. We show the rigidity spectra of the first and the second harmonics of the 27-day variation of the galactic cosmic ray intensity have similar time profiles. These spectra are hard (γ ≈ 0.5 ± 0.1) and soft (γ ≈ 1.1 ± 0.2) during solar maximum and minimum activity, respectively. We ascribe this to the alternation of the sizes of the modulation regions responsible for the 27-day variation of the galactic cosmic ray intensity in different epochs of solar activity. Especially, the average radial sizes of the modulation regions of the 27-day variation and the heliolatitudinal extension of the heliolongitudinal asymmetry are smaller during solar minimum than during solar maximum. We show also, that the temporal changes of the power rigidity spectra of the first and the second harmonics of the 27-day variation of the galactic cosmic ray intensity are in a negative correlation with the changes of the rigidity spectrum of the corresponding 11-year variation.     

Day-to-day variability of equatorial anomaly in GPS-TEC during low solar activity period

The ionospheric total electron content (TEC) in the northern hemispheric equatorial ionization anomaly (EIA) crest region is investigated by using dual-frequency signals of the Global Positioning System (GPS) acquired from Rajkot (Geog. Lat. 22.29°N, Geog. Long. 70.74°E; Geom. Lat. 14.21°N, Geom. Long. 144.90°E), India. The day-to-day variability of EIA characteristics is examined during low solar activity period (F10.7∼83 sfu). It is found that the daily maximum TEC at EIA crest exhibits a day-to-day and strong semi-annual variability. The seasonal anomaly and equinoctial asymmetry in TEC at EIA is found non-existent and weaker, respectively. We found a moderate and positive correlation of daily magnitude of crest, Ic with daily F10.7 and EUV fluxes with a correlation coefficient of 0.43 and 0.33, respectively indicating an existence of a short-term relation between TEC at EIA and the solar radiation even during low solar activity period. The correlation of daily Ic with Dst index is also moderate (r = −0.35), whereas no correlation is found with the daily Kp index (r = 0.14) respectively. We found that the magnitude of EIA crest is moderately correlated with solar flux in all seasons except winter where it is weakly related (0.27). The magnitude of EIA crest is also found highly related with EEJ strength in spring (r = 0.69) and summer (r = 0.65) than autumn (0.5) and winter (r = 0.47), though EEJ is stronger in autumn than spring.     

GRB host galaxies: An unbiased sample

We describe the current status and recent results from our Swift/VLT legacy survey, a VLT Large Programme aimed at characterizing the host galaxies of a homogeneously selected sub-sample of Swift   GRBs. The immediate goals are to determine the host luminosity function, study the effects of reddening, determine the fraction of Lyα$\mathrm{Ly}\mathrm{\alpha }$ emitters in the hosts, and obtain redshifts for targets without a reported one. The main effort so far has been the definition of a very carefully selected sample, obeying strict and well-defined criteria: 68 targets in total. Among the preliminary results is a large optical detection rate, the lack of extremely red objects (only one possible case in the sample) and an update of the Swift   GRB redshift distribution with 〈z〉∼2.0$〈z〉\sim 2.0$.