{'fraction_non_alphanumeric': 0.04295529894133286, 'fraction_numerical': 0.052196614187100594, 'mean_word_length': 4.125010260198637, 'pattern_counts': {'":': 0, '<': 7, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 7, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Optical and near-infrared follow-up observations of four Fermi/LAT GRBs : Redshifts, afterglows, energetics and host galaxies

Fermi can measure the spectral properties of gamma-ray bursts over a very large energy range and is opening a new window on the prompt emission of these energetic events. Localizations by the instruments on Fermi in combination with follow-up by Swift provide accurate positions for observations at longer wavelengths leading to the determination of redshifts, the true energy budget, host galaxy properties and facilitate comparison with pre-Fermi bursts. Multi-wavelength follow-up observations were performed on the afterglows of four bursts with high energy emission detected by Fermi/LAT : GRB090323, GRB090328, GRB090510 and GRB090902B. They were obtained in the optical/near-infrared bands with GROND mounted at the MPG/ESO 2.2m telescope and additionally of GRB090323 in the optical with the 2 m telescope in Tautenburg, Germany. Three of the events are classified as long bursts while GRB090510 is a well localized short GRB with GeV emission. In addition, host galaxies were detected for three of the four bursts. Spectroscopic follow-up was initiated with the VLT for GRB090328 and GRB090510. The afterglow observations in 7 bands are presented for all bursts and their host galaxies are investigated. Knowledge of the distance and the local dust extinction enables comparison of the afterglows of LAT-detected GRBs with the general sample. The spectroscopic redshifts of GRB090328 and GRB090510 were determined to be z=0.7354+/-0.0003 and z=0.903 +/- 0.001 and dust corrected star-formation rates of 4.8 Mdot yr^-1 and 0.60 M_dot yr^-1 were derived for their host galaxies, respectively. The afterglows of long bursts exhibit power-law decay indices alpha from less than 1 to ~2.3 and spectral indices (beta) values from 0.65 to ~1.2 which are fairly standard for GRB afterglows. Constraints are placed on the jet half opening angles of less than 2.1 deg to greater than 6.4 deg which allows limits to be placed on the beaming corrected energies. These range from less than 5x10^50 erg to the one of the highest values ever recorded, greater than 2.2x10^52 erg for GRB090902B, and are not consistent with a standard candle. The extremely energetic long Fermi bursts have optical afterglows which lie in the top half of the brightness distribution of all optical afterglows detected in the Swift era or even in the top 5% if incompleteness is considered. The properties of the host galaxies of these LAT detected bursts in terms of extinction, star formation rates and masses do not appear to differ from previous samples.


Introduction
The follow-up of gamma-ray bursts (GRBs) detected by the Swift satellite (Gehrels et al. 2004) has led to the determination of the distance scale for a large sample of bursts. The Burst Alert Telescope (BAT, Barthelmy et al. 2005) is sensitive in the energy range 15−150 keV and has good localization capabilities with typical uncertainties in the arcminute range. Rapid follow-up by Swift's narrow field instruments in the X-rays (XRT, Burrows et al. 2005) and optical/UV (UVOT, Roming et al. 2005) have lead to the arcsecond localizations required for ground-based observers and in turn to spectro- scopic redshift measurements of a large sample of GRBs (e.g., Fynbo et al. 2009) and investigation of their host galaxies (e.g., Savaglio et al. 2009;Perley et al. 2009b). To date distances to ∼ 200 GRB sources have been established with redshifts ranging from z = 0.0085 (GRB 980425: Tinney et al. 1998;Galama et al. 1998) to z ∼ 8.2 (GRB090423: Tanvir et al. 2009;Salvaterra et al. 2009). The BAT has a narrow spectral range and is not able to determine the spectral parameters of the prompt emission of GRBs over a broad energy range. Since the launch of the Fermi Gamma-Ray Space Telescope (Atwood et al. 1994;Michelson 1996) there is now the possibility to investigate the spectrum over seven decades in energy. These events can be localized by instruments from Fermi and followed up by Swift and ground-based obervatories, enabling investigation of their afterglows and host properties and facilitating comparison to a general sample. In some fortuitous cases instruments on the Fermi and Swift satellites may trigger on the same event with high-energy emission enabling both rapid follow-up and broadband prompt emission coverage (e.g. GRB 090510, see Section 2.3). The high-energy spectral properties of GRBs were investigated previously by instruments on-board the Compton Gammaray Observatory, however the redshifts of these events are unknown. The spectra of these GRBs were described by an extrapolation of the low energy spectra to energies >100 MeV (e.g., Hanlon et al. 1994), and in some cases an additional component at high-energies (González et al. 2003;Kaneko et al. 2008) and long-lived GeV emission (Hurley et al. 1995). Recently a photometric redshift of z = 1.8 +0.4 −0.3 was reported by Rossi et al. (2008) for GRB 080514B which was detected at energies up to 300 MeV by instruments on the AGILE satellite (Giuliani et al. 2008).
Fermi hosts two instruments that detect GRBs: the Gammaray Burst Monitor (GBM: Meegan et al. 2009) sensitive to photons from ∼8 keV to ∼40 MeV, and the Large Area Telescope with a spectral coverage from ∼ 30 MeV to ∼100 GeV (LAT: Atwood et al. 2009). Together they provide valuable information on the spectral properties of the prompt emission (e.g. Abdo et al. 2009a,b,c,d), the presence or absence of intrinsic spectral cut-offs or those due to the optical depth of the universe to high-energy γ-rays due to pair production on infrared diffuse Extragalactic Background Light (e.g., de Jager & Stecker 2002;Metcalfe et al. 2003;Kneiske et al. 2004;Kashlinsky et al. 2005;Stecker et al. 2006;Franceschini et al. 2008;Finke et al. 2009;Abdo et al. 2009c) and test for quantum gravity effects (e.g., Amelino-Camelia et al. 1998;Scargle et al. 2008;Mattingly 2005;Amelino-Camelia & Smolin 2009;Abdo et al. 2009d,a).
Among the most impressive of these events are GRB 080916C (Abdo et al. 2009d) and GRB 090510 (Abdo et al. 2009a). GRB 080916C is a long bright GRB for which the prompt emission spectrum could be fit over six decades in energy by the empirical Band function (see Abdo et al. 2009d, for detailed results). The burst was found to have a high photometric redshift of z = 4.35 ± 0.15 via ground-based follow up observations ). The distance information enabled the determination of the rest frame properties, energetics, and the placing of lower limits on the bulk Lorentz factor of the outflow ( 900 − 1100) (Abdo et al. 2009d;Greiner et al. 2009a). GRB 090510 is a short burst with a T 90 of 2.1 s and from which a 31 GeV photon was detected (Abdo et al. 2009a). The distance determination (Rau et al. 2009b) and high energy emission allow a limit to be placed on photon dispersion (Abdo et al. 2009a). Moreover, very recently a 33 GeV photon was detected from GRB 090902B (Abdo et al. 2009c) and this burst is at redshift z = 1.822 (Cucchiara et al. 2009). Clearly, the distance determination via optical spectroscopy is crucial to the interpretation of these events.
We report on the optical and near-infrared (NIR) observations of the afterglow of four bursts with high energy emission, GRB 090323, GRB 090328, GRB 090510 and GRB 090902B and compare them to the sample of GRB afterglows to date. Furthermore we also report on the spectroscopic redshift determination for GRB 090328 and GRB 090510. Throughout the paper, we adopt concordance ΛCDM cosmology (Ω M = 0.27, Ω Λ = 0.73, H 0 = 71 km/s/Mpc), and the convention that the flux density of the GRB afterglow can be described as F ν (t) ∝ ν −β t −α

GRB 090323
At 00:02:42.63 UT on 23 March 2009, Fermi GBM triggered and located the long burst GRB 090323 ). An excess of 5 σ was detected by the LAT and the burst was localized to R.A.(J2000), decl.(J2000)=190.69, +17.08 degrees with a 68% containment radius of 0.09 • . Furthermore, it was reported that emission continued up to a few kiloseconds post trigger . A triangulation by the Inter Planetary Network localized the burst to an arc which intersected the GBM and LAT positions (Hurley et al. 2009). Swift carried out a Target of Opportunity observation ∼ 19 hours post burst and a fading afterglow was found (Kennea et al. 2009a;). Updike et al. (2009b) observed the position of the Xray afterglow with the Gamma-ray Burst Optical/Near-infrared Detector (GROND) in 7 bands ∼ 27 hours after the burst and detected a source for which they reported a preliminary photometric redshift of 4.0 ± 0.3 (see below). Chornock et al. (2009) reported a spectroscopic redshift of z = 3.57 based on observations of the optical afterglow using the Gemini Multi-Object Spectrograph (GMOS) mounted on the Gemini South Telescope. The GRB was also detected in the radio band (Harrison et al. 2009;van der Horst 2009).
Using the reported spectral parameters (Bissaldi et al. 2010) and a redshift of z = 3.57 (Chornock et al. 2009) E γ,iso is 4.1 × 10 54 erg in the 1 keV to 10 MeV range, or 5.1 × 10 54 erg in the 1 keV to 10 GeV range with a restframe peak energy of 2.7 ± 0.2 MeV. The E γ,iso is one of the highest ever measured even rivals GRB 080916C for the most energetic GRB.
The field of GRB 090323 was observed simultaneously in 7 bands (g ′ r ′ i ′ z ′ JHK S ) by GROND, mounted at the 2.2 m MPG/ESO telescope at La Silla Observatory, Chile ) and using the Schmidt-camera at the Alfred-Jensch telescope of the Thüringer Landessternwarte Tautenburg (TLS), Germany. GROND and TLS observations started on 24th March 2009, 27 hours and 45 hours after the GBM trigger, respectively (Updike et al. 2009b;Kann et al. 2009c). Observations continued in the upcoming nights (Kann et al. 2009d,b) and the last GROND epoch was obtained ∼ 100 days after T 0 . The best optical position of the afterglow measured in the r ′ band against the SDSS catalogue is R.A.(J2000)=12:42:50.28, decl.(J2000)=+17:03:11.9 with systematic and statistical uncertainties of 0. ′′ 1 in each coordinate. A finding chart of the afterglow is shown in Fig. 1.  Fig. 1. GROND finding chart for the afterglow of GRB 090323. The main panel shows the field of GRB 090323 taken 3.5 days post burst. The field shown has 4 ′ × 4 ′ , where North is up and East is to the left. The afterglow position is indicated with two lines in each image. The inset in the top left shows a zoom into the afterglow region at the earliest epoch, where the afterglow was still bright. All images are i ′ band, except the host image (r ′ , top right inset), which is also enhanced in contrast for demonstration purposes.

GRB 090328
At 09:36:46 UT on 28 March 2009, Fermi GBM triggered and located a long burst GRB 090328. The event was significantly detected up to a few GeV at the 5 σ level by the Fermi/LAT and was localized to R.A.(J2000), decl.(J2000)=90.87, −41.95 with a 68% containment radius of 0.11 deg (McEnery et al. 2009b). The emission in the LAT was reported to last until about 900 seconds post trigger by Cutini et al. (2009). A Swift Target of Opportunity observation was initiated and XRT observations commenced ∼16 hours post burst (Kennea 2009) and X-ray and bright UV/optical afterglow candidates were subsequently reported by Kennea et al. (2009b) and Oates (2009). A spectrum of the optical afterglow was taken with the GMOS instrument mounted on the Gemini South Telescope (Cenko et al. 2009a) and the redshift was reported to be z = 0.736. Further photometric observations of the afterglow were reported in the optical (Allen et al. 2009) and optical/NIR with GROND (Updike et al. 2009a, see below) and additionally in the radio band ).
The time averaged spectrum (T 0 + 3.1 to T 0 + 29.7 s) of the prompt emission was reported by GBM to be best fit by a Band function with indices α = −0.93 ± 0.02 and β = −2.2 ± 0.1, a peak energy of E peak = 653 ± 45 keV with an event fluence of (9.5 ± 1.0) ×10 −5 erg/cm 2 in the 8 keV to 40 MeV band (Rau et al. 2009a). Bissaldi et al. (2010) report that the spectrum (T 0 + 0 to T 0 + 66.6 s) is best fit by a power-law with exponential cut off. The index of the power-law is −1.07 ± 0.02 and the peak energy is 744 +50 −47 keV. At a redshift of 0.736 , these spectral parameters correspond to E γ,iso in the 1 keV to 10 MeV range of The main panel shows the field of GRB 090328 taken 1.5 days post burst. The field shown has 4 ′ × 4 ′ , where North is up and East is to the left. The inset in the top left shows a zoom into the afterglow region at the same epoch, where the afterglow is indicated with two lines. The inset on the top right contains a 15 ′′ ×15 ′′ region and shows a host image taken 10 days post burst. The afterglow position is shown by a small white circle of radius 0. ′′ 3. All images are in the i ′ band.
Spectroscopy of the afterglow was performed on 2009 March 30.01 UT (∼ 1.6 days post burst) using the FOcal Reducer and low dispersion Spectrograph 1 (FORS1; Appenzeller et al. 1998) at the 8 m ESO-VLT UT2 telescope (Programme ID: 083.D-0903). Observations were obtained using the 600B and 600R grisms, providing a spectral resolution of roughly 5 Å full width at half maximum. Two 1800 s integrations were taken in each of the two grisms.

GRB 090510
At 00:23:00 UT on 10 May 2009 the Swift BAT triggered and located GRB 090510 ). The burst also triggered Fermi/GBM (Guiriec et al. 2009) and other instruments (Golenetskii et al. 2009;Ohmori et al. 2009;Longo et al. 2009;Giuliani et al. 2010). The duration (T 90 ∼ 2.1 sec) and the negligible spectral lags below ∼1 MeV of the event make it consistent with the class of short bursts (Abdo et al. 2009a). Significant emission was detected by the LAT (Ohno & Pelassa 2009;Omodei et al. 2009). The time integrated spectrum using Fermi GBM and LAT data is best fit by two spectral components, a Band function with the Epeak = 3.9 ± 0.3 MeV and a power-law at high energies. The total isotropic energy is (1.08 ± 0.06) × 10 53 erg (10 keV -30 GeV) (Abdo et al. 2009a). The event was also detected at 5σ significance at energies above 100 MeV by AGILE GRID Giuliani et al. 2010).
Swift slewed after 91 seconds and observations with the XRT and UVOT revealed an uncatalogued source ). The XRT light curve was reported to be fading (Grupe & Hoversten 2009) and the UVOT afterglow candidate was confirmed by Olofsson et al. (2009) using the Nordic Optical Telescope.  report that the afterglow of GRB 090510 is detected in almost all UVOT filters which implies that the redshift is less than about 1.5 ). See also De Pasquale et al. (2010).
GROND imaged the field of GRB 090510 starting 6.18 h after the Swift and Fermi triggers. At that point the position of the GRB was becoming visible above the pointing constraints of the telescope (Olivares et al. 2009b). Therefore GROND images taken in the first two hours suffer from high airmass and moderate image quality. Observations continued for 4 h until local twilight and were resumed in the following night for 1.5 h. In the first epoch, the afterglow is detected in the g ′ r ′ i ′ z ′ optical channels  Fig. 3.
Spectroscopy of the host galaxy of GRB 090510 was performed (∼ 2.3 days post burst) using FORS2 (Appenzeller et al. 1998) at the 8 m ESO-VLT UT1 telescope (Programme ID: 083.D-0903(A)). Observations were obtained using the 300I grism, providing a spectral resolution of ≈ 7 Å FWHM. Three 1800 s integrations were taken. The preliminary analysis and redshift z = 0.903 was reported by Rau et al. (2009b). We note that the redshift of this galaxy is consistent with the afterglow colours detected by UVOT ). Additionally, spectroscopy of the source located 5 ′′ south of the GRB host galaxy was performed ∼ 1.3 days post burst using FORS2 (300V grism), which yields a redshift consistent with that of the GRB host. Two 1800 s integrations were taken of this object.

GRB 090902B
At 11:05:08.31 UT on 2 September 2009 Fermi GBM triggered on a long, bright, hard burst GRB 090902B (trigger 273582310 / 090902462) (Bissaldi & Connaughton 2009). The burst was localized by Fermi/LAT to R.A.(J2000), decl.(J2000) = 265.00, 27.33 with a statistical uncertainty of only 0.06 degrees (de Palma et al. 2009b) and Target of Opportunity observations were initiated with the near-field instruments on Swift ∼12.5 hours post-trigger. A candidate X-ray afterglow within the LAT error circle was reported by Kennea & Stratta (2009) and subsequently confirmed to be fading ). Optical detections were reported by a number of observers (Swenson & Siegel 2009;Swenson & Stratta 2009;Perley et al. 2009a;Guidorzi et al. 2009;Pandey et al. 2009), in the NIR (Olivares et al. 2009a) and in the radio ). The afterglow redshift of z = 1.822 was obtained from an afterglow absorption spectrum by Cucchiara et al. (2009) using the GMOS spectrograph. Based on the spectral parameters (de Palma et al. 2009a), the isotropic energy released is E γ,iso = 2.26×10 54 in the 1 keV to 10 MeV range and 3.54 × 10 54 erg in the 1 keV to 10 GeV band. An excess of emission in addition to the Band function at both low energies 50 keV and above 100 MeV was reported by Abdo et al. (2009c).
GROND imaged the field of GRB 090902B starting ∼ 13 h after the Fermi trigger when the position of the GRB was becoming visible above the pointing constraints of the telescope. At the time of the first epoch observations, only the LAT localization with 3.5 arcmin uncertainty was available. The afterglow location was outside the narrower (5.4 ′ × 5.4 ′ ) optical field of view (FoV) but within the larger NIR FoV (10 ′ × 10 ′ ) (Olivares et al. 2009a). Observations continued in the following nights by GROND and an additional VLT/FORS2 R band observation was executed at ∼ 23 days post trigger (ESO DDT Proposal number : 283.D-5059). A finding chart of the field of GRB 090902B is shown in Fig. 4.

Analysis
GROND optical/NIR data were reduced and analyzed in standard manner using pyraf/IRAF tasks (Tody 1993). A detailed description of the reduction and relative photometry can be found in Krühler et al. (2008). In the afterglow-dominated regime we used PSF fitting techniques to derive the brightness of the optical transient. These were complemented by standard aperture photometry, which yields consistent results. The magnitudes of the GRB host galaxies were measured using aperture photometry with an appropriate aperture correction. In addition, SExtractor (Bertin & Arnouts 1996) was used to perform host photometry, again with consistent results to the reported aperture photometry.
Absolute photometry in g ′ r ′ i ′ z ′ -bands was obtained against the g ′ r ′ i ′ z ′ magnitudes of the Sloan Digital Sky Survey (SDSS DR7, Abazajian et al. 2009) for all bursts. SDSS field stars in the vicinity of GRB 090323 (Updike et al. 2009b) were used for GRBs 090323 and 090328. Both bursts were observed with GROND consecutively in 5 different nights, out of which 4 where photometric. The obtained calibration is consistent between all epochs. The field of GRB 090510 was calibrated against the primary Sloan Standards SA 114-750 and 114-650 and their SDSS field stars. Absolute photometry of the afterglow of GRB 090902B was obtained against a nearby SDSS field observed under photometric conditions. JHK S absolute photometry was measured against the magnitudes of selected stars in the field of GRBs 090323, 090328, 090510 and 090902B from the 2MASS catalogue (Skrutskie et al. 2006). This procedure results in a typical absolute accuracy of 0.04 mag in g ′ r ′ i ′ z ′ , 0.06 mag in JH and 0.08 in K S , which was added quadratically to the statistical error in the analysis of the broadband spectral energy distribution. TLS data of GRB 090323 was reduced in a standard manner with MIDAS and IRAF and calibrated against six SDSS stars in the field-of-view, where we derived R C magnitudes using the transformations of Lupton from 2005 1 . The standard stars were measured with SExtractor, and afterglow magnitudes were derived using aperture photometry. All magnitudes are corrected for the expected Galactic reddening in the direction of the bursts (Schlegel et al. 1998), which is E B−V = 0.03 mag for GRB 090323, E B−V = 0.06 mag for GRB 090328, E B−V = 0.02 mag for GRB 090510, and E B−V = 0.04 mag for GRB 090902B, respectively. All magnitudes are quoted in the AB system throughout this paper.
The FORS1+2 spectroscopy data were reduced with standard IRAF routines, and spectra were extracted using an optimal (variance-weighted) method. The observations were obtained with a long slit of 1. ′′ 0 width and corrected for slit losses due to finite slit width. Spectro-photometric calibrations were carried out using observations of the standard stars EG 274 and LTT 7379 (Hamuy et al. 1992) for the observations of GRB 090328 and GRB 090510, respectively.

Light curve
The GROND and TLS multicolor light curve of the optical/NIR afterglow starting at T 0 + 95 ks is presented in Fig. 5 and is well described with a single power-law plus a superimposed bump component and a constant host contribution dominating at later times (χ 2 = 161/138 d.o.f). All available optical/NIR data were fitted simultaneously, allowing for a color change in the bump and host dominated regime. The initial power-law decay slope is α = 1.90 ± 0.01, with no sign of an additional break at later times. The host magnitudes in the g ′ z ′ JHK S bands are not constrained by the data and have been fixed. Similarly, the flare is unconstrained in the H and K S bands. The host galaxy of GRB 090323 was detected in the r ′ and i ′ bands but the stellar mass is not constrained by the optical identification (magnitudes r ′ = 24.87 ± 0.15 and i ′ = 24.25 ± 0.18) obtained by GROND because observations probe the rest-frame UV, where the massto-light ratio can vary by a factor of more than 100.

Spectral energy distribution
The multi-band spectral energy distribution (SED) of the afterglow of GRB 090323 is shown in Fig. 6. It is well fit with a power-law of index β opt = 0.65 ± 0.13 and a small amount of Small Magellanic Cloud (SMC)-like dust extinction (Bouchet et al. 1985 . A pure power-law does not provide an acceptable fit (χ 2 = 26/5 d.o.f.), and Large Magellanic Cloud (LMC) (Fitzpatrick 1986) and Milky Way (MW) (Seaton 1979) dust models are strongly ruled out because of the absence of the 2200 Å feature, which would be in the z ′ band at this redshift. The flux difference between the extrapolation of the SED defined from i ′ to the K S and the g ′ and r ′ bands is nicely consistent with the spectroscopic redshift of z = 3.57. The final photo-z obtained with the calibrated GROND data is 3.44 +0.18 −0.16 .  Fig. 6. GROND broadband spectral energy distribution of the afterglows of GRB 090323, GRB 090328 and GRB 090902B. Note that the g ′ band points for GRB 090323 and GRB 090902B are overlapping.
In the standard fireball model (see e.g. Piran 2004; Mészáros 2006, for reviews) the spectral index, β of the afterglow is connected to the temporal index, α, via the closure relations. These depend on the type of circumburst medium, location of the characteristic frequencies in the synchrotron spectrum, and the evolutionary stage of the afterglow (e.g. Rees & Meszaros 1998;Sari et al. 1998;Chevalier & Li 2000;Zhang & Mészáros 2004;Panaitescu 2005;Panaitescu et al. 2006;Zhang et al. 2006;Racusin et al. 2009). However, afterglow light curves often display complex temporal behaviour with shallow decays, plateaus, rebrightenings, flares or smooth breaks such that limited sampling strongly affects the inferred light-curve slope (e.g. Lipkin et al. 2004;Nousek et al. 2006;Evans et al. 2009;Greiner et al. 2009b). Hence, here and later, we derive the spectral index from the optical to NIR SED, and use the temporal index from the light curve fitting and the closure relations to obtain constraints on the jet properties. In particular, for a constant spectral index and value of p, the expected light curve slopes significantly differ (∆α is between 0.5 and 1.3 depending on the model) between the pre and post jet break evolution, which facilitates a reliable discrimination between these two regimes.
For GRB 090323, the obtained spectral index β opt = 0.65 ± 0.13 is compatible with the observational frequency ν opt being between the typical synchrotron frequency ν m and the cooling frequency ν c for an ISM (density n = constant) or wind (n ∝ r −2 ) environment.
The power-law decay index of the light curve, α, is expected to be ∼1 in the pre-jet break, ISM case and is not compatible with the observations. A value of α ∼ 2.3 is obtained in a post-jet break evolution with significant lateral spreading of the ejecta, (both ISM and wind Chevalier & Li 2000)), and α ∼ 1.7 (ISM) or α ∼ 2.0 (wind) if lateral spreading of the ejecta is not considered (Panaitescu 2005). The latter case (post-jet break without spreading) is compatible within 1σ with the observed α while the case with spreading is compatible at the 2σ level. In the above cases the electron spectral index p can be derived from β = p−1 2 and a value of p = 2.3 ± 0.3 is obtained. The value of α for GRB 090323 indicates that the jet break occurred before the start of the observations. However it is not possible to distinguish between the ISM and wind environments. The post-jet evolution implies that the half opening angle of the jet is smaller than θ ISM jet 2.1 • in an ISM-type circumburst environment (Sari et al. 1999;Frail et al. 2001). The beaming corrected energy emitted in γ-rays E γ (1 keV to 10 GeV) of GRB 090323 is thus E γ 3.3 × 10 51 erg. In a wind environment (Chevalier & Li 2000) the opening angle of the jet would be θ wind jet 1.1 • and the beaming corrected energy would be E γ 1.0 × 10 51 erg. (1 keV to 10 GeV). The following values are assumed throughout the paper A * , n 1 , η 0 .2 = 1.

Light curve
The GROND multicolor light curve of the afterglow of GRB 090328 starting at T 0 +137 ks is presented in Fig. 7 and is remarkably similar to that of GRB 090323 (Fig. 5). It is well described with a single power-law of index α = 2.27 ± 0.04, which is typical of a post-jet break slope , plus a superimposed bump component and a constant host contribution dominating at later times (χ 2 =35/39 degrees of freedom). There is minor evidence for a slightly redder spectra in the bump component peaking at ∼350 ks post burst, but given the uncertainties in the measurement however, the spectral index is consistent with the afterglow slope at the 2σ level.

Spectral energy distribution
The broadband SED of the afterglow was constructed from data taken in the first epoch where the last epoch was used as a reference frame. Subtracting the host fluxes was performed using two methods: firstly the fluxes from the light curve fitting were used and directly subtracted from the first epoch magnitudes, where measurement errors were propagated accordingly. For the K S band, where the host is undetected, we used the mean flux between the H band detection and K S band upper limits with errors over the complete range. The effect of the host magnitude error on the afterglow SED, however, is negligible as its uncertainty is dominated by the absolute photometric accuracy. Secondly, the reference frame was subtracted from the 1st epoch image in each band, and magnitudes were derived using the difference image directly. Both methods resulted in consistent results. As the  second method tends to underestimate the true photon noise in the image, results from the direct flux difference with propagated errors are used in the following. The afterglow SED was fitted using a power-law modified by dust reddening as shown in Fig. 6. A pure power-law provides a reasonably good fit to the data with a χ 2 of 4.5/5 d.o.f. The resulting power-law index of β opt = 1.46 +0.07 −0.08 , however would be surprisingly red, both theoretically (Sari et al. 1998) and observationally (Kann et al. , 2009a if it were intrinsic to the afterglow. In fact, a small amount of extinction of A host V = 0.22 +0.06 −0.18 mag with a SMC type dust attenuation law provides a better fit (χ 2 = 2.6/4 d.o.f) and a more typical spectral power-law index of β SMC opt = 1.19 +0.21 −0.19 . Though statistically not necessarily required, it seems thus very likely that there is mild reddening by dust in the circumburst environment. Due to the redshift of z = 0.7354, the optical data obtained hardly probe the rest frame UV regime and do not probe the region of the 2200 Å bump. Therefore, LMC and MW dust models cannot be distinguished and return comparable results, with spectral indices of β LMC opt = 1.07 +0.29 −0.12 , β MW opt = 1.16 +0.24 −0.21 and A host V between 0 and 0.4 mag.
A spectral index of β opt ∼ 1.2 is only consistent with temporal index α = 2.27±0.04 in a post jet break evolution and the ISM or wind, slow cooling case in the spectral regime above the cooling frequency ν opt > ν c (Sari et al. 1998;Chevalier & Li 2000;Zhang & Mészáros 2004). In this case the electron energy index p is p = 2β = 2.4 +0.4 −0.4 as obtained from the spectral index. Hence, the observations indicate, that the jet break must have occurred before the first GROND observations. This yields an upper limit of θ ISM jet 5.5 • or θ wind jet 4.2 • (Sari et al. 1999;Frail et al. 2001), which constrains E γ (1 keV to 10 GeV) to E γ 4.6 × 10 50 erg for an ISM or E γ 2.7 × 10 50 erg for a wind type environment.

The host galaxy of GRB 090328
The optical spectrum of the afterglow of GRB 090328 is presented in Fig. 9  . This is about a factor of two higher than the mean < SFR >= 2.5 M ⊙ yr −1 derived for a sample 43 GRB hosts in the redshift interval 0 < z < 3.4 (Savaglio et al. 2009), and about one order of magnitude higher than in the LMC. The spectrum as been corrected for slit losses, however we note that the star formation rate quoted above is a fraction of the total and that the true SFR may be slightly higher.
The spectrum has a number of absorption lines at the same redshift of the [OII] line, associated with FeII, MnII, MgII, MgI and CaII. We derived column densities using the curve-ofgrowth analysis (Spitzer 1978), and cross-checked the results with the apparent optical-depth method (Savage & Sembach 1991). Identification, rest-frame equivalent width EW r and measured column densities are in Table 1. The effective Doppler parameter given by the best fit is, independently for FeII and MnII, b ≃ 78 km s −1 . We used this value to derive the column density for MgI and CaII. We estimate a very uncertain ironto-manganese relative abundance [Fe/Mn] = 0.08 +0.53 −0.37 . No intervening absorption systems have been identified. The redshift obtained from the optical spectrum is z = 0.7354 ± 0.0003.
The SED of the host galaxy of GRB 090328 is presented in Fig. 8 and was fit using hyperZ (Bolzonella et al. 2000) and models from Bruzual & Charlot (2003). The fitted galaxy spectral templates include elliptical, different types of early and late spirals, irregular and starburst galaxies at various ages and different star formation rates and solar metallicity. In addition, the spectral templates were modified by an extinction term, resembling the dust attenuation law of the SMC. The redshift of the galaxy was fixed at z = 0.7354, which is consistent with the apparent 4000 Å Balmer break being within the r ′ band.
The host galaxy of GRB 090328 is best fit with a starburst galaxy template, but all other templates provide acceptable fits. Similar to the afterglow SED, there is evidence for dust reddening from the host observations of the order of A host V = 0.5 +0.5 −0.3 in the starburst galaxy (although A host V is not well constrained ), with  Table 1. a best-fit age of the dominant stellar population of ∼50 Myr. By folding the restframe best-fit template spectrum with a typical B and K filter response, we derive absolute AB magnitudes of the host galaxy of GRB 090328 of M B = −20.67 and M K = −20.17, well within the samples of previous long GRB host galaxies (e.g., Christensen et al. 2004;Savaglio et al. 2009). The stellar mass of the host can be approximated from the best fitting template to M * ≃ 10 9.4 M ⊙ . When using the empirical relation derived by Savaglio et al. (2009) between K-band absolute magnitude and stellar mass: we derive M * ≃ 10 9.24 M ⊙ with a typical uncertainty of a factor 2, consistent with the previous estimate and similar to the stellar mass of the LMC. Using the dust-corrected star formation rate SFR = 4.8 M ⊙ yr −1 (see note above on slit loss) and the above stellar mass, the specific star formation rate per unit mass rate is 2.8 Gyr −1 again well within the sample of previous long burst host galaxies (Savaglio et al. 2009).
Using the galaxy distribution function in the corresponding redshift interval from z = 0.6 − 0.8 from Willmer et al. (2006) to compare the host's absolute magnitudes to the population of field galaxies, the luminosity of the GRB host is estimated to ∼ 0.46 − 0.50L * . We obtain a dust-corrected star formation rate per unit luminosity of SFR / L B ∼ 10 M ⊙ yr −1 L −1 * for this host galaxy. Hence, the specific star formation rate in the host galaxy of GRB 090328 is consistent with the values for the long bursts presented by Berger (2009).

Afterglow and host association
The properties of the afterglow of GRB 090510 are not well constrained due its faintness at the time of the GROND observation and the blending with the nearby galaxy. Fitting the afterglow SED with a reddened power-law results in an unconstrained power-law index, and a 3σ upper limit for the intrinsic extinction of A host V ≤ 1.1 mag, assuming an SMC-type extinction law. MW-or LMC-type extinction cannot be distinguished and result in comparable values for β opt and A V .
The putative host galaxy at a distance of 1. ′′ 2 , corresponding to a projected distance of 9.4 kpc from the afterglow to the brightest emission in the nearby galaxy, has r ′ -and z ′ -band magnitude of 23.4 and 22.3 mag, respectively and extends to the afterglow position in the GROND images of 0. ′′ 9 seeing (Fig. 3). According to object counts from the Hubble Ultra Deep Field (Beckwith et al. 2006) and the GOODS (Giavalisco et al. 2004) one would statistically expect less than 0.008 unrelated objects of similar brightness at this distance to the afterglow, corresponding to a chance probability of 0.8%. For comparison, the chance coincidence probability for the galaxies ∼5 ′′ south and ∼8 ′′ west of the afterglow are ∼7%, and ∼30% respectively.
In the merger of two compact objects, the progenitors may be kicked outside their host galaxy at scales of 1-100 kpc, with a distribution peaking at several kpc for different galaxies and merger conditions (Fryer et al. 1999). Most of the convincing host galaxy associations however are found for host positions which coincide or overlap with the afterglow (e.g., Fox et al. 2005;Levan et al. 2006a;Kann et al. 2008;Troja et al. 2008).  Table 2.
Short burst host samples show a typical offset of 0 ′′ −3 ′′ with clustering at small values (Berger et al. 2007;D'Avanzo et al. 2009) indicating that a large kick might not be a dominant process. Recently, based on HST observations of ten short burst hosts, Fong et al. (2010) report that the median offset of short GRBs from their host galaxies is 5 kpc and based on a larger sample of short bursts they report that 25% of bursts have projected offsets of less than 10 kpc. The host galaxy of GRB 090510 certainly falls within this distribution. Berger et al. (2007) argue that even fainter, undetected galaxies at low redshifts z 0.5 are very unlikely to host the dominant fraction of short GRBs. Also in the case of GRB 090510, direct observations of the optical afterglow place a strong constraint of z < 1.5 of the redshift ). Given that the observed properties of the closeby galaxy fit very well into previous samples of short burst host galaxies in terms of brightness and redshift (e.g., Berger et al. 2007;Savaglio et al. 2009), and the probability of chance association is very low, this galaxy can be confidently associated with the burst and is hence very likely the host of GRB 090510.
We further report that the redshift of the galaxy ∼5 ′′ south obtained using VLT/FORS is the same as that of the host galaxy based on the detection of [OII] and Hβ emission lines. Using a distance scale of 7.8 kpc/ ′′ this translates to a projected distance of ∼ 40 kpc and suggests that the host and this galaxy are part of a group or cluster association. Consequently in the unlikely event that this second galaxy were indeed the host galaxy, the distance to the source would therefore remain unchanged. The source to the east ∼ 8 ′′ away appears stellar and all the remaining nearby bright galaxies have chance coincidence probabilities of order ∼100%.
The [OII] and Hβ luminosities are L [OII] = (5.4 ± 0.6) × 10 40 erg s −1 and L Hβ = (2.1±0.6)×10 40 erg s −1 , respectively. The star formation rate, not corrected for dust extinction, derived from the [OII] luminosity is SFR = 0.30 M ⊙ yr −1 , consistent with the SFR = 0.26 M ⊙ yr −1 derived from the Hβ luminosity. The SFR conversions are those proposed by Savaglio et al. (2009) for GRB hosts. The derived value is one order of magnitude lower than the mean SFR measured in long GRB hosts (Savaglio et al. 2009) and consistent with relatively faint UV emission derived from the SED of the host (Fig. 8) and at the lower end of the distribution of SFRs of short GRBs (Berger 2009). Correcting for dust extinction yields a star formation rate of ∼ 0.60 M ⊙ yr −1 . We note as in the case of GRB 090328, the spectrum as been corrected for slit losses, however we note that the star formation rate quoted above is a fraction of the total and that the true SFR may be slightly higher.
The SED fitting for the host of GRB 090510 has been performed in a similar manner to that of GRB 090328 (Section 4.2.3). The redshift of the galaxy was fixed to the spectroscopic value of z = 0.903. Again all galaxy templates provide an acceptable fit, where the best fit is obtained with an elliptical galaxy. We note that the SFR is in the range of local star-forming elliptical galaxies (Huang & Gu 2009). From the best fit we derive AB absolute magnitudes M B = −19.91 and M K = −20.51, and slight evidence for a dust extinction reddening similar to the SMC with a visual extinction of A host V = 0.7 +0.2 −0.4 . The stellar mass of the host can be estimated directly from the galaxy templates to M * ≃ 10 9.6−10.2 M ⊙ . The dust-corrected star formation rate per unit mass is 0.04 to 0.15 Gyr −1 , significantly below that of the host galaxy of GRB 090328.
Using the galaxy distribution function in the corresponding redshift interval from z = 0.8 − 1 from Willmer et al. (2006), the luminosity of the GRB host is estimated to ∼ 0.26 − 0.30L * , and we obtain a dust-corrected star formation rate per unit luminosity of SFR / L B ∼ 2 M ⊙ yr −1 L −1 * for this host galaxy, consistent with the albeit wide range of SFR / L * values for the short bursts presented by Berger (2009).

Light curve
The early (t < T 0 + 1.10 Ms) GROND multi-color light curve of the optical/NIR afterglow of GRB 090902B is well fit by a single power-law of index α 1 = 0.94 ± 0.02. A late epoch of imaging obtained with the VLT 23 days after the burst indicates a break in the light curve at the 2σ level. No further followup observations could be performed due to the proximity of the GRB field to the Sun. Consequently the light curve coverage at late times is sparse and the evolution not fully conclusive. Fixing the late slope to a conventional post jet break slope of α 2 = 2.2 with a relatively sharp break as shown in Fig. 11 results in a break time of t ∼ 1.5 Ms (χ 2 = 43/35 d.o.f).

Spectral energy distribution
The optical/NIR SED of the afterglow of GRB 090902B is well described with a simple power-law template (Fig. 6). The spectral index is β opt = 0.79 +0.05 −0.19 , with 3σ upper limits of the intrinsic visual extinction of A host V 0.36 with an SMC/LMC type reddening and A host V 0.24 assuming a MW-like dust attenuation law. This is in agreement with the light curve slope of α ∼ 0.94 in the fireball model, which is typical for a pre-jet break evolution with an ISM surrounding and with an observational frequency between ν m and ν c . Following Sari et al. (1999) and Frail et al. (2001), the half opening angle for a break time of t ∼ 150 Ms would be 7.2 • , or setting t break 1.1 Ms as a conservative lower limit θ jet 6.4 • . Together with the energetics of the prompt emission, this implies a beaming corrected energy of E γ 2.2 × 10 52 ergs. We note that a wind environment would significantly relax the lower bounds on the opening angle and thus E γ , however it is not consistent with the expectations for α and β opt in the standard models. The power-law index p of the electron population is then p = 2.58 +0.10 −0.38 as obtained from the spectral index which is in reasonable agreement with the values estimated from the early LAT emission under  Fig. 11. GROND and VLT optical/NIR light curve of GRB 090902B. The data are fitted with a broken power-law. Due to the proximity of the source to the sun, the light curve coverage is sparse at later times, and the evolution is ambigous. There is however indication of a break at around 150 Ms post burst. Shown is the fit with a fixed late slope and relatively sharp break.
the assumption that it is generated by an external forward shock (Kumar & Barniol Duran 2009a).

Discussion
We have presented the optical and NIR follow-up observations of four bursts with high energy emission detected by Fermi/LAT and spectroscopic investigations of two of those events. Since the end of the commissioning phase Fermi/GBM has detected several hundred GRBs, many of which were in the field of view of the LAT instrument, however only fourteen of those have been also detected in by LAT (up to the end of January 2010). Follow-up of these bursts is of special interest due to the broadband coverage of the prompt emission (keV to GeV) which allows the spectral properties and energetics of these events to be constrained. It is noteworthy that the redshifts of bursts detected in the GeV regime by the LAT range from the relatively low-redshift GRB 090328 to the high-reshift of GRB 080916C and that the long bursts have high fluences and very high isotropic gamma-ray energies (although GRB 090328 has a smaller isotropic energy release due to its lower redshift). The properties of these events presented in this paper and another Fermi/LAT event, GRB 080916C, including their afterglows and energetics are summarized in Table 3. Three of the bursts presented in this paper belong to the class of long gamma-ray bursts while GRB 090510 is a short burst. The existence of two classes of gamma-ray bursts has been established for some time (e.g., Mazets et al. 1981;Norris et al. 1984;Hurley 1992;Kouveliotou et al. 1993). It is generally accepted that long GRBs have their origins in massive star progenitors because of their association with core-collapse supernovae (SNe) (e.g. Galama et al. 1998;Hjorth et al. 2003;Malesani et al. 2004;Zeh et al. 2004) and occurrence in starforming galaxies (Bloom et al. 2002) and in highly star-forming regions therein (Fruchter et al. 2006). For a recent review of gamma-ray bursts see Vedrenne & Atteia (2009).
The origin of short GRBs is still open, with mergers of compact objects being the leading concept (e.g., Eichler et al.  (Kann et al. , 2009a(Kann et al. , 2008. GROND and TLS data only are presented for the the long Fermi bursts while the early data points from GRB 090510 are from UVOT ) and include a point from Olofsson et al. (2009). All afterglows have been corrected for Galactic extinction and for host-galaxy contribution, where applicable. All afterglows have been additionally corrected for host-galaxy extinction, and have been shifted to a redshift of z = 1. Long GRBs are shown as light grey lines, whereas short GRBs are thicker grey lines with symbols. Squares are detections, and downward pointing triangles upper limits.
1989; Gehrels et al. 2005;Hjorth et al. 2005;Fox et al. 2005;Levan et al. 2006b). For a recent review of short bursts and their progenitors see Nakar (2007) and Lee & Ramirez-Ruiz (2007). In contrast to long GRBs, the afterglows of the short bursts are generally fainter (Kann et al. 2008;Nysewander et al. 2009), making it difficult to obtain an absorption line spectrum of the afterglow itself. Almost exclusively, the distances to short bursts are obtained by associations with and observations of the (putative) host galaxies (e.g., Berger 2009;Graham et al. 2009). See however Levan et al. (2007). One possible exception is the case of GRB 090426 at z = 2.609 (Levesque et al. 2010;Antonelli et al. 2009;Zhang et al. 2009) for which absorption line spectroscopy of the afterglow was successfully obtained, but may be considered on balance by the authors to be associated with the population of long gamma-ray bursts. The host galaxies of the short bursts were found initially to be mostly elliptical galaxies, e.g., GRB 050509B (Gehrels et al. 2005;Bloom et al. 2006) and GRB 050724 (Berger et al. 2005;Gorosabel et al. 2006) and that an old progenitor population was required. However, as the sample has since increased, all types of galaxies, from elliptical to star-forming, have been found to be associated with short bursts and the majority of short bursts are now found to reside in star forming galaxies (e.g., Berger et al. 2007;Berger 2009;Prochaska et al. 2006). A comparison of the short and long burst hosts by Berger (2009) revealed that the former have higher luminosities, and that it is unlikely that both samples are drawn from the same underlying galaxy distribution (p = 10 −3 ) and that the short burst hosts properties match well to field galaxies in the range z ∼ 0.1 − 1.

Afterglow properties and energetics
The bursts were followed up in the optical and NIR starting at times ranging from 6 hours to 1.6 days post trigger. In all cases a bright afterglow was detected and followed for many days, with the exception of GRB 090510 which was only detected on the first night post-trigger. The afterglow light curves of the three long bursts are well sampled in the seven GROND bands. The two March events are well fit by power-laws with indices consistent with post-jet break models constraining the opening angle to be before the start of GROND observations. In the case of GRB 090323 the opening angle (θ ISM jet ) is less than 2.1 • while for Table 3. Summary of the burst properties including the redshifts, isotropic gamma-ray energy (E γ,iso ), the power-law decay index of the afterglow light curve (α), the afterglow spectral index (β opt ), the template employed to constrain the host galaxy extinction (A host V ), limits on the half opening angle (θ ISM jet ) and on the beaming corrected gamma-ray energies (E γ ) and finally the star formation rates obtained for two of the host galaxies. · · · · · · 6 • 4.8 × 10 52 · · · GRB 090510 0.903 1.1 × 10 53 · · · · · · · · · ≤1.1 · · · · · · 0.30 a E γ,iso values are quoted in the energy range 1 keV -10 GeV, with the exception of GRB 080916C (10 keV-10 GeV: Abdo et al. 2009a) and GRB 090510 (10 keV -30 GeV : Abdo et al. 2009c).
b The opening angle and E γ are lower in a wind medium for GRB 090323 and 080916C. c Chornock et al. (2009). d LMC and MW templates fit the data equally well and it is not possible to distinguish between models. e Cucchiara et al. (2009). f A value of A host V 0.23 is obtained using the MW template. g Values of α, β opt and the jet break lower limit are taken from Greiner et al. (2009a). h The true SFR may be larger than measured because of slit losses.
GRB 090328 the observations began later and the opening angle (θ ISM jet ) is less than 5.5 • . The slope of the afterglow decay of GRB 090902B is flatter and consistent with a pre-jet break value, constraining the opening angle (θ ISM jet ) of the jet to be 6.4 • . Knowledge of the opening angle allows us to correct the isotropic energy and calculate limits on E γ , the beaming corrected energy release in γ-rays (Sari et al. 1999;Frail et al. 2001). The E γ of GRBs 090323 and 090328 is in a range compatible with the majority of GRBs observed to date. In the case of GRB 090902B, the late break places at least a lower limit on the opening angle and therefore a very large beaming corrected energy output of 2.2 × 10 52 erg is obtained. To date this is one of the largest values of E γ measured and indicates that the large E γ,iso values of Fermi/LAT GRBs are not always due to a narrow opening angle where we happen to be located in the cone. The value of E γ obtained for the three long bursts ranges from less than 5×10 50 erg to greater than 2×10 52 erg, spanning over a factor of at least 30. Values in excess of 10 52 erg have been reported for another Fermi/LAT detected burst ) and a number of Swift bursts ) and indicate that the distribution of E γ is broad and not compatible with a standard candle (Frail et al. 2001;Bloom et al. 2003). Furthermore these very luminous GRBs with high values of the redshift corrected E peak are not compatible with the E peak,i -E γ relationship as shown in Fig 13 which was obtained for GRBs with in general lower values of the peak energy (Ghirlanda et al. 2004;Campana et al. 2007). While being consistent with the Amati (E peak,i -E iso ) relation at the 2σ level (see also Amati et al. 2009), albeit with a large scatter (see also e.g. Butler et al. 2007), all four LAT bursts are inconsistent with the Ghirlanda relation (Ghirlanda et al. 2004) at 2σ using a standard approach. Two GRBs (080916C and 090902B) have too large beaming corrected E γ for their E peak,i , while E peak,i for GRBs 090323 and 090328 is too high to fit the Ghirlanda relation. In the earlier case, subtracting the energy in the additional power-law component of the prompt emission (Abdo et al. 2009c), a wind environment, or a two-component model could possibly relax the constraints on E γ directly (extra-component) or indirectly by reducing the lower limit on the opening angle (Wind-case), or by disconnecting the beaming of the prompt emission from the late afterglow (two-component model).
Although somewhat contrived, multi-component jets (e.g. Berger et al. 2003;Racusin et al. 2008;Krühler et al. 2009a) would allow the γ-ray emission to be tighter collimated than the late afterglow, and hence reduce the measure on the opening angle obtained from the afterglow light-curve to an upper limit. Similar arguments can not be invoked for GRBs 090323 and 090328, where the latter one in particular is a significant outlier from the E peak,i -E γ relationship (Fig. 13).

Host properties
Host galaxies were detected for two of the long bursts, GRBs 090323 and 090328, and the short burst GRB 090510. The host galaxy of GRB 090323 was detected in the r ′ and i ′ bands but further constraints could not be obtained due to the faintness of the source and the high redshift. The host galaxy of GRB 090328 is best fit with a burst galaxy template (although all other templates provide acceptable fits) and the absolute magnitudes are comparable with previous long GRB host galaxies (e.g., Christensen et al. 2004;Savaglio et al. 2009). A dust-corrected star-formation rate of 4.8 M ⊙ yr −1 , as well as the star formation rate per unit mass (2.8 Gyr −1 ) and per unit luminosity ( ∼ 10 M ⊙ yr −1 L −1 * ) derived from the afterglow spectrum are consistent with the long burst sample (Savaglio et al. 2009;Berger 2009). The host galaxy of GRB 090510 was detected by GROND and a spectrum was obtained by the VLT. A star formation rate of 0.30 M ⊙ Gyr −1 (not corrected for dust) was derived and the host galaxy is best fit by an elliptical template, however all templates provide acceptable fits. The dust-corrected star formation rate per unit luminosity of ∼ 2 M ⊙ yr −1 L −1 * for this host galaxy is consistent with the values for the short bursts presented by Berger (2009). The properties of the host galaxies of two bursts with high-energy emission detected by Fermi/LAT, Fig. 13. Amati and Ghirlanda relation challenged by Fermi/LAT bursts. The black crosses, the black dashed line, and the grey shaded areas represent the data and best-fit E peak,i -E iso relation and its 1σ, 2σ and 3σ scatter as compiled in Amati et al. 2008. The grey crosses, the grey dashed line, and the golden shaded areas represent the data and best-fit E peak,i -E γ relation and its 1σ, 2σ and 3σ scatter as compiled in Ghirlanda et al. 2007. Red points are the 5 Fermi/LAT bursts with respect to the Amati relation. The red circles are calculated in the bolometric energy range (1 keV to 10 MeV) and the red diamonds are calculated over the energy ranges quoted in Table 3. The four long LAT bursts are consistent with it at the 2σ level. Blue points with arrows are the limits on the four long Fermi/LAT bursts with respect to the Ghirlanda relation. The blue circles symbols indicate the data in the bolometric energy range (1 keV to 10 MeV) and the diamonds are again the values from Table 3. Over the bolometric energy range all four are inconsistent with the relation at 2σ level and over the wider energy range, three are inconsistent with it at the 3σ level and GRB 090323 is consistent. the long burst GRB 090328 and the short burst GRB 090510, are not exceptional with respect to pre-Fermi samples.

Optical afterglows of LAT-detected GRBs in comparison with the previous sample
With knowledge of the distance from the redshift determination and the local dust extinction from the SED fitting, we are able to compare the afterglows of our sample of LAT-detected GRBs (where we also include GRB 080916C, Greiner et al. 2009a) with a large sample of well observed GRB afterglows presented in Kann et al. (2006) (pre-Swift long GRBs), Kann et al. (2009a) (Swift-era long GRBs) and Kann et al. (2008)(Swift-era short GRBs). Using the method detailed in Kann et al. (2006), the afterglows are corrected for host-galaxy extinction and all shifted to a common redshift of z = 1, so they can be compared directly both in terms of temporal evolution as well as luminosity. The LAT-detected-GRB afterglows are shown in comparison to the total sample in Fig. 12.
In comparison to the sample of long GRB afterglows, the three afterglows presented in this work as well as GRB 080916C are seen to be diverse. The optical afterglow of GRB 090323, especially at early times, is one of the most luminous afterglows ever detected. The steep decay from discovery on gives it a more typical luminosity at later times. The afterglow of GRB 080916C  (Kann et al. , 2009a. The comparison sample comprises 50 Swift-era and 21 pre-Swift afterglows.
lies close to the mean of the sample distribution. At discovery, the afterglow of GRB 090328 has a similar luminosity, but again, a steep decay from the onset implies that it becomes much less luminous at later times. The afterglow of GRB 090902B has a flatter slope and the very late break causes it to be one of the most luminous afterglows known at late times.
Observationally, however, the four optical afterglows of the long Fermi/LAT bursts are exceptionally bright. The histogram of observed magnitudes from the sample in Fig. 12 is shown in Fig. 14 at a common time of 2 days after the trigger in the observers frame. At this time, all four bursts are brighter than the mean brightness (R C = 22.2) of the Swift-era comparison sample. Furthermore, the sample is not complete, comprising 50 out of a total of 370 bursts with X-ray afterglows as a rough indicator for observability between the launch of Swift and May 2009. Accounting for the observational bias, and assuming in the extreme case that the comparison sample is complete at the bright end of the distribution, then all four LAT afterglows would lie in the brightest 5%, with GRB 090328 being the brightest (top 2%, i.e. 8th brightest of Swift-era or 16th of all afterglows).
In comparison to the afterglows of other short GRBs (grey lines), it can be seen that the afterglow of GRB 090510 is among the most luminous (in agreement with its extreme prompt energy release, Abdo et al. 2009a, and the correlation between prompt energy release and afterglow luminosity, Kann et al. 2008;Nysewander et al. 2009;Gehrels et al. 2008), at least at early times. After just a few hours, it goes over into a steep, probably post-jet break decay, and at one day, its luminosity is comparable to the few other short GRB afterglows detected at this time. The upper limit derived from the second GROND epoch is not restrictive.
The host galaxy extinction found in all three long GRB afterglows lies within the range of A V ≈ 0.2 which is typical for bright optical/NIR afterglows (Kann et al. 2009a), with no extinction detected in the case of GRB 080916C .
The intrinsic spectral slopes β opt are also diverse, with that of GRB 090328 being a typical value for an afterglow with a cooling break ν c redwards of the optical/NIR regime, while the obtained optical/NIR data for GRBs 090323 and 090902B very probably probe the spectral region of ν m < ν opt < ν c .

Bumps in the afterglow light curves
Variability superimposed on to a power-law decay as seen in GRBs 090323 and 090328 (Figs. 5 and 7) is frequently observed in well sampled optical afterglow light curves (e.g., Lazzati et al. 2002;Greiner et al. 2003;Lipkin et al. 2004;Updike et al. 2008;Perley et al. 2008;Krühler et al. 2009b). Given the low amplitude, smooth structure, time, and roughly equal color to the power-law component, the observed variability in the light curve of GRB 090323 is most likely related to a process intrinsic to the generic forward shock or the circumburst medium, and not to late inner engine activity. The most commonly inferred mechanism for producing variability as observed in the optical afterglow of GRB 090323 is additional energy injected into the afterglow via refreshed shocks (Rees & Meszaros 1998), although an internal origin, i.e. an optical flare related to late central engine activity can not be ruled out. In GRB 090328 there is minor evidence for a slightly redder spectra in the bump component peaking at ∼350 ks post burst, which would support a different origin of the observed variability than forward shock emission. Given the uncertainties in the measurement however, the spectral index is consistent with the afterglow slope at the 2σ level.
Various models are emerging which explain delayed GeV with respect to the keV -MeV emission in Fermi bursts, among others an external shock model (e.g. Ghisellini et al. 2009;Kumar & Barniol Duran 2009b) or a hadronic model ). Razzaque et al. (2009) demonstrate that synchrotron radiation from cosmic-ray protons accelerated in a GRB, delayed by the proton synchrotron cooling timescale in a jet of magnetically-dominated shocked plasma moving at highly relativistic speeds, explains the delayed GeV emission as observed in many Fermi/LAT bursts. A second generation prompt electron synchrotron component from attenuated proton synchrotron radiation makes enhanced soft X-ray to MeV gammaray emission, which is, however, too weak to be detected. In this scenario the implicit assumption is made that the outflow impacts either a previously ejected shell or a pre-existing uniform density medium with an extent of r ∼ 6 × 10 16 (Γ/1000) 2 (δt/10s) cm. For a 3 second delay in GRB 090323/28 this corresponds to a light travel time of 8-10 days, surprisingly similar to the timing of the bump in the optical/NIR light curve as observed with GROND. It is conceivable that the delayed GeV emission and the late, likely achromatic bumps in the afterglow light curve of GRBs could be related. One possible scenario would be Compton scattering of this proton synchrotron emission at the shell material, thus creating an achromatic bump with a powerlaw of photon index -2 (or even steeper) depending on the efficiency of the scattering. We note however, that density enhancements in the circumburst medium alone can not account for significant rebrightnings or strong bumps observed in optical afterglow light curves (Nakar & Granot 2007). To reproduce the observed variability in optical afterglow light curves, a second process is required, possibly the variation of the microphysical shock parameters along the shock condition as suggested by Kong et al. (2010).

Conclusion
We have presented the multi-colour afterglow observations of four bursts with high energy emission including three long and one short burst. In addition, we present spectroscopic observations and redshift determinations of two of these bursts. It is now possible for the first time to combine the high energy information from the prompt emission with the afterglow, host galaxy and redshifts of these sources. Follow-up of GRBs detected by Fermi with the seven band imager, GROND, allows simultaneous determination of the temporal, α, and spectral, β opt , indices of the burst afterglows. The long bursts exhibit power-law decay indices (α) from less than 1 to ∼2.3 and spectral indices (β opt ) from 0.65 to ∼1.2 which are fairly standard for GRB afterglows. Moreover, an estimation of the jet break time is vital to determine the broadband properties and energetics of these events. The redshifts of the long bursts span the range from 0.7354 to 3.57 and the beaming corrected energies differ by a factor of at least ∼30 with a value of E γ 2.2 × 10 52 erg for GRB 090902B and are not compatible with a standard candle. Interestingly, the higher redshift bursts detected by the LAT are very luminous in the keV-MeV region of the prompt spectrum, have high E γ,iso and E γ and also exceptionally bright afterglows. The lower redshift GRB 090328 is the exception in this respect, with a more standard E γ,iso as compared to Swift GRBs. The host galaxies of GRB 090328 and GRB 090510 fit well within the distributions of the long and short burst host galaxies respectively. Future photometric and spectroscopic follow-up of these rare Fermi bursts with high-energy emission is crucial to further investigate the energetics and host galaxies of a larger sample.