starburst galaxies are defined as galaxies that are currently undergoing a period of intense star - formation , mainly within ( but not exclusive to ) the nuclear region , at a rate that can not be maintained over their lifetime @xcite . the radio emission from such galaxies
is not dominated by an active supermassive black hole , but by thermal and non - thermal emission that traces the star - forming activity ( e.g. @xcite ) .
the thermal emission is the result of free - free ( or thermal bremsstrahlung ) radiation from gas ionised by hot , young stars , while the non - thermal emission is from synchrotron radiation due to the acceleration of electrons to relativistic speeds by the core - collapse of massive stars to form type  ii or type  ibc supernovae @xcite .
interaction of the supernova with the local circumstellar medium generates a hot , shocked region that produces a prompt emission typically observed as a radio supernova ( rsn ) @xcite . as the rsn expand into the surrounding interstellar medium , they evolve into supernova remnants ( snrs ) , forming shells that emit non - thermal synchrotron radiation that can be visible over many years @xcite .
studies of snrs within nearby starbust galaxies can provide important information on the astrophysical processes occurring within these galaxies .
obscuration by dust and gas limit observations at shorter wavelengths , while radio wavelengths are largely unaffected .
in addition , high - resolution , wide - field radio interferometry offers an opportunity to detect and resolve individual snrs ( e.g. @xcite ) .
such observations over multiple epochs allow monitoring of the evolution of existing snrs and the detection of new rsne and/or snrs , enabling the supernova and star formation history of the galaxy to be reconstructed .
such studies have been successfully applied to nearby starburst galaxies such as : m82 @xcite ; arp  220 @xcite ; arp  299 @xcite ; ngc  4945 @xcite ; and ngc  253 @xcite .
these studies have provided a wealth of information concerning rsne and the evolution of snrs , the star - formation and supernova rates of the host galaxies , and the interstellar medium of the starburst region .
for example , almost three decades of observations have identified @xmath5100 compact sources in m82 @xcite , where most have been resolved into parsec - scale shell - like structures @xcite .
in addition , these multi - epoch observations have been instrumental in measuring the expansion speeds of the resolved snrs ( 2,000 - 11,000  kms@xmath2 ) , and with giving a direct estimate of the supernova rate ( 0.09  yr@xmath2 ) in m82 . the discovery of new sources such as sn2008iz @xcite and the radio transient 43.78 + 59.3 @xcite would not have been possible without regular radio monitoring of m82 .
in addition , by combining multi - wavelength , high resolution radio observations of 46 compact sources in m82 , @xcite were able to identify low - frequency spectral turnovers , due to free - free absorption by the surrounding ionized medium .
as one of the nearest star - forming galaxies , ngc  253 has been extensively studied from gamma - rays to radio wavelengths .
recently , as part of the advanced camera for surveys ( acs ) nearby galaxy survey treasury ( angst ) , @xcite estimated the distance to ngc  253 as 3.47  @xmath6  0.24  mpc , 3.46  @xmath6  0.07  mpc and 3.40  @xmath6  0.09  mpc .
the weighted average of the @xcite distance estimates , 3.44  @xmath6  0.13  mpc , is used as the distance to ngc  253 in this paper .    the first high - resolution , wide - field radio interferometric observations of ngc  253 were made by ( * ? ? ?
* hereafter th85 ) with the very large array ( vla ) at 15  ghz ( 2  cm ) . with a resolution of @xmath7 , th85 discovered nine compact sources ( designated th1 to th9 ) within the central nuclear region . following the results of th85 , ulvestad and antonucci conducted almost a decade of multi - frequency radio (
1.8 , 5 , 8.3 , 15 and 23  ghz ) observations of the nuclear star - forming region of ngc  253 with the vla , that culminated in a series of papers @xcite ; they identified 64 individual compact sources , including the original nine from th85 .
spectral index measurements were obtained for the 17 brightest sources between frequency pairs of 5/15 and 8.3/23 ghz .
almost half of the 17 sources were identified with thermal hii regions while the remaining sources were taken to be associated with snrs .
over the course of their survey , no new rsne or snrs appeared and the radio flux densities of the detected snrs were found to be stable ( * ? ? ?
* hereafter ua97 ) .
this led to an estimate for the upper limit on the supernova rate of 0.3  yr@xmath2 by ua97 .
in an attempt to resolve the low frequency radio emission in the inner 300  pc , @xcite conducted the first wide - field vlbi observations of ngc  253 .
the observations were carried out with the australian long baseline array ( lba ) at 1.4  ghz and matched the angular resolution of the ua97 23  ghz vla observations . while the observations detected only two sources ( th7 and th9 ) , @xcite showed that the radio emission at low frequencies is absorbed by ionized gas with a free - free optical depth range at 1  ghz of @xmath8  2.5 to @xmath9  8 .
this result was consistent with observations between 0.33 - 1.5 ghz with the vla @xcite and radio - recombination line ( rrl ) modelling of the nuclear region of ngc  253 @xcite .
motivated by these indicators of free - free absorption in ngc  253 , ( * ? ? ?
* hereafter lt06 ) began a program to observe the nuclear region of this galaxy with the lba .
their observations , conducted in 2004 at 2.3  ghz , were higher in sensitivity and resolution than @xcite and covered the region observed by ua97 .
they identified six compact sources , which were also seen in the higher frequency observations with the vla ( ua97 ) .
one of the sources , 5.48 - 43.3 was also resolved into a shell - like structure approximately 90  mas ( 1.7  pc ) in diameter . combining the lt06 data with the multi - wavelength radio data from ua97
, the spectra of 20 compact sources in the nuclear region of ngc  253 were found to be consistent with the free - free absorption interpretation of @xcite .
the results indicated that while the free - free optical depth is highest towards the supposed nucleus , it varies significantly ( @xmath10  1 to @xmath9  20 ) throughout the nuclear region , implying variations in the gas density .
of the 20 sources , eight were found to have flat spectral indices , indicative of thermal h  ii regions , while the remaining sources were taken to be associated with snrs due to their steep spectral indices at high frequencies . with no new sources detected in ngc  253 over almost two decades
, lt06 developed a monte - carlo method based on the work of @xcite , to estimate the upper limit on the supernova rate .
their model took into consideration improved distance measurements , a median free - free opacity , and the sensitivity limits of six observations over a 17  yr period .
a value of 2.4  yr@xmath2 was derived for the upper limit on the supernova rate .
this high value suggests that there may be a large number of undetected rsne , with observations over that period only detecting the rare , bright events .
detecting weaker or short - lived snrs would provide tighter constraints on the supernova rate .
however , this can only be done through frequent , high sensitivity , high resolution observations of ngc  253 .    in this paper
, we present the results of multi - epoch , wide - field vlbi observations conducted at 2.3  ghz with the lba , of the nuclear region of ngc  253 .  [ sec : obsall ] describes the observations , data analysis methods and the sources detected , including cross identifications with previous observations and investigations of possible flux density variations .
free - free absorption modelling of the spectra is described in  [ sec : cs ] .
the morphology of the resolved snr , 5.48 - 43.3 , is examined in detail in  [ sec:548 ] .
an improvement to the monte - carlo method of lt06 is presented with new estimates of the upper limits on supernova and star - formation rates within the inner 300  pc region of ngc  253 in  [ sec : sr ] @xmath11 [ sec : sfr ] .
the results are summarised in  [ sec : sum ] .
[ cols="<,^,^,^",options="header " , ]     ngc  253 was observed at 2.3  ghz with the australian long baseline array ( lba ) , at three epochs as described in  [ tab : table1 ] .
the observations were carried out with : the 64 m parkes ( pa ) antenna of the australia telescope national facility ( atnf ) ; the atnf australia telescope compact array ( atca ) 22 m antennas were used for the 2006 @xmath11 2008 epochs , while 4 @xmath12 22 m antennas were used for the 2007 epoch .
] , used as a phased array ; the atnf mopra ( mp ) 22 m antenna ; the university of tasmania s 26 m antenna near hobart ( ho ) ; and the university of tasmania s 30 m antenna near ceduna ( cd ) .
in addition , the 70 m nasa deep space network antenna at tidbinbilla ( tid ) was used for the 2008 epoch .
the data for each observation were recorded from 4 @xmath12 16 mhz bands ( digitally filtered 2-bit samples ) in the frequency ranges given in  [ tab : table1 ] .
all bands had dual circular polarization at parkes , atca , and mopra , with right circular polarization only at the remaining antennas .
[ tab : table1 ] lists the parameters associated with the lba observations . during each observation , 3  minute scans of ngc  253 ( centred on : @xmath13 = 00@xmath1447@xmath15 33.178@xmath16 ; @xmath17 = @xmath025@xmath1817@xmath19 17@xmath20.060 [ j2000.0 ] ) were scheduled , alternating with 3 minute scans of a nearby phase reference calibration source , pks j0038 - 2459 ( @xmath13 = 00@xmath14 38@xmath15 14.735@xmath16 ; @xmath17 = @xmath024@xmath18 59@xmath19 02@xmath20.235 [ j2000.0 ] ) , located 2.13 degrees from the target . the recorded data for all epochs were correlated using the difx software correlator @xcite , with an integration time of 2 seconds and 64 frequency channels across each 16 mhz band ( channel widths of 0.25 mhz ) .
the @xmath21 coverage for the 2007 and 2008 epochs are shown in  [ fig : uvcov ] .
the initial data reduction and calibration were performed using the data reduction package aips .
prior to calibration , flagging of data during times at which each of the antennas were known to be slewing and time ranges that contained known bad data , were carried out via application of flag files and information provided by the observing logs .
data from the first 30 seconds of each scan from baselines involving atca or parkes were flagged , in order to eliminate known corruption of the data at the start of each scan at these two telescopes .    during correlation , nominal ( constant ) system temperatures ( in jansky ) for each antenna
were applied to the correlation coefficients .
the nominal calibration was refined by application of antenna system temperatures ( in kelvin ) measured during the observation , along with the gain ( in janskys per kelvin ) for each antenna .
further refinements to the amplitude calibration were derived from simultaneously recorded atca observations of pks j0038 - 2459 , which is a bright , compact radio source , unresolved at this frequency on the lba baselines .
since the lba is nt resolving out any extended structure , the atca and lba baselines will measure the same flux density .
thus , the flux density measured at the atca can be used to check and refine the amplitude calibration for the lba data .
suitable simultaneous recorded atca data were obtained for all but the 2008 observation .
the nrao s very long baseline array ( vlba ) calibrator survey routinely observes pks j0038@xmath02459 at 2.3  ghz , but did not in 2008 .
[ fig : j0038lightcurve ] plots the 2.3  ghz light curve of pks j0038@xmath02459 for eight epochs from 1997 to 2010 .
the error bars are @xmath610@xmath22 of the measured flux densities .
a weighted non - linear least - squares fit was used to interpolate to the flux density of pks j0038@xmath02459 at the time of the 2008 lba observation .
the fit gave a flux density of 337  mjy , with a 1@xmath23 error of 19  mjy at the time of the 2008 observation .
this flux density was used to refine the amplitude calibration for the 2008 lba data , by defining a model with this flux density in difmap @xcite .
the difmap task gscale was used to determine corrections to the uncalibrated amplitudes of pks j0038@xmath02459 from the model .
the square point plots the final imaged flux density following amplitude refinement for the 2008 epoch ( square ) .
the arrows indicate the dates of the lba observations . ]    following amplitude calibration refinement , global fringe - fitting solutions were determined for pks j0038@xmath02459 ( aips task fring ) with a three minute solution interval , finding independent solutions for each of the 16 mhz bands .
the delay and phase solutions were examined and , following editing of bad solutions , applied to pks j0038@xmath02459 .
the pks j0038@xmath02459 data were exported to difmap , where the data were vector - averaged over 30  s , flagged of bad data and imaged using standard imaging techniques ( deconvolution and self - calibration of both phase and amplitude ) .
the resulting images of pks j0038 - 2459 for the 2006 , 2007 and 2008 epochs show a highly compact source , with no significant structure on these baselines at this frequency to a dynamic range of 400 . the final calibration solutions ( phase and amplitude ) of pks j0038@xmath02459 were exported via the difmap task , cordump @xcite to an aips - compatible solutions table .
the solutions table was then transferred to aips and applied to all sources in the dataset .
the pks j0038 - 2459 data were also used to derive a bandpass calibration via the aips task bpass which was applied to the ngc  253 data .
the edge channels of each band were flagged from the data set ( 2 channels from both the lower and upper edge of each 16 channel band ) .
the final calibration solutions were applied to both pks j0038@xmath02459 and ngc  253 and the visibility data exported as fits files .
the loss of amplitude due to bandwidth smearing at the field s edge ( radius of @xmath24 from the phase centre ) is @xmath25 on the longest baseline ( hobart - ceduna ) .
thus , to allow imaging of the inner 300  pc at a resolution of @xmath26  milli - arcseconds ( mas ) , the ngc  253 dataset was not averaged in time or frequency . to facilitate imaging in difmap ,
the frequency channels were converted into intermediate frequencies ( ifs ) .
this conversion allowed difmap to treat the frequency channels independently in the @xmath21-plane rather than averaging them together , thus avoiding any further bandwidth smearing effects during the imaging process .
ngc  253 imaging was initially performed with a reduced resolution by excluding data from the hobart and ceduna antennas ( the longest baselines ) .
the imaging parameters were chosen to closely match the imaging parameters of lt06 , with cellsize of 11  mas and application of natural weighting to minimise noise at the expense of resolution .
[ fig : wf ] displays the resulting contour images of ngc  253 for the three epochs .
the 1@xmath23 noise measurement for the three images are listed in  [ tab : table1 ] .
the 1@xmath23 noise levels at the 2006 and 2007 epochs are 2 - 3 times higher than at the 2008 epoch , which can be attributed to the absence of the 70  m tidbinbilla antenna in 2006/07 . however , the 1@xmath23 noise of the 2008 image is 28@xmath22 lower than that of the lt06 ( 2004 epoch ) image due to increased bandwidth .
[ fig:5.48 ] shows the images of the snr 5.48 - 43.3 for the three epochs .
the images were made with the full array ( i.e. including data from both hobart and ceduna ) , with natural weighting and a cell size of 4  mas .      * 15c & & & & & & & & & & & & & + @xmath13 ( j2000.0 ) & @xmath17 ( j2000.0 ) & & & & & + ( 00@xmath1447@xmath27 & ( @xmath02517 ) & ua97 & th85 & lt06@xmath28 & & @xmath29 & @xmath30 & @xmath31 & @xmath32 & & @xmath33 & @xmath34 & @xmath35 & @xmath36 + [ 1pt ] 32.876 & 21.383 & 5.48@xmath043.3 & th9 & y & & 8.6 @xmath6 0.9 & 13.5@xmath6 1.6 & 16.0 @xmath6 1.9 & 7.8 @xmath6 01.2 & & 32.0 @xmath6 3.2 & 22.8 @xmath6 2.4 & 30.2 @xmath6 3.4 & 25.1 @xmath6 9.9 + [ 1pt ] 33.012 & 19.425 & 5.62@xmath041.3 & th7 & y & & 1.7 @xmath6 0.2 & 2.7 @xmath6 0.6 & 2.4@xmath6 0.5 & 1.3 @xmath6 0.3 & & 8.8 @xmath6 0.9 & 6.3@xmath6 0.8 & 5.7 @xmath6 0.8 & 6.8 @xmath6 2.7 + [ 1pt ] 33.115 & 18.039 & 5.72@xmath040.1 & th6 & y & & 1.5 @xmath6 0.2 & @xmath32.0 & @xmath32.2 & 1.0 @xmath60.2 & & 3.8 @xmath6 0.4 & @xmath32.0 & @xmath32.2 & 3.0 @xmath6 1.2 + [ 1pt ] 33.185 & 19.084 & 5.78@xmath041.1 & .... & n & & @xmath31.2 & @xmath32.0 & @xmath32.2 & 0.9 @xmath60.2 & & @xmath31.2 & @xmath32.0 & @xmath32.2 & 4.6 @xmath6 1.8 + 33.177 & 17.830 & 5.79@xmath039.7 & th3 & y & & 1.8 @xmath6 0.2 & @xmath32.0 & @xmath32.2 & 1.0 @xmath60.2 & & 6.2 @xmath6 0.6 & @xmath32.0 & @xmath32.2 & 3.2 @xmath6 1.3 + [ 1pt ] 33.181 & 17.148 & 5.79@xmath039.0 & th2 & y & & 2.4 @xmath6 0.2 & 4.1 @xmath6 0.6 & 4.9 @xmath6 0.7 & 2.6 @xmath60.4 & & 5.7 @xmath6 0.6 & 7.8 @xmath6 1.0 & 7.3 @xmath6 0.9 & 5.9 @xmath6 2.3 + [ 1pt ] +   +   +   +   +    a number of detected sources are clearly visible in  [ fig : wf ] .
we used the source extraction software , blobcat @xcite to identify and measure the flux densities of compact sources in the three images .
seven sources were detected above 5@xmath23 by blobcat in the 2008 epoch , while only the three brightest sources ( 5.48 - 43.3 , 5.62 - 41.3 and 5.79 - 39.0 ) were found via blobcat for the less sensitive 2006 and 2007 epochs ( see  [ tab : lbafluxes ] ) .
included in  [ tab : lbafluxes ] are the flux densities for the compact sources detected in the 2004 epoch by lt06 .
only three of these sources were detected in all four 2.3  ghz lba epochs ( i.e. 2004 , 2006 , 2007 and 2008 ) , while three were detected in only the 2004 and 2008 epochs , with one possible additional source detected solely in the 2008 image .
the flux density of the sources detected by lt06 were also estimated via blobcat and found to agree with lt06 s values within measurement errors . for comparison , the published source flux densities of lt06 are included in  [ tab : lbafluxes ] .
while stacking of the 2004 and 2008 images improved image noise by a factor of @xmath37 , the resulting image did not reveal any additional sources .
errors in the flux - density estimates have two different origins .
firstly , there is the overall uncertainty in the absolute flux density scale for the lba , estimated as @xmath38 @xcite .
secondly , there are the errors in determining the flux density by blobcat , which ranges from 5 - 10@xmath22 .
both terms are added in quadrature for the 2006 and 2007 observations . for the 2008 observation
, we determine a third contribution to the error , which is described in  [ sec : fluxvar ] .
table  [ tab : lbafluxes ] lists cross identifications with sources detected at 15  ghz and 23  ghz by ua97 and at 2.3  ghz with the lba by lt06 . the mean position difference between our 2008 and lt06 s 2004 j2000.0 source positions
was found to be 8  mas , with a standard deviation of 5  mas , which is @xmath39 of the beam - width .
the mean position difference between our 2008 j2000.0 source position and the b1950.0 source positions of the 23 ghz sources from ua97 is 87  mas , with a standard deviation of 51  mas .      )
are plotted on the y - axis , while the observation date is listed on the x - axis .
5@xmath23 upper limits are represented by arrows for epochs without detections .
the error - bars are the absolute flux density scaling and measurement errors ( see  [ sec : flux ] ) added in quadrature . ]
lower integrated flux densities are recorded for all sources except for 5.79@xmath039.0 in the 2008 epoch compared to the 2004 epoch ( see  [ tab : lbafluxes ] and  [ fig : lcurve ] ) .
since it was necessary to interpolate the flux density of a variable calibrator to obtain a value for the 2008 epoch ( see  [ sec : calib ] ) it would not be surprising if there were an error on the flux density scale . given the covariance of the flux densities between 2006 and 2008 ( and possibly for the other epochs too ) this would appear to be the case .
adjusting empirically for the difference in flux density between 2006 and 2008 by taking the mean of each , we find that we must increase the absolute flux density calibration of the 2008 epoch by 20@xmath22 .
after this adjustment , no source appears significantly variable with the exception of 5.805 - 38.92 , whose decrease in flux density would still be a 9@xmath23 error ( assuming of course that the distribution of the errors is gaussian ) .    to investigate whether the flux density decrease is genuine we used the aips task uvmod to simulate datasets with the same @xmath21-coverage and noise level as the 2004 and 2008 lba datasets . in uvmod
the parameter factor was set to zero to allow the injection of the fake sources at the same positions as the real ones .
the fake sources were injected with the properties ( i.e. size , position angle , integrated flux density and position ) of the sources detected by lt06 .
the resulting datasets were exported and imaged in difmap following the procedure described in  [ sec : calib ] .
the integrated flux densities of the injected sources were measured with blobcat and plotted in  [ fig : uvmod ] ( left plot ) against the model flux densities .
if the flux density decrease in the observed dataset were due to only source variability , the measured ( or recovered ) flux densities would be equal to the injected flux densities , represented by the dashed lines .
the deviation of the uvmod flux densities from the dashed lines follows a similar trend to the observed flux densities ( right plot of  [ fig : uvmod ] ) .
thus it is possible the flux density decrease results mainly from systematic effects such as deconvolution errors ( e.g. clean bias @xcite ) and high side - lobes .
the effect of clean bias is known to produce a systematic underestimate of flux densities @xcite by redistributing flux density from sources to noise peaks in the image during deconvolution .
the magnitude of the flux density redistributed is generally independent of flux density , and thus the fractional error is largest for weak sources .
this effect may have a stronger effect in the 2008 image ( as compared to the 2004 image ) due to the presence of high side - lobes .
the high side - lobes may have resulted from the combination of poor @xmath21-coverage and an increased sensitivity due to a higher bandwidth .
the rms ] difference between the 2008 recovered flux densities and the model flux densities is 37@xmath22 , which is added in quadrature with the absolute flux density scaling and measurement errors for the 2008 epoch ( see  [ sec : flux ] ) and listed in  [ tab : lbafluxes ] .
the resulting flux density decrease of 5.805@xmath038.92 is then reduced to a significance of 4.5@xmath23 .
while this variability may be due to unaccounted for calibration errors , we can not rule out contributions from intrinsic variability .
l|cccc|cccc source & @xmath40 & @xmath41 & @xmath42 & type@xmath43 & @xmath44 & @xmath45 & @xmath46 & type + 4.81 - 43.6 & 19.92@xmath612.81 & @xmath01.46@xmath60.62 & 9.16@xmath65.16 & s & 24.1@xmath616.4 & @xmath01.53@xmath60.73 & 11.59@xmath66.84 & s + 5.48 - 43.3 & 105.6@xmath60.3 & @xmath00.77@xmath60.01 & 3.30@xmath60.02 & s & 76.4@xmath60.13 & @xmath00.63@xmath60.01 & 2.89@xmath60.01 & s + 5.49 - 42.2 & 2.31@xmath60.53 & @xmath00.33@xmath60.15 & 2.19@xmath61.65 & t & 2.47@xmath60.84 & @xmath00.35@xmath60.24 & 4.12@xmath63.11 & t + 5.54 - 42.2 & 1.54@xmath60.18 & 0.33@xmath60.06 & 2.54@xmath61.66 & t & 1.68@xmath60.30 & 0.3@xmath60.09 & 4.21@xmath62.60 & t + 5.59 - 41.6 & 0.72@xmath60.04 & 0.67@xmath60.03 & 0.29@xmath61.12 & t & 0.77@xmath60.11 & 0.65@xmath60.07 & 2.23@xmath62.30 & t + 5.62 - 41.3 & 16.11@xmath60.04 & @xmath00.31@xmath60.01 & 1.82@xmath60.02 & t & 12.25@xmath60.07 & @xmath00.21@xmath60.01 & 1.49@xmath60.02 & t + 5.65 - 40.7 & 2.71@xmath60.73 & @xmath00.25@xmath60.17 & 3.46@xmath62.28 & t & 2.89@xmath61.02 & @xmath00.27@xmath60.24 & 5.40@xmath63.66 & t + 5.72 - 40.1 & 8.25@xmath60.37 & @xmath00.02@xmath60.02 & 4.12@xmath60.53 & t & 9.03@xmath60.94 & @xmath00.05@xmath60.04 & 5.70@xmath61.23 & t + 5.73 - 39.5 & 2.07@xmath60.47 & 0.09@xmath60.12 & 3.59@xmath62.33 & t & 2.21@xmath60.65 & 0.07@xmath60.17 & 5.52@xmath63.58 & t + 5.75 - 41.8 & 76.97@xmath65.12 & @xmath01.16@xmath60.04 & 17.46@xmath60.64 & s & 86.16@xmath65.80 & @xmath01.20@xmath60.04 & 18.98@xmath60.66 & s + 5.78 - 39.4 & 121.18@xmath62.49 & @xmath00.81@xmath60.01 & 22.67@xmath60.51 & s & 129.41@xmath63.35 & @xmath00.83@xmath60.01 & 24.61@xmath60.67 & s + 5.79 - 39.0 & 150.79@xmath60.59 & @xmath00.46@xmath60.01 & 17.58@xmath60.08 & s & 149.60@xmath62.32 & @xmath00.46@xmath60.01 & 17.45@xmath60.26 & s + 5.79 - 39.7 & 4899.0@xmath62737.6 & @xmath02.58@xmath60.93 & 25.15@xmath610.52 & s & 5244.7@xmath62959.9 & @xmath02.6@xmath61.14 & 30.09@xmath613.9 & s + 5.78 - 41.1 & 75.97@xmath65.53 & @xmath01.38@xmath60.04 & 17.09@xmath60.60 & s & 33.03@xmath64.71 & @xmath01.07@xmath60.08 & 6.20@xmath60.88 & s + 5.805 - 38.92 & 87.74@xmath625.72 & @xmath01.23@xmath60.15 & 8.26@xmath61.68 & s & 97.7@xmath649.61 & @xmath01.27@xmath60.55 & 15.24@xmath67.14 & s + 5.87 - 40.1 & 70.27@xmath636.82 & @xmath01.50@xmath60.51 & 16.00@xmath66.63 & s & 84.0@xmath616.12 & @xmath01.57@xmath60.11 & 18.17@xmath61.33 & s + 5.90 - 37.4 & 2.30@xmath60.35 & 0.35@xmath60.07 & 4.72@xmath62.30 & t & 2.48@xmath60.48 & 0.32@xmath60.09 & 6.20@xmath63.05 & t + 5.95 - 37.7 & 5.08@xmath61.94 & @xmath00.80@xmath60.30 & 4.45@xmath62.67 & s & 6.15@xmath62.89 & @xmath00.87@xmath60.42 & 6.89@xmath64.37 & s + 6.00 - 37.0 & 82.41@xmath68.22 & @xmath01.32@xmath60.05 & 17.28@xmath60.80 & s & 94.1@xmath69.4 & @xmath01.37@xmath60.05 & 18.98@xmath60.82 & s + 6.40 - 37.1 & 25.72@xmath612.40 & @xmath01.10@xmath60.47 & 12.35@xmath65.91 & s & 31.17@xmath615.94 & @xmath01.17@xmath60.55 & 14.78@xmath67.40 & s +   +   +   +
a downturn in the spectra of the compact sources in ngc  253 at low frequencies was observed previously by @xcite and lt06 .
different mechanisms were explored by lt06 to explain the observed effect including ; a simple power law , a power law spectrum with free - free absorption by a screen of ionized gas , and a self - absorbed bremsstrahlung spectrum .
lt06 demonstrated that the spectra of the sources were consistent with a free - free absorbed power - law spectrum , given by the following equation :    @xmath47    where @xmath48    and @xmath13 is the optically thin intrinsic spectral index , @xmath49 is the free - free optical depth at 1  ghz , and @xmath50 is the intrinsic flux density of the source at 1  ghz .
similar spectra have been obtained for compact sources in other nearby starburst galaxies , such as m82 @xcite , arp  220 @xcite , and ngc  4945 @xcite .
we therefore adopt this model and the method described by lt06 .    to model the free - free absorption towards the compact sources within ngc  253 we compile multi - wavelength radio flux density measurements for these sources from the literature ( listed in  [ tab : allfluxes ] ) , including : 1.4  ghz lba data @xcite ; 2.3  ghz lba data ( lt06 and this paper )
; 23  ghz , 15 ghz , 8.3  ghz and 5  ghz vla data ( ua97 ) .
@xcite detected 10 ua97 sources at 23  ghz with the vla ( a configuration ) in 2004 .
they reported higher flux densities for all sources compared to ua97 , at the same frequency .
only the source 5.79@xmath039.0 ( the assumed core th2 ) was similar in both epochs within a 10@xmath22 error .
the authors attributed the higher flux densities in their data to a larger beam ( @xmath51 ) . to keep the spectral modelling consistent with lt06 , and
maintain a comparable beam size with the lba 2.3 ghz images , the flux densities from @xcite were not used .
[ fig : ffspectra ] presents multi - epoch monitoring of the free - free spectra of 20 compact sources in ngc  253 .
the free - free spectra obtained by lt06 are plotted as dashed lines , while the solid lines are the spectra derived using the 2008 2.3  ghz lba data . for sources not detected in the 2008 observation , an upper limit of 0.9  ( 5@xmath23 ) was adopted .
table  [ tab : ffparams ] lists the fitted free parameters @xmath13 , @xmath49 , and @xmath50 for each compact source .
included for reference are the values of the parameters for the 2004 data .
differences in the free - free spectra and derived parameters between 2004 and 2008 , while potentially large , are still within our errors , suggesting no evidence for changes in the free - free absorbing medium between epochs .
these results further confirm the free - free absorption interpretation of @xcite and lt06 .    from the radio spectra
, lt06 deduced that eight of the 20 sources were thermally dominated h  ii regions ( indicated with a flat intrinsic spectrum , @xmath520.4 ) and the remaining 12 as consistent with steep intrinsic power - law spectra typical of rsne or snrs .
we find no difference in the source classification in our study .
the flux density for the sources with spectra consistent with rsne or snrs , ranges from 2 - 30 mjy at 5  ghz .
the relationship between the rise time and peak luminosity for core - collapse supernovae by @xcite predicts a rise time @xmath53110 days . a typical
rsne follows a power law decay that begins after this fast rise , and would be detectable over many years .
however , no such decay has been observed for these compact sources over two decades , suggesting that they are probably very old snrs .
the brightest source in ngc  253 at 2.3  ghz , 5.48 - 43.3 , is the only fully resolved snr in the starburst galaxy with the full lba array ( i.e. including the long baselines of cd and ho ) . in  [ fig:548v2 ]
we present the two high sensitivity vlbi images of 5.48 - 43.3 ( the 2004 and 2008 epochs ) , restored with a beam size of 13@xmath1215 mas and beam position angle of 58 .
the contours were chosen to represent identical surface brightness levels in both images .    as first noted by lt06 , 5.48@xmath043.3 appears to be a shell - type snr , with a diameter of @xmath570 - 90 mas ( 1.4 - 1.8 pc ) . at both epochs ,
the structure of 5.48@xmath043.3 is dominated by the eastern lobe , which has a higher flux density than the western lobe by a factor of 3 .
this may be the result of interactions with a denser interstellar medium in the direction of the eastern lobe .    comparing the epochs
, we notice that there are several differences in the lobes .
such apparent changes in the small scale structures are possibly due to ambiguities caused by the combination of structural evolution and image fidelity limitations , resulting from incomplete sampling of the @xmath21-plane .
the effects of @xmath21-plane sampling on the appearance and evolution of complex small - scale structures , within spherically symmetric shell - like sources , are well documented by @xcite ; they demonstrated that sparse @xmath21-plane sampling , and the non - uniqueness of deconvolution , can add complex azimuthal structure to a radio brightness distribution that is in reality morphologically simple .
one striking feature of the @xcite results is the departure from a spherically symmetric shell to a two - lobed structure ( similar to sn  1987a @xcite and 5.48@xmath043.3 ) when the @xmath21-plane coverage lacks intermediate baselines , similar to the @xmath21 coverage of both the 2004 ( see  2 of lt06 ) and the 2008 (  [ fig : uvcov ] ) observations .
thus , it is possible that 5.48 - 43.3 possesses a spherically symmetric ( or slightly elliptical ) morphology that is not recovered due to low sensitivity associated with sparse sampling of the @xmath21-plane .
this effect is noticeable in the epochs without the intermediate baselines to tidbinbilla ( 2006 and 2007 ) , with the non - detection of the weaker , western lobe .      with multiple high resolution observations at the same frequency
, it may be possible to determine the expansion speed of 5.48 - 43.3 . in order to carry out such measurements
, we use the aips task iring to measure radial profiles ( averaged in azimuth ) of the source at each epoch .
this method has been successfully implemented to determine the expansion speed of resolved snrs in m82 @xcite .
owing to differences in the small - scale structure between the two epochs ( as discussed in the previous section ) , determination of the geometrical centre of 5.48 - 43.3 is difficult . to account for any positional offset between the epochs ,
the position of the peak surface brightness in the 2008 image was aligned to the position of the 2004 peak surface brightness , using the aips task ogeom .
a common geometrical centre for both images was estimated by visual inspection , and indicated by the crosses in  [ fig:548v2 ] .
using this position as the reference point , the radial profiles were obtained by measuring the integrated flux density within a series of 4  mas thick annuli , and are plotted in  [ fig : iring ] .
no discernible expansion between the two epochs can be seen in  [ fig : iring ] . at the smaller radii
, higher integrated flux densities were recovered for the 2004 image .
this is possibly due to differences in the short spacings of the @xmath21 coverages and/or deeper cleaning in the 2008 image .
if we consider an upper limit of the expansion speed on @xmath54  km  s@xmath2 @xcite for snrs in ngc  253 , an expansion of @xmath52.2  mas , @xmath55 of the 2.3  ghz lba beam , is expected for 5.48@xmath043.3 over 4.3 years between the 2004 and 2008 lba observations . to measure the expansion of 5.48@xmath043.3 , multiple observations separated by @xmath910 years with the current lba baselines , or observations with longer baselines ( @xmath93000  km ) separated by @xmath95 years
would be needed .
a lower limit on the supernova rate can be estimated based on the number of detected snrs , their size , and an assumed expansion rate ( lt06 ) . based upon the method by @xcite and an assumed radial expansion rate of @xmath1  km  s@xmath2 , lt06 estimated the lower limit of the supernova rate to be @xmath56  km  s@xmath57 .
our results from  [ sec:548 ] are consistent with this as the minimum supernova rate .
no new sources have been detected in ngc  253 after two decades of multi - wavelength , high resolution radio observations .
the consequences of the non - detections for estimates of the supernova rate were first investigated by ( * ? ? ? * hereafter ua91 ) .
they assumed a hypothetical population of rsne whose flux densities at 5  ghz peak 100 days after the optical maxima and decay as the @xmath00.7 power of time @xcite .
the rsne population was assumed to have peak luminosities uniformly distributed between 5 and 20 times that of cas  a , and an ngc  253 distance of 2.5  mpc ( th85 ) was also assumed .
ua91 found that @xmath52/3 of the rsne that occurred during the 18 month period between two 5  ghz epochs should be detectable above the second epoch s sensitivity limit . by assuming the events are poisson distributed in time , ua91 determined with 95@xmath22 confidence an upper limit to the supernova rate of 3.0  yr@xmath2 in ngc253 .
this model was used to estimate the supernova rate upper limit for two subsequent 5  ghz epochs , @xmath51.4  yr@xmath2 for a third epoch 2.5 years later @xcite and @xmath50.3  yr@xmath2 for a fourth epoch , 4.0 years later ( ua97 ) .    with new distance measurements for ngc  253 ( 3.94@xmath60.5  mpc , @xcite ) and two additional high resolution observations ,
lt06 developed a new model based on the principles of ua91 and using new data : a 5  ghz vla observation by @xcite ; and a 2.3  ghz lba observation by lt06 .
the rsne at 2.3  ghz were assumed to peak 200 days after the optical maxima at the same flux density as at 5  ghz ( see @xcite ) .
lt06 also considered the effects of free - free absorption by adjusting the flux densities of the rsne with  , assuming a median value of @xmath49 = 6 .
the number of rsne that occur between epochs was determined via a poisson - distributed random number , given a specified supernova rate .
the radio luminosities of the rsne were allowed to evolve using the parametric equation of @xcite into supernovae remnants ( snrs ) and then , given the time and sensitivity of the observation , a test was made to determine if each snr could be detected .
a monte carlo simulation was used to drive the model ( @xmath5 10000 iterations ) , where the proportion of snrs detected , @xmath58 at the end of each epoch was obtained .
the simulation was seeded with an initial supernova rate of 0.1  yr@xmath2 producing a confidence level for the @xmath58 at each epoch . using linear interpolation , the supernova rate required to drive the simulation to a 95@xmath22 confidence limit was determined and the simulation was repeated until the confidence limit was achieved .
the resulting supernova rate was then used to seed the next epoch and the process was continued until all epochs were processed .    using this model
, lt06 found that the increased distance , combined with the effects of free - free absorption , decreased @xmath58 which resulted in an upper limit of the supernova rate at the end of the final epoch of 2.4  yr@xmath2 .
a summary of the model parameters used by @xcite and lt06 to derive upper limits on the supernova rate are listed in  [ tab : snrmodelparams ] .      with additional high resolution radio observations of ngc  253 since lt06 , improved distance measurements to ngc  253 , and availability of a more realistic rsne luminosity function @xcite , it is worth revisiting the estimates of the upper limit on the supernova rate .
the additional ngc  253 observations are the three new 2.3  ghz lba observations in this paper , plus 23  ghz vla observations by @xcite and 15  ghz vla observation by @xcite .
[ fig : timeline ] gives a time - line of all observations of ngc  253 used to determine the supernova rate . to further constrain the upper limit on the supernova rate with these observations ,
we have developed a new model based on the principles of @xcite and lt06 , as well as incorporating suggested improvements by lt06 .
15  ghz ( mga05 ) , vla 23  ghz ( bct09 ) and three lba 2.3  ghz ( this paper ) observations .
note : mga05 - @xcite ; bct09 - @xcite ]    llccccccc name & galaxy & type & distance@xmath43 & detected?@xmath59 & s@xmath60 @xmath28 & l@xmath61 @xmath62 & l@xmath63 @xmath62 & ref .
+ & & & ( mpc ) & & ( mjy ) & & & + sn1983n & m83 & ibc & 4.5 & d & 40.10 & 9.72 & 132.45 & 1 + sn1984l & ngc 4991 & ibc & 21.7 & d & 4.59 & 25.86 & 352.54 & 2 + sn1990b & ngc 4568 & ibc & 18.2 & d & 1.26 & 4.99 & 68.08 & 3 + sn1994i & m51 & ibc & 8.0 & d & 16.20 & 12.50 & 170.38 & 4 + sn1998bw & eso184-g082 & ibc & 37.3 & d & 37.40 & 622.59 & 8487.24 & 5 + sn2000c & ngc 2415 & ibc & 59.0 & e & 0.19 & 7.79 & 106.18 & 6 + sn2001b & ic391 & ibc & 24.0 & d & 2.56 & 17.64 & 240.51 & 6 + sn2001ci & ngc 3079 & ibc & 17.0 & d & 1.06 & 3.67 & 50.06 & 6 + sn2002ap & ngc 628 & ibc & 7.3 & d & 0.45 & 0.29 & 3.94 & 7 + sn2002cj & eso582-g5 & ibc & 106.0 & d & 0.22 & 29.58 & 403.19 & 6 + sn2002dg & anon . &
ibc & 215.0 & e & 0.09 & 50.88 & 693.65 & 6 + sn2003l & ngc 3506 & ibc & 92.0 & d & 2.56 & 258.75 & 3527.31 & 8 + sn2003bg & mcg-05 - 10 - 15 & ibc & 19.6 & d & 52.81 & 242.74 & 3309.07 & 9 + sn2004cc & ngc 4568 & ibc & 18.0 & e & 3.26 & 12.64 & 172.28 & 10 + sn2004dk & ngc 6118 & ibc & 23.0 & d & 1.80 & 11.39 & 155.31 & 10 + sn2004gq & ngc 1832 & ibc & 26.0 & d & 5.75 & 46.51 & 634.00 & 10 + sn2007bg & anon . &
ibc & 152.0 & d & 1.10 & 304.91 & 4156.62 & 11 + sn2007gr & ngc 1058 & ibc & 7.7 & d & 1.00 & 0.71 & 9.67 & 12 + sn2008d & ngc 2770 & ibc & 27.0 & e & 3.00 & 26.17 & 356.72 & 13 + sn2009bb & ngc 3278 & ibc & 40.0 & d & 12.00 & 229.73 & 3131.69 & 14 + sn2010as & ngc 6000 & ibc & 27.3 & d & 1.40 & 12.48 & 170.19 & 15 + ptf11qcj & sdss@xmath64 & ib & 124.0 & d & 5.43 & 998.98 & 13618.22 & 16 + sn1970 g & ngc 5457 & ii & 7.2 & d & 2.50 & 1.55 & 21.14 & 17 + sn1974e & ngc 4038 & ii ?
& 20.9 & d & 1.24 & 6.48 & 88.35 & 18 + sn1978k & ngc 1313 & ii & 4.1 & d & 518.00 & 104.19 & 1420.28 & 18 + sn1979c & ngc 4321 & iil & 22.0 & d & 12.40 & 71.81 & 978.91 & 18 + sn1980k & ngc 6946 & iil & 7.0 & d & 2.45 & 1.44 & 19.58 & 18 + sn1981k & ngc 4258 & ii & 6.6 & d & 5.15 & 2.68 & 36.59 & 18 + sn1982aa & ngc 6052 & ii & 66.0 & d & 19.10 & 995.48 & 13570.58 & 18,19 + sn1985l & ngc 5033 & ii & 12.9 & d & 1.56 & 3.11 & 42.34 & 20 + sn1986e & ngc 4302 & ii & 16.8 & d & 0.33 & 1.11 & 15.19 & 21 + sn1986j & ngc 891 & iin & 12.0 & d & 135.00 & 232.60 & 3170.83 & 18 + sn1988z & mcg+03 - 28 - 022 & iip & 94.5 & d & 1.85 & 197.67 & 2694.71 & 22 + sn1993j & m81 & iib & 3.6 & d & 96.90 & 15.28 & 208.26 & 23 + sn1995n & mcg-02 - 38 - 17 & iin & 24.0 & e & 9.00 & 62.03 & 845.56 & 24 + sn1997eg & ngc 5012 & ii & 40.0 & e & 10.00 & 191.44 & 2609.74 & 25 + sn1998s & ngc 3877 & ii & 17.0 & e & 1.04 & 3.60 & 49.02 & 26 + sn1999em & ngc 1637 & ii & 7.8 & e & 0.30 & 0.22 & 2.98 & 26 + sn2000 ft & ngc 7469 & ii ?
& 70.0 & d & 1.76 & 103.30 & 1408.25 & 27 + sn2001ig & ngc 7424 & iib & 11.5 & d & 21.90 & 34.65 & 472.41 & 28 + sn2001gd & ngc 5033 & iib & 20.0 & d & 7.96 & 38.10 & 519.34 & 29 + sn2004dj & ngc 2403 & iip & 3.6 & e & 1.90 & 0.29 & 4.02 & 30 + sn2004et & ngc 6946 & ii & 5.5 & e & 2.50 & 0.90 & 12.34 & 31 + sn2006jd & ugc 4179 & iib & 76.8 & d & 2.22 & 156.39 & 2131.91 & 32 + sn2008ax & ngc 4490 & iip & 8.1 & d & 4.00 & 3.13 & 42.70 & 33 + sn2008iz & m82 & ii & 3.2 & e & 180.00 & 22.05 & 300.64 & 34 + sn2011cb & pgc 69707 & iib & 29.2 & d & 0.73 & 7.46 & 101.73 & 35 + sn2011dh & m51 & iip & 8.0 & d & 7.80 & 6.02 & 82.04 & 36 + sn2011ei & ngc 6925 & iib & 28.5 & d & 0.56 & 5.48 & 74.72 & 37 + a25 & arp220 & iip & 44.8 & e & 0.30 & 7.20 & 98.21 & 38 + a27 & arp220 & iip ? & 44.8 & e & 0.75 & 18.01 & 245.52 & 38 + & & & & & & & & +   +   +   +   +   +   +   +   +   +   +   +   +   +   +   +   +   +    llccccccc name & galaxy & type & distance & s@xmath60 & l@xmath60 & l@xmath63 & p - value@xmath43 & ref .
+ & & & ( mpc ) & ( mjy ) & & & & + sn1980i & ngc 4374 & ibc & 22.0 & @xmath30.50 & @xmath32.90 & @xmath339.47 & 0.27 & 1,2 + sn1980n & ngc 1316 & ibc & 33.0 & @xmath30.60 & @xmath37.82 & @xmath3106.58 & 0.39 & 1 + sn1981b & ngc 4536 & ibc & 33.0 & @xmath30.30 & @xmath33.91 & @xmath353.29 & 0.30 & 1 + sn1999ex & ic 5179 & ibc & 54.0 & @xmath30.07 & @xmath32.48 & @xmath333.77 & 0.25 & 3 + sn2000cr & ngc 5395 & ibc & 54.0 & @xmath30.04 & @xmath31.47 & @xmath319.98 & 0.20 & 3 + sn2000ew & ngc 3810 & ibc & 15.0 & @xmath30.07 & @xmath30.18 & @xmath32.46 & 0.06 & 3 + sn2000fn & ngc 2526 & ibc & 72.0 & @xmath30.05 & @xmath32.79 & @xmath338.05 & 0.26 & 3 + sn2001 m & ngc 3240 & ibc & 56.0 & @xmath30.03 & @xmath31.05 & @xmath314.32 & 0.17 & 3 + sn2001ai & ngc 5278 & ibc & 118.0 & @xmath30.05 & @xmath38.16 & @xmath3111.29 & 0.39 & 3 + sn2001bb & ic 4319 & ibc & 72.0 & @xmath30.06 & @xmath33.72 & @xmath350.73 & 0.29 & 3 + sn2001ef & ic 381 & ibc & 38.0 & @xmath30.12 & @xmath31.99 & @xmath327.09 & 0.23 & 3 + sn2001ej & ugc 3829 & ibc & 63.0 & @xmath30.04 & @xmath32.09 & @xmath328.48 & 0.23 & 3 + sn2001is & ngc 1961 & ibc & 60.0 & @xmath30.07 & @xmath33.02 & @xmath341.10 & 0.27 & 3 + sn2002j & ngc 3464 & ibc & 58.0 & @xmath30.09 & @xmath33.42 & @xmath346.64 & 0.28 & 3 + sn2002bl & ugc 5499 & ibc & 74.0 & @xmath30.05 & @xmath32.95 & @xmath340.19 & 0.27 & 3 + sn2002bm & mcg013219 & ibc & 85.0 & @xmath30.07 & @xmath35.71 & @xmath377.78 & 0.35 & 3 + sn2002cg & ugc 10415 & ibc & 150.0 & @xmath30.07 & @xmath319.92 & @xmath3271.58 & 0.51 & 3 + sn2002cp & ngc 3074 & ibc & 80.0 & @xmath30.04 & @xmath32.68 & @xmath336.54 & 0.26 & 3 + sn2002dn & ic 5145 & ibc & 115.0 & @xmath30.04 & @xmath36.80 & @xmath392.76 & 0.37 & 3 + sn2002ge & ngc 7400 & ibc & 47.0 & @xmath30.13 & @xmath33.44 & @xmath346.84 & 0.29 & 3 + sn2002gy & ugc 2701 & ibc & 114.0 & @xmath30.04 & @xmath36.06 & @xmath382.67 & 0.35 & 3 + sn2002hf & mcg05320 & ibc & 88.0 & @xmath30.04 & @xmath33.71 & @xmath350.52 & 0.29 & 3 + sn2002hn & ngc 2532 & ibc & 82.0 & @xmath30.04 & @xmath33.54 & @xmath348.26 & 0.29 & 3 + sn2002ho & ngc 4210 & ibc & 42.0 & @xmath30.05 & @xmath31.14 & @xmath315.54 & 0.17 & 3 + sn2002hy & ngc 3464 & ibc & 58.0 & @xmath30.07 & @xmath32.98 & @xmath340.60 & 0.27 & 3 + sn2002hz & ugc 12044 & ibc & 85.0 & @xmath30.03 & @xmath32.51 & @xmath334.18 & 0.25 & 3 + sn2002ji & ngc 3655 & ibc & 23.0 & @xmath30.04 & @xmath30.27 & @xmath33.71 & 0.08 & 3 + sn2002jj & ic 340 & ibc & 66.0 & @xmath30.03 & @xmath31.72 & @xmath323.45 & 0.21 & 2 + sn2002jp & ngc 3313 & ibc & 58.0 & @xmath30.04 & @xmath31.53 & @xmath320.85 & 0.20 & 3 + sn2002jz & ugc 2984 & ibc & 24.0 & @xmath30.03 & @xmath30.17 & @xmath32.35 & 0.06 & 3 + sn2003jd & mcg0159021 & ibc & 81.0 & @xmath30.05 & @xmath34.08 & @xmath355.65 & 0.31 & 4 + sn2011hp & ngc 4219 & ibc & 23.7 & @xmath30.20 & @xmath31.34 & @xmath318.32 & 0.19 & 5 + sn2005cs & m51 & ii & 8.0 & @xmath30.37 & @xmath30.29 & @xmath33.89 & 0.08 & 6 + sn2009mk & pgc 474 & iib & 18.3 & @xmath30.15 & @xmath30.60 & @xmath38.19 & 0.12 & 7 + & & & & & & & &
+   +   +   +   +   +   +      at 5  ghz and 2.3  ghz our hypothetical rsn is the same as that used by ua91 and lt06 ( see previous section and  [ tab : snrmodelparams ] ) . at 8.3  ghz , 15  ghz and 23  ghz
we assume that the rsn luminosity peaks 100 , 70 and 50 days after the optical maximum at the same luminosity as at 5  ghz .
this assumption closely agrees with the light curves of sn  1993j @xcite .    in the models of ua91 and lt06 , the distribution of
the 5  ghz peak luminosities of rsne was assumed to be uniformly distributed between 5 and 20 times the luminosity of cas  a. this assumption follows @xcite , who came to this conclusion after comparing the peak flux densities at 5  ghz of 16 type ii supernovae ( three detections and 13 upper limits ) . since the study by @xcite
there have been many more new detections of type ii as well as type ibc supernovae in the radio .
thus it is worth re - investigating the 5  ghz peak luminosity distribution of these objects .    ) .
the upside down triangles are supernovae whose peak luminosities have not been detected ( 2 type ii and 32 type ibc , see  [ tab : snrupperlims ] ) .
the dashed lines is the average 3@xmath23 detection limit of the vla from @xcite . ]
table  [ tab : snrflux ] lists 51 core - collapse ( 29 type  ii @xmath11 22 type  ibc ) supernovae whose 5  ghz peak flux density have either been detected or estimated via spectral index or extrapolation using the models by @xcite or @xcite .
this table is represented graphically in  [ fig : snedist ] , where the type ibc ( grey symbols ) and type ii ( black symbols ) luminosity ratio to cas  a s at 5  ghz is plotted on the y - axis and their distances on the x - axis .
included in the plots are the 5  ghz peak luminosity upper limits for 34 core - collapse supernovae not detected in the radio (  [ tab : snrupperlims ] ) , that follow the average 3@xmath23 detection limit of the vla ( the dashed lines , from @xcite ) . using the detected core - collapse supernovae listed in  [ tab : snrflux ] , we derived a more accurate luminosity distribution at 5  ghz shown in  [ fig : snehist ] .
the data are divided into bins of size @xmath65log@xmath66(l@xmath63 is the rsne 5  ghz peak luminosity ratio to cas  a s ] ) = 1 .
the dashed line presents the best fit gaussian with mean luminosity ratio of 245.5 , standard deviation of 1.27 and @xmath67 = 0.018 .
this distribution is similar to that presented by @xcite .
however , this analysis ignores the upper limits (  [ tab : snrupperlims ] ) entirely , and so we should consider whether these offer any significant additional information . while their true luminosities are unknown , we need to determine if , given the upper limits , these supernovae are distributed differently to  [ fig : snehist ] . to examine this
we calculate the probability ( or p - value ) , using the cumulative distribution function derived from the gaussian distribution in  [ fig : snehist ] , that a supernova will have a luminosity less than the upper limits listed in  [ tab : snrupperlims ] .
the p - values (  [ tab : snrupperlims ] ) , although on average are lower than 0.5 , are probably consistent with the chosen luminosity function given that the upper limits may not have been observed at a time corresponding to when the peak flux density actually occurs . however , detailed work on testing whether the upper limits in  [ tab : snrupperlims ] are consistent with the chosen model are beyond the scope of this paper .
( l@xmath63 ) . the x - axis mirrors the y - axis of  [ fig : snedist ] and
is given in  [ tab : snrflux ] .
the data are binned to @xmath65log@xmath66(l@xmath63 ) = 1 .
the dashed curve shows a non - linear least - squares fitted gaussian to the underlying data ( black x s ) .
the resulting fit : mean log@xmath66(l@xmath63)= 2.39 , @xmath23 = 1.27 and @xmath67 = 0.018 . ]
@xcite showed that the spectral index between 6  cm and 20  cm evolves with time , for a number of rsne .
this has also been shown to be true for sn  1993j for wavelength pairs of 1.2/2  cm , 2/3.6  cm and 3.6/6  cm @xcite .
thus , it is expected that 2.3/5  ghz , 5/8.3  ghz , 5/15  ghz and 5/23  ghz spectral indices will display similar variations with time . to account for this
, we have included light curves that best describe the spectral indices between 6 @xmath11 20  cm from sn  1980k @xcite and between 6 @xmath11 3.6  cm , 6 @xmath11 2  cm , 6 @xmath11 1.2  cm from sn  1993j @xcite .
[ tab : snrate ] lists the results of applying the improved model parameters described above and listed in column 4 of  [ tab : snrmodelparams ] , to the monte - carlo simulation of lt06 ( see  3.3.2 ) .
we find that the proportion of snrs detected , @xmath58 at the end of each epoch are relatively consistent , with minor variations between 0.66 and 0.84 .
the major changes occur during epochs where the sensitivity varies by almost an order of magnitude .
the consistently high @xmath58 resulted in an upper limit on the supernova rate of @xmath68 = 0.2  yr@xmath2 at the end of the final epoch .
the results of the simulation indicate that the radio observations are detecting 60@xmath22 - 80@xmath22 of all snrs in ngc  253 , thus pointing to the existence of a small population of undetected snrs .
moreover , the fact that no rsne or snrs has been observed in ngc  253 during 21 years of radio observation suggests a low rate of rsne production in ngc  253 .
cccccc    epoch & time@xmath43 & @xmath69 & sensitivity@xmath28 & @xmath58 @xmath62 & @xmath68 @xmath70 + & ( yr ) & ( ghz ) & ( mjy ) & - & ( yr@xmath2 ) + 1 & ... & 5.0 & ... & ... & + 2 & 1.5 & 5.0 & 3.0 & 0.746 & @xmath32.71 + 3 & 2.5 & 5.0 & 3.0 & 0.713 & @xmath31.04 + 4 & 4.0 & 5.0 & 3.0 & 0.683 & @xmath30.53 + 5 & 4.0 & 8.3 & 0.3 & 0.803 & @xmath30.34 + 6 & 0.25 & 15 & 0.6 & 0.838 & @xmath30.33 + 7 & 4.5 & 2.3 & 1.2 & 0.697 & @xmath30.25 + 8 & 0.5 & 23 & 2.0 & 0.661 & @xmath30.24 + 9 & 1.7 & 2.3 & 2.0 & 0.676 & @xmath30.22 + 10 & 0.8 & 2.3 & 2.2 & 0.689 & @xmath30.21 + 11 & 1.0 & 2.3 & 0.9 & 0.791 & @xmath30.20 +   +   +   +   +   +   +   +    the final result for @xmath68 agrees with estimates determined from near infrared ( nir ) observations of the [ feii ] line ( 0.24  yr@xmath2 ) by @xcite .
observations of the [ feii ] line have been shown to be a strong tracer of shocks associated with snrs @xcite .
dust grains in the interstellar medium containing fe atoms are destroyed via the shocks .
the process releases the fe atoms , which are then ionized by the interstellar radiation field . in the post shock region ,
fe@xmath71 is excited by electron collisions , causing it to emit at nir wavelengths .
however , there are possible contributions to the [ feii ] line from shocks due to other processes such as mergers @xcite . thus , it is possible that the nir [ feii ] observations are overestimating the supernova rate upper limit .
the major uncertainty in our model arises from the rsne peak luminosity distribution .
first , the detection of rsne has primarily been through follow - up observations of optical supernovae , resulting in a small and incomplete sample , whose luminosity depends on the sensitivity of the survey .
second , there is the possibility that many rsne are not visible in the optical due to extinction and we may be looking at a slightly different sub - population further uncertainty comes from the assumed rise time of the rsne at the different frequencies . in this model , for simplicity , constant rise times ( the time to reach maximum radio luminosities ) were assumed for the hypothetical rsne at the different frequencies .
however , this has been observed to be variable for different rsne ( e.g. @xcite ) .
additionally , the relation between the rise time and peak luminosity has also been observed to be variable @xcite and does not necessarily follow a one - to - one relationship as assumed in the model presented
.    further refinement to the supernova rate of ngc  253 could be obtained from either a large unbiased sample of the peak flux densities of rsne , or direct observations of the rsne population in ngc  253 .
both possibilities can be provided by the _ _ square kilometre array _ _ , ( ska ) and new generation telescopes , such as the karl g. jansky vla ( jvla ) . the sensitivity of both instruments may provide the opportunity to detect a population of weak snrs that are faint due to either free - free absorption or intrinsic properties .
the ska is expected to have a type ii supernova detection rate of @xmath5620 yr@xmath2  deg@xmath72 out to a redshift @xmath73  5 @xcite .
the jvla and precursors to the ska may provide a detection rate of @xmath5160 yr@xmath2  deg@xmath72 out to redshift @xmath73  3 @xcite .      with new upper limits on the supernova rate
, we can provide a new estimate for the star formation rate ( sfr ) in ngc  253 .
following @xcite , the sfr can be determined directly from the supernova rate ( and vice versa ) . by assuming that all stars with mass @xmath74 eventually form supernovae
, the relation between sfr and the supernova rate is given by :    @xmath75 = 24.4\left[\dfrac{\nu_{sn}}{\mathrm{yr}^{-1}}\right]\ ] ]    using the new supernova rate upper limits of @xmath76  yr@xmath2 (  [ sec : sr ] ) , we find @xmath4  yr@xmath2 .
we can compare this upper limit to estimates from independent methods .
if one assumes that the main contribution to the observed radio luminosity in a starburst are from rsne , snrs and thermal h  ii regions , then the sfr of that region is directly proportional to its radio luminosity at the observed radio wavelength and can be calculated using the relation given by @xcite and @xcite .
@xcite observed a total flux density of 0.59 jy for ngc  253 at 24  ghz . using the relation by @xcite and the data from @xcite we obtain @xmath77 = 9.3@xmath78  yr@xmath2 .
additional observations of ngc  253 at 23  ghz ( total flux density of 0.56  jy ) by @xcite give a similar @xmath77 = 8.7@xmath78  yr@xmath2 for the @xmath5300  pc nuclear region , using the same method .    the @xmath77 can also be determined from far infrared ( fir ) emission using equation 26 from @xcite . from fir luminosity measurements with _ _ iras _ _ @xcite and scaling for the new distance 3.44  mpc , the sfr is 1.3 - 2.0 @xmath78  yr@xmath2 for the inner @xmath5300  pc nuclear region and 2.6 - 3.2 @xmath78  yr@xmath2 for the entire galaxy .
fir luminosity measurements of the inner @xmath5350  pc nuclear region by _ _
_ @xcite lead to an estimate of 2.7@xmath78  yr@xmath2 .
it is encouraging that the @xmath77 calculated using different estimators are in broad agreement with each other .
however , the differences point to the existence of systematic or intrinsic differences between the different estimators .
the supernova rate derived from the fir luminosity ( via  20 from @xcite ) is 0.05 - 0.13  yr@xmath2and is similar to the lower limit of the supernova rate deduced from expansion rates ( lt06 ) .
we have presented the results of multi - epoch observations of the southern starburst galaxy , ngc  253 with the lba at 2.3  ghz . the results presented here are complementary to previous 2.3  ghz lba observations by @xcite .
we find the following results :    1 .
seven compact sources were detected in the highest sensitivity observation ( the 2008 epoch ) .
all sources were identified with higher - frequency vla observations ( ua97 ) , while six were identified with a 2.3  ghz lba detections by lt06 .
the three brightest sources were also detected in the lower sensitivity observations ( 2006 and 2007 ) .
the shell - like snr , 5.48 - 43.3 , was successfully imaged with the highest resolution in all three epochs .
the structure in the 2008 image shows a double - lobed morphology which is dominated by the eastern lobe .
the weaker , western lobe was not detected in the 2006 and 2007 observations , which is attributed to lower sensitivity and fidelity , due the lack of the 70  m tidbinbilla antenna .
3 .   observed differences in the small - scale structure of 5.48@xmath043.3 between the 2004 ( lt06 ) and the 2008 ( this paper ) observations can be explained by deconvolution errors associated with sparse @xmath21-plane sampling .
no discernible expansion has been observed for 5.48@xmath043.3 between the 2004 and 2008 epochs .
the spectra of the 20 compact sources in ngc  253 from lt06 are compared with the spectra obtained with the new 2.3  ghz lba flux density measurements ( 2008 epoch ) . as with lt06
, we find that the spectra fit a free - free absorption model , with little difference in the spectra between the epochs .
our results show no change in the classification of the compact sources by lt06 .
12 of the 20 sources have steep spectra associated with snrs , while the remaining eight have flat intrinsic power - law spectra ( @xmath79 ) , indicative of h  ii regions .
we derive an improved rsne peak luminosity distribution at 5  ghz using data from the literature for 51 core - collapse supernovae ( 29 type  ii @xmath11 22 type  ibc ) . 7 .
we estimate a value for the upper limit of the supernova rate of 0.2  yr@xmath2 for the inner 300  pc of ngc  253 , using an improved model that is based upon the principles of ua91 and lt06 .
this result was found to be consistent with estimates from near infrared ( nir ) observations .
the results of our model indicate that the radio observations are detecting 60@xmath22 - 80@xmath22 of all snrs in ngc  253 , thus pointing to the existence of a small population of undetected snrs .
moreover , no rsne or snrs has been observed in ngc  253 during 21 years of radio observation suggests a low rate of rsne production in ngc  253 .
a new upper limit to the star formation rate of @xmath4  yr@xmath2 is estimated directly from the supernova rate limits for the inner 300  pc region of the galaxy , which is consistent with independent estimates .
the international centre for radio astronomy research is a joint venture between curtin university and the university of western australia , funded by the state government of western australia and the joint venture partners .
the long baseline array is part of the australia telescope which is funded by the commonwealth of australia for operation as a national facility managed by csiro .
sjt is a western australian premier s research fellow , funded by the state government of western australia .
the centre for all - sky astrophysics ( caastro ) is an australian research council centre of excellence , funded by grant ce110001020 .
we wish to thank paola castangia @xmath11 andreas brunthaler for use of their vla data , cormac reynolds for assisting with the calibration of the lba data and cath trott for the helpful discussions on statistics .
this research has made use of material from the bordeaux vlbi image database ( bvid ) and the the nasa / ipac extragalactic database ( ned ) which is operated by the jet propulsion laboratory , california institute of technology , under contract with the national aeronautics and space administration .
we thank the anonymous referee for comments that improved this paper .


today string theory is a promising candidate for the quantum theory of gravity .
since it requires more than four dimensions for its internal consistency , there has been growing interest in higher dimensional solutions and , in particular , in higher dimensional black holes ( see e.g. @xcite ) .
the higher - dimensional generalizations of the stationary axisymmetric kerr black holes were found by myers and perry @xcite , who anticipated already the existence of higher - dimensional black rings . in 2001
emparan and reall @xcite then found such black rings in five dimensions .
they rotate in the direction of the ring and possess a horizon topology of @xmath0 . the phase diagram of these black rings together with the myers - perry black holes showed , that the uniqueness of 4-dimensional vacuum black holes does not generalize to higher dimensions .    in myers - perry black hole
spacetimes the geodesic equations are separable @xcite .
however , the geodesic equations of five dimensional black rings do not seem to be separable , in general . nevertheless , it is possible to separate the equations of motion on the rotational axis ( which is actually a plane ) , in the equatorial plane and in the case @xmath1 ( which is only possible in the ergosphere ) @xcite .
hoskisson @xcite studied the geodesic motion of a singly spinning black ring and discussed the separability of the hamilton - jacobi equation , and provided numerical solutions .
he analyzed numerically the motion on the rotational axis and the equatorial plane in detail .
some aspects of nullgeodesics in the equatorial plane were studied by elvang , emparan and virmani @xcite .
the pomerasky - senkov doubly spinning black ring @xcite was studied by durkee @xcite .
he showed that it is possible to separate the hamilton - jacobi equation of the doubly spinning black ring in the case @xmath1 and analyzed the effective potential on the two axes ( planes ) of rotation and in the case @xmath1 .
the zero energy nullgeodesics of the singly spinning dipole black ring were analyzed by armas @xcite . in @xcite igata ,
ishihara and takamori concentrated on stable bound orbits in the singly spinning black ring spacetime , which they found numerically on and near the rotational axis .
so far the equations of motion for test particles in black ring spacetimes were only solved numerically , and no analytic solutions have been given .
the first to solve the geodesic equations in a black hole spacetime analytically was hagihara @xcite .
he presented the solution of the geodesic equation for test particles in the schwarzschild spacetime in terms of the elliptic weierstra @xmath2 function .    when adding the cosmological constant to the schwarzschild metric one encounters hyperelliptic curves in the geodesic equations .
the equations of motion in schwarzschild-(anti ) de sitter spacetimes in four dimensions were solved in @xcite .
also analytical solutions of the geodesic equations in higher dimensional myers - perry spacetime @xcite as well as in higher dimensional schwarzschild , schwarzschild-(anti ) de sitter , reissner - nordstrm and reissner - nordstrm-(anti ) de sitter @xcite were found .
the mathematical method is based on the jacobi inversion problem .
the solution can be found if the problem is restricted to the theta - divisor , the set of zeros of the theta function .
enolski , pronine and richter developed this method in 2003 to solve the problem of the double pendulum @xcite .
in this paper we present analytical solutions of the equations of motion of a singly spinning black ring . in the case @xmath1 and in the equatorial plane the equations
are of elliptic type , however , on the rotational axis the geodesic equations are of hyperelliptic type .
the singly spinning black ring solution can be written in the form @xcite @xmath3 .\ ] ] the metric is given in toroidal coordinates ( see [ pic : ringcoord ] ) where @xmath4 , @xmath5 and @xmath6 .
@xmath7 and @xmath8 are @xmath9-periodic .
the metric functions are @xmath10     and @xmath8 .
solid circles correspond to @xmath11 and dashed circles corredspond to @xmath12 .
@xcite , width=453 ]    the parameters @xmath13 and @xmath14 describe the shape , mass and angular momentum of the ring .
@xmath13 lies in the range @xmath15 to ensure the black ring is balanced .
a spacelike curvature singularity is located at @xmath16 .
the metric has a coordinate singularity at @xmath17 , so the event horizon lies at @xmath18 .
the ergosphere of the singly spinning black ring is determined by @xmath19 , which is at @xmath20 .
since we have @xmath21 , it is clear that an ergoregion does exist .
the topology of the event horizon and the ergosphere is @xmath22 .
+ the inverse metric is @xmath23 \ , .\ ] ] the singly spinning black ring metric and its hamiltonian @xmath24 do not depend on the coordinates @xmath25 , @xmath7 and @xmath8 , so we have three conserved momenta @xmath26 with the associated killing vector fields @xmath27 , @xmath28 and @xmath29 . a dot denotes the derivative with respect to an affine parameter @xmath30 .
@xmath31 @xmath32 is the energy , @xmath33 and @xmath34 are the angular momenta in @xmath7- and @xmath8-direction .
the conjugate momenta in @xmath35- and @xmath36-direction are : @xmath37 to obtain the equations of motion for a particle in the singly spinning black ring spacetime we need the hamilton - jacobi equation : @xmath38 we already have three constants of motion ( @xmath32 , @xmath33 and @xmath34 ) and the mass shell condition @xmath39 gives us a fourth , so we can make the ansatz @xmath40 inserting this ansatz into ( [ eqn : hamilton - jacobi ] ) gives @xmath41 \ , .
\label{eqn : hjd - sring}\ ] ] the hamilton - jacobi equation does not seem to be separable in general . however , it is possible to separate the equation in the special case @xmath1 .
these zero energy null geodesics are only realisable in the ergoregion .
we can also obtain equations of motion for geodesics on the @xmath7- and @xmath8-axis by setting @xmath42 ( @xmath7-axis ) or @xmath43 ( @xmath8-axis ) .
the plane @xmath42 which is called the @xmath7-axis , is the equatorial plane of the black ring .
the plane @xmath43 which is called the @xmath8-axis corresponds to the rotational axis of the singly spinning black ring .
+ in the next sections we will study these special cases and solve the corresponding equations of motion analytically .
for @xmath1 it is possible to separate the hamilton - jacobi equation : @xmath44 with a separation constant @xmath45 , the equation ( [ eqn : ham - jac - sing ] ) splits into two : @xmath46 so that @xmath47 using @xmath48 and ( [ eqn : t - impuls])-([eqn : y - impuls ] ) the separated hamilton - jacobi equation gives the equations of motion : @xmath49 where we have introduced the mino - time @xcite @xmath50 .      equation ( [ eqn : sing - x - gleichung ] ) and ( [ eqn : sing - y - gleichung ] ) can be written as @xmath51 @xmath52 and @xmath53 can be regarded as effective potentials ( see @xcite ) . to get real solutions for the @xmath35- and @xmath36-equation the effective potentials have to be negative .
so @xmath54 is required , because @xmath55 for @xmath15 and @xmath56 .
the zeros of the effective potentials and hence @xmath57 and @xmath58 mark the turning points of the motion of light or a test particle ( in this case we only have light since @xmath59 ) . a good way to determine the number of zeros are parametric diagrams .
figure [ pic : sing - parameter ] shows a parametric @xmath33-@xmath13-@xmath45-diagram for @xmath52 .
it turns out that @xmath60 has two real zeros in the allowed range of @xmath35 or none .
one zero is always negative while the other can be positive or negative .
if @xmath60 has two zeros the @xmath35-motion takes place between these two values , if @xmath60 has a single zero at @xmath61 the @xmath35-motion stays constant .
@xmath62 and accordingly @xmath53 determine the type of the orbit . if @xmath63 and @xmath64 , @xmath62 has no real zeros .
otherwise @xmath62 has always one real zero in the allowed range of @xmath36 .
that means the only possible orbit is a terminating orbit ( to ) , where light crosses the horizon and falls into the singularity .
see figure [ pic : sing - potential ] for examples of the effective potentials .
-@xmath13-@xmath33-diagram for the singly spinning black ring in the case @xmath1 . inside the structure
@xmath60 has two real zeros , outside the structure it has no real zeros in the allowed range of @xmath35.,width=302 ]      equation ( [ eqn : sing - x - gleichung ] ) can be written as @xmath65 where @xmath57 is a polynomial of third order with the coefficients @xmath66 the substitution @xmath67 transforms the polynomial into the standard weierstra form @xmath68 where @xmath69 equation ( [ eqn : sing - weierstrass - form ] ) is of elliptic type and is solved by the weierstra elliptic function @xcite @xmath70 where @xmath71 and @xmath72
. then the solution of ( [ eqn : sing - x - gleichung ] ) takes the form @xmath73 .\ ] ]      equation ( [ eqn : sing - y - gleichung ] ) can be written as @xmath74 where @xmath58 is a polynomial of third order with the coefficients @xmath75 the problem can be solved analogously to the @xmath35-equation . here the solution is @xmath76 , \ ] ] where @xmath77 and @xmath78 .
using ( [ eqn : sing - x - gleichung ] ) the equation ( [ eqn : sing - phi - gleichung ] ) becomes @xmath79 @xmath80 we substitute @xmath81 to transform @xmath60 into the weierstra form @xmath82 ( see ( [ eqn : sing - weierstrass - form ] ) ) : @xmath83 now @xmath84 has the zeros @xmath85 and @xmath86 .
we next apply a partial fractions decomposition upon equation ( [ eqn : sing - ix ] ) : @xmath87 @xmath88 are constants which arise from the partial fractions decomposition and depend on the parameters of the metric and the test particle
. then we substitute @xmath89 with @xmath90 .
equation ( [ eqn : sing - ix - partial ] ) now simplifies to @xmath91 with @xmath92 and @xmath93 .    after solving the integrals of the third kind ( see e.g. @xcite ) , the final solution reads @xmath94 with @xmath95 .
the index @xmath35 refers to the weierstra - functions with respect to the parameters @xmath96 and @xmath97 .      using ( [ eqn : sing - y - gleichung ] ) equation ( [ eqn : sing - psi - gleichung ] ) becomes @xmath98 @xmath99 the @xmath8-equation
can be solved analogously to the @xmath7-equation . with @xmath100 , @xmath93 and @xmath95
the solution is @xmath101 @xmath102 are constants which arise from the partial fractions decomposition and depend on the parameters of the metric and the test particle .
the index @xmath36 refers to the weierstra - functions with respect to the parameters @xmath103 and @xmath104 .      using ( [ eqn
: sing - y - gleichung ] ) we can write ( [ eqn : sing - t - gleichung ] ) as @xmath105 @xmath106 the @xmath25-equation can be solved in an analogous way to the @xmath7- and @xmath8-equation .
we substitute @xmath107 .
the integral ( [ eqn : sing - em0-tint ] ) is of the third kind and has the poles @xmath108 and @xmath109 ( with respect to @xmath110 ) .
then we apply a partial fractions decomposition , where the constants @xmath111 arise and we substitute again @xmath112 .
after the solution of the occurring elliptic integrals of the third kind , the solution of ( [ eqn : sing - t - gleichung ] ) yields @xmath113      in the ergosphere of a singly spinning black ring only tos with @xmath1 are possible .
figure [ pic : sing - to ] shows a to plotted in the @xmath35-@xmath36-plane in cartesian coordinates ( @xmath114 , @xmath115 ) .    to change from the ring
coordinates to the polar coordinates ( @xmath116 , @xmath117 ) the transformation @xmath118 is used .
then conventional cartesian coordinates take the form @xmath119 ( see @xcite or @xcite ) .
the singularity of the black ring is at @xmath120 , @xmath121 .    , @xmath122 , @xmath123 and @xmath124 @xmath114-@xmath115-plot of a to for the singly spinning black ring in the case @xmath1 .
the black dashed circles are the event horizon and the red dotted circles denote the ergosphere.,width=377 ]
the surface @xmath43 is the axis of rotation of the singly spinning black ring . here
the hamilton - jacobi equation depends on the coordinate @xmath35 only .
we set @xmath43 , @xmath125 and @xmath126 in the hamilton - jacobi equation ( [ eqn : hjd - sring ] ) : @xmath127 \ , .\ ] ] this can be rearranged to @xmath128 - \frac{h(x)}{g^2(x)}\phi^2 : = x_s .\ ] ] then we have @xmath129 now we set the partial derivatives of @xmath130 with respect to the constants @xmath131 , @xmath32 and @xmath33 to zero in order to obtain the equations of motion . with the mino - time @xcite @xmath132 the equations of motion take the form @xmath133 where @xmath60 is a polynomial of fifth order .      from ( [ eqn : sing - psi - x - gleichung ] ) we can read off the effective potential consisting of two parts @xmath134 and @xmath135 ( to be consistent with the effective potential on the equatorial plane later on ) :
@xmath136 since @xmath60 can be written as @xmath137 the effective potential takes the form @xmath138 where @xmath139 and @xmath140 .
figure [ pic : s - psi - orbits1 ] shows the effective potential for the motion on the @xmath8 axis .
@xmath141 is plotted in red ( solid ) while @xmath142 is plotted in blue ( dotted ) .
the grey area between the two parts of the potential is a forbidden zone where no motion is possible because @xmath60 becomes negative there . @xmath141 and @xmath142 are symmetric ( @xmath143 ) and meet at the horizon .
@xmath144 is always a zero of @xmath60 , but since the point @xmath144 , @xmath43 corresponds to infinity in cartesian coordinates it is not a real turning point of the test particle . if @xmath145 then @xmath60 has the zeros @xmath144 and @xmath146 and possibly a third zero between -1 and + 1 .
the coordinate range of @xmath35 ( @xmath147 ) only covers the space from infinity ( -1 ) to the center of the black ring ( + 1 ) . since there is no potential barrier at @xmath146 for @xmath145 , light and test particles with the right amount of energy cross the center of the black ring and continue their orbit at the other side of the black ring .    if @xmath145 then none or one turning point exists . if @xmath148 there is a potential barrier which prevents the geodesics from reaching @xmath146 and one or two turning points exist . for larger @xmath13 and @xmath149
the potential has local extrema which lead to three turning points .
possible orbits are bound orbits ( bo ) , where light or test particles circle the black ring , and escape orbits ( eo ) , where light or test particles approach the black ring , turn around at a certain point and escape the gravitational field .
there are five different types of orbits ( see table [ tab : s - psi - typen - orbits1 ] ) .    *
type a : + @xmath60 has no zero in the range @xmath150 . eos without a turning point exist .
the orbit crosses the equatorial plane ( @xmath146 ) and reaches infinity ( @xmath144 and @xmath43 ) . *
type b : + @xmath60 has one zero in the range @xmath150 .
bos with a turning point on each side of the ring exist , so that the orbit crosses the equatorial plane ( @xmath146 ) . * type c : + @xmath60 has one zero in the range @xmath150 .
eos with a turning point exist . *
type d : + @xmath60 has two zeros in the range @xmath150 .
bos which do not cross the equatorial plane exist . *
type e : + @xmath60 has three zeros in the range @xmath150 .
bos which do not cross the equatorial plane and eos exist .
a & 0 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( -2.5,0)(3,0 )    & eo + b & 1 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( -1.5,0)(3,0 )    & bo + c & 1 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( -2.5,0)(2,0 )    & eo + d & 2 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( -1,0)(2,0 )    & bo + e & 3 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( -0,0)(2,0 ) ( -2.5,0)(-1.0,0 )    & eo , bo +      equation ( [ eqn : sing - psi - x - gleichung ] ) can be written as @xmath151 where @xmath60 is a polynomial of fifth order with the coefficients @xmath1522\lambda e^2}{(1-\lambda)^2 } - 2r^2\lambda^2 \left ( \frac{(1+\lambda^2)}{(1-\lambda)^2}e^2 -m^2 \right ) \nonumber\\   a_{x,3 } & = & 2\phi^2\lambda + \frac{[2r^2\lambda^2-r^2(1+\lambda^2)]2\lambda e^2}{(1-\lambda)^2 } + ( -2r^2\lambda - r^2(1+\lambda^2)\lambda ) \left ( \frac{(1+\lambda^2)}{(1-\lambda)^2}e^2 -m^2
\right ) \nonumber\\   a_{x,2 } & = & -4 \phi^2 \lambda-\phi^2 ( 1+\lambda^2 ) + \frac{[2r^2\lambda^2+r^2(1+\lambda^2)]2\lambda e^2}{(1-\lambda)^2 } + ( 2r^2\lambda - r^2(1+\lambda^2 ) ) \left ( \frac{(1+\lambda^2)}{(1-\lambda)^2}e^2 -m^2 \right ) \nonumber\\   a_{x,1 } & = & -2\phi^2\lambda-2\phi^2(1+\lambda^2 ) + \frac{2r^2(1+\lambda^2)\lambda e^2}{(1-\lambda)^2 } + ( 2r^2\lambda + r^2(1+\lambda^2)\lambda   ) \left ( \frac{(1+\lambda^2)}{(1-\lambda)^2}e^2 -m^2 \right ) \nonumber\\   a_{x,0 } & = & -\phi^2(1+\lambda^2 ) + r^2(1+\lambda^2)\left ( \frac{(1+\lambda^2)}{(1-\lambda)^2}e^2 -m^2 \right ) .
\label{eqn : sing - xkoeff}\end{aligned}\ ] ] a separation of variables gives the hyperelliptic integral @xmath153 a canonical basis of holomorphic ( @xmath154 ) and meromorphic ( @xmath155 ) differentials associated with the hyperelliptic curve @xmath156 is given by ( see @xcite or @xcite ) @xmath157 @xmath158 furthermore we introduce the holomorphic and meromorphic period matrices @xmath159 and @xmath160 : @xmath161 with @xmath162 , where @xmath163 is the canonical basis of closed paths
. the normalized holomorphic differentials are @xmath164 the solution of equation ( [ eqn : sing - psi - x - quadrgleichung ] ) is extensively discussed in @xcite , and is given by the derivatives @xmath165 of the kleinian sigma function @xmath166 : @xmath167 where @xmath168 the constant @xmath169 depends on @xmath170 and @xmath171 only .
@xmath172 is defined by the vanishing condition of the kleinian sigma function @xmath173 so that @xmath174 is an element of the theta divisor @xmath175 ( the set of zeros of the theta function ) where @xmath176 is the vector of riemann constants and @xmath30 is the riemann period matrix defined as @xmath177 .      with
( [ eqn : sing - psi - x - gleichung ] ) equation ( [ eqn : sing - psi - phi - gleichung ] ) yields @xmath178 or @xmath179 the integral ( [ eqn : psi - phiint ] ) has poles at @xmath180 and @xmath181 .
now we apply a partial fractions decomposition upon ( [ eqn : psi - phiint ] ) : @xmath182 where @xmath183 are constants which arise from the partial fractions decomposition .
the differentials in the equation above are of the third kind and can be solved with the help of the following equation ( see @xcite ) . @xmath184\\ & + \ln \frac{\sigma\left(\int_{\infty}^p \mathrm{d}\boldsymbol{u}- \int_{(e_2,0)}^{(z , w ) } \mathrm{d}\boldsymbol{u } - \boldsymbol{k}_\infty   \right)}{\sigma\left(\int_{\infty}^p \mathrm{d}\boldsymbol{u}+ \int_{(e_2,0)}^{(z , w ) } \mathrm{d}\boldsymbol{u } - \boldsymbol{k}_\infty \right ) } - \mathrm{ln } \frac{\sigma\left(\int_{\infty}^{p ' } \mathrm{d}\boldsymbol{u } - \int_{(e_2,0)}^{(z , w ) } \mathrm{d}\boldsymbol{u } - \boldsymbol{k}_\infty   \right)}{\sigma\left(\int_{\infty}^{p ' } \mathrm{d}\boldsymbol{u } + \int_{(e_2,0)}^{(z , w ) } \mathrm{d}\boldsymbol{u } - \boldsymbol{k}_\infty   \right)}.\end{split } \label{main1 - 2}\end{aligned}\ ] ] @xmath185 and @xmath186 are points on the hyperelliptic curve , @xmath187 is a pole , @xmath188 and @xmath156 . the zeros @xmath189 of @xmath190 are the branch points of the curve @xmath191 .
@xmath192 is the vector of the holomorphic differentials of the first kind @xmath193 with @xmath194 .
@xmath195 and @xmath196 are kleinian functions and @xmath197 is the vector of riemann constants .    the vector @xmath198 identified with each branch point @xmath189 is defined as @xcite @xmath199 with the vectors @xmath200 and @xmath201 whose entries @xmath202 , @xmath203 are @xmath204 or @xmath205 for every @xmath206 , @xmath207 .
the matrix @xmath208 = \left [   \begin{array}{c }    \boldsymbol{\varepsilon}_i'\\    \boldsymbol{\varepsilon}_i   \end{array } \right]\ ] ] is called the characteristic of a branch point @xmath189 .
+ the @xmath209th component ( in this case genus @xmath210 ) of the vector @xmath211 is @xmath212 and for @xmath213 we have @xmath214 the @xmath215 are elementary symmetric functions of order @xmath216 built on @xmath217 branch points @xmath218 : @xmath219 , @xmath220 , etc .
+ then the solution of the @xmath7-equation reads @xmath221   \end{split}\ ] ] where @xmath222 and @xmath223 .      with ( [ eqn : sing - psi - x - gleichung ] ) equation ( [ eqn : sing - psi - t - gleichung ] ) yields @xmath224 or @xmath225 next we apply a partial fractions decomposition upon ( [ eqn : psi - tint ] ) where the constants @xmath226 arise .
@xmath227 first we will solve the holomorphic integrals @xmath228 .
we introduce a variable @xmath229 so that @xmath230 .
the inversion of this integral yields @xmath231 ( see section [ sec : sing - x - solution ] or @xcite ) , where @xmath232 the function @xmath233 can be found from the condition @xmath234 .
equation ( [ eqn : sing - t - pbz ] ) now reads @xmath235 \ , .\ ] ] the remaining differential is of the third kind .
its solution is presented in ( [ main1 - 2 ] ) .
then the solution of the @xmath25-equation ( [ eqn : sing - psi - t - gleichung ] ) is @xmath236 + k_2(v - v_0 ) + k_3f_1(v - v_0 ) \right\rbrace \\ \end{split}\ ] ] where @xmath237 .      on the rotational axis of the singly spinning black ring bound orbits and escape orbits are possible .
the orbits either move around the ring or move directly through the center of the ring . figure [ pic : sing - psi - bo ] and [ pic : sing - psi - bo2 ] show bound orbits and figure [ pic : sing - psi - eo ] and [ pic : sing - psi - eo2 ] show escape orbits in @xmath114-@xmath115-coordinates ( see section [ sec : em0-orbit ] ) and the corresponding solution @xmath238 . since the orbits on the rotational axis are presented by lines in the @xmath35-@xmath36-plane , we will also show them in the @xmath35-@xmath7-plane where we will use the coordinates @xmath239 and @xmath7
.    one can think of ring coordinates as two pairs of polar coordinates @xmath240 where @xmath241 ( see @xcite ) .
if @xmath8 is constant , the horizon of the black ring consists of two @xmath242 spheres .
if we look at the rotational axis where @xmath43 , the coordinates @xmath243 and @xmath244 describe the plane between these two spheres , so the horizon can not be seen in this plane . if @xmath7 is constant , the horizon has @xmath245 topology .
so if @xmath246 the coordinates @xmath247 and @xmath248 describe the equatorial plane `` as seen from above '' .
+ the bound orbits of type b and the escape orbits of type a , which move through the center of the black ring , are lines in every plane since here both angles @xmath7 and @xmath8 are constant .
so in that case we only show the @xmath114-@xmath115-plot .
the surface @xmath42 is the equatorial plane of the black ring , which is divided into two parts .
the first part @xmath146 is the plane enclosed by the ring ( or more precisely : enclosed by the singularity ) , which we will refer to as `` inside '' the ring .
the second part @xmath144 describes the equatorial plane around the black ring ( or more precisely : around the singularity ) , which we will refer to as `` outside '' the ring .
if we set @xmath42 , @xmath249 and @xmath250 in the hamilton - jacobi equation ( [ eqn : hjd - sring ] ) , it depends on the coordinate @xmath36 only : @xmath251 \ , .\ ] ] this can be rearranged to @xmath252 - \frac{h(y)}{g^2(y)}(\psi + \omega_\psi e)^2 : = y_s \ , .\ ] ] then we have @xmath253 now we set the derivatives of @xmath130 with respect to the constants @xmath131 , @xmath32 and @xmath34 to zero in order to obtain the equations of motion . with the mino - time @xcite @xmath254 and the relation @xmath255 the equations of motion
take the form @xmath256 - \frac{(\pm 1-y)^2h(y)}{(1\pm\lambda)^4 } [ h(y)\psi + \omega_\psi]^2 \right\rbrace ^{1/2 }   \nonumber\\        & : = & \sqrt{y(y ) } \ , , \label{eqn : sing - phi - y - gleichung}\\   \frac{\mathrm{d}\psi}{\mathrm{d}\gamma } & = & -\frac{(\psi + \omega_\psi e)(\pm 1-y)h(y)}{(1\pm\lambda)^2g(y ) } \ , , \label{eqn : sing - phi - psi - gleichung}\\   \frac{\mathrm{d}t}{\mathrm{d}\gamma } & = & \frac{r^2e}{(\pm 1-y)h(y ) } + \frac{(\omega_\psi\psi+\omega_\psi^2e)(\pm 1-y)h(y)}{(1\pm\lambda)^2g(y ) } \label{eqn : sing - phi - t - gleichung } \ , .\end{aligned}\ ] ] it might not be obvious at first glance , but @xmath62 is a polynomial of third order in @xmath36 and therefore the equations of motion are of elliptic type .
from ( [ eqn : sing - phi - y - gleichung ] ) we can read off an effective potential consisting of the two parts @xmath257 and @xmath258 : @xmath259 since @xmath62 can be written as @xmath260 the effective potential takes the form @xmath261 @xmath262 the two cases @xmath146 ( geodesics inside the ring ) and @xmath144 ( geodesics outside the ring ) have to be discussed separately .
let us first take a look at the motion on the surface outside the black ring .
here we have @xmath144 .
figure [ pic : sing - phi - orbits1 ] shows the effective potential @xmath53 for different values of the parameters . @xmath263 and @xmath264 meet at the horizon .
mainly the angular momentum @xmath34 defines the shape of the effective potential . for @xmath125
the potential is symmetric and @xmath62 has none or one zero . if @xmath265 the potential is no longer symmetric and if @xmath266 is large enough up to two zeros of @xmath62 are possible .
possible orbits are terminating orbits ( to ) with or without a turning point , where light or test particles cross the horizon and fall into the singularity , and escape orbits ( eo ) , where light or test particles aproach the black ring , turn around at a certain point and escape the gravitational field .
the zero of @xmath62 of a to can lie directly on the event horizon .
there are three different types of orbits ( see table [ tab : sing - phi - typen - orbits1 ] ) .
* type a : + @xmath62 has no zeros and only tos exist . *
type b : + @xmath62 has one zero and only tos exist . in a special case
the zero of @xmath62 lies on the horizon . *
type c : + @xmath62 has two zeros .
tos and eos exist . in a special case the zero of @xmath62
lies on the horizon .
a & 0 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( 0.5,-0.2)(0.5,0.2 ) ( -2.5,0)(3,0 )    & to + b & 1 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( 0.5,-0.2)(0.5,0.2 ) ( -2.5,0)(1,0 )    & to + b@xmath267 & 1 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( 0.5,-0.2)(0.5,0.2 ) ( -2.5,0)(0.5,0 )    & to + c & 2 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( 0.5,-0.2)(0.5,0.2 ) ( -2.5,0)(1,0 ) ( 2.0,0)(3,0 )    & to , eo + c@xmath267 & 2 &    ( -2.5,-0.2)(3,0.2)(-2.5,0)(3,0 ) ( 0.5,-0.2)(0.5,0.2 ) ( -2.5,0)(0.5,0 ) ( 2.0,0)(3,0 )    & to , eo +      the effective potential for geodesics on the surface enclosed by the black ring ( @xmath146 ) is shown in figure [ pic : sing - phi - orbits2 ] . again , if we have @xmath125 the potential is symmetric and @xmath62 has none or one zeros . in the case
@xmath265 a potential barrier appears which prevents test particles and light from reaching @xmath43 .
then @xmath62 has always a single zero in the allowed range of @xmath36 .
the higher @xmath266 the higher the energy where @xmath263 and @xmath264 meet .
note that @xmath146 , @xmath43 is the location of the center of the black ring .
possible orbits are terminating orbits where @xmath62 has one zero ( type b ) or none zero ( type a ) .
for the type of orbits see previous section . in a special case
the zero of @xmath62 lies on the horizon .
( type b@xmath267 ) .    , @xmath268 , @xmath122 und @xmath145 effective potentials @xmath257 ( red , solid ) and @xmath258 ( blue , dotted ) on the equatorial plane inside the ring .
the grey area is a forbidden zone , where no motion is possible .
the horizon is marked by a vertical dashed line .
green dashed lines represent energies and green points mark the turning points .
possible orbits are of type a , b and b@xmath267.,width=226 ]      equation ( [ eqn : sing - phi - y - gleichung ] ) can be solved analogously to ( [ eqn : sing - x - gleichung ] ) . ( [ eqn : sing - phi - y - gleichung ] ) can be written as @xmath269 where @xmath270 the solution is ( see section [ sec : ergo - xsol ] ) @xmath271 , \ ] ] where @xmath272 and @xmath273 .
the coefficients @xmath274 and @xmath275 of the polynomial in the weierstra form are @xmath276      with ( [ eqn : sing - phi - y - gleichung ] ) equation ( [ eqn : sing - phi - psi - gleichung ] ) yields @xmath277 or @xmath278 this can be rewritten as @xmath279 this equation can be solved analogously to the @xmath7- and @xmath8-equation for nullgeodesics in the ergosphere ( see section [ sec : ergo - phisol ] and [ sec : ergo - psisol ] ) . with @xmath280 , @xmath93 and @xmath95
the solution is @xmath281 @xmath102 are constants which arise from the partial fractions decomposition and depend on the parameters of the metric and the test particle .      with (
[ eqn : sing - phi - y - gleichung ] ) equation ( [ eqn : sing - phi - t - gleichung ] ) yields @xmath282}{(1\pm\lambda)^2(1-y)(1+\lambda y)h(y ) } \right ) \frac{\mathrm{d}y}{\sqrt{y(y)}}\ ] ] or @xmath283}{(1\pm\lambda)^2(1-y')(1+\lambda y')h(y ' ) } \right)\ , \frac{\mathrm{d}y'}{\sqrt{y(y ' ) } } \ , .
\label{eqn : phi - tint}\ ] ] this equation can be solved analogously to the @xmath7- and @xmath8-equation for nullgeodesics in the ergosphere ( see section [ sec : ergo - phisol ] and [ sec : ergo - psisol ] ) . with @xmath280 , @xmath93 and @xmath284
the solution is @xmath285 @xmath111 are constants which arise from the partial fractions decomposition and depend on the parameters of the metric and the test particle .      on the equatorial plane around the singly spinning black ring terminating orbits and
escape orbits are possible .
figure [ pic : sing - phiout - eo ] and [ pic : sing - phiout - to ] show some orbits in @xmath114-@xmath115-coordinates ( see section [ sec : em0-orbit ] ) and the corresponding solution @xmath286 .
also the @xmath36-@xmath8-plane is shown in the coordinates @xmath247 and @xmath248 ( see section [ sec : psi - axis - orbits ] ) .
an escape orbit is depicted in figure [ pic : sing - phiout - eo ] .
figure [ pic : sing - phiout - to ] shows a terminating orbit which starts at its turning point and then falls into the singularity .
the frame dragging effect can be seen in figure [ pic : framedragg ] .
once the particle enters the ergosphere it is dragged along by the rotation of the black ring . if the angular momentum of the particle and the black ring have opposite signs , the particle changes its direction when approaching the ergosphere .
+    on the equatorial plane enclosed by the black ring only terminating orbits are possible .
figure [ pic : sing - phiin - to ] shows a terminating orbit which starts at the center of the black ring and then falls into the singularity .
in this paper we presented the analytical solutions of the geodesic equations of the singly spinning black ring for special cases .
since the hamilton - jacobi equation seems not to be separable in general , we had to concentrate on the nullgeodesics in the ergosphere ( @xmath1 ) , geodesics on the rotational axis ( @xmath43 ) and geodesics on the equatorial plane ( @xmath246 ) .
+ we discussed the general structure of the orbits and gave a complete classification of their types .    in the ergosphere
there is just one possible orbit , where light crosses the event horizon and falls inevitably into the singularity ( terminating orbit ) .
the @xmath35-motion bounces back and forth between two values or stays constant at @xmath61 , while @xmath36 ranges from a turning point to @xmath287 . on the rotational axis
@xmath36 is constant , so here the @xmath35-motion determines the type of orbit .
we found escape orbits and bound orbits , the latter were also shown numerically by igata et al .
@xcite . on the equatorial plane we found terminating orbits and
escape orbits .
+ the separability of the hamilton - jacobi equation is a coordinate related phenomenon , so one might think of a coordinate system in which it would be possible to separate the hamilton - jacobi equation in general .
but recently igata , ishihara and takamori found evidence of chaotic motion in the singly spinning black ring spacetime using the poincar map @xcite . from that one
could conclude that it is not possible to separate the hamilton - jacobi equation in any coordinate system .
+ besides the singly spinning black ring , one could consider black rings with two angular momenta ( doubly spinning black ring @xcite ) or add charge to the black ring @xcite . also a supersymmetric black ring solution was found @xcite , @xcite .
the methods shown in this paper can be applied to ( charged ) doubly spinning black rings as well as to supersymmetric black rings .
this will be done in future work .
we would like to thank victor enolski , norman grlebeck and volker perlick for helpful discussions .
we gratefully acknowledge support by the dfg , in particular , also within the dfg research training group 1620 `` models of gravity '' .
r.  emparan and h.  s.  reall , living rev .
rel .   * 11 * , 6 ( 2008 ) [ arxiv:0801.3471 [ hep - th ] ] .
r.  c.  myers and m.  j.  perry , annals phys .   *
172 * , 304 ( 1986 ) .
r.  emparan and h.  s.  reall , phys .
lett .   * 88 * , 101101 ( 2002 ) [ hep - th/0110260 ] .
d.  kubiznak and v.  p.  frolov , class .
* 24 * , f1 ( 2007 ) [ gr - qc/0610144 ] .
d.  n.  page , d.  kubiznak , m.  vasudevan and p.  krtous , phys .
lett .   * 98 * , 061102 ( 2007 ) [ hep - th/0611083 ] .
v.  p.  frolov , p.  krtous and d.  kubiznak , jhep * 0702 * , 005 ( 2007 ) [ hep - th/0611245 ] .
j.  hoskisson , phys .
d * 78 * , 064039 ( 2008 ) [ arxiv:0705.0117 [ hep - th ] ] .
m.  durkee , class .
* 26 * , 085016 ( 2009 ) [ arxiv:0812.0235 [ gr - qc ] ] .
h.  elvang , r.  emparan and a.  virmani , jhep * 0612 * , 074 ( 2006 ) [ hep - th/0608076 ] .
a.  a.  pomeransky and r.  a.  senkov , hep - th/0612005 .
j.  armas , class .
* 28 * , 235014 ( 2011 ) [ arxiv:1011.5618 [ hep - th ] ] .
t.  igata , h.  ishihara and y.  takamori , phys .
d * 82 * , 101501 ( 2010 ) [ arxiv:1006.3129 [ hep - th ] ] . y. hagihara , em jpn . j. astron .
* 8 * , 67 ( 1931 )    e.  hackmann and c.  lmmerzahl , phys .
d * 78 * , 024035 ( 2008 ) .
v.  z.  enolski , e.  hackmann , v.  kagramanova , j.  kunz and c.  lmmerzahl , j.  geom .
phys .   * 61 * , 899 ( 2011 ) [ arxiv:1011.6459 [ gr - qc ] ] .
e.  hackmann , v.  kagramanova , j.  kunz and c.  lmmerzahl , phys .
d * 78 * , 124018 ( 2008 ) [ erratum - ibid .   * 79 * , 029901 ( 2009 ) ] [ arxiv:0812.2428 [ gr - qc ] ] . v. z. enolski , m. pronine , p.h .
richter , _ j. nonlinear sci .
* 13 * , pp .
157 - 174 ( 2003 )    r.  emparan and h.  s.  reall , class .
* 23 * , r169 ( 2006 ) [ hep - th/0608012 ] . y.  mino , phys .  rev .
d * 67 * , 084027 ( 2003 ) [ gr - qc/0302075 ] .
a. i. markushevich , _ theory of functions of a complex variable _
( prentice - hall , englewood cliffs , nj , 1967 ) , vol .
v.  enolski , b.  hartmann , v.  kagramanova , j.  kunz , c.  lmmerzahl and p.  sirimachan , journal of mathematical physics * 53 * , 012504 ( 2012 ) [ arxiv:1106.2408 [ gr - qc ] ] .
y. k. lim , e. teo , _ geodesics through the black ring _ , report , university of singapore , http://staff.science.nus.edu.sg/  phylyk / downloads / reports / blackrings_report.pdf ( 2008 )    t.  igata , h.  ishihara and y.  takamori , phys .
d * 83 * , 047501 ( 2011 ) [ arxiv:1012.5725 [ hep - th ] ] .
h.  elvang , phys .
d * 68 * , 124016 ( 2003 ) [ hep - th/0305247 ] .
j.  hoskisson , phys .
d * 79 * , 104022 ( 2009 ) [ arxiv:0808.3000 [ hep - th ] ] .
d.  v.  galtsov and n.  g.  scherbluk , phys .
d * 81 * , 044028 ( 2010 ) [ arxiv:0912.2771 [ hep - th ] ] .
h.  elvang , r.  emparan , d.  mateos and h.  s.  reall , phys .
lett .   * 93 * , 211302 ( 2004 ) [ hep - th/0407065 ]
. h.  elvang , r.  emparan , d.  mateos and h.  s.  reall , phys .
d * 71 * , 024033 ( 2005 ) [ hep - th/0408120 ] .