We show how the distances to binary black holes measured in gravitational wave observations with ground-based interferometers can be used to constrain the redshift-distance relation and, thereby, measure the Hubble constant ($H_0$). Gravitational wave observations of stellar-mass binary black holes are not expected to be accompanied by any electro-magnetic event that may help in accessing their redshifts. We address this deficiency by using an optical catalog to get the distribution of galaxies in redshift. Assuming that the clustering of the binaries is correlated with that of the galaxies, we propose using that correlation to measure $H_0$. We show that employing this method on simulated data obtained for second-generation networks comprising at least three detectors, e.g., advanced LIGO - advanced VIRGO network, one can measure $H_0$ with an accuracy of $\sim$8% with detection of a reference population of 25 binaries, each with black holes of mass 10$M_\odot$. As expected, with third-generation detectors like the Einstein telescope (ET), which will measure distances much more accurately and to greater depths, one can obtain better estimates for $H_0$. Specifically, we show that with 25 observations, ET can constrain $H_0$ to an accuracy of $\sim$7%. This method can also be used to estimate other cosmological parameters like the matter density $\Omega_m$ and the dark energy equation of state.

We correct for the use of electron densities instead of total gas density in the pressure fits which were used to derive the local free-fall times (t$_{\rm ff}$) in the original paper.

Low surface brightness (LSB) galaxies are dominated by dark matter halo from the innermost radii; hence they are ideal candidates to investigate the influence of dark matter on different dynamical aspects of spiral galaxies. Here, we study the effect of dark matter halo on grand-design, m = 2, spiral modes in a galactic disk, treated as a collisionless system, by carrying out a global modal analysis within the WKB approximation. First, we study a superthin, LSB galaxy UGC 7321 and show that it does not support discrete global spiral modes when modeled as a disk-alone system or as a disk plus dark matter system. Even a moderate increase in the stellar central surface density does not yield any global spiral modes. This naturally explains the observed lack of strong large-scale spiral structure in LSBs. An earlier work (Ghosh, Saini & Jog, 2016) where the galactic disk was treated as a fluid system for simplicity had shown that the dominant halo could not arrest global modes. We found that this difference arises due to the different dispersion relation used in the two cases and which plays a crucial role in the search for global spiral modes. Thus the correct treatment of stars as a collisionless system as done here results in the suppression of global spiral modes, in agreement with the observations. We performed a similar modal analysis for the Galaxy, and found that the dark matter halo has a negligible effect on large-scale spiral structure.

We present a modal analysis of instabilities of counter-rotating, self-gravitating collisionless stellar discs, using the recently introduced modified WKB formulation of spiral density waves for collisionless systems (Gulati \& Saini). The discs are assumed to be axisymmetric and in coplanar orbits around a massive object at the common center of the discs. The mass in both discs is assumed to be much smaller than the mass of the central object. For each disc, the disc particles are assumed to be in near circular orbits. The two discs are coupled to each other gravitationally. The perturbed dynamics of the discs evolves on the order of the precession time scale of the discs, which is much longer than the Keplerian time scale. We present results for the azimuthal wave number $m=1$ and $m=2$, for the full range of disc mass ratio between the prograde and retrograde discs. The eigenspectra are in general complex, therefore all eigenmodes are unstable. Eigenfunctions are radially more compact for $m = 1$ as compared to $m = 2$. Pattern speed of eigenmodes is always prograde with respect to the more massive disc. The growth rate of unstable modes increases with increasing mass fraction in the retrograde disc, and decreases with $m$; therefore $m=1$ instability is likely to play the dominant role in the dynamics of such systems.

The short--wave asymptotics (WKB) of spiral density waves in self-gravitating stellar discs is well suited for the study of the dynamics of tightly--wound wavepackets. But the textbook WKB theory is not well adapted to the study of the linear eigenmodes in a collisionless self-gravitating disc because of the transcendental nature of the dispersion relation. We present a modified WKB of spiral density waves, for collisionless discs in the epicyclic limit, in which the perturbed gravitational potential is related to the perturbed surface density by the Poisson integral in Kalnaj's logarithmic spiral form. An integral equation is obtained for the surface density perturbation, which is seen to also reduce to the standard WKB dispersion relation. We specialize to a low mass (or Keplerian) self-gravitating disc around a massive black hole, and derive an integral equation governing the eigenspectra and eigenfunctions of slow precessional modes. For a prograde disc, the integral kernel turns out be real and symmetric, implying that all slow modes are stable. We apply the slow mode integral equation to two unperturbed disc profiles, the Jalali--Tremaine annular discs, and the Kuzmin disc. We determine eigenvalues and eigenfunctions for both $m = 1$ and $m = 2$ slow modes for these profiles and discuss their properties. Our results compare well with those of Jalali--Tremaine.

Low surface brightness (LSB) galaxies form a major class of galaxies, and are characterized by low disc surface density and low star formation rate. These are known to be dominated by dark matter halo from the innermost regions. Here we study the role of dark matter halo on the grand-design, $m=2$, spiral modes in a galactic disc by carrying out a global mode analysis in the WKB approximation. The Bohr-Sommerfeld quantization rule is used to determine how many discrete global spiral modes are permitted. First a typical superthin LSB galaxy, UGC 7321 is studied by taking only the galactic disc, modelled as fluid; and then the disc embedded in a dark matter halo. We find that both cases permit the existence of global spiral modes. This is in contrast to earlier results where the inclusion of dark matter halo was shown to nearly fully suppress local, swing-amplified spiral features. Although technically global modes are permitted in the fluid model as shown here, we argue that due to lack of tidal interactions, these are not triggered in LSB galaxies. For comparison, we carried out a similar analysis for the Galaxy, for which the dark matter halo does not dominate in the inner regions. We show that here too the dark matter halo has little effect, hence the disc embedded in a halo is also able to support global modes. The derived pattern speed of the global mode agrees fairly well with the observed value for the Galaxy.

We have examined the effect of slow growth of a central black hole on spherical galaxies that obey S\'ersic or $R^{1/m}$ surface-brightness profiles. During such growth the actions of each stellar orbit are conserved, which allows us to compute the final distribution function if we assume that the initial distribution function is isotropic. We find that black-hole growth leads to a central cusp or ``excess light', in which the surface brightness varies with radius as $R^{-1.3}$ (with a weak dependence on S\'ersic index $m$), the line-of-sight velocity dispersion varies as $R^{-1/2}$, and the velocity anisotropy is $\beta\simeq -0.24$ to $-0.28$ depending on $m$. The excess stellar mass in the cusp scales approximately linearly with the black-hole mass, and is typically 0.5--0.85 times the black-hole mass. This process may strongly influence the structure of nuclear star clusters if they contain black holes.

Future space-based gravity wave experiments such as the Big Bang Observatory (BBO), with their excellent projected, one sigma angular resolution, will measure the luminosity distance to a large number of gravity wave (GW) sources to high precision, and the redshift of the single galaxies in the narrow solid angles towards the sources will provide the redshifts of the gravity wave sources. One sigma BBO beams contain the actual source only in 68 per cent cases; the beams that do not contain the source may contain a spurious single galaxy, leading to misidentification. To increase the probability of the source falling within the beam, larger beams have to be considered, decreasing the chances of finding single galaxies in the beams. Saini, Sethi and Sahni (2010) argued, largely analytically, that identifying even a small number of GW source galaxies furnishes a rough distance-redshift relation, which could be used to further resolve sources that have multiple objects in the angular beam. In this work we further develop this idea by introducing a self-calibrating iterative scheme which works in conjunction with Monte-Carlo simulations to determine the luminosity distance to GW sources with progressively greater accuracy. This iterative scheme allows one to determine the equation of state of dark energy to within an accuracy of a few percent for a gravity wave experiment possessing a beam width an order of magnitude larger than BBO (and therefore having a far poorer angular resolution). This is achieved with no prior information about the nature of dark energy from other data sets such as SN Ia, BAO, CMB etc.

We study the linear $m=1$ counter-rotating instability in a two-component, nearly Keplerian disc. Our goal is to understand these \emph{slow} modes in discs orbiting massive black holes in galactic nuclei. They are of interest not only because they are of large spatial scale--and can hence dominate observations--but also because they can be growing modes that are readily excited by accretion events. Self-gravity being nonlocal, the eigenvalue problem results in a pair of coupled integral equations, which we derive for a two-component softened gravity disc. We solve this integral eigenvalue problem numerically for various values of mass fraction in the counter-rotating component. The eigenvalues are in general complex, being real only in the absence of the counter-rotating component, or imaginary when both components have identical surface density profiles. Our main results are as follows: (i) the pattern speed appears to be non negative, with the growth (or damping) rate being larger for larger values of the pattern speed; (ii) for a given value of the pattern speed, the growth (or damping) rate increases as the mass in the counter-rotating component increases; (iii) the number of nodes of the eigenfunctions decreases with increasing pattern speed and growth rate. Observations of lopsided brightness distributions would then be dominated by modes with the least number of nodes, which also possess the largest pattern speeds and growth rates.

Two new statistics, namely $\Delta_\chi^2$ and $\Delta_\chi$, based on extreme value theory, were derived in \cite{gupta08,gupta10}. We use these statistics to study direction dependence in the HST key project data which provides the most precise measurement of the Hubble constant. We also study the non-Gaussianity in this data set using these statistics. Our results for $\Delta_\chi^2$ show that the significance of direction dependent systematics is restricted to well below one $\sigma$ confidence limit, however, presence of non-Gaussian features is subtle. On the other hand $\Delta_\chi$ statistic, which is more sensitive to direction dependence, shows direction dependence systematics to be at slightly higher confidence level, and the presence of non-Gaussian features at a level similar to the $\Delta_\chi^2$ statistic.

We study the accretion of modified Chaplygin gas upon different types of black hole. Modified Chaplygin gas is one of the best candidates for a combined model of dark matter and dark energy. In addition, from a field theoretical point of view the modified Chaplygin gas model is equivalent to that of a scalar field having a self-interacting potential. We formulate the equations related to both spherical accretion and disc accretion, and respective winds. The corresponding numerical solutions of the flow, particularly of velocity, are presented and are analyzed. We show that the accretion-wind system of modified Chaplygin gas dramatically alters the wind solutions, producing faster winds, upon changes in physical parameters, while accretion solutions qualitatively remain unaffected. This implies that modified Chaplygin gas is more prone to produce outflow which is the natural consequence of the dark energy into the system.

By observing mergers of compact objects, future gravity wave experiments would measure the luminosity distance to a large number of sources to a high precision but not their redshifts. Given the directional sensitivity of an experiment, a fraction of such sources (gold plated -- GP) can be identified optically as single objects in the direction of the source. We show that if an approximate distance-redshift relation is known then it is possible to statistically resolve those sources that have multiple galaxies in the beam. We study the feasibility of using gold plated sources to iteratively resolve the unresolved sources, obtain the self-calibrated best possible distance-redshift relation and provide an analytical expression for the accuracy achievable. We derive lower limit on the total number of sources that is needed to achieve this accuracy through self-calibration. We show that this limit depends exponentially on the beam width and give estimates for various experimental parameters representative of future gravitational wave experiments DECIGO and BBO.

We revise and extend the extreme value statistic, introduced in \cite{gup08}, to study directional dependence in the high redshift supernova data; arising either from departures from the cosmological principle or due to direction dependent statistical systematics in e data. We introduce a likelihood function that analytically marginalises over the Hubble constant, and use it to extend our previous statistic. We also introduce a new statistic that is sensitive to direction dependence arising from living off-centre inside a large void, as well as previously mentioned reasons for anisotropy. We show that for large data sets this statistic has a limiting form that can be computed analytically. We apply our statistics to the gold data sets from \cite{rie04} and \cite{rie07}, as in our previous work. Our revision and extension of previous statistic shows that 1) the effect of marginalsing over Hubble constant instead of using its best fit value has only a marginal effect on our results. However, correction of errors in our previous work reduce the level of non-Gaussianity in the 2004 gold data that was found in our earlier work. The revised results for the 2007 gold data show that the data is consistent with isotropy and Gaussianity. Our second statistic confirms these results.

Thin accretion discs around massive compact objects can support slow pressure modes of oscillations in the linear regime that have azimuthal wavenumber $m=1$. We consider finite, flat discs composed of barotropic fluid for various surface density profiles and demonstrate--through WKB analysis and numerical solution of the eigenvalue problem--that these modes are stable and have spatial scales comparable to the size of the disc. We show that the eigenvalue equation can be mapped to a Schr\"odinger-like equation. Analysis of this equation shows that all eigenmodes have discrete spectra. We find that all the models we have considered support negative frequency eigenmodes; however, the positive eigenfrequency modes are only present in power law discs, albeit for physically uninteresting values of the power law index $\beta$ and barotropic index $\gamma$.

We present the simplest model that permits a largely analytical exploration of the m=1 counter-rotating instability in a "hot" nearly Keplerian disc of collisionless self-gravitating matter. The model consists of a two-component softened gravity disc, whose linear modes are analysed using WKB. The modes are slow in the sense that their (complex) frequency is smaller than the Keplerian orbital frequency by a factor which is of order the ratio of the disc mass to the mass of the central object. Very simple analytical expressions are derived for the precession frequencies and growth rates of local modes; it is shown that a nearly Keplerian disc must be unrealistically hot to avoid an overstability. Global modes are constructed for the case of zero net rotation.

The emission from neutral hydrogen (HI) clouds in the post-reionization era (z < 6), too faint to be individually detected, is present as a diffuse background in all low frequency radio observations below 1420 MHz. The angular and frequency fluctuations of this radiation (~ 1 mK) is an important future probe of the large scale structures in the Universe. We show that such observations are a very effective probe of the background cosmological model and the perturbed Universe. In our study we focus on the possibility of determining the redshift space distortion parameter, coordinate distance and its derivative with redshift. Using reasonable estimates for the observational uncertainties and configurations representative of the ongoing and upcoming radio interferometers, we predict parameter estimation at a precision comparable with supernova Ia observations and galaxy redshift surveys, across a wide range in redshift that is only partially accessed by other probes. Future HI observations of the post-reionization era present a new technique, complementing several existing one, to probe the expansion history and to elucidate the nature of the dark energy.

The most detailed constraints on the accelerating expansion of the universe and details of nature of dark energy are derived from the high redshift supernova data, assuming that the errors in the measurements are Gaussian in nature. There is a possibility that there are direction dependent systematics in the data, either due to uncorrected, known physical processes or because there are tiny departures from the cosmological principle, making the universe slightly anisotropic. To investigate this possibility we introduce a statistic based on the extreme value theory and apply it to the gold data set from Riess et al. (2004). Our analysis indicates a systematic, direction dependent non-gaussianity at about one sigma level.

We investigate the case of a homogeneous tachyon field coupled to gravity in a spatially flat Friedman-Robertson-Walker spacetime. Assuming the field evolution to be exponentially decaying with time we solve the field equations and show that, under certain conditions, the scale factor represents an accelerating universe, following a phase of decelerated expansion. We make use of a model of dark energy (with p=-\rho) and dark matter (p=0) where a single scalar field (tachyon) governs the dynamics of both the dark components. We show that this model fits the current supernova data as well as the canonical \LambdaCDM model. We give the bounds on the parameters allowed by the current data.

We reconstruct the equation of state $w(z)$ of dark energy (DE) using a recently released data set containing 172 type Ia supernovae without assuming the prior $w(z) \geq -1$ (in contrast to previous studies). We find that dark energy evolves rapidly and metamorphoses from dust-like behaviour at high $z$ ($w \simeq 0$ at $z \sim 1$) to a strongly negative equation of state at present ($w \lleq -1$ at $z \simeq 0$). Dark energy metamorphosis appears to be a robust phenomenon which manifests for a large variety of SNe data samples provided one does not invoke the weak energy prior $\rho + p \geq 0$. Invoking this prior considerably weakens the rate of growth of $w(z)$. These results demonstrate that dark energy with an evolving equation of state provides a compelling alternative to a cosmological constant if data are analysed in a prior-free manner and the weak energy condition is not imposed by hand.

In a recent paper (astro-ph/0311364) Alam et al argued that the SNe data of Tonry et al 2003 and Barris et al 2003 appear to favour DE which evolves in time, provided no other priors are invoked. (The effect of invoking priors such as the age of the Universe, the values of $H_0$ and $\Omega_m$ and CMB/LSS observations could modify this conclusion, as demonstrated in astro-ph/0403687 and other recent papers.) The approach adopted by Alam et al to reconstruct the properties of DE was severely (and, as we shall show below -- unfairly) criticized by J\"onsson et al in astro-ph/0404468. In this paper we re-examine the parametrisation used in astro-ph/0311364 and show that, contrary to the claims of J\"onsson et al, the results obtained from this reconstruction are robust and therefore representative of the true nature of dark energy.

The coming few years are likely to witness a dramatic increase in high quality Sn data as current surveys add more high redshift supernovae to their inventory and as newer and deeper supernova experiments become operational. Given the current variety in dark energy models and the expected improvement in observational data, an accurate and versatile diagnostic of dark energy is the need of the hour. This paper examines the Statefinder diagnostic in the light of the proposed SNAP satellite which is expected to observe about 2000 supernovae per year. We show that the Statefinder is versatile enough to differentiate between dark energy models as varied as the cosmological constant on the one hand, and quintessence, the Chaplygin gas and braneworld models, on the other. Using SNAP data, the Statefinder can distinguish a cosmological constant ($w=-1$) from quintessence models with $w \geq -0.9$ and Chaplygin gas models with $\kappa \leq 15$ at the $3\sigma$ level if the value of $\om$ is known exactly. The Statefinder gives reasonable results even when the value of $\om$ is known to only $\sim 20%$ accuracy. In this case, marginalizing over $\om$ and assuming a fiducial LCDM model allows us to rule out quintessence with $w \geq -0.85$ and the Chaplygin gas with $\kappa \leq 7$ (both at $3\sigma$). These constraints can be made even tighter if we use the Statefinders in conjunction with the deceleration parameter. The Statefinder is very sensitive to the total pressure exerted by all forms of matter and radiation in the universe. It can therefore differentiate between dark energy models at moderately high redshifts of $z \lleq 10$.

We apply the Bayesian concept of `evidence' to reveal systematically the nature of dark energy from present and future supernova luminosity distance measurements. We express the unknown dark energy equation of state w(z) as a low order polynomial in redshift and use evidence to find the polynomial order, thereby establishing the minimum order required by the data. We apply this method to the current supernova data, and with a prior -1 < w(z) < 1 and Omega_m = 0.3 \pm 0.05, obtain a large probability of 91% for the cosmological constant model, with the remaining 9% assigned to the two more complex models tested. We also investigate the use of evidence for future supernova data sets such as distances obtainable from surveys like the Supernova Acceleration Probe (SNAP). Given a low uncertainty on the present day matter density we find that, if the underlying dark energy model is only modestly evolving, then a constant $w(z)$ fit is sufficient. However, if the evolution of the dark energy equation of state to linear order is larger than |w_1| ~ 0.5, then the evolution can be established with statistical significance. For models where we can assume the prior -1 < w(z) < 1, the correct polynomial order can be established even for modestly evolving equations of state.

We show that the distance measures (such as the luminosity and angular diameter distances) are linear functionals of the equation-of-state function w(z) of the dark energy to a fair degree of accuracy in the regimes of interest. That is, the distance measures can be expressed as a sum of (i) a constant and (ii) an integral of a weighting function multiplied by the equation of state parameter w(z). The existence of such an accurate linear response approximation has several important implications: (a) Fitting a constant w model to the data drawn from an evolving model has a simple interpretation as a weighted average of w(z). (b) Any polynomial (or other expansion coefficients can also be expressed as weighted sums of the true w(z). (c) A replacement for the commonly used heuristic equation for the effective w(z), as determined by the CMB, can be derived and the result is found to be quite close to the heuristic expression commonly used. (d) The reconstruction of w(z) by Huterer et al. (2002) can be expressed as a matrix inversion. In each case the limitations of the linear response approximation are explored and found to be surprisingly small.

We construct a simple, regularized estimator for the dark energy equation of state by using the recently introduced linear response approximation. We show that even a simple regularization substantially improves the performance of the free form fitting approach. The use of linear response approximation allows an analytic construction of maximum likelihood estimator, in a convenient and easy to use matrix form. We show that in principle, such regularized free form fitting can give us an unbiased estimate of the functional form of the equation of state of dark energy. We show the efficacy of this approach on a simulated SNAP class data, but it is easy to generalize this method to include other cosmological tests. We provide a possible explanation for the sweet spots seen in other reconstruction methods.

We introduce a new cosmological diagnostic pair $\lbrace r,s\rbrace$ called Statefinder. The Statefinder is dimensionless and, like the Hubble and deceleration parameters $H(z)$ and $q(z)$, is constructed from the scale factor of the Universe and its derivatives only. The parameter $r(z)$ forms the next step in the hierarchy of geometrical cosmological parameters used to study the Universe after $H$ and $q$, while the parameter $s(z)$ is a linear combination of $q$ and $r$ chosen in such a way that it does not depend upon the dark energy density $\Omega_X(z)$. The Statefinder pair $\lbrace r,s\rbrace$ is algebraically related to the the dark energy pressure-to-energy ratio $w=p/\epsilon$ and its time derivative, and sheds light on the nature of dark energy/quintessence. Its properties allow to usefully differentiate between different forms of dark energy with constant and variable $w$, including a cosmological constant ($w = -1$). The Statefinder pair can be determined to very good accuracy from a SNAP type experiment.

We develop an algorithm for the reconstruction of the two-dimensional mass distribution of a gravitational lens from the observable distortion of background galaxies. From the measured reduced shear, the lens mapping is obtained, from which a mass distribution is derived. This is unlike other methods where the convergence ("kappa") is directly obtained. We show that this method works best for sub-critical lenses, but can be applied to a critical lens away from the critical lines. For finite fields the usual mass-sheet degeneracy is shown to exist in this method as well. We show that the algorithm reproduces the mass distribution within acceptable limits when applied to simulated noisy data.

We investigate the possibility of detecting HI emission from gravitationally lensed HI clouds (akin to damped Lyman-$\alpha$ clouds) at high redshift by carrying out deep radio observations in the fields of known cluster lenses. Such observations will be possible with present radio telescopes only if the lens substantially magnifies the flux of the HI emission. While at present this holds the only possibility of detecting the HI emission from such clouds, it has the disadvantage of being restricted to clouds that lie very close to the caustics of the lens. We find that observations at a detection threshold of 50 micro Jy at 320 MHz (possible with the GMRT) have a greater than 20% probability of detecting an HI cloud in the field of a cluster, provided the clouds have HI masses in the range 5 X 10^8 M_{\odot} < M_{HI} < 2.5 X 10^{10} M_{\odot}. The probability of detecting a cloud increases if they have larger HI masses, except in the cases where the number of HI clouds in the cluster field becomes very small. The probability of a detection at 610 MHz and 233 MHz is comparable to that at 320 MHz, though a definitive statement is difficult owing to uncertainties in the HI content at the redshifts corresponding to these frequencies. Observations at a detection threshold of 2 micro Jy (possible in the future with the SKA) are expected to detect a few HI clouds in the field of every cluster provided the clouds have HI masses in the range 2 X 10^7 M_{\odot} < M_{HI} < 10^9 M_{\odot}. Even if such observations do not result in the detection of HI clouds, they will be able to put useful constraints on the HI content of the clouds.

We examine the behavior of an anisotropic brane-world in the presence of inflationary scalar fields. We show that, contrary to naive expectations, a large anisotropy does not adversely affect inflation. On the contrary, a large initial anisotropy introduces more damping into the scalar field equation of motion, resulting in greater inflation. The rapid decay of anisotropy in the brane-world significantly increases the class of initial conditions from which the observed universe could have originated. This generalizes a similar result in general relativity. A unique feature of Bianchi I brane-world cosmology appears to be that for scalar fields with a large kinetic term the initial expansion of the universe is quasi-isotropic. The universe grows more anisotropic during an intermediate transient regime until anisotropy finally disappears during inflationary expansion.

Observations of high-redshift supernovae indicate that the universe is accelerating. Here we present a {\em model-independent} method for estimating the form of the potential $V(\phi)$ of the scalar field driving this acceleration, and the associated equation of state $w_\phi$. Our method is based on a versatile analytical form for the luminosity distance $D_L$, optimized to fit observed distances to distant supernovae and differentiated to yield $V(\phi)$ and $w_\phi$. Our results favor $w_\phi\simeq -1$ at the present epoch, steadily increasing with redshift. A cosmological constant is consistent with our results.

Supernovae in distant galaxies that are gravitationally lensed by foreground galaxy clusters make excellent cosmological candles for measuring quantities like the density of the Universe in its various components and the Hubble constant. Distant supernovae will be more easily detectable since foreground cluster lenses would magnify such supernovae by up to 3-4 magnitudes. We show that in the case of the lens cluster Abell 2218, the detectability of high-redshift supernovae is significantly enhanced due to the lensing effects of the cluster. Since lensed supernovae will remain point images even when their host galaxies are stretched into arcs, the signal-to-noise ratio for their observation will be further enhanced, typically by an order of magnitude. We recommend monitoring well-modelled clusters with several known arclets for the detection of cosmologically useful SNe around z=1 and beyond.