dwarf galaxies are thought to dominate the cosmic scenery in terms of number density @xcite and therefore form an important morphological class of galaxies .
several investigators have studied the mass distributions of these galaxies ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* and references therein ) .
@xcite was the first to study the dark matter ( dm ) properties of a large representative sample of nearby dwarf galaxies as part of the westerbork  survey of spiral and irregular galaxies ( whisp , * ? ? ? * ) .
more recently , a morphologically diverse sample of nearby galaxies , including dwarf systems , was observed as part of the  nearby  galaxy  survey ( things , * ? ? ?
these cumulative efforts have led to an understanding that dm plays an important role in the dynamics and evolution of dwarf galaxies .
@xcite found that none of the observed rotation curves in his sample of late - type dwarf galaxies decline at outer radii .
this was further confirmation that these systems have a large dark mass fraction .
dwarf galaxies therefore serve as ideal candidates to test theories of dm .
observational studies of the mass distributions of these dwarf systems have shown their dm halos to be well - modelled by a pseudo iso - thermal sphere with an approximately constant - density core @xcite .
theoretically , the study of the properties of dm halos in a hierarchical clustering context is carried out mainly in the form of cold dark matter ( cdm ) numerical simulations ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the equilibrium density profile , @xmath7 , of simulated dm halos is found to vary with radius , @xmath8 , as @xmath9 with @xmath10 @xcite , thereby predicting extremely steep inner density profiles .
this discrepancy between the observed shapes of dm halos and the predictions from numerical simulations has become known as the `` cusp / core '' problem .
attempts at reconciling observations with theoretical predictions depend crucially on accurate dynamical analyses of nearby galaxies . to test dm halo models
one wants galaxies that are as dark - matter - dominated as possible . in this respect ,
a potentially useful candidate for dm studies is ngc  2915 .
this galaxy is classified as a nearby ( @xmath11  3.78  mpc , * ? ? ?
* ) blue compact dwarf according to its optical appearance , yet has the  morphology of a late - type spiral .
detailed imaging by @xcite showed that the optical appearance of ngc  2915 is dominated by two main stellar populations : a compact blue population , which is the location of on - going high - mass star formation ; and a more diffuse , older red population .
what makes this galaxy so interesting is the fact that , when observed at 21  cm , the stellar core of ngc  2915 is seen to be completely embedded in a huge  disk extending out to @xmath12  22 @xmath1-band scale lengths @xcite .
furthermore , this gas disk has well - defined spiral structure that is not seen in the stellar disk .
no significant star - formation is observed in the outer parts of this gas disk .
only a few faint hii regions were detected by @xcite after carrying out deep h@xmath13 imaging of the disk .
the  disk of ngc  2915 serves as an ideal tracer of the gravitational potential out to radii of @xmath12  10  kpc , far beyond the visible radial extent of the stellar disk .
due to the apparent lack of stars in its outer parts , the dynamics of the  disk are an almost direct tracer of the dm distribution .
@xcite targeted ngc  2915 for dm studies and found the galaxy to be dm - dominated at nearly all radii with a total mass to @xmath1-band light ratio m@xmath14/l@xmath3  @xmath15  76  /l@xmath16 ( assuming a distance of 5.1 mpc ) .
they also found the galaxy to have a dense ( @xmath17  pc@xmath5 ) and compact ( @xmath18  1  kpc ) dm core .
several investigators have attempted to explain the observed  distribution of ngc  2915 .
@xcite studied the dynamics of the central  region and the spiral pattern of the outer disk . for both
they calculated a common , slow pattern speed of @xmath19  which they associated with the figure rotation of a tri - axial dm halo .
they also proposed that some dm is distributed in the disk of ngc  2915 , thereby making it gravitationally unstable to the formation of the observed spiral structure .
@xcite , using hydrodynamical simulations , further explored the ideas of @xcite .
they showed that the observed spiral structure can be accounted for by either an unseen bar or a rotating tri - axial dm halo .
however , the mass of the required bar , @xmath20  , is very large in comparison to the total stellar mass , thereby making the nature of such a bar problematic .
the required pattern speed of the tri - axial dm halo is significantly larger than those from numerical simulations .
@xcite disfavoured the external perturber scenario .
they also found that while a heavy disk is able to account for the main features of the observed  morphology , it fails to match the observed gas dynamics .
a satisfactory explanation for the various morphological and kinematic features of the  in ngc  2915 is therefore still lacking .    in this paper
we investigate the  dynamics and dm distribution of ngc  2915 .
we use data from new  synthesis observations of ngc  2915 carried out using the australian telescope compact array as part of the southern hemisphere extension of the  nearby galaxy survey .
our  observations are significantly deeper and have better spatial resolution than any other  observations of ngc  2915 .
these data are complemented by high - quality 3.6  infrared observations of the stellar disk , carried out as part of the _ spitzer _ infrared nearby galaxies survey ( sings , * ? ? ?
the  observations are presented in sec .
[ data_cube ] while the  data products appear in sec .
[ data_products ] .
this paper focuses on the regular gas dynamics of ngc  2915 . in sec .
[ vrot ] a rotation curve out to @xmath12  9.3  kpc is determined . the far - infrared observations of the stellar disk are combined with our new  observations in sec .  [ mass_modeling ] to produce a mass model that determines the contributions at various radii of the stars , gas and dm .
we use two different parameterisations of the dm halo to reconstruct the observed rotation curve , one of which is the nfw halo @xcite favoured by numerical simulations while the other is the observationally motivated pseudo - isothermal sphere . finally , in sec .
[ conclusions ] , we summarise our results and present our conclusions .
ngc 2915 was observed between 23  october  2006 and 2  june  2007 ( project number c  1629 ) with six different atca configurations using all six antennas .
a single run consisted of a primary calibrator observation , regular secondary calibrator observations and source observations .
the calibrator sources pks  193463 and pks  0823500 were used as primary and secondary calibrators respectively .
each of the ew-352 , 750d , 1.5b and 1.5c runs was approximately 12 hours long while 24 hours were spent in the 6a configuration .
besides these data , archival data were also incorporated .
@xcite determined the hi diameter of ngc 2915 to be @xmath21  0.32@xmath22 .
the hi extent of ngc 2915 therefore falls well within the @xmath21  0.54@xmath22 field - of - view of the atca dishes when observing at 21  cm . for our observations , the telescope pointing centres
were set to the optical centre of ngc 2915 ( @xmath23  =  09@xmath24  26@xmath25  11.5@xmath26 , @xmath27  =  -76@xmath22
37@xmath28  35@xmath29 ) with no mosaicking required .
the correlator was set to use 512 channels with a bandwidth of 8  mhz , centered at 1418  mhz .
the resulting velocity range is @xmath30  to 1353  with an approximate channel spacing of 3.2  .
atca_table ] provides a summary of all our observing setups .
* comments on columns : * column 1 : atca configuration used ; column 2 : date of observation ( ut ) ; column 3/4 : start / end of observations ( ut ) ; column 5 : time on source .
the miriad software package @xcite was used to take the raw @xmath31 data from the correlator through to the image analysis stage .
all velocities are measured relative to the barycentre rest frame .
the first five and last five correlator channels were flagged .
the data were then split into primary calibrator , secondary calibrator , and source subsets .
the required corrections to the antenna gains , delay terms and bandpass shapes were determined .
the time - varying phases and antenna gains were calibrated based on observations of the secondary calibrator .
the calibrated source data were continuum - subtracted by fitting and subtracting a first - order polynomial to the line - free channels which , for our observations , were channels @xmath32 and @xmath33 .
having isolated the hi signal , image cubes were produced from the continuum - subtracted @xmath31 data .
data cubes were produced using natural weighting as well as robust weighting with a brigg s visibility weighting robustness parameter of 0.2 .
each of the dirty images was deconvolved using a steer clean algorithm @xcite to produce an output map of clean components .
channels were either cleaned down to a flux cut - off of 2.5 times the typical r.m.s .
of the flux in a line - free channel or for 50 000 iterations , whichever condition was met first .
the typical number of iterations was @xmath21 20000  30000 .
after the deconvolution process , each of the clean components was convolved with a gaussian approximation of the dirty beam .
the beam s full width at half maximum was 17@xmath34 and @xmath35 for the naturally - weighted ( na ) and robust - weighted ( rw ) data cubes respectively while the channel width was set to @xmath36  for both .
no hanning smoothing was applied to either data cube . for each cube ,
the noise in a line - free channel is gaussian distributed with a mean of @xmath37  mjy  beam@xmath38 and a standard deviation of @xmath39  0.60  mjy  beam@xmath38 for the na cube , while @xmath40  mjy  beam@xmath38 and @xmath41  mjy  beam@xmath38 for the rw cube .
all hi moment maps were generated using the groningen image processing system ( gipsy , @xcite ) together with the na  data cube .
the entire  data cube was smoothed , using a gaussian convolving function , to a resolution of 34@xmath42 .
a flux cut - off was applied to the smoothed cube at 2.5@xmath43 . for each channel of the smoothed data cube ,
any flux that was clearly not spatially connected to , or associated with the galaxy emission was flagged and removed .
the blotted smoothed cube was then applied as a mask to the original @xmath44 cube .
the hi total intensity and 2nd - order moment maps were constructed from the zeroth- and second - order moments of the data cube respectively . a map in which each pixel is equal to @xmath45 , where @xmath46 is the number of un - blanked channels for a given line profile that would contribute to the moment maps , was created .
such a construction was permitted since no hanning smoothing was applied to the  data cube . dividing each pixel in the total intensity map by @xmath45 produced a dimensionless signal - to - noise ( s / n ) map .
the average intensity of all the pixels in the total intensity map that had a corresponding signal - to - noise ratio in the range 2.75 @xmath47 s / n @xmath47 3.25 was used as an intensity cut - off for the total intensity map .
this removed , in a robust manner , pixels from the outer disk that had low s / n values .
this new  total intensity map was then used as a mask on other hi maps .
the method of using gauss - hermite polynomials to parameterise line profiles of galaxies was introduced by @xcite .
while gaussians are reasonable approximations to many realistic line profiles , higher order gauss - hermite polynomials are able to better capture profile asymmetries and non - gaussian deviations .
the hermite method has already been successfully implemented by various authors @xcite .
we fitted a third - order gauss - hermite polynomial to each na  data cube line profile to generate a so - called gauss - hermite  velocity field .
three filters were used simultaneously when fitting the profiles : 1 ) profiles with fitted peak fluxes below 2.5@xmath43 were excluded ; 2 ) profiles with a fitted line width less than the channel width were excluded ; and 3 ) fitted profile peaks had to be within the velocity range of the data cube .
a grey - scale representation of the @xmath48 resolution channel maps is shown in fig . [ channels ] .
@xmath49  mjy in all channels .
structures are visible at large and small scales . at large scales
, the emission clearly exhibits the usual pattern of a rotating disk .
the  disk of ngc 2915 exhibits clear spiral structure ( c.f .  total intensity map , fig . [ mom0 ] ) . on small scales , the gas in the central regions forms two dominant over - densities in nearly all of the channels , with fainter emission in - between them .              a global  profile , extracted from the blotted na  data cube by summing the  emission in each velocity channel , is presented in fig .
[ global_hi_profiles ] . using a synthesis array to estimate
the total  flux of a galaxy can result in an underestimate due to missing short baselines .
the shortest baseline used for our observations was 31  m ( as part of the ew352 antenna configuration ) , implying that the observations are less sensitive to structures larger than @xmath21 23@xmath28 .
@xcite , however , estimated the hi extent of ngc  2915 to be @xmath21 19@xmath28 .
the amount of hi missed in the global hi profile is therefore expected to be small . in fig .
[ global_hi_profiles ] we compare our determination of the global  profile for ngc  2915 to the hipass global  profile from @xcite .
the hipass data are calibrated with 8@xmath28 spatial binning and 13 km s@xmath38 channel width . due to the large hipass beam , as well as the fact that the parkes telescope is a single dish telescope
, the amount of flux missed in the hipass observations should be negligible .
it is clear from fig .
[ global_hi_profiles ] that no significant amount of flux has been missed by our observations .
the resulting global  profile is double - horned and asymmetric . from this profile
we derived the profile widths , the systemic velocity , and the total hi mass .
the profile widths were calculated as the differences between the high and low velocities of the galaxy at 20@xmath50 and 50@xmath50 of the peak flux density ( v@xmath51 , v@xmath52 , v@xmath53 , v@xmath54 respectively ) .
the profile widths were not corrected for inclination .
two systemic velocity estimates were determined : 1 ) from the global profile , at the midpoint of the two velocities at the 50@xmath50 level , 2 ) from the equation @xmath55 finally , the total hi mass was calculated as @xmath56 where @xmath57 is the distance to the galaxy in mpc and @xmath58 is the total hi line flux in jy  km  s@xmath38 .
this determination of the total hi mass assumes that the hi is optically thin . throughout this paper
we adopt the @xmath59 mpc distance determination of @xcite .
table  [ hi_global_profile_table ] lists our determinations of the above mentioned quantities .
assuming a distance of @xmath60  mpc , @xcite measured a total  mass of @xmath61  for ngc  2915 . using a distance of @xmath62  mpc ,
their total mass estimate reduces to @xmath63  , consistent with our estimate of @xmath64  .
[ hi_global_profile_table ] +    * comments on columns : * column 1 : total hi mass ; column 2 : velocity width at 50@xmath50 of the peak flux ; column 3 : velocity width at 20@xmath50 of the peak flux ; column 4 : systemic velocity from w@xmath65 midpoint ; column 5 : systemic velocity from eqn .
[ vsysest ] .    .
]      the full  total intensity map extracted from the na  data cube is shown in the upper panel of fig .
the edge of the  disk ( the outer - most contour level ) lies at @xmath12  510@xmath29 and therefore extends out to @xmath12  20 @xmath1-band scale lengths . for comparative purposes the irac 3.6  image of the stellar disk
is shown in fig .
[ irac_image ] . to determine the photometric centre we fitted ellipses to three 3.6  flux density annuli near the edge of the old stellar disk .
a 3.6  flux density contour at a level of 1.2  mjy  ster@xmath38 , the average flux density of the second annulus , is shown in fig .
[ irac_image ] .
the average flux densities of the other two annuli are 0.8 and 1.7  mjy  ster@xmath38 .
the average centre position of the ellipses fitted to the flux density annuli , @xmath23  =  09@xmath66  26@xmath67  12.611@xmath68 , @xmath27  =  @xmath69  37@xmath28  37.80@xmath29 , was used to estimate the position of the photometric centre ( shown as a cross in fig .
[ irac_image ] ) . for comparison ,
the 2mass determination of the centre of the system is @xmath23  =  09@xmath66  26@xmath67  11.53@xmath68 , @xmath27  =  @xmath69  37@xmath28  34.80@xmath29 @xcite . to facilitate length - scale comparisons between ir and  images , an ellipse with semi - major axis length @xmath70 , roughly equal to the @xmath71-band @xmath72 radius determined by @xcite ,
is included in the figures . in their study of the  content of 73 late - type dwarf galaxies
, @xcite found that the ratio of the  extent to the optical diameter , defined as 6.4 disk scale lengths , is 1.8 @xmath73 0.8 on average .
the same ratio for ngc  2915 is @xmath74 , close to double that value .
( white ) to 0.2  jy  beam@xmath38 ( black ) .
contour levels run from 10    60  mjy  beam@xmath38 in steps of 10  mjy  beam@xmath38 and from 80    230  mjy  beam@xmath38 in steps of 20  mjy  beam@xmath38 .
the position of the photometric centre estimated using the 3.6  image is marked with a cross .
the hatched circles in the lower corners represent the half power beam width of the synthesised beam .
the black and white ellipses have a semi - major axis length @xmath70 , an inclination @xmath75 , a position angle @xmath76 and a centre position equal to our determination of the photometric centre .
the solid black lines represent the position - velocity slices shown in fig .
[ pv_slices ] . ]
( white ) to 2.0  mjy  ster@xmath38  ( black ) .
the single white contour is at a level of 1.2  mjy  ster@xmath38 and represents the average flux of one of three flux annuli used to estimate the position of the photometric centre .
the estimated photometric centre is marked with a cross .
the black ellipse is the same as that shown in fig .
[ mom0 ] ]    our high - resolution  observations resolve the central regions of ngc  2915 into two distinct  concentrations ( fig .
[ mom0 ] , lower panel ) .
their centres are separated by @xmath21  60@xmath29 ( 1.1  kpc ) and have a combined mass of @xmath77  .
also visible is a plume - like  feature to the north - west of these central  concentrations with a mass of @xmath78  .
the observed spiral structure of the outer disk is asymmetric : the arm beginning at the north - western central  concentration can be traced over @xmath21 200@xmath22 in azimuthal angle while the arm beginning at the south - eastern  concentration can be traced for only @xmath21 120@xmath22 . from the hi total intensity map an
azimuthally averaged , inclination - corrected hi surface density profile was constructed .
each ring has a width of @xmath79 , a position angle ( measured anti - clockwise from north to the receding major axis ) of @xmath80 and an inclination of @xmath81 .
the profile is shown in fig .
[ surf_dens_prof ] .    in an attempt to better constrain the  distribution of ngc  2915
, an  total intensity map was extracted from the rw  data cube . for the sake of brevity
this map is not shown here , yet it suffices to say that the higher spatial resolution confirms the above - mentioned nature of the central  concentrations .
however , the very low s / n ratio in the outer parts of this map yields it unsuitable for the study of low column density  features , with only the central - most  features standing out about the noise .
this map does not allow for a more accurate study of the  distribution and is therefore not discussed further .          the third - order gauss - hermite velocity field extracted from the na  data cube
is shown in fig .
the blank pixels in this map are those that were rejected by the fitting routine according to the criteria mentioned in sec .
[ moment_maps ] .
@xcite showed that the gauss - hermite parameterisation of a line profile gives a robust estimate of the peak - associated velocity of the profile .
inner contours running approximately parallel to the kinematic minor axis are caused by @xmath82 rotation while the outer contours running radially away from the centre are due to a constant rotation speed .
the overall `` s - shaped '' distortion of the contours is indicative of a kinematic warp within the disk .
small wiggles along the outer contours are caused by streaming gas motions .
the sharp kinks at inner radii suggest the presence of non - circular velocity components within the gas . finally , differences in the shapes of the iso - velocity contours on the receding and approaching halves of the galaxy suggest a certain degree of kinematic lopsidedness within the galaxy .
a third - order gauss - hermite velocity field was also extracted from the rw  data cube . due to the high noise levels in the rw  cube
, most line profiles were rejected by the fitting filters and hence not fitted .
the general characteristics of the fitted peaks of the remaining profiles are very similar to those of the na velocity field , differing on average by @xmath83  . rather than showing this  velocity field here
, a tilted ring model is fitted to it in sec .
[ tr_output ] to demonstrate that the resulting rotation curve is very similar to the one derived using the na third - order gauss - hermite velocity field .
motivated by this result is the decision to use only the na velocity field for our kinematic analyses .     with the thick contour marking the systemic velocity at 465 km s@xmath38 derived by fitting a tilted ring model to the velocity field .
white and black contours represent receding and approaching halves of the galaxy respectively .
the hatched circles in the lower corners represent the half power beam width of the synthesised beam .
the ellipse is the same as that shown in fig .
[ irac_image ] . ]      the second - order  moment map extracted from the na  data cube is shown in fig .
the highest second - order moments are found near the stellar core of the galaxy .
the radial profile of the second - order moment map is shown in fig .
[ disp_prof ] .
if the second - order moments are interpreted as a measure of the gas velocity dispersion , then the average outer disk dispersion of @xmath84 10 km s@xmath38 is close to the constant gas velocity dispersion value of @xmath85  found by @xcite for a sample of things dwarf galaxies .    .
]          position - velocity slices taken though the centre of the na  data cube betray the presence of broad and multi - component line profiles at inner and outer radii .
figure [ pv_slices ] shows two position - velocity slices , the positions of which are shown in the  total intensity map ( fig .
[ mom0 ] , top panel ) .
the first slice was taken at a position angle of 322 , through both of the central  concentrations .
this slice suggests that the galaxy contains a fast - rotating gas component near its centre , clearly manifesting itself as a sharp velocity spike with a peak velocity of @xmath86 550 .
line profiles near the centre of this position - velocity slice show significant departures from gaussianity .
this position - velocity slice also reveals asymmetric line profiles on the receding side of the galaxy , seen in the upper - left quadrant ( @xmath87  , angular offset @xmath88 ) as an ``  beard '' that is lagging in velocity relative to the main disk .
the beard emission spatially coincides with the plume - like  feature seen in the  total intensity map .
the estimated mass of this morphological feature , @xmath78  , is 5.6@xmath50 of total  mass .
one interpretation of this kinematically anomalous gas component is that it is associated with pristine gas being accreted from the nearby inter - galactic space onto the outer disk of ngc 2915 .
if confirmed , this would be the first evidence of such an event in a dwarf system like ngc  2915 .
the second position - velocity slice was extracted almost parallel to the minor axis , through the south eastern central  concentration which contains many double  component profiles .
line profiles in this central region are split by @xmath89  on average , implying that the sharp rise seen toward the centre of the second - order  moment map ( fig .
[ mom2 ] ) over - estimates the true  velocity dispersions of the inner disk . since the disturbed  line profiles are located within @xmath90 of the centre of ngc 2915 , it is plausible that the central gas dynamics of ngc  2915 are largely dictated by the stellar winds of its young stellar population . in a forthcoming paper ( elson et al .
2010 , in prep . )
we present detailed modeling results of the central  dynamics .
( top ) and 2 ) a position angle of 220@xmath22 ( bottom ) .
the slice positions are shown in the  total intensity map ( fig .
[ mom0 ] ) .
the thickness of the slices is that of a single pixel in the data cube ( 2.5@xmath29 ) and the velocity resolution is 3.5 km s@xmath38 .
the intensity scale runs from 3@xmath43 to ( white ) 40@xmath43 ( black ) .
contours start at 2@xmath43 in steps of 3@xmath43 .
the dashed horizontal line marks the systemic velocity at @xmath91  . ]
the hi line data allow us to study the hi gas dynamics . in this section
we derive a rotation curve for ngc  2915 .
the standard method of producing a rotation curve involves fitting a tilted ring model to a velocity field @xcite .
ngc  2915 s intermediate inclination of @xmath9255@xmath22 makes it a well - suited candidate for such an analysis .
the method involves modeling the galaxy as a set of concentric rings within which the gas orbits about the kinematic centre .
each ring has a set of defining parameters : central coordinates @xmath93 and @xmath94 , inclination @xmath95 , systemic velocity @xmath96 , position angle @xmath97 ( in the sky plane ) , and the rotation velocity @xmath98 of the material within the ring . in the case that expansion velocities within the disk are ignored , the standard algorithm carries out a least squares fit to @xmath99 where @xmath100 is the line - of - sight velocity , @xmath101 and @xmath102 are rectangular coordinates on the sky and @xmath103 is the angle from the major axis in the galaxy plane .
@xmath103 is related to the position angle @xmath97 in the sky plane by @xmath104 where @xmath8 is the radius of the ring in the galaxy plane .
the gipsy task rotcur @xcite was used to fit tilted ring models to the third - order gauss - hermite hi velocity field extracted from the na  data cube .
rings width of @xmath105 were used to ensure that adjacent rings were largely independent of one another . when fitting eqn .
[ vlos ] to the data , each datum was weighted by @xmath106 so that points closer to the major axes held more weight .
all points within 10@xmath22 of either side of the minor axes were excluded from the fit .
both sides of the galaxy were used for the fitting .
each tilted ring should be centred on the dynamical centre of the galaxy , assuming it is constant , which is usually estimated by the least - squares fitting algorithm together with the other tilted ring parameters . in the case of ngc  2915 , however , the kinematic centre determinations from rotcur varied with radius in a non - systematic manner .
we therefore used the photometric centre of the 3.6  emission as the kinematic centre .
using the photometric and kinematic centres interchangeably in this manner is reasonable if one assumes that the stars lie at the bottom of the gravitational potential .
we carried out a single rotcur iteration with all parameters allowed to vary freely .
the resulting @xmath93 and @xmath94 values are shown as filled black circles in panels a and b respectively of fig .
[ tr_results ] . these fitted @xmath93 and
@xmath94 positions deviate on average by 8.3@xmath29  @xmath73  @xmath107 and 7.2 @xmath29  @xmath73 0.25@xmath29 ( 152  pc and @xmath21 132  pc ) respectively from the photometric centre , thereby placing them well within a single half power beam width .
this average deviation is consistent with the results of @xcite who found that approximately 50@xmath50 of their kinematic centres derived for @xmath21  1000 individual tilted rings fitted to the things  velocity fields differed by less than a beam width from their best centre estimates . in the case of the things sample
, one beam width corresponds to @xmath21  200  pc on average . throughout the tilted ring fitting procedure , @xmath93 and @xmath94
were therefore each fixed to the position of the photometric centre ( solid lines in panels a and b of fig .
[ tr_results ] )        the remaining tilted ring parameters ( @xmath96 , @xmath97 , @xmath95 , @xmath98 ) were derived iteratively .
we first allowed all parameters to vary freely .
this gave a general feel for the parameter behaviour with @xmath8 .
little scatter was seen in the radial run of @xmath96 , which was then fixed to the average value of 465  in subsequent iterations ( compare with the estimate of @xmath108  derived from the global  profile using eqn .  [ vsysest ] ) .
the @xmath97 , @xmath95 , and v@xmath109 parameters were allowed to vary freely in the subsequent iterations .
the @xmath97 and @xmath95 runs showed relatively little scatter except for a clear change from outer to inner disk beginning at @xmath71 @xmath21 240@xmath29 where the inclination value changed from @xmath92  55@xmath22 to @xmath92  75@xmath22 .
our kinematically derived position angles were consistent with the position angles of ellipses fitted to flux contours in the hi total intensity map .
several more iterations with various combinations of fixed and free parameters were carried out in order to check the stability of the radial parameter profiles .
the sharp rise in @xmath95 for r @xmath110 240@xmath29 always remained with little scatter .    .
error bars represent the r.m.s .
spread of the velocities in a given ring . ]        for the final iteration , the @xmath111 , @xmath112 and @xmath96 parameters were fixed to values that were constant with @xmath8 while the inclination and position angle radial profiles were smoothed and then fixed .
the smoothed versions of the profiles are shown as solid curves in panel a , b , c , d , e in fig .
[ tr_results ] while the filled circles represent the parameter values just before they were smoothed and fixed . smoothing the @xmath97 and @xmath95 radial distributions
did not affect the final rotation curve in any noticeable manner .
as @xcite point out , the smoothing only affects the point - to - point scatter of the profiles and in no way affects the resolution of the radial distributions . to derive the final rotation curve , only @xmath98
was allowed to vary .
one of the main sources of uncertainty in the derivation of the inner rotation curve is the rise in inclination of @xmath21  20  towards inner radii preferred by the tilted ring models .
it is not clear whether this rise is due to a genuine kinematic warp , or whether it is merely an artefact of the non - gaussian line profiles near the galaxy s centre .
we therefore derived two final tilted ring models , one with and one without the rise in inclination ( solid and dashed curves in panel d of fig .
[ tr_results ] respectively ) .
both models used the same @xmath93 , y@xmath113 , @xmath96 and @xmath97 radial profiles .
we hereafter refer to these two tilted ring models as models si ( steep inclination ) and ci ( constant inclination ) . in panel f of fig .
[ tr_results ] , the rotation curves of both models are presented . in each of the panels a  -  e in fig .
[ tr_results ] , the error bars represent the formal least squares errors from the rotcur fitting routine . in panel f ,
the error bars represent the r.m.s .
spread of the velocities in a single ring .
the final tilted ring models fitted to the third - order gauss - hermite velocity field are presented in fig .
[ tr_results ] . while the stellar disk seen in the 3.6  image ( fig .  [ irac_image ] ) extends as far as @xmath114 , the inner portion of the rotation curve rises as @xmath82 out to @xmath90 . beyond this radius
the rotation velocity remains almost constant . for model
si , the average rotation velocity for @xmath115 is @xmath116  @xmath73  1.6  . despite using different inclination profiles to derive the two final rotation curves , the maximum absolute difference between them is @xmath117  at @xmath118 .
the average absolute difference for @xmath119 is @xmath120 .
the resulting two rotation curves are thus very similar to one another .
for each tilted ring model we constructed separate rotation curves for the approaching and receding sides of the galaxy . in the case of model
si ( fig .
[ 2vrots ] ) , the rotation curve differs significantly between the two sides of the galaxy . within the steeply rising portion
, differences of @xmath21  20  are observed as well as at @xmath121 .
the same is true in the case of model  ci .
@xcite state that galaxies with asymmetric global profiles often have rotation curves that are more slowly rising on one side of the galaxy than the other .
this statement seems true in the case of ngc  2915 .    to check that the two tilted ring models are consistent with the data , we used the smoothed radial parameter profiles to project the rotation curves onto various position - velocity slices extracted from the  data cube .
the overlays for tilted ring model si are presented in fig  [ pv_overlays ] .
the projected rotation curves of the two tilted ring models are almost identical to one another with model si fitting well the high - intensity regions of the position - velocity slices . given a set of fixed radial profiles , the rotcur routine will adjust the rotation curve so that the _ line - of - sight _ velocities best match the data .
thus although the circular rotation curves ( i.e. panel f of fig .  [ tr_results ] ) differ , the good agreement between the projected rotation curves is not surprising .    using the smoothed , fixed radial parameter profiles of model si and ring widths of @xmath122 , a tilted ring model was also fitted to the third - order gauss - hermite velocity field extracted from the rw  cube .
the resulting rotation curve is compared to the rotation curve of model si in fig .  [ 10as_17as_vrots ] .
the figure demonstrates that for @xmath123 , the derived rotation curves are very similar to one another . beyond @xmath124
the filling factor of the rw  velocity field becomes too low for meaningful derivations of the rotation velocity .
for these reasons , only the rotation curves of tilted ring models ci and si are used as mass model inputs in sec .
[ mass_modeling ] .
dwarf galaxies such as ngc 2915 serve as useful probes of dark matter ( dm ) . in this section
we derive two mass models for ngc 2915 , one for each of the tilted ring models ci and si .
the gravitational potential within a galaxy is determined by the combined gravitational potentials of all the mass components .
the derived 3.6  surface brightness profile of ngc 2915 shows little deviation from an exponential light distribution at any radius .
we therefore concluded that the galaxy does not contain a significant central stellar bulge .
the total mass was treated as the sum of the masses of the stellar and gas disks as well as the dm halo . since @xmath125 at a given @xmath71 , the rotation curves of the individual mass components , when summed in quadrature , will yield the square of total rotation curve , @xmath126 : @xmath127 where @xmath128 , @xmath129 , and @xmath130 are the contributions to the total rotation curve of the gas , the stars and the dm respectively .
@xmath126 should match the observed rotation curve .
@xmath131 is the stellar mass - to - light ratio used to convert the observed distribution of star light to a stellar mass distribution while @xmath132 is the scaling factor used for all mass models to take into account the contribution of helium to the rotation curve .
we used the gipsy task rotmod @xcite to convert the observed mass distributions into rotation curves .
each of the fitted tilted ring models ( sec .
[ tr_output ] ) was used together with the  total intensity map to yield an  surface density profile .
these profiles were used as rotmod input which , assuming an infinitely thin gas disk , yielded two different  rotation curves .
ellipses fitted to 3.6  surface brightness isophotes suggest a constant inclination of @xmath133 for the old stellar population . using this inclination and a position angle of 306 , a single surface density profile for this mass component
was derived .
this surface density profile was extrapolated out to @xmath134 to match the radial extent of the observed rotation curve , and then converted into a stellar rotation curve assuming a @xmath135h@xmath136 vertical distribution @xcite and @xmath137 , yielding @xmath138  kpc .    in principle the dm density distribution of a suitably rotating system can be determined knowing the contributions to the observed rotation curve from the gas and the stars .
the  rotation curve is directly related to the observed  distribution , assuming that the  is optically thin .
the situation is not as simple for the stellar component of a galaxy , however , since a stellar mass - to - light ratio is needed to link the stellar light and mass distributions . determining the stellar mass - to - light ratio from the rotation curve alone is problematic . given an observed rotation curve , equally well - fitted mass models can be obtained for a range of stellar mass - to - light ratios ( e.g. * ? ? ?
* ; * ? ? ?
this degeneracy therefore leads to uncertainties in stellar mass - to - light ratios determined solely by rotation curves .
@xcite used spectroscopic and photometric spiral galaxy evolution models to show that a reasonable correlation exists between the optical and infrared stellar mass - to - light ratios and the colours of integrated stellar populations .
they further showed that the relative trends between model stellar mass - to - light ratios and colours are fairly insensitive to uncertainties in stellar population and galaxy evolution modeling .
@xcite used stellar population synthesis models ( together with a scaled salpeter 1955 imf ) to determine a relationship between the stellar mass - to - light ratios in the 3.6  and @xmath139 bands ( @xmath140 and @xmath141 respectively ) : @xmath142 from @xcite , the relation between the @xmath143 colours and the @xmath139-band mass - to - light ratio is @xmath144 where @xmath145 and @xmath146 for varying colours and metallicities are provided in their table 4 . by combining these two equations
we can estimate the stellar mass - to - light ratio at 3.6  .
we use @xmath147 and @xmath148 for @xmath143 colours as used by @xcite for ic 2574 , the justification being that ic  2574 is also a low metallicity , star - forming dwarf galaxy with @xmath149 , much like ngc 2915 .
@xmath150 for ngc  2915 @xcite thereby yielding ( from eqn .  [ k - colour ] ) @xmath151 . substituting this value into eqn .
[ emp_ml ] yields @xmath152 .
@xcite , applying the @xcite relation between @xmath143 colour and @xmath141 and assuming a @xcite imf , found @xmath153 for the @xmath154band stellar mass - to - light ratios of all of the things galaxies in their sample .
the way in which dark matter is distributed on galactic length - scales has been a topic of much debate . over the past decade , the constant increase in computational power has allowed detailed @xmath46-body simulations of the clustering properties of cold dark matter particles to be carried out .
one of the first large - scale studies was carried out by @xcite in which they studied the equilibrium density profiles of dark matter halos in a hierarchial galaxy clustering scenario .
their results suggested the existence of a `` universal dark matter density profile '' @xcite which is independent of clustering scale , mass or size as well as the power spectrum of initial fluctuations .
their parameterisation of this universal profile is @xmath155 where @xmath156 is a measure of the density of the universe at the time of collapse of the dm halo , @xmath157 is the characteristic scale radius and @xmath158 is the critical density required for closure
. this profile scales as @xmath159 , with @xmath160 and @xmath161 at large and small radii respectively , thereby predicting extremely steep inner density profiles of dm halos , known as central cusps .
the aquarius project @xcite recently completed a set of high - resolution numerical simulations to study the formation and structure of galaxy - sized dm halos in the standard @xmath162cdm cosmology .
@xcite reported that the spherically averaged density profile of the six different - sized simulated galaxy halos becomes progressively shallower inwards and , at the innermost resolved radius , has a logarithmic slope @xmath163 .
these predictions of cuspy inner dm density profiles can be tested against real , dark - matter - dominated galaxies .
if cuspy nfw halos are found , constraints can be put on halo concentrations and cosmological parameters .
the rotation curve resulting from an nfw density profile is @xmath164 where @xmath165 is the circular rotation speed at @xmath166 , the radius at which the density of the dm halo equals 200 times the critical density , @xmath167 ; @xmath168 is the radius in units of the virial radius and @xmath169 is the `` concentration '' parameter of the halo .
low surface brightness dwarf galaxies are thought to be dark - matter - dominated down to small radii @xcite and can therefore be more directly compared to the cold dark matter simulations . when this is done , it is often the case that the density profiles of simulated dm halos are too steep to fit the rotation curves of dwarf galaxies @xcite .
the distribution of the dm in these dwarfs is instead consistent with approximately constant inner dark matter densities .
a particular dark matter density profile that describes the dm distribution of dwarf systems well is that of the pseudo - isothermal sphere @xcite : @xmath170 where @xmath171 is the central dm density and @xmath172 is the core radius .
this parameterisation , which has no particular physical justification , has a constant density core .
the rotation curve resulting from such a density profile is @xmath173\right)^{1/2}.\ ] ]      we used the gipsy task rotmas to construct mass models for ngc  2915 .
the task subtracts from the observed rotation curve the scaled rotation curves of the stellar and gas disks .
it then fits to the residuals a rotation curve corresponding to a particular parameterisation of the dm halo .
when matching the total and observed rotation curves , the points were uniformly weighted to ensure that each point in the rotation curve contributed equally to the fit .
this essentially reduced the fitting procedure to a non - weighted least - squares fitting method .
as previously mentioned , each of the rotation curves from the two tilted ring models fitted to the  velocity field were used as mass modeling input . for each of these mass models we carried out fits for a pseudo - isothermal sphere and an nfw halo .
upper and lower limits for the stellar mass - to - light ratio were determined by scaling the contribution of the stellar rotation curve to the total rotation curve . under the so - called `` maximum disk assumption ''
, the stellar rotation curve is scaled up to contribute maximally to the observed rotation curve at inner radii .
such a scaling usually allows the stellar disk to explain most of the inner rotation curve @xcite . under the `` minimum disk assumption ''
, the contribution of the stellar disk is made as low as possible ( often zero ) while still allowing a good fit of the total rotation curve to the observed rotation curve .
each of the halo fits was divided into 3 sub - cases : a fit 1 ) with @xmath131 fixed to our predetermined value of 0.3 ; 2 ) using a maximum disk assumption ; and 3 ) using a minimum disk assumption . for the minimum disk cases ,
@xmath131 was fixed to zero ( no contribution from the stellar disk to the observed rotation curve ) . for the maximum disk cases ,
regardless of the dm halo parameterisation , @xmath131 was fixed to 0.64 and 0.44 for mass models si and ci respectively . in all models ,
@xmath174 was fixed to 1.37 .
thus , six mass models were fitted to each of the 2 rotation curves derived from the tilted ring modeling .
* comments on rows : * row  1 :  reduced @xmath175 goodness - of - fit statistic ; row  2 :  r.m.s . of the difference between the observed and total rotation curves ; row  3 :  stellar mass - to - light ratio ; row  4 :  central density for pseudo - isothermal sphere ; row  5 :  core radius for pseudo - isothermal sphere , row  6 :  circular rotation speed at virial radius for nfw halo ; row  7 :  concentration parameter for nfw halo ; row  8 :  dynamical mass ; row  9 :  to total mass ratio ; row  10 :  dark to total mass ratio ; row  11 :  dark to luminous mass ratio ; row  12 :  total mass to @xmath1-band light ratio .
+ * comments on columns : * for each model , columns 1 and 4 show the mass modelling results for the case in which @xmath131 is fixed to the predetermined value of 0.3 .
columns 2 and 5 show the results for the maximum disk case while columns 3 and 6 are for the minimum disk case .
the fitted parameters for all of the mass models are summarised in table [ mass_modeling_results ] .
the results for mass model si are also presented in fig .
[ rotmas_tr_a_plot ] .
this figure shows that ngc 2915 is dark - matter - dominated at nearly all radii ( depending on modeling assumptions ) .
we calculate mass ratios as well as the total mass to @xmath1-band light ratio at the last measured point of each fitted total rotation curve . for mass models si and ci , dm constitutes , on average , 92.8@xmath50 of the dynamical mass .
this dark - matter - dominance is confirmed by the total - mass - to - light ratio of ngc 2915 which is 141 /@xmath176 and 129 /@xmath176 on average for mass models si and ci respectively .
these estimates are almost double that of the 76  /@xmath176 upper limit set by @xcite thereby making ngc  2915 one of the most dark - matter - dominated late - type dwarf galaxies known .
the second main mass - modeling result that fig .
[ rotmas_tr_a_plot ] demonstrates is that , for either parameterisation of the dm halo , the steeply rising portion of the observed rotation curve is poorly matched by the total rotation curve .
it is usually the stellar rotation curve that is used to match the inner observed rotation curve @xcite , yet no such scaling is useful in the case of ngc 2915 due to the intrinsically different shapes of the stellar and observed rotation curves .
the stellar disk is clearly contained well within the steeply rising portion of the observed rotation curve which continues to rise as @xmath82 out to @xmath90 .
the best - fitting mass models are those in which the stellar rotation curve contributes zero to the observed total rotation curve .
figure [ rotmas_tr_a_plot ] shows that when the fitting routine attempts to fit the inner observed rotation curve with the dm rotation curve , the result is an over - estimation of the observed rotation velocities for @xmath177 kpc ( 110@xmath29 ) and an under - estimation for radii 2  kpc  @xmath178  5.5  kpc .
the poor matches between the total and observed inner rotation curves do not allow a particular dm halo parameterisation to be confidently ruled out or confirmed .
furthermore , the observed non - gaussian line profiles near the centre of ngc 2915 ( sec . [ line_profiles ] ) lead to uncertainties in the shape of the inner portion of the observed rotation curve .
these , in turn , lead to associated uncertainties in the shapes of the inner portions of the fitted dm halos . assuming that the stars are less susceptible than the gas to the effects of non - cricular motions , spectroscopic observations of the stars would have to be carried out to better constrain the inner rotation curve .
what is the explanation for the discrepancy between the observed and total rotation curves ?
have we failed to include a significant mass component during the decomposition process ?
a stellar bulge can not be the answer since our derived 3.6  surface brightness profile shows no significant deviation from an exponential light distribution .
a significant mass of molecular gas at inner radii is a possibility , yet co emission in low - mass galaxies is usually weak or not detected @xcite .
however , the conversion factor used to convert co luminosity to molecular gas mass is largely uncertain .
if we consider the mass - modeling results as they appear in table [ mass_modeling_results ] then the pseudo - isothermal sphere allows for the best match to the observed rotation curve .
the fitted pseudo - isothermal sphere , under the minimum disk assumption , has a core density and a core radius of @xmath179  @xmath180  pc@xmath5 and @xmath181  kpc for model si .
these results are similar to those of @xcite who , for their favoured mass model , estimated @xmath182  @xmath180  pc@xmath5 and @xmath183  kpc .
we also fitted the mass models with @xmath131 as a free parameter yet always found the best - fitting value to be negative , consistent with the mass modeling results of @xcite .    in fig .
[ thingsmm ] we compare our derived pseudo - isothermal halo parameters for ngc 2915 to those derived for the things sample from @xcite .
the infrared stellar mass - to - light ratios used in their mass models were derived from the 3.6  images in combination with stellar population synthesis arguments .
our models for the minimum disk and @xmath184 cases for both models si and ci are shown in fig . [ thingsmm ] .
this comparison suggests that ngc 2915 has a central dm core that is very compact and dense compared to other late - type systems .
cases in which the dm halo is modelled as a pseudo - isothermal sphere , compared to the mass modeling results of the things sample of @xcite .
the black - filled and open triangles correspond to the minimum disk results for models si and ci respectively while the black - filled and open squares correspond to the @xmath184 results for models si and ci respectively .
the mass modeling results from things are shown as open circles for dwarf galaxies and filled circles for all other galaxies . ]
we have obtained new deep , high - resolution  synthesis observations of ngc  2915 and have used them to carry a detailed study of the  distribution and dynamics of this galaxy .
the galaxy has very different optical and  properties .
optically , it is similar to other dwarf galaxies , consisting of a small stellar disk .
surrounding the stellar disk , however , is a huge  disk extending out to galactocentric radii greater than @xmath185 with well - defined spiral structure in its outer parts .
we estimated the total mass of the  disk to be @xmath186 .
our high resolution , high sensitivity observations clearly resolve the inner  distribution of ngc  2915 into two  concentrations separated by @xmath21 1.1 kpc .
this observation discredits the claim by previous investigators that the central gas disk of the system is made up of a massive  bar .
we have shown the shapes of  line profiles to be significantly non - gaussian in the vicinity of these central  over - densities , being split by @xmath21  30  on average over the spatial extent of the lower  concentration .
these profiles are indicative of non - circular velocity components within the gas near the centre of ngc  2915 . in our forthcoming paper ( elson et al . , 2010 , in prep .
) we investigate in further detail the central gas dynamics of this system .
linked to these non - circular velocity components , we have provided speculative evidence for the possibility of cold gas accretion from the inter - galactic medium onto the outer disk of ngc  2915 which , if confirmed , would have significant implications for the evolutionary history of this system .
these complicated dynamics lead to the new picture that ngc  2915 is not a simple , isolated , low - surface - brightness galaxy but rather an evolving system with a complex interplay between its various mass components .
the main focus of this paper has been the regular  dynamics of ngc  2915 .
we have fitted a tilted ring model to the  velocity field .
the best - fitting model was one in which the  disk is severely warped , with the inner and outer disks inclined at @xmath92 70 - 75  and @xmath187 55  respectively .
this result is , however , made uncertain by the non - gaussian line profiles near the centre of the galaxy .
we therefore also fitted a model with an almost constant inclination of @xmath92  55  for the entire  disk .
both models yielded a rotation curve typical of a late - type spiral with @xmath82 rotation out to @xmath90 and constant velocity thereafter .
rotation curves separately derived for the approaching and receding sides of the galaxy differed significantly within their steeply rising portions , providing the clear evidence that ngc  2915 is a kinematically lopsided galaxy .
our results show that these rotation curves can not be accounted for by the stellar or gas disks of ngc  2915 . by using the observed rotation curves as mass model inputs
, we found the galaxy to be dark - matter - dominated at nearly all radii with the stellar disk unable to account for the @xmath82 portion of the rotation curve .
the best - fitting mass model is one in which the stellar disk does not contribute to the observed rotation curve and in which the inferred dark matter halo is parameterised as a pseudo - isothermal sphere .
this model has a central core density of @xmath188  pc@xmath5 and a core radius of @xmath189  kpc .
the non - circular gas motions at inner radii lead to uncertainties in the fitted dark matter halo parameters . these , together with the fact that neither dm halo parameterisation allows for an accurate match of the observed rotation curve at inner radii , do not allow a particular dm halo parameterisation to be confidently ruled out or confirmed . further observations and modeling
are required to more accurately determine the inner rotation curve of ngc  2915 and to therefore more confidently discriminate between various parameterisations of the dm halo .
it is clear , however , that ngc  2915 has a very dense and compact dark matter core . at the last measured point on the rotation curve
, this galaxy has a total - mass - to - light ratio of @xmath190  140  /l@xmath16 , almost twice that of the upper limit placed my @xcite , thereby making it one of the most dark - matter - dominated galaxies known .
the work of ece is based upon research generously supported by the south african ska project .
ece would like to thank prof .
renzo sancisi for his enthusiasm and endless patience in discussing and helping to analyse the  data .
all authors acknowledge funding recieved from the south african national research foundation .
the work of wjgdb is based upon research supported by the south african research chairs initiative of the department of science and technology and the national research foundation .
the australian telescope compact array is part of the australia telescope which is funded by the commonwealth of australia for operation as a national facility managed by csiro .
this work is based [ in part ] on observations made with the _ spitzer _ space telescope , which is operated by the jet propulsion laboratory , california institute of technology under a contract with nasa .
finally , all authors thank the referee , claude carignan , for constructive comments that improved the quality of the paper .
r.  c. , bendo g. , engelbracht c. , gordon k. , li a. , rieke g.  h. , rieke m.  j. , smith j.  d. , armus l. , helou g. , jarrett t.  h. , roussel h. , calzetti d. , leitherer c. , malhotra s. , meyer m. , regan m.  w. , dale d.  a. , draine .
2003 , in bulletin of the american astronomical society vol .
35 of bulletin of the american astronomical society , sings : the sirtf nearby galaxies survey .
pp 1351+    b.  s. , staveley - smith l. , kilborn v.  a. , ryder s.  d. , kraan - korteweg r.  c. , ryan - weber e.  v. , ekers r.  d. , jerjen h. , henning p.  a. , putman m.  e. , zwaan m.  a. , de blok w.  j.  g. , calabretta m.  r. , 2004 , , 128 , 16                  g.  r. , freeman k.  c. , bland - hawthorn j. , jones m.  r. , knezek p.  m. , allen r.  j. , 1999 , in bulletin of the american astronomical society vol .  31 of bulletin of the american astronomical society , star formation in the dark disk of ngc 2915 .
pp 828+                      r.  j. , teuben p.  j. , wright m.  c.  h. , 1995 , in shaw r.  a. , payne h.  e. , hayes j.  j.  e. , eds , astronomical data analysis software and systems iv vol .
77 of astronomical society of the pacific conference series , a retrospective view of miriad .
pp 433+      m.  f. , cutri r.  m. , stiening r. , weinberg m.  d. , schneider s. , carpenter j.  m. , beichman c. , capps r. , chester t. , elias j. , huchra j. , liebert j. , lonsdale c. , monet d.  g. , price s. , 2006 , , 131 , 1163      v. , white s.  d.  m. , jenkins a. , frenk c.  s. , yoshida n. , gao l. , navarro j. , thacker r. , croton d. , helly j. , peacock j.  a. , cole s. , thomas p. , couchman h. , evrard a. , colberg j. , pearce f. , 2005 , , 435 , 629                  j.  m. , terlouw j.  p. , begeman k.  g. , zwitser w. , roelfsema p.  r. , 1992 , in worrall d.  m. , biemesderfer c. , barnes j. , eds , astronomical data analysis software and systems i vol .
25 of astronomical society of the pacific conference series , the groningen image processing system , gipsy .