two of the most enigmatic types of high redshift galaxy are the extremely red objects ( eros , eg . * ? ? ? * ; * ? ? ?
* ) , and the ly@xmath1 blobs ( eg . * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
eros are defined as having extremely red observed - frame optical / near - ir colors , and may be dusty starburst galaxies or agn ( eg . * ? ? ?
* ; * ? ? ?
* ) , or spheroidal galaxies ( eg . * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the ly@xmath1 blobs are characterized by 100 kpc scale low surface brightness ly@xmath1 nebulae , with total ly@xmath1 luminosities @xmath6 .
all ly@xmath1 blobs identified to date appear to lie in proto - cluster environments .
it has been proposed that these blobs are a form of agn , that they are associated with cooling flows @xcite or that they are generated by galaxy superwinds @xcite .
can a galaxy be both an ero and a ly@xmath1 blob ?
many high redshift radio galaxies contain ero nuclei embedded within dense gaseous regions , which are often associated with extended diffuse ly@xmath1 emission ( eg .
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
in this paper we present a detailed study of a _ radio quiet _
ly@xmath1 blob that contains two compact eros .
@xcite identified 2142@xmath04420 b1 , which is a luminous ly@xmath1 and emitting source at redshift @xmath7 .
the enormous ly@xmath1 luminosity of this source comes from an extended @xmath8 kpc low surface brightness nebula .
ground - based near - ir observations showed an unresolved , luminous , extremely red source lying within this blob .
the whole system appears to lie within a cluster of lyman - limit qso absorption - line systems @xcite .
the optical properties of this source closely resemble those of many high redshift radio galaxies , but @xcite obtained an upper limit on its radio flux of 0.27 mjy at 1344 and 2378 mhz .
we describe our observations in section  [ obs ] , and show the results in section  [ results ] .
we then discuss the nature of the ero ( section  [ continuum ] ) and of the ly@xmath1 blob surrounding it ( section  [ lines ] ) . in section  [ conclusions ]
we conclude that this object is probably a pair of young merging elliptical galaxies , though a compact dusty starburst model is hard to completely exclude .
we further conclude that the ly@xmath1 nebula surrounding them is either photoionized by an agn or shock excited by the energy of gravitational infall or of a superwind .
we assume a flat universe with @xmath9 , @xmath10 and @xmath11 . at @xmath7 , given this cosmology and redshift , one arcsecond corresponds to @xmath12 proper kpc , and the luminosity distance is @xmath13 .
our previous observations of this object can be found in @xcite and @xcite .
the observation log is in table  [ obstab ] .
hst imaging was used to determine the morphology and colors of b1 .
ground - based imaging measured its integrated colors and diffuse ly@xmath1 emission .
slitless spectroscopy and fabry - perot imaging was used to constrain the velocity structure .
in addition to our own observations , we searched the rosat all sky survey @xcite for data on this object .
there is only a 308 sec exposure toward the field , giving a 90% upper limit x - ray count rate of 0.0075 counts sec@xmath14 . assuming a galactic column density of @xmath15 @xmath16 along this line of sight ( which is the weighted mean of the four nearest 21 cm measurements ) and a power law spectrum with an energy index of 1.0 , we obtain a flux @xmath17 e@xmath18 photons @xmath16 sec@xmath14 kev@xmath14 where e is the energy in kev .    lccclc hst / wfpc2 & f410 m & @xmath19 & 67200 & 1999 aug 27sep 8 & + hst / wfpc2 & f450w & @xmath20 & 19400 & 1999 aug 31sep 2 & + hst / wfpc2 & f702w & @xmath21 & 14400 & 1999 aug 23 & + hst / nicmos & f110w & @xmath22 & 5120 & 1998 jun 27 & + hst / nicmos & f160w & @xmath23 & 10112 & 1998 oct 30 & + hst / nicmos & f164n & @xmath24 & 10236 & 1997 oct 30 & + hst / stis & clear & & 3120 & 1998 oct 2527 & + hst / stis & g430l & & 17000 & 1998 oct 2527 & + ntt / sofi & @xmath25 & @xmath26 & 3540 & 1999 aug 16 & +
ntt / sofi & @xmath27 & @xmath28 & 4740 & 1999 aug 16 & +
ntt / sofi & @xmath29 & @xmath30 & 3480 & 1999 aug 16 & + ctio / mosaic & 4107 & @xmath31 & 18000 & 1999 aug 7 & + ctio / mosaic & @xmath32 & @xmath33 & 9000 & 1999 aug 8 & + ctio / mosaic & @xmath34 & @xmath35 & 7200 & 1999 aug 78 & + ctio / mosaic & @xmath36 & @xmath37 & 7800 & 1999 aug 78 & + ctio / mosaic & @xmath38 & @xmath39 & 2400 & 1999 aug 8 & + ctio / mosaic & @xmath40 & @xmath41 & 4800 & 1999 aug 8 & + aat / ttf & 4110 & @xmath42 & 5400 & 1999 sep 11 & + aat / ttf & 4114 ? & @xmath43 & 5400 & 1999 sep 11 & + aat / ttf & 4118 & @xmath44 & 9000 & 1999 sep 12 & +      the cluster region was observed for 38 orbits with the wide field and planetary camera ii ( wfpc2 , * ? ? ?
* ) on the hubble space telescope .
galaxy b1 was placed on the wf3 chip .
exposures were dithered on a sub - pixel grid , to allow reconstruction of a better sampled image .
broad - band images were obtained through the f450w and f702w filters .
relatively narrow - band images were taken through the f410 m filter , which matches the wavelength of ly@xmath1 at the cluster wavelength .
the f410 m images suffer from a very low surface brightness , time - varying mottling .
we hypothesise that this is caused by scattered earth - light .
the images were pipeline processed , then cosmic rays were removed and the frames co - added using the drizzle algorithm @xcite to minimise undersampling .
the space telescope imaging spectrograph ( stis , * ? ? ?
* ) observations ( pairs of slitless spectra and images ) were taken at two different orientation angles separated by 11.5 deg , in an attempt to separate slitless spectra from objects aligned along the same row for a given angle .
for both the direct clear - aperture and the dispersed images , the ccd was binned @xmath45 , producing a plate scale of @xmath46 arcsec / pixel , and spectral resolution @xmath47   for point sources .
there were 20 exposures each for the direct and dispersed images , with dithering of @xmath48 arcsec between each exposure in order to minimize the effects of cosmic rays and hot pixels .
the direct and dispersed images were combined using the stis gto team software calstis @xcite .
the spectra were extracted using both the stis_extract @xcite and slwidget @xcite software .
results from both were consistent .      the near infrared camera and multi object spectrometer ( nicmos , * ? ? ? * ) observations were made with camera 3 , through the f110w , f160w and f164n filters .
all exposures were dithered for each filter over 4 ( or 8 in the case of f160w ) positions spanning a @xmath45 arcsec box , to minimize the effects of bad pixels .
the f160w and f164n observations are affected by vignetting from the field offset mirror , and were taken before the focus became optimal .
the nicmos data were reduced using the same methods as for sofi ( section  [ irim ] ) .
the data were also reduced using nicred @xcite ; there is little difference between final images reduced by the two methods .
a comparison of broadband sofi and nicmos data show the photometry to be consistent between the two , with the exception of a vignetted band along one end of the chip in the f160w and f164w images .      near - ir photometry was obtained with the sofi camera on the eso new technology telescope ( ntt , * ? ? ?
conditions were photometric , but seeing was @xmath49 1.1  and the telescope was subject to considerable wind shake .
the data were reduced using a modified version of iraf scripts written by peter mcgregor .
nicmos standards @xcite were used .
we obtained broad- and narrow - band ( ly@xmath1 ) imaging of this field with the mosaic camera @xcite , at the prime focus of the cerro tololo inter - american observatory ( ctio ) 4 m blanco telescope .
the narrow - band images were taken through a specially purchased filter , which working in the f/2.8 beam has a bandpass of 54 , centered at 4107 .
conditions were mostly photometric , with seeing @xmath50 .
our anglo - australian telescope imaging was obtained with the taurus tunable filter ( ttf , * ? ? ?
* ) , which is a fabry - perot etalon system .
conditions were photometric on the first night , and the typical image quality was @xmath51 .
a mit - lincoln labs ccd was used behind the ttf .
the ttf was used with a full - width at half maximum ( fwhm ) spectral resolution of between 3.35 and 4.4 .
images were obtained at central wavelengths of 4110 , 4114 , and 4118 .    an electronics fault
was subsequently found to induce significant wavelength drifts on timescales of a few hours .
the 4110 and 4118  observations were taken within four hours of a calibration , but the 4114  observations were taken between four and eight hours after the calibration .
thus the wavelength calibration of all three images are somewhat uncertain , particularly that of the 4114  image .
we obtained a long - slit spectrum of b1 with the low dispersion survey spectrograph ( ldss , * ? ? ? * ) .
these observations are described in @xcite
our images of b1 are shown in figures [ wfpc_frame ] & [ zoom_frame ] .
the appearance of b1 varies dramatically as a function of wavelength and resolution . in the observed - frame red and near - ir
, it breaks up into two compact components : b1a and b1b ( fig  [ zoom_frame ] ) , separated by 0.8 . in this respect
it is similar to another strong , high redshift ly@xmath1 source , the coup - foure galaxy @xcite .
a third continuum source , b1c , appears as chain of diffuse knots extending about 1  to the south from b1a and b1b .    in ly@xmath1 , the picture is different again .
when viewed with wfpc2 , we see a knot of ly@xmath1 emission from b1b , but no emission from b1a .
diffuse patches of low surface brightness ly@xmath1 emission are found @xmath52  to the south and north - east of b1b .
the ground - based ctio image has much greater sensitivity to low surface brightness ly@xmath1 emission , and confirms that it extends at least 10   north - east of b1a and b1b , as reported by @xcite and @xcite .
the non - detection of this emission in the wfpc2 image indicates that it is truly diffuse , and not coming from a series of compact sources separated by less than the ground - based resolution limit ( as was the case for the ly@xmath1 nebula observed by @xcite ) .    the ly@xmath1 contribution to the blue continuum images , and the blue continuum contribution to the narrow - band ly@xmath1 images , were removed .
b1 has a total ly@xmath1 flux of about @xmath53 : a more exact number is hard to calculate as the diffuse flux fades gradually into the sky ( this number is for a rectangular aperture 10  long and 5  wide , aligned along the major axis of the ly@xmath1 flux ) .
b1b has a ly@xmath1 flux of @xmath54 , which is only @xmath55% of the total .
the remaining flux is diffuse : the average ly@xmath1 surface brightness of the region within 2  of b1b is @xmath56 .
4further east , the average surface brightness drops to @xmath57 .
continuum fluxes ( table  [ magtab ] ) were calculated for the three components of b1 .
small circular apertures were used : 0.3  in radius for b1a and b1b , and 0.65  for b1c .
these small apertures will miss much of the flux of these components , but were chosen to avoid contamination from the other components .
no aperture corrections were attempted .
a total flux for b1 was also measured , using a 2radius aperture .
error estimates were made using apertures placed randomly on sky regions .
a copy of the image was made with all the sources detected by the sextractor @xcite package masked out .
we randomly placed 1000 apertures of identical radius on the image .
any aperture that fell on a masked region was rejected , and the procedure continued until 1000 clean measurements had been obtained .
the standard deviation of these measurements is our sky error estimate .
a poisson noise model was combined with the sky error estimate to calculate object flux errors .
the ly@xmath1 contribution to the f450w flux has been removed , but the probable emission - line contributions to f110w , f160w and @xmath29 have not been removed : these contributions are discussed in section  [ linesub ] .
finally , the fluxes were corrected for the local milky way extinction ( @xmath58 ) , as estimated by @xcite .
the extinction curve of @xcite was used .
integrated near - ir magnitudes for b1 were measured from the ntt image , using sextractor with 4.2  radius circular apertures .
we measure @xmath59 , @xmath60 and @xmath61 .
the measured coordinates of the components of b1 are shown in table  [ coords ] .
lc b1a & 21:42:27.45@xmath044:20:28.69 + b1b & 21:42:27.51@xmath044:20:28.99 + b1c & 21:42:27.46@xmath044:20:30.2 +    lcccc f450w ( @xmath34 ) & @xmath62 & @xmath63 & @xmath64 & @xmath65 + f702w ( @xmath38 ) & @xmath66 & @xmath67 & @xmath68 & @xmath69 + f110w ( @xmath25 ) & @xmath70 & @xmath71 & @xmath72 & @xmath73 + f160w ( @xmath27 ) & @xmath74 & @xmath75 & @xmath76 & @xmath77 +      the near - ir filters are subject to possible emission - line contamination at this redshift , z=2.38 .
[ ] ( 3727   ) is redshifted into the @xmath25-band , h@xmath78 and the [ ] ( 4959 and 5007   ) doublet are shifted into the @xmath27-band , and h@xmath1 and [ ] ( 6583   ) are shifted into the @xmath79 band .
these lines could in principle be very strong ( eg . * ? ? ?
* ) . in this section
, we estimate their strength , and conclude that it is probably small .    as described in @xcite , we imaged the field with a narrow - band filter of wavelength coverage @xmath80 m , which includes both h@xmath1 and [ ] . using a 5  radius aperture , we determine that these lines contribute @xmath81% of the integrated @xmath29 flux of b1 .
the f164n nicmos image covers the wavelength of redshifted h@xmath78 .
no significant narrow - band excess was detected : we can place a @xmath82 upper limit on the fraction of the integrated @xmath27-band light coming from h@xmath78 in b1 of 3% .
we have no direct measurements of the other line fluxes , so we estimate their importance by making assumptions about their ratios with measured balmer lines .
typical line ratios were taken from @xcite , @xcite and @xcite , for a variety of assumed chemical compositions and ionization sources , and including both narrow- and broad - line agn ratios .
we conclude that line emission probably contributes only @xmath83% at most of the total continuum flux of b1 in @xmath27 and @xmath29 . in @xmath25 ,
the contribution is probably @xmath84% .
furthermore , only about 10%  of the line flux will likely come from within the photometric aperture of b1b , and @xmath83%  from within the aperture of the b1c .
the remainder of the optical emission lines will come from the region of the diffuse blob .
as b1a shows no ly@xmath1 emission , its rest - frame optical line flux is also probably small .
the line corrections will thus probably be minimal .
we have not , therefore , applied any emission - line corrections to the near - ir fluxes .
the ly@xmath1 line is generally extremely optically thick and often strongly self absorbed , making it a very unreliable tracer of dynamics .
it is , however , the only sufficiently strong line available to us , and can be used to place some limits on the dynamics of the ly@xmath1 nebula surrounding b1 .
the velocity dispersion of the gas around b1 ( measured at the peak of the surface brightness , as seen in ground - based images ) is @xmath85 @xcite .
ly@xmath1 is redshifted with respect to , indicating probable self absorption .
the deep multislit spectroscopy of @xcite marginally spatially resolves this emission along the slitlet .
fig  [ longslit ] shows the spectra extracted from different ccd columns : notice the extra component at @xmath86 in the southernmost spectrum ( and another possible component at @xmath87 ) .
the diffuse ly@xmath1 structure differs significantly between the three ttf images .
this confirms the existence of multiple velocity components separated by @xmath88 .
it is hard to say more , given the uncertainty in the wavelength calibration of the ttf data ( section  [ ttf ] ) .
the stis spectrum detected the ly@xmath1 flux from b1b .
a flux of @xmath89 was measured , consistent with the value derived from the wfpc2 imaging .
the line was not significantly spectrally resolved ; an upper limit on the velocity dispersion of @xmath90 can be placed .
the spectral energy distributions ( seds ) of the three components of b1 are shown in fig  [ seds ] .
b1a and b1b are quite red , while the sed of b1c is very blue .
b1c s colors are similar to those of lyman - break galaxies ( eg . * ? ? ?
. a power - law fit to its colors gives a slope of the form @xmath91 , which is close to the median slope of lyman break galaxies ( 1.5 , * ? ? ?
our ntt photometry demonstrates that while b1 is red at wavelengths shortward of @xmath27 , it has a relatively blue @xmath92 color .
this confirms the result of @xcite .
b1 s sed thus peaks in the @xmath27-band . as we discussed in section  [ linesub
] , this is unlikely to be an artifact of emission - line contamination .
we modelled the seds of b1 using the 1997 version of the spectral synthesis models of @xcite .
models both with continuous uniform star formation and models in which all the stars were formed in an instantaneous burst ( a simple stellar population ) were used .
all models had salpeter stellar initial mass functions , with no mass cut - off .
dust was modelled using the empirical absorption curves of @xcite .
our predicted colors are compared with the observations in figures [ brhmodel ] and [ jhkmodel ] .
the red components , b1a and b1b , can be fit by two models :    * an unreddened burst model of age @xmath93 myr ( 1@xmath94 limits ) and stellar mass @xmath95 .
this naturally fits the @xmath27-band peak ( fig  [ jhkmodel ] ) and the photometry from @xmath96 thru @xmath29 .
it greatly underpredicts the observed f450w flux of b1a .
we do not consider this to be a problem : a star formation rate of only @xmath97 will produce ob stars sufficient to produce the observed blue flux . * a dusty starburst model .
the best fit is a model with a continuous star formation rate of @xmath98 , an extinction of @xmath99 mag , and an age of 500 myr
. a slightly worse fit ( but still acceptable at the @xmath100 level ) is younger ( 100 myr ) , dustier ( @xmath101 mag ) and has a star formation rate of @xmath102 .
the dusty starburst model is hard to reconcile with the strong ly@xmath1 emission of b1b , as ly@xmath1 is resonantly scattered and thus very strongly absorbed by even small amounts of dust .
b1c is much bluer , and can be fit either by an unreddened instantaneous burst model of age @xmath8 myr and mass @xmath103 , or by a slightly dusty continuous star formation model , with a star formation rate of @xmath104 and @xmath105 .      in this section ,
we show that b1a and b1b are spatially extended .
they are , however , very compact objects : too compact to be easily modelled as spiral galaxy disks .
they appear to show colour gradients : bluer in their central regions .    in figures  [ rpsf ] and [ jpsf ]
we compare the radial profiles of b1a and b1b with the relevant point spread functions ( psfs ) measured from a bright star in the same image . b1a and b1b are marginally resolved in the f702w ( rest - frame 2100  ) and f110w ( rest - frame 3250  ) images .
the f160w ( rest - frame 4730 ) image was taken with nicmos out of focus , and it does not significantly resolve b1a and b1b .
the observed psf does not vary significantly across the images , and is consistent with analytic predictions .
although b1a and b1b are clearly resolved , both are very compact , with the observed surface brightness dropping to half its peak value within about 0.1  ( @xmath52 kpc ) .
this places a strong constraint on the size of the galaxies .
consider a disk galaxy , with a typical radial surface brightness profile of the form @xmath106 @xcite .
the surface brightness of such a disk drops to half its peak value at @xmath107 ( fig  [ newfig ] ) .
thus our observations require that @xmath108 kpc .
only @xmath109% of modern galaxy disks are this compact @xcite . now consider a spheroidal galaxy , with a radial surface brightness profile of the form @xmath110\}\ ] ] @xcite .
this profile is much more sharply peaked than the exponential disk profile , and its surface brightness drops to half its peak value at @xmath111 ( fig  [ newfig ] ) .
thus @xmath112 kpc .
over 50% of modern elliptical galaxies have effective radii this small @xcite .
more detailed modelling confirms these results .
a wide variety of model galaxy profiles were convolved with psfs and fit by @xmath113 minimisation to the observed radial profiles .
the modelling confirms that b1a and b1b are significantly extended .
both exponential disk models and spheroidal ( de vaucouleurs ) models give acceptable fits to the data , as do models with a point source embedded in a fainter halo .
pure exponential disk models only give acceptable fits if @xmath114 ( @xmath115 pc ) , which would place b1a and b1b in the most compact 1% of present - day disk galaxies .
spheroidal models give best fits for effective ( half - light ) radii @xmath116 of around @xmath117  ( @xmath118 kpc ) , which are typical of modern elliptical galaxies .
models which combine a nuclear spheroid or point source with a disk can give acceptable fits , but only if the disk is unusually small and/or contributes only a small fraction ( @xmath119 10% ) of the light .
the compact radial profiles of b1a and b1b are typical of those of high redshift lyman - break galaxies @xcite .
the inferred surface brightnesses in the rest - frame uv are comparable to those of faint lyman - break galaxies , but the rest - frame optical surface brightnesses are an order of magnitude higher . as @xcite note ,
these compact radial profiles are consistent with present - day elliptical galaxy profiles or bulges , but are much more compact than present - day disks .
curiously , both b1a and b1b appear to be significantly larger at longer wavelengths .
this can be seen in figures  [ rpsf ] and [ jpsf ] : despite the wider psf for the f110w filter , both components are clearly more extended than as seen through the f702w filter .
our modelling confirms this result : even after convolution with the relevant psfs , no single model gives a good fit to either b1a or b1b at both f702w and f110w .
the central @xmath52 kpc of both components is @xmath120 magnitudes bluer in f702w@xmath0f110w than the outer regions .
the change in size as a function of wavelength could be explained by a colour gradient within a single spheroidal component .
it could also be explained by a two component model : a compact blue component that dominates at rest - frame 2100  and a more diffuse ( but still compact ) red component that dominates at 3250  and beyond . this red component , if disk - like , must have @xmath121 , which would still place it in the most compact 1% of low redshift disks .
what are the continuum components of b1 ?
b1c , with its blue colors , appears relatively straightforward : it is a region of extended moderate star formation , obscured by little dust .
this star formation may be triggered by the interaction of b1a and b1b .
the nature of b1a and b1b is less straightforward .
we consider four models in turn .      could b1a and b1b be young elliptical galaxies ?
this hypothesis fits the data well .    *
the radial profiles of b1a and b1b are well fit with de vaucouleurs profiles , and the inferred radii are typical of present - day elliptical galaxies . * the red colors are most easily explained by a @xmath122 stellar population which completed its star formation about 750 myr before we observe it
. it would thus be similar to modern e@xmath123a galaxies . * b1 lies in an overdense region of the early universe @xcite , which may be the ancestor of a galaxy cluster .
* the velocity field around b1 , if virial , implies masses of @xmath124 .    if b1a and b1b are elliptical galaxies , how do their surface brightnesses and sizes fit on the kormendy relation @xcite ?
the f160w filter corresponds quite closely to the rest - frame @xmath34-band .
the radial profile fitting suggests that both components have effective radii in the observed - frame near - ir of @xmath118kpc . integrating a de vaucouleurs profile over our photometric aperture
, we can convert our observed @xmath125 magnitudes into rest - frame @xmath34-band surface brightnesses at the effective radius .
were these components at low redshift , their surface brightnesses at the effective radius would be @xmath126 .
this is roughly an order of magnitude higher than the kormendy relation would predict for an elliptical galaxy of this radius .
this is consistent with passive evolution from our unreddened burst model ( section  [ sed ] ) .    if b1a and b1b are elliptical galaxies , why do they apparently show color gradients ?
the uv emission from the central @xmath52 kpc could be caused by a starburst , a hole in the dust or ( for b1b ) an agn .
alternatively , some chemodynamical models predict color gradients in young elliptical galaxies ( eg . * ? ? ?
* ; * ? ? ?
a dusty agn can certainly have colors as red as b1a and b1b ( eg . * ? ? ?
the red components of b1a and b1b are however spatially extended , and hence can not be produced by an agn .
could the observed color gradient , however , be caused by the superposition of a compact blue agn on a more extended , lower surface brightness red galaxy ?    even in the relatively blue central regions , the continuum slope is redder than @xmath127 , while radio - quiet qsos have typical continuum slopes of @xmath128 @xcite .
any central agn would thus have to be reddened by dust with @xmath129 .
this would make it hard for the ly@xmath1 emission from b1b to escape .
in addition , the agn in b1a and b1b would need remarkably equal luminosities and dust extinctions to explain the very similar observed colors and radial profiles of the two components .
this fine tuning seems implausible .
furthermore , the high measured surface brightness was measured in the f160w band , to which the central blue component does not significantly contribute .
we therefore conclude that while agn may be present , they are unlikely to be responsible for the observed colors .
sensitive hard x - ray or far - ir observations will be required to determine whether a dusty agn is present .
the more extended red component of b1a and b1b can be fit by a very compact exponential disk .
could these two objects then be compact spiral galaxies ?
the more compact bluer central component could be caused by agns or nuclear starbursts .    with sufficient dust
, a spiral galaxy can certainly appear as red as these two components ( section  [ sedmodel ] ) .
two lines of reasoning however oppose this hypothesis .    1 .
even the extended red components are extremely compact by the standards of modern disk galaxies . this could , however , be explained if disk galaxies form from the inside outwards .
this model requires that the outer ( red ) components of b1a and b1b be greatly reddened .
this would make it difficult for the observed uv light and ly@xmath1 emission from the bluer central components to escape
.      could b1a and b1b be the bulges of two disk galaxies ?
the disks themselves may not have formed , or might be too low in surface brightness to be detected by the hubble space telescope .
two lines of reasoning oppose this hypothesis .
firstly , the red colors and @xmath27-band peak of b1 are seen in our ground - based photometry with large ( @xmath130 ) apertures .
thus any extended component must have a similar stellar population to that inferred here .
secondly , the inferred masses and sizes for b1a and b1b are far larger than those typical of the bulges of disk galaxies .
b1 has a total ly@xmath1 luminosity of @xmath131 , spread over a region at least @xmath132 kpc in size .
the emission consists of a number of diffuse , discrete components with relative velocities of @xmath88 .
what can produce such a luminous , diffuse , fast moving ly@xmath1 nebula ?
the bulk motions of the ly@xmath1 nebula suggest that fast shocks must be present .
we will show that such shocks are quite capable of producing the observed ly@xmath1 luminosity .
the puzzle then becomes one of explaining the origins of the bulk motions .
it has also been suggested that photoionzation by an agn , or a cooling flow may power such nebulae .
we discuss these possibilities , which are hard to exclude .
if we assume that the ly@xmath1 velocities are representative of the gas , and not an artifact of the high optical depth in ly@xmath1 , then the dynamical timescale ( crossing time ) of the nebula is @xmath133 years .
this is interestingly similar to the inferred age of the stellar population ( section  [ sed ] ) .    what is the density of the nebula ?
the sight - line to background qso 2139@xmath04434 passes 20  from b1 ( though it does not pass through the ly@xmath1 nebula ) .
the qso spectrum shows a strong absorption - line system at the redshift of b1 , with a column density @xmath134 @xcite .
if we assume that the column density through the ly@xmath1 nebula is at least as great as this , and that the nebula is @xmath83 kpc thick ( ie . as thick as it is wide ) along the line of sight ( both big assumptions ) , then its density must be at least @xmath135 . if , however , the density were this low , the gas would be fully ionized by the metagalactic uv background at this redshift .
if the nebula is @xmath83 kpc thick , it must have a density of @xmath136 for its recombination rate to balance the photoionization from the uv background @xcite .
the diffuse star formation of b1c , taking place within part of the ly@xmath1 nebula , also suggests that densities may be higher . in reality , of course , the cloud may be highly inhomogeneous , with dense neutral ` bullets ' embedded in low density ionized gas .
if we take this lower limit on the density of the nebula ( @xmath137 ) , and assume that the nebula is 100 kpc long , 30 kpc wide and 10 kpc thick , then its hydrogen mass is @xmath138 .
whenever there are blobs of gas moving at supersonic relative speeds , as are observed in the nebula , shocks are inevitable .
we thus know that shocks are present .
could the emission from these shocks produce the observed ly@xmath1 emission ?
we used the fast shock models computed with the mappings  iii code to calculate the line emission from the shocks .
mappings  iii is an updated version of the mappings  ii code @xcite and uses improved radiation continuum and hydrogenic calculations to calculate the line and continuum radiation from shocks .
the new version of the code enabled the calculation of a grid of shocks for a range of metallicities , up to 1000 km / s , although a small number of low metallicity models above @xmath139 failed to converge after several global iterations .
a fiducial model with a velocity of @xmath140 was used here , having achieved a converged solution over a range of metallicities .    if the shocks are radiative , and if we assume , for example , a density of @xmath141 , then the ly@xmath1 emissivity will be @xmath142 , almost independent of the metallicity of the gas .
a quarter of this comes from the shock itself , the rest coming from photoionization of the assumed neutral gas in front of the shock ( the precursor region ) .
if the shock is face - on to our sight - line , this would correspond to an observed surface brightness of @xmath143 : nearly two orders of magnitude greater than the observed ly@xmath1 surface brightness ( section  [ results ] ) .
much of this may be self absorbed , but if even 1% escapes , it can explain our observations .
the predicted ly@xmath1/ ratio is also consistent with the observations @xcite .
how does this vary with the assumed density of the gas ? if the shocks are radiative , the emissivity will scale roughly with the density .
thus even with our lower limit on the density @xmath144 ( section  [ lines ] ) , shocks can produce all the ly@xmath1 we observe .    would these shocks be radiative ? for our lower limit on the density ,
@xmath145 shocks become radiative on a timescale of @xmath133 years , which is comparable to the dynamical timescale .
thus it seems likely that at least some fraction of the shocks will be radiative .
thus shocks are probably present , and can easily produce the ly@xmath1 emission we observe . a number of factors could , however , suppress the ly@xmath1 emission from shocks .    1 .
the high optical depth in ly@xmath1 may have caused us to overestimate the velocity dispersion of the ly@xmath1 nebula .
the emission from shocks is a strong function of their velocity : @xmath146 shocks emit two orders of magnitude less ly@xmath1 per unit area than our fiducial model .
2 .   the density of the ly@xmath1 nebula may be close to our lower limit of @xmath144 , in which case the ly@xmath1 emission , after the probable self absorption losses , could be an order of magnitude below the observed ly@xmath1 surface brightness .
3 .   the ly@xmath1 nebula could consist of small dense ` bullets ' of neutral gas , moving at high speeds through a hot , highly ionized , low density medium .
the bowshocks in the ionized medium would not be radiative , while the reverse shocks within the ` bullets ' would be slow , and hence have low surface brightnesses ( as well as having a small area ) .
if the shock model is correct , the bulk gas motions explain the ly@xmath1 luminosity .
what then could be the energy source for the bulk gas motions ?
if we divide the total kinetic energy of the nebula ( assuming a mass of @xmath147 ) by the ly@xmath1 luminosity , we can derive an upper limit on the damping time , which is @xmath148 years .
this is less than the dynamical timescale and stellar age of b1 , suggesting that continued energy input is required .
note that the ly@xmath1 we observe may only be a small fraction of the ly@xmath1 emitted , which in turn will only be a small fraction ( @xmath52% ) of the total energy dissipation in the shocks . on the other hand
, the gas mass may be greater than we are assuming .
if continued energy input is required , what can the source be ?
perhaps the most natural energy source is the gravitational potential energy of b1 .
if we assume that b1 has a total mass of @xmath149 , then the gravitational potential energy released in its formation would be @xmath150 ergs , which is more than sufficient to power the ly@xmath1 luminosity for the dynamical timescale .
such a mass for b1 would imply virial velocities in the ly@xmath1 nebula of @xmath151 , which are comparable to the bulk gas motions observed .
this energy could be released by the merger of b1a and b1b , or by a continuing mass infall rate of @xmath152 .
@xcite point out that chemically primordial gas at a temperature of @xmath153 k has few cooling mechanisms other than ly@xmath1 emission .
thus if most of the gas around b1 is primordial , a ly@xmath1 luminosity comparable to the gravitational binding energy ( section  [ energy ] ) of the system must be radiated .
alternatively , @xcite suggest that galactic superwinds driven by a starburst could power extended ly@xmath1 nebulae .
this would naturally explain the similarity between the dynamical timescale and the inferred age of the stellar population in b1a and b1b .
the maximum plausible amount of energy which star formation can inject into the intergalactic medium is @xmath154 ( eg .
as b1a and b1b have a combined stellar mass of @xmath155 , the total energy liberated would be @xmath156 ergs , which could in principle drive the observed ly@xmath1 luminosity for @xmath157 years .
the ly@xmath1 emission in high redshift radio galaxies has sometimes been ascribed to mechanical energy deposition by a radio jet ( eg . * ? ? ?
could such a jet be present in b1 ?
we do not detect any radio emission from b1 @xcite .
our @xmath82 upper limit of 0.23 mjy at 2.4 ghz corresponds to a luminosity limit of @xmath158 at rest - frame 1.4 ghz ( assuming a typical radio galaxy radio spectrum of the form @xmath159 ) .
typical ratios of the total jet mechanical power to the monochromatic radio flux at 1.4 ghz are @xmath160@xmath161 @xcite .
thus we can place an upper limit of @xmath162 on the energy injection from a radio jet .
this is sufficient to power the observed ly@xmath1 luminosity only if @xmath163% of the mechanical energy is converted into observable ly@xmath1 emission , which seems implausible .
deeper radio observations should test this hypothesis .      if the shock emission is suppressed , is there an alternative mechanism capable of producing the ly@xmath1 emission ?
stellar ionisation ( by hot young stars ) is incapable of producing the observed emission .
photoionisation by the uv continuum radiation from a qso can , however , produce both the ly@xmath1 and , as discussed by @xcite .
we do not see strong uv continuum emission from any component of b1 , so such an agn would have to be concealed from our sight - line . @xcite
detected strong sub - mm emission from a ly@xmath1 blob at @xmath164 , suggesting that a dusty agn or starburst was present , but we have no evidence for such a source in b1 .
extended emission - line nebulosities are associated with many low redshift cooling flow clusters , such as ngc 1275 .
as all high - redshift ly@xmath1 blobs seem to lie in cluster environments , could the ly@xmath1 luminosity be associated with a cooling flow ( eg . * ? ? ?
@xcite suggest that ly@xmath1 nebulae are caused by gas cooling from higher temperatures , in such a cooling flow . if each cooling hydrogen atom produces one ly@xmath1 photon , however , we require that @xmath165 of hydrogen be cooling onto b1
. this would build the observed baryonic mass of b1 in @xmath166 years .
this is far smaller than either the dynamical time or the age of the stellar population in b1a and b1b , so the simple cooling model seems implausible .
note that the nebular emission - line luminosities of low redshift cooling - flow clusters are also orders of magnitude too great to be explained by this simple recombination mechanism ( eg . * ? ? ?
the cause of this discrepancy is still controversial even at low redshifts .    on purely energetic grounds , could the ly@xmath1 luminosity be powered by the thermal energy in the nebula ?
if we divide the thermal energy of @xmath167 of hydrogen at @xmath168 k by the observed ly@xmath1 luminosity , we get a timescale of @xmath148 years : shorter than the dynamical timescale or the stellar age .
@xmath167 may however be an underestimate of the true gaseous mass of the ly@xmath1 nebula , particularly if the neutral gas we observe is immersed in a much more massive cloud of hot , highly ionised x - ray emitting gas .
so what is b1 ? on balance , we conclude that b1a and b1b are probably young elliptical galaxies , perhaps analogues of e@xmath123a galaxies today .
both their red colors and compact sizes are naturally explained by this model .
an extremely dusty starburst in the central regions of a disk galaxy can however , with some fine tuning , reproduce the colors and sizes .
the bluer colors of the centres of both components , and the ly@xmath1 emission from the core of b1b are , however , arguments against the presence of large amounts of dust in these systems .
near - ir h - band spectroscopy on an 8m - class telescope should resolve this ambiguity , by detecting the spectral signatures of a balmer/4000break .
b1c appears to be a region of relatively unobscured star formation , perhaps induced by the interaction of b1a and b1b .
if b1a and b1b are elliptical galaxies , their surface brightnesses are an order of magnitude greater than those that the kormendy relation would predict for galaxies of this size , and they show strong radial color gradients , with blue cores . the very existence of massive elliptical galaxies at this redshift , if confirmed , is interesting ( eg .
* and references therein ) .
these galaxies would give us a chance to study young elliptical galaxies for the first time .
our observations would pose several puzzles :    * why is the surface brightness of these galaxies so great ? is this simply due to passive evolution ? * why do both galaxies show strong , but very similar color gradients ?
* is the remarkable similarity and proximity of the two galaxies coincidental , or is it telling us something about how cluster elliptical galaxies form ? * if these galaxies are so massive and red at @xmath169 , what would they have looked like at @xmath170 when they were forming ?
the nature of the ly@xmath1 nebula is less well constrained by our observations .
the nebula could be excited by photoionization from an agn , by whatever powers low redshift cooling flow nebulae , or by fast shocks powered by either superwinds or gravitational potential energy .
spatially resolved spectroscopy of other emission lines will help discriminate between these models .
we wish to thank joss bland - hawthorn and catherine drake for their assistance at with the ttf observations , and all the staff at ctio for heroic efforts to get the mosaic camera working for us
. we are grateful to mike dopita , megan donahue and martijn de kool and ian george for helpful discussions .
hsiao - wen chen was instrumental in getting the slitless stis spectrum extraction software working .
david hogg s cosmology crib sheet ( astro - ph/9905116 ) was very helpful to us .