A&A 405, 859-866 (2003)
DOI: 10.1051/0004-6361:20030648
C. Gunnarsson - A. Goobar
Department of Physics, Stockholm University, AlbaNova, 106 91 Stockholm, Sweden
Received 21 November 2002 / Accepted 24 April 2003
Abstract
We investigate the potential of using massive clusters as
gravitational telescopes for searches and studies of supernovae of
Type Ia and Type II
in optical and near-infrared bands at central wavelengths in the interval
0.8-1.25 m.
Using high-redshift supernova rates derived from the measured star
formation rate,
we find the most interesting effects for the detection of
core-collapse SNe in
searches at limiting magnitudes
mag, where
the total detection rate could be significantly enhanced and the
number of detectable events is
considerable even in a small field. For shallower searches,
24 mag,
a net gain factor of up to 3 in the discovery rate could be obtained, and yet
a much larger factor for very high source redshifts. For programs
such as the GOODS/ACS transient survey, the discovery rate of supernovae
beyond
could be significantly increased if the observations were done
in the direction of massive clusters.
For extremely deep observations,
mag, or for very
bright SNe (e.g. Type Ia) the
competing effect of field reduction by lensing dominates, and fewer supernovae
are likely to be discovered behind foreground clusters.
Key words: cosmology: gravitational lensing - cosmology: distance scale - galaxies: clusters: general - stars: supernovae: general
Massive galaxy clusters, used as gravitational telescopes (GTs), are extremely useful tools for the studies of faint high redshift galaxies in wavelengths ranging from optical to submillimetre as demonstrated e.g. in Hu et al. (2002); Biviano et al. (2003); Smail et al. (2002); Ellis et al. (2001); Lemoine-Busserole et al. (2002). In this work we explore the potential of GTs for magnifying high redshift supernovae (SNe) in optical and near-infrared (NIR) wavelengths, thereby increasing the chances of detection. A competing effect related to the use of GTs is due to the spreading of the field by the lens, analogous to what happens when looking through a magnifying glass: a smaller, although magnified, portion of the field is actually observed. The effect is sometimes referred to as amplification bias. For any specific field-of-view (FOV) it is not obvious a priori which of the two effects dominates when looking for distant supernovae. The net gain depends upon 1) the lens and source parameters: the mass distribution of the lensing cluster, the intrinsic rate and brightness for core collapse and Type Ia supernovae as a function of redshift. 2) The observational set-up: the limiting search magnitude and the choice of wavelength band. In this paper we consider several scenarios relevant for supernova searches.
Studying supernova (SN) rates at the highest possible redshifts provides critical information for the understanding of cosmic star formation rate. Lensed SNe at high redshifts can in principle also be useful as distance indicators, provided the magnification is known to high precision. In addition, strongly lensed, multiply imaged supernovae may also be used to constrain cosmological parameters through the time delay measurements of the SN images (Goobar et al. 2002a).
Similar work on the feasibility of galaxy clusters as GTs was carried out by Sullivan et al. (2000). They focused on data sets through the HST I814W-band with limiting magnitudes between 26.0 and 27.0. We extend this further by also investigating the Z- and J-bands where we find larger effects than in I-band. This is not surprising as the redder filters are sensitive to higher-redshift sources which in turn are fainter. In addition, we quantify the limiting search magnitudes and source brightness where GTs enhance or deplete the discovery rates. Also (Gal-Yam et al. 2002) conducted an I-band search for SNe in HST fields with galaxy clusters. Earlier feasibility studies were also carried out by Saini et al. (2000). They concentrate their study to the lens cluster Abell 2218 where the assumed background sources were supernovae Type Ia and IIL and constrained their analysis to SN searches at shorter wavelengths.
Our study thus generalises previous work and is more applicable to instruments currently being used for high-z SN searches, e.g. the Z- band search with the ACS camera on HST by the GOODS Treasury Team (Riess et al. 2002; Giavalisco et al. 2002), or future missions such as JWST (formerly NGST). HST/ACS has a considerable depth but a relatively small field-of-view compared to optical cameras used for ground based SN searches. However, the small FOV makes a good match to the high amplification region of massive clusters at the moderate redshifts considered here.
Throughout the paper we use natural units in which c=G=1.
We have also adopted the "concordance'' cosmology with
and h=0.7, where
.
All magnitudes are in the Vega system.
In almost all lensing situations of practical astrophysical interest,
the deflection of the light takes place in a very small
fraction of the total light path. This justifies the common
approximation that all deflection occurs in a plane called the
lens plane. So to compute the lensing properties of a halo we
need to project its mass density onto this plane. This projected mass
density we denote by .
Define the 2-D vector
as the image position (or impact
parameter) in the
lens plane and
as an arbitrary length scale in this plane which can be chosen
suitably to simplify the appearance of the equations. Furthermore we
introduce the corresponding quantities in the source plane:
and
.
and
below are angular diameter distances between observer and
source, observer and lens and lens and source
respectively.
is the source position.
With these definitions, the lens equation in dimensionless form can be
written as
![]() |
Figure 1:
Schematic picture of the lensing situation.
![]() |
Open with DEXTER |
![]() |
Figure 2:
Number of SNe detected per year in redshift bins of
![]() ![]() ![]() ![]() |
Open with DEXTER |
The number of SNe that can be detected is found by integrating the
number density of supernovae over the volume swept out by the FOV.
Naturally, we count all SNe in this volume
that are bright enough to be observed,
either intrinsically or after having been magnified by lensing, see
Fig. 1. The trumpet shaped shaded area beyond
indicates where the SNe are magnified above the detection
threshold,
.
This figure is to be compared to
Fig. 2, which is an example of a typical set of
parameters for Type Ia
SNe and a NFW lens and shows the number of detected SNe, N, per year
and redshift bin
.
The shaded areas (volumes) denoted 1-4 in
Fig. 1 also have counterparts here. Up to the lens
redshift,
(at 0.2 here), there should be no difference between the
cases with and without the lens (area 1). Beyond the lens, there is a
region where the spreading of the field diminishes the area compared
to the case with no lens (area 2), from
to
.
Beyond
there is a small area where the visible
field still is the limiting factor to the number count (area
3). Finally, the magnification is limiting in area 4. Under specific
circumstances, e.g. if the lens is very massive or at low redshifts
and the FOV is
small, another effect might become important. As indicated above,
there are no primary images inside the Einstein radius. In all our
considered models, primary images are
much more
numerous than secondary or tertiary, and
in principle a
situation can arise where
only
giving secondary and tertiary images at this source redshift. This
will lead to a decrease in the number
count when many sources are affected and the effect can be
seen in Fig. 9.
As multiple imaging is possible we have taken this into account by counting secondary and tertiary images that are bright enough as events independent of the primary image, resulting in a slightly higher rate. As indicated above this is a very small effect in all cases considered here.
![]() |
Figure 3:
Supernova rates versus redshift in bins of
![]() |
Open with DEXTER |
We have concentrated our study on Type Ia and Type II supernovae. The Type II SNe are divided into three subclasses: IIn, IIp and IIL, with abundances of 0.02, 0.3 and 0.3 relative to the total core-collapse (cc) SN rate (IIn+IIp+IIL+Ibc+87a-like) following the predictions in Dahlén & Fransson (1999), see Fig. 3. In order to study the power of galaxy clusters as GTs we have neglected the intrinsic brightness dispersion of the supernovae as well as extinction. Thus, the magnitude dispersion considered in this work stems only from gravitational lensing. For comparison, the intrinsic properties (from Richardson et al. 2002) of SNe are listed in Table 1.
We have also taken the conservative viewpoint that the SNe are not discovered at maximum by adding 0.5 mag to the absolute magnitudes in the table. This has the same effect as decreasing the limiting magnitude of the telescope by 0.5 mag.
The simple identification of areas 1-4 in Figs. 1 and
2 is not possible in the case of our Type II SN
plots. This is due to the fact that they have different luminosities
and will be seen out to different
(defined in Fig. 1). This instead gives
rise to three peaks (at most) in the
vs. zplots for Type II SNe as they show the sum of all three types considered.
Table 1: Supernova properties, assuming discovery at maximum luminosity in B-band. Type Ibc and 87a-like events are not considered but are listed for completeness.
We use two different models for the clusters, the singular isothermal sphere (SIS) and the Navarro-Frenk-White profile (NFW) obtained in N-body simulations (Navarro et al. 1997) on which we will focus. The large lens masses considered here seem to be better fitted by the NFW profile (Li & Ostriker 2002) but we will see in Sect. 4 that they do not differ much for the cases we consider.
There is also an ongoing debate on whether the dark matter halos are centrally cuspy or possess a roughly constant density core. The issue is not settled but there is some evidence of cuspy central regions on cluster scales which is what we have assumed in this paper (see e.g. Dahle et al. 2002).
For a description of the lensing properties of the SIS, see e.g. Schneider et al. (1992, , Ch. 8.1.4). The SIS only has primary and secondary images as the tertiary image always is infinitely faint.
Gravitational lensing by NFW halos is well described in
Wright & Brainerd (2000) and we refer to that paper for further reading.
To find the parameters for a specific cluster mass and redshift we use the
Fortran 77 code charden.f available on
the homepage of Julio Navarro.
The NFW profile produces both primary, secondary and tertiary images.
To meaningfully compare the results for SIS and NFW we have assumed
,
the mass within the virial radius, to be equal for both of them.
The virial overdensity used to obtain the virial radius is computed
using the results found in Bryan & Norman (1998) which is a fit to the
calculation in (Eke et al. 1996). However, the Einstein radii of a NFW and a SIS halo
of equal virial mass will be quite different. For a
halo at
,
for source redshifts between 1.5 and 4. This difference
in ER does not
affect the results strongly as is seen e.g. in Fig. 11.
In Fig. 4 we have plotted the source redshift dependence
of the Einstein radius, within which strong lensing effects occur,
for two different
NFW cluster redshifts. As can be seen, the ER has a
weaker dependence on
for higher redshifts.
![]() |
Figure 4:
Redshift dependence of the Einstein radius, within which
strong lensing effects occur, for two different cluster redshifts,
![]() ![]() |
Open with DEXTER |
The main reason for studying gravitationally magnified supernovae is that it might enable us to explore supernova explosions otherwise too faint to be detected by the search instruments currently used or being planned. Because these SNe are predominantly at very high redshifts we focus on the use of the largest wavelength optical filters and NIR for their discovery. We concentrate on the feasibility to discover magnified SNe in the I-, Z- and J-bands.
We are assuming a circular 16 sq. arcmin field centered on the
galaxy cluster. The considered solid angle is thus comparable to the
foreseen FOV of NIRCam/JWST and
slightly larger than the
dimensions of the Advanced Camera for Surveys (ACS,
),
on HST where SN searches are currently conducted in Z-band by e.g. the GOODS
team
.
Ground based NIR cameras on 8-m class telescopes with deep imaging
capabilities
in J-band like ISAAC/VLT (
),
NIRI/Gemini (
)
or CISCO/Subaru (
)
are also
contained within the maximum FOV considered. We also briefly consider
searches in I-band where
most of the efforts of the SCP and High-Z teams are concentrated, although with significantly
larger FOV. K-corrections were calculated using the SNOC package (Goobar et al. 2002b).
![]() |
Figure 5: Redshift reach as a function of limiting (Z-band) magnitude in the supernova search. For comparison, the ACS/GOODS transient survey has a limiting magnitude around Z=25.5. |
Open with DEXTER |
![]() |
Figure 6:
Number of Type Ia SNe detected per year and
![]() ![]() |
Open with DEXTER |
We have investigated limiting magnitudes ranging from 22nd to 30th
magnitude, reaching up
to a limiting redshift
for some of the SN types considered.
Primarily, we study the properties of quite heavy clusters of virial
mass
as GTs but we also briefly consider a cluster of mass
.
Cluster redshifts between
and
were assumed.
First, we consider the case of observational programs such as the ongoing ACS/GOODS transient survey, comparing the potential gain in the discovery rate if the observations were done looking towards a heavy galaxy cluster. These search observations are done in Z-band, with an approximate limiting magnitude of 25.5 mag.
![]() |
Figure 7:
Number of Type IIn+p+L SNe detected per year and
![]() ![]() |
Open with DEXTER |
![]() |
Figure 8:
Gain factor of Type IIn+p+L SNe detected
vs.
![]() ![]() |
Open with DEXTER |
If nothing else is stated, the lens model we use is a NFW halo.
Figures 6 and 7 show the
expected number of discoveries in a 16 sq. arcmin field for
Type Ia and Type II SNe respectively. In the Ia case, although the search
becomes significantly deeper, the effect of spreading the field
dominates, giving
an overall decrease in the number count, relative to the case with no
lens, of more than 30% for the assumed z-dependence of the rates. The
qualitative
effect proves to be roughly independent of the
lens redshift even though it is less severe for higher cluster
redshifts. For the fainter Type II SNe, the effect is reversed at this limiting
magnitude and more supernovae could be detected at redshifts far beyond the
cutoff without the lens. The quantitative gain can be seen in
Fig. 8 where the gain factor,
,
is shown
as a function of limiting magnitude. At
,
a factor of 1.5 to 1.6 is obtained depending on
the cluster redshift.
Since the simulations were quite
time-consuming, some time was saved by not computing too many points
in each plot. This gives the artificial non-smooth behavior seen in
the figure. The dip
at
for
originates from the fact that
for Type IIL SNe
i.e. slightly less than the lens redshift implying that all SNe of this
kind beyond
the cluster are missed in our "standard candle''
treatment. Increasing
pushes
above 0.5. Thus, a sudden rise is expected since
the Type IIL SNe behind the lens that are sufficiently magnified are
seen. If these SNe were not present, the curve would have followed the
qualitative behavior of the
curve, i.e. decreasing
after
.
Next, we consider the possibility of doing the searches at near-IR wavelengths, e.g. in J-band. In this discussion we ignore the observational difficulties of such a program, e.g. the rather long exposure times required to overcome the high sky noise due to atmospheric emission if a ground based search were to be attempted. Searching with NICMOS on HST would be an alternative. A limiting factor is the small size of the HgCdTe arrays.
To see what effect the FOV and different cluster masses have on the
simulations we plot the gain factor vs.
(defined in Fig. 1) for two different cluster masses
and
in
Fig. 9.
The aforementioned effect that for a small FOV it is possible to end up
inside the Einstein radius where no primary images of sources at some
specific distance are found can be seen for small
for both clusters. Although
depends on the source
redshift, the dependence is rather weak for the higher redshift
range as already noted. Between z=1.5 and z=4 the Einstein radius
changes by about
20% and is
for the heavier of the two clusters as
can be seen in Fig. 4.
![]() |
Figure 9:
Gain factor of Type IIn+p+L SNe detected
vs.
![]() ![]() ![]() |
Open with DEXTER |
Figure 10
shows a gain in
the number count discernible by comparing the areas under the
graphs with and without the lens. This is seen quantitatively in
Fig. 11 as a factor 1.25 at
for both lens redshifts. Figures 10
and 11
also makes it quite clear that
the difference between the NFW and SIS models is small. Again we
see the dip at
in Fig. 11. To continue
the comparison; we see a quite large increase in the gain factor at
(max. at 23.5 to be precise) for the
cluster,
and therefore we display in Fig. 12 the distribution of
these Type IIs in redshift space for
showing a clear
and large displacement of the distribution towards higher redshifts
compared to the case with no lens. This would seem particularly interesting
as J=24 is clearly within reach for 8-m class ground based telescopes.
However, Fig. 13 reveals that the absolute number of SNe
for this limiting magnitude is quite low. One important thing to note
though is that the absolute numbers are uncertain since they depend
strongly on the supernova rates which in turn are quite uncertain, especially
for sources at higher redshifts.
![]() |
Figure 10:
Number of Type IIn+p+L SNe detected per year and
![]() |
Open with DEXTER |
![]() |
Figure 11:
Gain factor of Type IIn+p+L SNe detected
vs.
![]() ![]() |
Open with DEXTER |
![]() |
Figure 12:
Number of Type IIn+p+L SNe detected per year and
![]() |
Open with DEXTER |
![]() |
Figure 13:
Number of Type IIn+p+L SNe detected per year vs.
![]() ![]() |
Open with DEXTER |
If we now study similar plots for Type Ia SNe we see a huge increase
in the gain factor at small
22-23 in Fig. 14 for the
cluster. Unfortunately this is not so
exciting when comparing with Fig. 15 where
we again see that the absolute numbers are quite low. However, in the
region
the situation looks a
bit more interesting for the nearer cluster.
![]() |
Figure 14:
Gain factor of Type Ia SNe detected vs.
![]() ![]() |
Open with DEXTER |
![]() |
Figure 15:
Number of Type Ia SNe detected per year vs.
![]() ![]() |
Open with DEXTER |
Finally, we briefly compare the results shown so far with the expectations from an I-band SN search, the currently favoured search filter for ground based high-z supernova searches. Figure 16 shows the expected benefit of heavy clusters as GTs at two redshifts for I- and J-band. The gain in the I-band is generally smaller in the interesting region where the gain is greater than one except maybe for a high redshift cluster. However, clusters this heavy at such high redshifts are very rare (Rosati et al. 2002). Again, the bumpy behaviour of the z=1.0 J-band graph is due to the effect described in Sect. 4.1 for Type IIL SNe.
![]() |
Figure 16:
Gain factor of Type IIn+p+L SNe detected
vs.
![]() ![]() |
Open with DEXTER |
As stated above we have treated all SNe as standard candles. This
perhaps crude assumption when applied to Type II SNe might appear plausible
because one can assume approximately equally many SNe brighter than average
as fainter for each SN type. This could then be assumed to average
out. However, the effect of adding a dispersion is not trivial
since it will affect the lensed and the unlensed case differently. The
magnification will boost the faint end tail,
thereby adding from this tail supernovae in area 2 of
Fig. 1 to the number count. The effect on the unlensed
case will be to add some SNe beyond and remove some SNe
prior to
but with a larger volume for the ones beyond.
The effect will also depend on the differential SN rate, especially
around
.
A more careful investigation should of course
include the SN brightness dispersion.
A more robust justification is found by studying the
magnification as a function of the impact parameter, see
Fig. 17. For a
NFW cluster
at
and a
source at z=2.0 it is seen that within 1
of the cluster
core, the magnification exceeds one magnitude, approximately the intrinsic
dispersion of the Type II SNe considered, see Table 1. Thus,
within this radius the
magnification dominates the magnitude dispersion for this
lens-source configuration. The situation does not change much for
sources in the most interesting regions where the SN rate is
high. However for lighter clusters we see that the approximation
becomes worse.
Similarly, extinction losses have been neglected. Thus, the net benefit
of gravitational lensing is somewhat underestimated, as GTs would
certainly enhance the detectability of extincted supernovae.
![]() |
Figure 17: Magnification as a function of the impact parameter in the lens plane for a source at z=2. The lens redshift is 0.2. |
Open with DEXTER |
In Fig. 17 we can also identify the regions where secondary and tertiary images appear. Tertiary images show up in the innermost regions up to the first big dip. Secondary images are located between the two dips and outside the second dip (the locus of the Einstein radius) there are only primary images. Although not obvious from Fig. 17, in all our considered models primary images are much more abundant for larger fields-of-view and roughly make up about 95% of the images. However, the fraction of secondary and tertiary images is larger for lower redshift clusters (and also for SIS clusters which are more efficient image-splitters) since we then look at the more central parts of the cluster.
We have only studied spherically symmetric clusters which of course is an oversimplification since clusters undoubtedly have substructure in the form of cluster member galaxies. To see in more detail how a specific cluster can affect the searches, a better way is to follow the example in Saini et al. (2000) and use a more sophisticated lens model.
The adopted cluster masses in this work,
,
are consistent with virial and lensing mass estimates of X-ray luminous clusters at
(see Irgens et al. 2002, , and references therein).
Finally, since we have left out Type Ibc (and 87a-like) SNe we have underestimated the total number of SNe that would be detectable.
Observations pointing behind massive heavy galaxy clusters
at redshifts
z=0.2-0.5 may be a useful way
to enhance the detectability of very high-z supernovae.
For the considered supernova rates, the most interesting
effects show up for core-collapse SNe in
searches at limiting magnitudes
mag where
the detection rate
could be significantly enhanced and where the absolute number of
detectable events is considerable. However, the SN rates at high-zand thus also the
absolute numbers are quite uncertain. The relative gain, which is less
sensitive to the SN rates, peaks at lower
limiting magnitudes of
24 mag, where the number of detections could
increase by up to a factor 3 and by a much
larger factor at the highest redshifts. For programs
such as the GOODS/ACS transient survey, the discovery of supernovae
beyond
may be significantly increased with the aid of
GTs. This technique may be used to probe the cosmic star formation
rate at very high redshifts where there is significant controversy
(see e.g. Lanzetta et al. 2002) although this requires good knowledge of
the cluster
properties such as mass and density profile.
When interpreting the measured rates of distant SNe as well as their apparent
brightness, special care must be taken so that magnification effects
from GTs are
properly accounted for. The potential bias increases with the central
wavelength
of the filter used for the supernova search, due to the enhancement of
acceptance for very redshifted sources.
For extremely deep
observations
mag or for very bright SNe (e.g. Type Ia) the
competing effect of field reduction dominates and fewer supernovae
are likely to be discovered behind foreground clusters.
However, the most important effect is that
the detection rate at high redshift could be significantly enhanced
when there is an intervening cluster acting as a gravitational telescope.
Acknowledgements
The authors would like to thank Anne Green, Edvard Mörtsell and Joakim Edsjö for helpful comments and suggestions, Håkon Dahle for cluster information and Julio Navarro for providing the NFW code on his homepage. CG would like to thank the Swedish Research Council for financial support. AG is a Royal Swedish Academy Research Fellow supported by a grant from the Knut and Alice Wallenberg Foundation.