the importance of atomic diffusion @xcite inside stars is now well established : not only can it modify the atmospheric abundances , as observed in the so - called `` chemically peculiar stars '' , but it can also have strong implications for the stellar internal structure . the variations with depth of the radiative accelerations on specific elements can lead to their accumulation or depletion in various layers inside the stars . in this framework , the accumulation of iron ( together with less important metals like nickel ) may have strong consequences , which are different according to the stellar type .    in main sequence a stars
, it may lead to the creation of an extra convective zone due to an iron - induced increase of opacity @xcite .
if overshooting leads to complete mixing between this extra convective zone and the one due to helium , helium settling may be slowed down in the outer stellar layers .
such an effect could explain the existence of evolved oscillating am stars @xcite .
the occurrence of an iron convective zone may also help explaining the excitation of g - modes in @xmath2 doradus stars , through the `` convective blocking '' process , as suggested by @xcite .
furthermore , in some cases iron accumulation can trigger stellar pulsations through the iron - induced @xmath0-mechanism .
this process has been invoked to explain main sequence b - type pulsators like spb and @xmath3 cephei stars @xcite .
@xcite predicted that the extreme horizontal branch stars referred to as sdb stars could also be destabilized due to iron - induced @xmath0-mechanism .
the oscillations of these stars were independently observed ( e.g. * ? ? ? * ) .
the qualitative success of the iron accumulation theory is impressive .
however when quantitatively comparing models and observations , some difficulties remain .
we will further discuss these difficulties in section [ discussion ] .    in previous stellar modelling done in this framework ,
the question of the stability of a heavy iron layer lying above layers with smaller mean molecular weights was not addressed .
this represents a serious problem , because of the existence of local inverse @xmath1-gradients , which are highly unstable against thermohaline ( double diffusive ) convection ( see * ? ? ?
the induced mixing has to be taken into account in the computations , otherwise the results which are obtained do not correctly represent the real situation found inside the stars .
the importance of thermohaline convection has been invoked in a few astrophysical cases .
the first case which was pointed out concerned the @xmath4he burning regions in giant stars @xcite .
situations with @xmath5he enhancement in outer stellar layers were also discussed in several studies .
this enhancement could be due to mass transfer @xcite or to helium diffusion in a stellar wind @xcite .
the observations of main sequence helium - rich stars , with effective temperatures around 20  000k , show that helium is enhanced in average by a factor of two .
this demonstrates that thermohaline mixing is important , otherwise the helium overabundance would be much larger .
however , thermohaline convection mixing can not be extremely efficient since this would lead to no helium overabundance . in this case , the reason for this lack of efficiency may be due to the presence of a magnetic field .
the same kind of processes may occur inside roap stars , as discussed by @xcite .
thermohaline convection has also been studied for the case of accretion of metal - rich matter onto main sequence stars .
this situation may occur for exoplanet host stars , if they accrete hydrogen poor matter in the process of planetary formation .
because of thermohaline convection , and contrarily to preliminary conclusions found in early papers @xcite , the accreted matter does not remain inside the convective zone , but diffuses downwards @xcite . another case concerns carbon enhanced main - sequence stars ( cemps ) , which are assumed to have suffered some accretion of material coming from an agb companion .
this heavy matter falls down inside the star due to thermohaline convection @xcite .
owing to the stabilizing @xmath1-gradient induced by helium diffusion , some of it can remain in the outer layers and may account for the observed abundances @xcite .    in this paper , we first give a general discussion of thermohaline convection in stars ( section [ thermoconv ] ) and show how this physical process can be introduced in the stellar evolution computations ( section [ computations ] ) .
then results are presented for heavy elements accumulation inside a - type stars .
we mainly concentrate on the evolution of a 1.7 @xmath6 star , but also give results for 1.5 @xmath6 and 1.9 @xmath6 stars ( section [ results ] ) . for these computations , the tgec code ( cf . * ? ? ?
* ; * ? ? ?
* ) was modified to include the computations of radiative accelerations using the svp ( standing for single valued parameters ) approximation proposed by @xcite .
this modified version of the tgec code was first used to compute stellar evolution models ( without introducing thermohaline convection @xcite ) which were compared to the results obtained by the montreal group .
these computations were repeated while including thermohaline convection below the iron accumulation layer . due to the stabilizing effect of helium settling , which fortunately occurs in the same stellar region ,
iron can still accumulate , but the abundance profiles are strongly modified .
we also discuss the potential important consequences of this process for other types of stars in section [ discussion ] .
a short conclusion will then follow .
thermohaline convection is a well - known process in oceanography : salted water layers lying above fresh water ones are rapidly mixed downwards even in the presence of stabilizing temperature gradients , due to the different diffusivities of heat and salt .
when a warm salted blob falls down in cool fresh water , the heat diffuses out more quickly than the salt .
the blob goes on falling due to its weight until it mixes with the surroundings .
this leads to the formation of so - called `` salt fingers '' .
a similar kind of convection occurs in stellar radiative zones in the presence of positive @xmath1-gradients , that is @xmath1-values increasing upwards .
an``iron layer '' lying above hydrogen rich medium will lead to `` iron fingers '' , because when a blob falls down heat diffuses between the blob and its surroundings more quickly than iron ions .
with reference to oceanography , this process is called `` thermohaline '' convection .
it is also sometimes referred to as double diffusive convection , because it is directly related to the presence of two `` scalars '' which diffuse on very different time scales .
another kind of double diffusive convection occurs in the reversed case , in the presence of destabilizing temperature gradients and stabilizing @xmath1-gradients like , for example , at the edge of convective cores .
this process , quite different from thermohaline convection , is generally referred to as `` semi convection '' .
it will not be considered here .    in the case of positive @xmath1-gradients ,
the medium can become dynamically unstable if ( ledoux criterion ) : @xmath7 where @xmath8 ln @xmath9 ln @xmath10 and @xmath11 ln @xmath9 ln @xmath12 .
when @xmath13 vanishes , marginal stability is achieved . for positive or null values of @xmath13 ,
the medium is stable against dynamical convection .
then thermohaline convection occurs on a time scale which is typically of a few thousand years , long compared to the dynamical time scales but short compared to stellar lifetimes .
fingers may form if the following condition is verified : @xmath14 with @xmath15 = @xmath16 where @xmath17 and @xmath18 are the thermal and molecular diffusion coefficients while @xmath19 and @xmath20 are the corresponding time scales .    in stars
the value of the @xmath21 ratio is typically on the order of @xmath22 if one assumes that @xmath18 is the `` atomic '' diffusion coefficient .
however , it can become larger when the shear flow instabilities which induce mixing between the edges of the fingers and the surroundings are taken into account ( see * ? ? ?
the effects of thermohaline convection as a mixing process in stars are far from trivial .
many detailed studies related to the occurence of this physical process in water have been published .
for example , @xcite gave precise comparisons between numerical simulations and laboratory experiments .
however the stellar case is different since mixing occurs in a compressible fluid .
two different parametrisation recipes , given respectively by @xcite and @xcite can differ by up to two orders of magnitude : this illustrates that treating thermohaline convection as a simple diffusion process may lead to inaccurate results .
the basic problem here concerns the vertical shear flow instability which occurs between the fingers and the inter - fingers medium .
this instability leads to local turbulence which increases the mixing at the edge of the fingers .
consequently , a process of self destruction appears for the blobs , so that fingers eventually reach a regime where they can not form anymore : this effect is taken into account in the kippenhahn et al .
procedure .    in the present paper , we take into account the fact that the thermohaline mixing time scales are short compared to the main sequence evolution time scales , and that mixing stops when the positive @xmath1-gradients vanish . in our computations ( see discussion below ) , we model rapid mixing so that the stellar medium always adjusts to keep either flat @xmath1-profiles or @xmath1-values increasing downwards . as helium settles downwards , it creates a stabilizing @xmath1-gradient which in turn may allow some heavy elements accumulation .
we computed stellar models to describe the evolution of a - type stars while using the toulouse - geneva evolution code ( tgec ) .
the code is described in detail in @xcite and @xcite , except for the computations of opacities and atomic diffusion which have undergone major improvements @xcite .      in the present computations ,
important improvements were introduced in the tgec code concerning the opacities and radiative accelerations .
the opacities were computed using the opcd v3.3 codes and data @xcite .
it allows to compute self consistent rosseland opacities taking into account the detailed composition of the chemical mixture .
the opacities were recalculated at each time step .    to more realistically treat atomic diffusion in the tgec code ,
the radiative accelerations on c , n , o , ca and fe were included following the improved version @xcite of the semi - analytical prescription proposed by @xcite .
this method allows for very fast computations of radiative accelerations with a reasonable accuracy .
the svp approximation may be implemented in existing codes in a simple way .
there are much less data to process than for detailed radiative acceleration calculation , because complete and detailed monochromatic opacites for each ion are not needed .
a new grid of svp - parameters , well fitted to the stellar mass range considered in this work , was computed following the procedure described in @xcite .    the diffusion of h
, he and the five metals mentioned above were included in the stellar models .
the abundance of the other elements were assumed to be constant throughout the star .
the radiative accelerations of these five heavy elements as well as the detailed abundances for all seven elements were computed at each diffusion time step , chosen to be about 50 times smaller than the evolution time step .      for our modelling
, the convective regions were instantaneously homogenised .
the hi and heii convective regions were supposed close enough to be connected by overshooting and mixed together . on the other hand , the iron convective region which may appear in some models in much deeper regions ,
was supposed disconnected from the surface convective zone .
the models were evolved from pre - main sequence up to hydrogen core exhaustion .
they were assumed homogeneous on the pre - main sequence and atomic diffusion was introduced at the beginning of the main sequence .
we first computed models without thermohaline convection . to avoid the appearance of steep and unrealistic abundance gradients at the transition between radiative and convective regions , we introduced mild mixing at the bottom of each convective zone .
this mixing was modeled as a diffusion process through a coefficient d@xmath23 introduced in the chemical transport equation : @xmath24 @xmath25 represents the mean concentration of a species i , @xmath26 is the density , r is the radius , d@xmath23 is given by : @xmath27 d@xmath28 and r@xmath28 are respectively the value of d@xmath23 and the value of the radius at the base of the convective zone , @xmath29 is the half width of the mixing region .
d@xmath28 and @xmath29 are free parameters , choosen to produce a mild mixing on a small extent .
the value of d@xmath28 is taken equal to @xmath30@xmath31.s@xmath32 and @xmath29 is fixed to 0.5% of the stellar radius below the surface convective region and to 0.05% of the stellar radius below the iron convective zone .      as discussed in the previous sections ,
thermohaline convection arises in regions where the mean molecular weight increases towards the surface .
it leads to mixing on short time scales , compared to those of evolution , and it disappears when the mean molecular weight gradient vanishes . to simulate this phenomenon , we found that the most efficient and simple way was to model the mixing as usual , provided that the chosen diffusion coefficient is large enough to flatten the @xmath1-profile any time a destabilizing @xmath1-gradient is created ( at each time step ) .
the diffusion coefficient proposed by @xcite is convenient in this framework : @xmath33 where @xmath34 is the pressure scale height , a the radiative pressure , c the speed of light , t the temperature , @xmath35 the adiabatic gradient , @xmath36 the real temperature gradient , @xmath37 the specific heat at constant pressure , @xmath0 the opacity and @xmath1 the mean molecular weight of the material .
we used this prescription in our computations and checked that the @xmath1-gradient remained always null or negative at each time step in all our models .
we computed models for 3 different masses : 1.5 , 1.7 and 1.9m@xmath38 .
each set of models includes a first series of computations without thermohaline convection and a second series including it .
the initial metal mixture used in these computations is the solar mixture presented by @xcite with the following values as h and he initial mass fractions : @xmath39 and @xmath40 .
the initial iron mass fraction was consequently @xmath41 .
figure [ diaghr ] presents the evolutionary tracks for the 3 sets of models .
inclusion of thermohaline convection does not modify these tracks in a sensible way : the two curves coincide .
the crosses show the positions of the models which are discussed below .     to [ dlnmu ] .
those situated on the 1.7 m@xmath42 track locate models at 62 myrs , 86 myrs , 299 myrs , 403 myrs , 588 myrs , 788 myrs , 988 myrs , 1188 myrs and 1388 myrs . on the 1.5 and 1.9 m@xmath42 track , only two models are presented , respectively at 838 myrs and 632 myrs.,width=302 ]      for the discussion here , we specifically focus on the case of 1.7m@xmath38 models .
figures [ fe ] to [ dlnmu ] compare the results obtained with ( solid lines ) and without ( dashed lines ) thermohaline mixing , at various epochs during the h burning phase ( crosses in fig .
[ diaghr ] ) .
figure [ fe ] displays the ratio of the iron mass fraction at various times to that of its initial value .
figure [ logk ] shows the opacity profiles at these same evolution times .
figure [ conv ] presents the differences between the radiative and adiabatic gradients , which allows to locate the convective regions .
figure [ mu ] displays the mean molecular weight profiles and figure [ dlnmu ] displays the mean molecular weight gradient , for the interval @xmath43 .
the models reach the main sequence with chemically homogeneous envelopes .
then , due to atomic diffusion , helium sinks and creates a stable ( negative ) @xmath1-gradient below the convective envelope , whereas the radiative acceleration on iron induces an increase of its abundance in the upper stellar layers , and more specifically in the iron - peak elements opacity bump centered at @xmath44 . for the 1.7m@xmath38 case , the models with and without thermohaline mixing begin to diverge at about 400 myrs .    at the age of 62 myrs , in both cases the iron accumulation reaches a factor of 7 below the surface , which does not alter the opacity , and a factor 4 in the opacity bump region , which leads locally to a small opacity increase .
consequently , the radiative gradient ( @xmath45 ) increases , but it remains smaller than the adiabatic gradient so that the corresponding region is still convectively stable . at this time ( 62 myrs )
the model only has a single surface convective zone , due to the ionization of hydrogen . as can be seen in figure [ conv ] , from the close - to - zero value of @xmath46 around @xmath47 , the heii convective zone has just disappeared , due to he settling . at this point ,
the iron accumulation is not large enough to significantly alter the @xmath1-profile . here
, the effect of the he - induced @xmath1-gradient is indeed much more important than the iron contribution . during stellar evolution , as the atomic diffusion proceeds , the region where the helium - induced stable @xmath1-gradient is the steepest slowly deepens inside the star .
as soon as it occurs below the iron opacity bump , the effects of the iron accumulation on the @xmath1-profiles become important .
at 86 myrs , iron accumulation exceeds a factor of 19 and affects the @xmath1-profile .
it is responsible for the change in the slope of the @xmath1-profile at @xmath48 .
the iron enrichment leads to a significant opacity increase and therefore to an increase of the radiative gradient which now exceeds the adiabatic one . in both cases , with or without thermohaline
mixing , a new very narrow convective region appears around @xmath44 .    at 299 myrs ,
this iron convective region becomes significant and leads to chemical homogenization between @xmath44 and 5.4 as seen in the @xmath1- and @xmath49-profiles shown .
the homogenization of the iron convective zone leads to an expansion of the iron accumulation region and reduces the iron maximum value from a factor of 19 to a factor of 12.5 . from the beginning of the main - sequence until 299 myrs ,
the @xmath1-gradient remains null or negative in most models , except below the surface convective zone where it sometimes takes a slightly positive value , due to iron enrichment in the outer regions .
after 400myrs , the @xmath1-gradient induced by he - settling is then located below the iron accumulation region and the effects of iron accumulation become visible . from now on ,
models including thermohaline convection below the iron accumulation layers significantly differ from those without thermohaline convection .    during the subsequent stellar evolution
, radiative levitation leads to iron accumulation in the iron opacity bump region . in models without thermohaline convection ,
this iron accumulation drastically increases with time , reaching a factor of 20 at 400 myrs and up to a factor of 95 at 1388 myrs ( fig .
[ fe ] , dashed lines ) .
the consequences on the opacity profiles can be seen in figure [ logk ] , and those on the radiative gradients in figure [ conv ] .
it also leads to a new convective zone which persists during most of the main sequence lifetime .
iron accumulation in these specific layers leads to a spectacular increase of the @xmath1-values ( fig .
[ mu ] , dashed lines ) .
after 400 myrs , `` @xmath1-bumps '' develop with a rapid increase slightly above @xmath44 and a steep decrease below @xmath50 , which is highly unstable .
the corresponding @xmath1-gradients are shown in figure [ dlnmu ] .    in the models including thermohaline convection induced by these unstable @xmath1-gradients , iron accumulation in the opacity bump region
is drastically reduced and never exceeds a factor of 15 .
it is never completely suppressed however , due to the stabilizing effect of helium settling .
the consequences on the opacity profiles , radiative gradients , @xmath1-values and their gradients can be seen in figure [ logk ] to [ dlnmu ] , solid lines .
the evolution of the internal stellar structure in the two cases studied here is sensibly different .
contrarily to the models without thermohaline mixing , those which include it do not have persistent iron convective zones , although marginal convection appears and disappears several times during stellar evolution .
actually , in models including thermohaline mixing , the iron abundance in the opacity bump always stays close to the critical value for the onset of convection .
when iron accumulation exceeds this critical value , convective mixing occurs and reduces the iron accumulation , making this region stable again . note that when this occasional convective zone appears , it has a smaller extent than the iron bump itself , so that only part of the iron enriched zone is homogenized
this is the reason why a small iron peak remains around @xmath44 in the 58 myrs , 788 myrs and 988 myrs models .
the results obtained for 1.5 and 1.9m@xmath38 models are qualitatively similar to those of the 1.7m@xmath38 models . to illustrate the efficiency of thermohaline convection in these models , we present in figure [ compar3 ] the molecular weight and the iron profiles obtained with and without thermohaline mixing for one model for each mass , located at the middle of the main - sequence ( cf .
[ diaghr ] ) . the radiative acceleration on iron increases with the stellar mass , so that , when no thermohaline convection is taken into account , the iron accumulation is larger in more massive models .
the consequent `` @xmath1-bump '' is also more prominent , as can be seen in figure [ compar3 ] ( dashed lines ) . in the results obtained including thermohaline mixing ( solid lines )
, we can see that iron accumulation is again drastically reduced . in the 1.9m@xmath38 model
, it does not exceed a factor of 10 whereas it goes up to 150 when thermohaline mixing is neglected .
figure [ zconv ] presents the position of the boundaries of the convective zones inside the models with and without thermohaline mixing for the three masses studied here .
the crosses locate the boundaries of the iron convective zone during phases with mixing episodes . when the convective zone settles durably , its boundaries are represented with solid lines .
for the three stellar masses considered , an iron convective zone resulting from a diffusion - induced iron accumulation rapidly appears after the beginning of the main sequence . in models without thermohaline mixing ( fig .
[ zconv ] , upper panels ) , a first phase occurs during which convective episodes are followed by convectively stable periods .
when convection appears in the opacity bump region , the iron peak decreases due to convective mixing .
consequently the opacity also decreases and the radiative gradient becomes smaller than the adiabatic gradient .
the region becomes radiatively stable again . in the 1.5m@xmath38 models ,
this succession of convective and radiative episodes goes on until the end of the main sequence . in the 1.7 and 1.9m@xmath38 models , after a few hundred megayears , iron accumulation has become large enough ( more than a factor of 15 ) , so that convective mixing can not reduce the iron abundance enough to stop convection .
the iron convective zone then persists during the rest of the main sequence lifetime .    in models including thermohaline convection
, the situation is different . here
again , there is a first phase during which mixing episodes alternate with convectively stable periods .
afterwards the onset of thermohaline mixing strongly modifies the convective ( un)stability of the iron opacity bump region .
iron accumulation is reduced and leads , most of the time , to a radiative gradient slightly below the adiabatic gradient .
when , during stellar evolution , the accumulation of iron exceeds , somewhere in the opacity bump region , the critical value for the onset of convection , the induced mixing rapidly reduces the iron abundance and the radiative stability is restored . as a result
the iron convective zone does not persist in these models .
thermohaline convection is expected to take place inside stars any time inverted @xmath1-gradients are built .
this happens in all types of stars where heavy element accumulation occurs in specific stellar layers , either due to radiative levitation or to nuclear reactions during stellar evolution .
it also occurs in case of accretion of hydrogen poor material onto stellar surfaces , due to planet swallowing in early phases , or to accretion of matter from a companion .
here we presented how introducing thermohaline convection modifies the results obtained for the accumulation of iron due to diffusion inside a - type stars .
some other cases where this process has to be introduced were already discussed in previous papers ( e.g. @xcite ) . before giving our general conclusions in the next section , we first propose a discussion of three other cases where thermohaline convection should be introduced and for which we expect important consequences .
the spb ( slowly pulsating b stars ) and the @xmath3 cephei stars are main sequence b type pulsators .
the spb stars oscillate with periods between 0.5 to 5 days , while the @xmath3 cephei stars possess 2 to 8 hours - oscillations .
the former are interpreted as high - order g - modes , the latter as low order p and/or g - modes . in the early nineties , the new opal opacity determinations @xcite allowed to solve the long standing enigma of b type pulsators .
rosseland mean opacities computated from these tables were substancially increased in the z - bump region compared to the previously computed opacities . this opacity enhancement ( as large as a factor of 3 ) provided a boost for the @xmath0-mechanism potentially activated in the metal opacity bump region and
could then explain the excitation of pulsations in @xmath3 cephei and spb theoretical standard models @xcite .
however , as observational capabilities progress , b type pulsators offer new challenges to the theorists .
some of them are now found in low metallicity environments ( e.g. * ? ? ? * and references therein ) .
a large number of pulsation modes are indeed detected in b stars , with frequency spectra sometimes revealing new discrepancies between theory and observations .
as an example , the two @xmath3 cephei stars , @xmath51 eri and 12 lac ( * ? ? ? * ; * ? ? ?
* and references therein ) present frequency spectra which can not be reproduced with b type stars standard models ( e.g. for @xmath51 eri , * ? ? ?
@xcite showed that an `` ad hoc '' factor 4 enhancement of the iron group elements , in the z - bump region could help reproducing the observed pulsation spectrum of @xmath51 eri .
the combined effects of radiative levitation and gravitional settling provides a natural explanation for this local metal enrichment and could then explain the observation of more than 60 @xmath3 cephei stars in low metallicity environments .
following the results of @xcite , @xcite presented preliminary evolutionary models of @xmath3 cephei stars evolved from the zams up to 1 myr : their evolutionary models did not include consistent diffusion computations but they included , at each timestep , a rough estimation of the diffusion - induced iron accumulation in the z - bump region .
after one million years , their inferred iron enrichment is increased by a factor of 12 in the iron opacity bump .
the stability analysis of these models showed that such iron enrichments makes the theoretical frequency spectra denser , which could help reducing the discrepancies between asteroseismic observations and theoretical stability analysis .
following the same idea , @xcite carried out a parametric study of the effects of local iron enhancement on the stability of spb and @xmath3 cephei stars . the parametric iron accumulation profile introduced in their models
is described by a gaussian function centered at @xmath44 , calibrated to be consistent with the iron accumulation profile deduced by @xcite in their a and f star models . as @xcite
, @xcite showed that such accumulation profiles widen the theoretical spb and @xmath3 cephei instability strips and increases the number of excited modes .    the computations by @xcite and @xcite reduced the discrepancies between theoretical and observed pulsations .
however they did not discuss the credibility of the iron profiles introduced in their models , nor the stability of the induced @xmath1-gradients .    on the other hand , @xcite analysed the effects of uncertainties in the opacity computations on the excitation of pulsation modes in b type stars .
they compared models computed with opal opacities @xcite with models computed with the recently updated op opacity tables @xcite .
the op updated data lead to an enhancement of 18% of the opacity in the z - bump .
@xcite showed that these differences considerably affect the theoretical instability strips of @xmath3 cephei and spb stars .
using the op opacities drastically increases the number of excited modes and significantly extends the theoretical instability strips .
these new opacity computations could help solve , at least partly , the discrepancies between theoretical predictions and observations in @xmath3 cephei , spb and hybrid stars .
up to now , no evolutionary models have been computed for b - type stars including consistent diffusion computations and the new op opacities .
however in the light of @xcite results , we expect the iron accumulation derived by @xcite to recouncile observations and theory ( by a factor of 4 ) to be overestimated . in this context
the iron enrichment inferred by @xcite ( by a factor of 12 ) for a 1 myr - model should also be too large , which suggests the presence of a `` forgotten '' process ( such as the thermohaline convection ) which could reduce the diffusion - induced iron accumulation .
contrarily to spb stars , subdwarf b stars ( sdb stars ) are evolved , compact objects with low - mass cores and outer h - rich envelopes .
they are experiencing the core helium burning phase on the extended horizontal branch .
they all present chemical peculiarities attributed to gravitational settling and radiative levitation in the presence of weak stellar winds .
subdwarf b stars host two groups of pulsators : the first one , composed of hot sdb stars , is characterised by rapid ocillations , with periods ranging from 80 to 600s , caused by low - order , low - degree p - modes .
the second group , composed of the coolest sdb stars , oscillates with periods in the 2000 - 9000s range , due to high - order , low - degree g - modes .
both types of pulsators are driven by the @xmath0-mechanism acting in the iron - peak element opacity bump .
the relation between the envelope metal content and the pulsation driving mechanism was first established by @xcite for the short - period b subdwarf pulsators .
they showed that a solar metal content in the envelope of their models is unable to excite modes but that models with envelopes containing an overabundance of metals can globally destabilize modes such as observed .
they argued that this metal enrichment needs only to occur in the driving region itself and not necessarily in the whole envelope as assumed in their models .
a natural explanation for the required local metal enhancement is the radiative levitation process , which is expected to act efficiently in the radiative envelopes of subdwarf b stars . in these stars , as in a and f stars , the competition between gravitational settling and radiative aceleration is expected to produce local accumulations of heavy elements . @xcite
presented more sophisticated and realistic models of sdb stars , the so - called `` second generation models '' , including non uniform iron abundances .
these envelope models rely on the assumption that a state of diffusive equilibrium is reached between radiative levitation and gravitational settling operating on iron ( disregarding other potentially competing processes ) , leading to stable non uniform iron profiles .
@xcite demonstrated that such local iron accumulations in the z - bump region may produce efficient excitation for pulsation modes .
similar models were successfully used for detailed asteroseismic studies of several sdb stars : they globally nicely reproduced the global properties of the pulsating sdb stars and were able to closely reproduce their observed periods @xcite .
as emphasized by @xcite , these `` second generation models '' however suffer from several shortcomings .
the theoretical instability strip predicted by these models is larger than the observed one and the observed period ranges in individual stars are also usually narrower than predicted .
these discrepancies are likely due to an overestimate of the iron abundance in the z - bump region .    in the `` second generation models '' of charpinet
the opacities are computed using the opal opacity tables @xcite .
as discussed previously for spb and @xmath3 cephei stars , the use of the most recent op opacities @xcite may also have a direct impact on the properties of the sdb pulsation driving mechanism .
similarly to other types of stars @xcite , using op opacities instead of opal data , would increase the opacity in the z - bump , resulting in a more efficient pulsation driving mechanism .    here
again we find evidences that the iron accumulation needed to account for the observations must be lower than obtained by the atomic diffusion alone . introducing thermohaline convection could certainly help reducing the discrepancies .
the @xmath2 doradus stars are a - f main sequence pulsators .
their oscillations are interpreted as high - order g - modes with periods ranging from 0.35 to 3 days .
they are driven by a flux blocking mechanism at the base of their convective envelope , as shown by @xcite with frozen convection models and by @xcite and @xcite with time - dependent convection models .
it has been demonstrated that the `` convection blocking '' of radiation can drive high order g - modes only in stellar models for which the temperature at the bottom of the convective envelope , ranges between @xmath52k ( which corresponds nearly to the z - bump region central temperature ) and @xmath53k .
the depth of the convective envelope plays a crucial role in the driving mechanism proposed by @xcite to explain the @xmath2 doradus pulsations .    as a consequence ,
the seismic properties of @xmath2 doradus stars are quite sensitive to atomic diffusion including radiative levitation because of the way they influence the convective zones .
as emphasized by @xcite , the diffusion effects on convection are of various kinds .
the he settling which leads to h - enrichment in the external stellar layers induces an opacity increase which affects the depth of the surface convective zone . on the other hand ,
as discussed in previous paragraphs , radiative levitation produces reservoirs of iron - peak elements which may induce the appearance of an iron convective zone near @xmath52k .
consequently the driving mechanism of @xmath2 doradus stars is expected to be sensitive to the thermohaline convection which occurs below the heavy elements accumulation zone .
in the present paper , we tested the effect of thermohaline convection on the amplitude of iron accumulation which occurs inside a - f type stars due to atomic diffusion .
we computed two series of models for three stellar masses : 1.5 , 1.7 and 1.9m@xmath38 . in both series , atomic diffusion was consistently computed for the elements h , he , c , n , o , ca and fe .
the computations included the radiative accelerations on c , n , o , ca and fe , by using the semi - analytical prescription proposed by @xcite and improved by @xcite .
no thermohaline convection was introduced in the first series of models .
we verified that our results were similar to those obtained by @xcite and @xcite , with a large accumulation of iron in the iron peak elements opacity bump region .
thermohaline convection was then introduced for the second series of models .
considering the short time scales of this mixing process compared to the evolution time scales , and the fact that mixing occurs until the unstable @xmath1-gradients disappear , we introduced it in such a way that the @xmath1-gradients were kept close to zero or negative all along the evolutionary tracks . in practice , we used an extra diffusion coefficient , as proposed by @xcite , which is large enough to reach this purpose .
we presented in detail the results obtained for the 1.7m@xmath38 case , during the main sequence phase .
we also showed some results for the two other stellar masses considered here .
the effect of thermohaline convection is very efficient in strongly reducing the heavy element accumulation in the opacity bump region .
whereas the iron abundance can increase by up to factors larger than 100 when thermohaline convection is neglected , it never increases by more than a factor of 15 when it is taken into account .
it is important to note that some iron accumulation always remains , even in the presence of thermohaline convection .
the reason why this happens is directly related to helium settling , which leads to a stabilizing @xmath1-gradient .
this stabilizing effect allows in turn some heavy material abundance increase up to the point where the @xmath1-profile becomes flat .    as a consequence ,
the cases where iron accumulation has been invoked to account for specific observations in stars remain .
indeed , in many cases the observations would be better explained with a smaller iron accumulation than that given by atomic diffusion alone .
all these cases ( am stars , spb stas , sdb stars , @xmath3 cephei stars and @xmath2 doradus stars as discussed in section 5 ) will have to be studied individually for a detailed comparison with the observations .    in the present computations , nickel was not included .
although less important than iron , it may significantly contribute to the opacity and will have to be included in future work .
no mass loss nor extra turbulence were taken into account . including such macroscopic motion may reduce the heavy element accumulation even without thermohaline convection
. however , thermohaline convection is a physical process which must occur in any case every time the @xmath1-gradient is unstable .
the final accumulation of elements thus strongly depends on this process , even in the presence of other macroscopic motion .    in this paper
, we have demonstrated the importance of thermohaline convection which has to be taken into account in the computation of the abundance variations induced by atomic diffusion .
precise applications to specific stellar cases still have to be undertaken .
they will lead to important consequences on the abundance anomalies observed in stars , as well as their oscillating properties through several processes like @xmath0-mechanism , convective blocking , etc .
the possible applications of thermohaline convection in stellar astrophysics are vast and will certainly be a topic of interest in years to come .
alecian , g. , & leblanc , f. 2002 , , 332 , 891 ausseloos , m. , scuflaire , r. , thoul , a. , & aerts , c. 2004 , , 355 , 352 badnell , n. r. , bautista , m. a. , butler , k. , delahaye , f. , mendoza , c. , palmeri , p. , zeippen , c. j. , & seaton , m. j. 2005 , , 360 , 458 balmforth , n. j. , cunha , m. s. , dolez , n. , gough , d. o. , & vauclair , s. 2001 , , 323 , 362 bourge , p .- o . , & alecian , g. 2006 , asp conference series , vol .
349 , 201 brassard , p. , billres , m. , charpinet s. , liebert , j. , & saffer , r. a. 2001 , , 563 , l1013 charbonnel , c. , zahn , j.p . , 2007 , , 476 , l29 charpinet s. , fontaine , g. , brassard , p. , & dorman , b. 1996 , , 471 , l103 charpinet s. , fontaine , g. , brassard , p. , chayer , p. , rogers , f. j. , iglesias , c. a. , & dorman , b. 1997 , , 483 , l123 charpinet , s. , fontaine , g. , & brassard , p. 2001
, , 113 , 775 charpinet s. , fontaine , g. , brassard , p. , green , e. m. , & chayer , p. 2005a
, , 437 , 575 charpinet s. , fontaine , g. , brassard , p. , billres , m. , green , e. m. , & chayer , p. 2005b
, , 443 , 251 charpinet s. , silvotti , r. , bonanno , a. , fontaine , g. , brassard , p. , chayer , p. , green , e. m. , bergeron , p. , et al .
2006 , , 459 , 565 charpinet s. , van grootel , v. , reese , d. , fontaine , g. , green , e. m. , brassard , p. , & chayer , p. 2008
, , 489 , 377 charpinet s. , fontaine , g. , & brassard , p. 2009
, , 493 , 595 cox , a. n. , morgan , siobahn , m. , rogers , f. j. , iglesias , & c. a. 1992 , , 393 272 dupret , m .- a . ,
grigahcne , a. , garrido , r. , gabriel , m. , & scuflaire , r. 2004 , , 414 , l17 dupret , m .- a . , grigahcne , a. , garrido , r. , gabriel , m. , & scuflaire , r. 2004 , , 435 , 927 dziembowski , w. a. , & pamyatnykh , a. a. 1993a , , 262 , 204 dziembowski , w. a. , moskalik , p. , & pamyatnykh , a. a. 1993b , , 265 588 gargett , a. , ruddick , b. , eds .
2003 , double - diffusion in oceanography ( oxford : pergamon ) , 381 grevesse , n. , & noels , a. 1993 , in origin and evolution of the elements , eds kubono , s. and kajino , t. , 14 guzik , j. a. , kaye , a. b. , bradley , p. a. , cox , a. n. , & neuforge , c. 2000 , , 542 , l57 handler , g. , jerzykiewicz , m. , rodrguez , e. , uytterhoeven , k. , amado , p. j. , dorokhova , t. n. , dorokhov , n. i. , poretti , e. , sareyan , j .- p .
, parrao , l. , and 21 coauthors , 2006 , , 365 , 327 hui - bon - hoa , a. 2008 , , 316 , 55 iglesias , c. a. , & rogers , f. j. 1996 , , 464 , 943 jeffery , c. s. , & saio , h. 2007 , , 378 , 379 jerzykiewicz , m. , handler , g. , shobbrook , r. r. , pigulski , a. , medupe , r. , mokgwetsi , t. , tlhagwane , p. , & rodrguez , e. 2005 , , 360 , 619 kippenhahn , r. , ruschenplatt , g. , & thomas , h .- c . 1980 , , 91 , 175 kiriakidis , m. , el eid , m. f. , & glatzel , w. 1992 , , 255 , 1 koaczkowski , z. , pigulski , a. , soszyski , i. , udalski , a. , kubiak , m. , szymaski , m. , ebru , k. , pietrzyski , g. , woniak , p. r. , szewczyk , o. , & wyrzykowski ,  .
2006 , memorie della societa astronomica italiana , vol .
77 , 336 leblanc , f. , & alecian , g. 2004 , , 352 , 1329 michaud , g. 1970 , , 160 , 641 miglio , a. , montalbn , j. , & dupret , m .- a .
2007 , , 375 , 21 miglio , a. , bourge , p .- o . ,
montalbn , j. , & dupret , m .- a .
2007 , coast , 150 , 209 montalbn , j. , miglio , a. , & thado , s. 2007 , coast , 150 , 137 moskalik , p. , & dziembowski , w. a. 1992 , , 256 , l5 odonoghue , d. , lynas - gray , a. e. , kilkenny , d. , stobie , r. s. , & koen , c. 1997 , , 285 , 657 pamyatnykh , a. a. , handler , g. , & dziembowski , w. a. 2004 , , 350 , 1022 richard , o. , michaud , g. , & richer , j. 2001 , , 558 , 377 richard , o. , thado , s. , & vauclair , s. 2004 , soph , 220 , 243 richer , j. , michaud , g. , & turcotte , s. 2000 , , 529 , 338 rogers , f. j. , & iglesias , c. a. 1992 , , 79 , 507 santos , n.c . , israelian , g. , & mayor , m. 2001 , , 373 , 1019 seaton , m. j. 2005 , , 362 , l1 stancliffe , r.j .
2009 , , 394 , 1051 stancliffe , r.j . ,
glebbeek , e. , izzard , r.g . , & pols , o.r .
2007 , , 464 , 57 stancliffe , r.j .
, & glebbeek , e. 2008 , , 389 , 1828 stankov , a. , & handler , g. 2005 , , 158 , 193 stothers , r. , & simon , n.r .
1969 , , 157 , 673 thado , s , alecian , g. , & leblanc , f. 2009 , in prep .
thompson , i.b .
, ivans , i.i , bisterzo , s. , sneden , c. , gallino , r. , vauclair , s. , burley , g.s . ,
shectman , s.a .
, & preston , g.w .
2008 , , 677 , 556 turcotte , s. , richer , j. , michaud , g. , & christensen - dalsgaard , j. 2000 , , 360 , 603 ulrich , r. k. 1972 , , 172 , 165 vauclair , s. 1975 , , 45 , 233 vauclair , s. 2004 , , 605 , 874


it is well established that a number of early - type galaxies contain small amounts of cold interstellar gas and dust ( eg .
knapp _ et al . _  1989
; colbert , mulchaey & zabludoff 2001 ; sadler _ et al . _
2000 ) presumably obtained through mergers / accretion or stellar mass loss .
what is not well known is the ultimate fate of this cold gas and dust .
van gorkom et al . (
1989 ) estimate that several hundredths of a solar mass per year may be falling onto and feeding the nuclear black hole in some early - type galaxies .
it is also possible that the cold gas may be destroyed due to the presence of a hostile environment . on the other hand ,
if the neutral gas does not get stripped , destroyed , or fall into a central agn , it may settle into a disk and form stars .
it has been argued by okuda _
et al . _  ( 2005 ) that gas disks in co - rich early - types are stabilized by their large dynamical masses , and so star formation is much less likely to occur in these galaxies than in late - type galaxies . however ,
many authors ( e.g. kauffman 1996 ; scorza & bender 1995 , 1998 ) agree that if a significant episode of star formation does occur in an early - type galaxy , then it may be transformed into something resembling a late - type galaxy .
indeed , many early - types are known to contain small stellar disks at their centers ( scorza & bender 1998 ; emsellem _ et al . _
the possibility of star formation in early - type galaxies is important since it challenges the long held belief that early - type galaxies are passive evolvers .
recently , kaviraj _
et al._(2006 ) used the galaxy evolution explorer ( galex ) to observe a large sample of low redshift early - type galaxies in the near ultraviolet ( nuv ) .
they combine these observations with optical data from the sloan digital sky survey ( sdss ) and find that at least @xmath030 percent of the sample galaxies have uv to optical colors consistent with some recent ( @xmath4gyr ) star formation .
their findings are in agreement with semi - analytical @xmath5cdm hierarchical merger models which predict that low redshift early - types will contain @xmath01 to 3 percent of their total stellar mass in stars less than 1 gyr old .
they also determine that the observed blue colors in the star forming fraction of their sample can not be produced by assuming a `` monolithic '' model of evolution .
furthermore , their results show that star formation is equally likely to occur in either high density or field environments .
although these results strongly suggest that a large fraction of early - type galaxies are currently forming stars , no conclusions can be drawn as to the actual star formation rates ( which could show some environmental dependence ) .
whether a galaxy is forming stars is sometimes traced by its placement on the fir - radio correlation known to exist in normal spiral galaxies ( young & scoville 1991 ; condon 1992 ) .
the radio continuum emission in these galaxies arises in hii regions , type ii supernovae and the synchrotron emission from the electrons accelerated by those supernovae .
dust heated by these same massive stars is thought to produce the fir emission from spiral galaxies .
studies of interstellar medium of early - type galaxies have shown that their cold gas is physically similar to that of spirals ( rupen 1997 ; young & lo 1996 ; young 2002 ) .
based on this , it is easy to imagine that star formation might be traced in early - types in the same manner as normal spiral galaxies .
however , many early - types contain agn , and so they may have multiple sources of radio continuum .
clearly , the total radio fluxes combined with a knowledge of the radio morphology are needed in order to disentangle these sources of radio emission .
wrobel & heeschen ( 1988 ) were the first to suggest that several early - type galaxies might be undergoing significant star formation , based on both placement on the fir - radio correlation and the fact that their 6 cm continuum detections are extended on scales of 10 kpc or more .
cotton and condon ( 1998 ) find support for this idea in the radio and fir flux ratios of e / s0s detected in the nrao vla sky survey ( nvss ) . about a quarter of them
have flux ratios consistent with star formation activity .
the other three quarters are presumably powered by agn .    while studies of this kind can rule out the presence of a luminous agn , they can not definitively prove that the radio continuum is associated with recent star formation for several reasons .
firstly , early - type galaxies are fundamentally different from spirals , with different stellar populations and different stellar dynamics , hence different types of energy input into the ism .
secondly , there is evidence to suggest that post asymptotic giant branch ( agb ) stars can provide enough uv radiation to heat the dust and ionize the gas and account for the observed ir luminosities ( ferrari _ et al . _  1999 ) . in other words ,
perhaps both the ir and the radio continuum could be traced back to an old stellar population .
one way to determine if the radio continuum emission is associated with star formation is to compare high resolution continuum images with molecular gas maps .
since molecular gas is the raw material for star formation , it constrains the location and the amount of star formation in the galaxy .
if the radio continuum emission is extended and coincident with the molecular gas , that is strong evidence for star formation activity . on the other hand ,
if the continuum emission is well described by a point source , whereas the co is extended , then the radio continuum is more likely associated with an agn or perhaps a nuclear starburst than with normal disk - type star formation . if the radio emission is extended but in a very different location from the molecular gas , the situation is similar : the radio continuum is more likely to be powered by an agn or the evolved stellar population .
up until now there has been very little knowledge of the co morphology in early - type galaxies .
recently , studies by young ( 2000 , 2001 , 2002 , 2005 ) have produced the largest sample of high resolution molecular gas maps of early - type galaxies to date .    in this paper
we present a high resolution search for radio continuum emission from young s ( 2002 ) seven galaxy sample plus two dwarf early - types .
the radio maps are then compared to high resolution co maps .
this information is then combined with radio morphology and fir - radio correlation placement to trace recent star formation activity .
we also evaluate the toomre q parameter in order to address the question of gravitational stability .
this project allows us the opportunity to independently corroborate previous optical and radio studies that suggest star formation activity in early - type galaxies and to probe the origin of the radio - fir correlation in those galaxies .
the layout of the paper is as follows : a sample description is given in section 2 .
existing and new observations are discussed in section 3 and 4 , respectively .
the results of the new observations are given in section 5 .
radio continuum emission mechanisms are discussed in section 6 .
the origin of fir emission is discussed in section 7 . in sections 9 and 10
, we discuss the implications of the results . in section 11
we present the paper conclusions .
the iras satellite detected many early - type galaxies at 60 @xmath6 m and 100 @xmath6 m ( knapp et al . 1989 ) .
the brightest ones ( 100 @xmath6 m flux @xmath7 1 jy ) have been searched for co using single dish telescopes ( wiklind , combes , and henkel 1995 ( hence forth wch95 ) ; lees et al .
1991 ; knapp & rupen 1996 ) .
young ( 2002 ) recently completed a high resolution survey of the molecular gas from seven bright co detections of wch95 , using the berkeley illinois maryland association ( bima ) and the owens valley radio observatory ( ovro ) interferometers .
of the seven early - types observed , five were detected and resolved in co. the molecular gas in these five galaxies is located in very regular , symmetric rotating disk with diameters ranging from 2 to 12 kpc .
we chose the same seven - galaxy sample in which to image the radio continuum .
we removed ngc 4649 from young s ( 2002 ) original sample since it has recently been determined that the fir from this galaxy is actually associated with its spiral companion ngc 4647 ( temi _ et al . _
2004 ) , and also no co was detected by young ( 2002 ) . two dwarf elliptical companions of m31 , ngc 205 and ngc 185 ,
are added to the sample in order to explore the possibility of star formation in very low mass early - type galaxies .
these two galaxies have been imaged and detected in co using the bima millimeter interferometer ( young 2000 , 2001 ) .
all of the galaxies in our final sample have an r@xmath8 profile or `` a consistent classification as e or e / so in several catalogs '' ( wch95 ) . in some cases ,
these classifications are based on photographic evidence .
recently , ccd imaging suggests that one of our sample galaxies may have a more complicated morphology ( ngc 5666 , donzelli & davoust 2003 ) so we prefer to refer to the galaxies in our sample as `` early - type '' galaxies .
the sample galaxies have a variety of properties and are found in a variety of environments .
some of the galaxies are isolated ; some are in loose groups , and some are deep within the virgo cluster .
some have very smooth , regular morphologies , whereas others are classified as `` peculiar '' because of dust lanes , stellar shells and ripples , etc .
the sample galaxies have optical luminosities that range from 10@xmath9 to 10@xmath10 solar luminosities .
table 1 contains ephemeris information for our galaxy sample .
vla data with matching resolution and signal to noise at 20 cm to that of young s ( 2000 , 2001 , 2002 ) co maps exist for only three of the sample galaxies ( ngc 3656 , ngc 7468 , and ngc 4476 ) .
ngc 4476 suffers from confusion in the nvss but is located within a deep multiconfiguration ( d , c , and b ) vla map of m87 at 20 centimeters ( f. owen , priv . comm . 2001 ) .
we obtained new vla 1.4 ghz c configuration data for the galaxies that are not detected in the nvss . for the purposes of spatial comparisons with young s ( 2002 ) co maps ,
new 20 cm observations were also obtained for ngc 807 , ngc 5666 , and ugc 1503 using the vla b array . we also re - observed ngc 3656 using the vla b array in order to produce a map with similar noise level to the observations just mentioned .
both the c and b array observations were made using two if pairs , each consisting of 50 mhz each in left and right circular polarizations centered on two different frequencies .
this gives a total bandwidth of 200 mhz .
the requested rms noise for these observations is 0.01 mjy beam@xmath2 .
the main parameters of the observed fields are listed in table 3 .
the data were edited , calibrated , and imaged using standard tasks in the nrao astronomical image processing system ( aips ; see van moorsel , kemball , & greisen 1996 for an up - to - date description ) . in the case of ngc 807
it was necessary to combine the c and b array data to minimize confusion due to the bright sidelobes of a nearby bright source .
the integrated radio fluxes , position angles , major axes , and minor axes for our detected sources ( see table 4 ) were determined using the miriad software package ( sault , teuben , & wright 1995 ) .
position angles and deconvolved source sizes were determined by fitting elliptical gaussians to the radio emission .
integrated fluxes were summed in polygonal regions and in concentric ellipses .
great care was taken to insure that no confusing sources contributed to the integrated flux .
a comparison of the co and radio is discussed for each galaxy below :    * ngc 7468 * - ngc 7468 is a member of a small group of galaxies .
it has been characterized as a very peculiar galaxy because of the presence of elliptical isophotes surrounding a nucleus that has been resolved into at least two sources of emission aligned along the major axis of the galaxy ( deeg , duric , & brinks 1997 ) .
another extended optical feature connected with this galaxy is visible approximately 30   to the southwest ( nordgren _ et al . _  1995 ) .
et al . _  ( 1990 ) consider this galaxy to be a possible polar - ring galaxy , with these structures being part of the polar ring .
recent broad band optical and h@xmath11 images reveal the presence of several hii regions coincident with the two central optical components as well as h@xmath11 emission associated with the polar ring ( cair@xmath12s _ et al . _  2001
; paz , madore , & pevunova 2003 ) .
duprie & schneider ( 1996 ) detect a large hi disk ( d@xmath0 11@xmath13 ) using the 300 meter arecibo telescope .
ngc 7468 has been observed in co at several single dish telescopes with varying results .
this galaxy has not been detected in the 2 - 1 line of co ( lees _ et al . _
1991 ) , but there are several tentative detections from the iram 30 meter telescope in the 1 - 0 line ( wch95 ; sage , welch & young 2006 ) .
interferometric observations at 3 mm produced a non - detection ( young 2002 )
. however , an upper limit ( derived from summing over the same spatial scale and velocity area as that of the single dish results ) produces a flux consistent with that measured in the most recent single dish observations .
this is highly suggestive that the single dish detections may in fact be real , and that there is a high probability that ngc 7468 contains some molecular gas .
ngc 7468 already has sufficiently high resolution data in the vla archive ( cox & sparke 2004 ) . in both sets of archive data
the image quality suffers due to confusion from a bright source to the north .
the poor image quality is more pronounced in the 20 cm c array data and is perhaps due to its poorer uv coverage .
the emission in this image is elongated in the same manner as the beam .
both images show that the radio emission in this galaxy is resolved into two large regions ( see figure 1 ) .
the northern blob is much larger , and contains more than 80@xmath14 of the total flux .
the radio emission is coincident with the two central features observed in the aforementioned optical and h@xmath11 images .
no radio continuum is detected along the polar ring .
both sets of data show some low level radio emission that appears to extend outside the eastern part of the galaxy .
the image quality is too poor to say whether this emission is real or just noise .
perhaps deeper vla observations can alleviate this difficulty in the future .
* ngc 3656 * - this is a field early - type that shows indications of a gravitational interaction sometime in the near past ( balcells _ et al . _  2001 ) .
there are several small spherical objects , presumed to be dwarf galaxies , surrounding ngc 3656 ( wch95 ) .
young ( 2002 ) favors the explanation that the molecular gas was acquired during a merger event .
both the co and radio emission peak at the optical center , are elongated in the same position angle and have similar minor axis widths .
however , the radio emission is less extended in radius ( by about 10 ) than the co in the northern parts of the galaxy .
the radio continuum also does not follow the co extension curving off the south east side of the molecular disk .
the extension probably traces gas that has not yet reached dynamical equilibrium and so has not settled into the disk ; we suspect that if the gas is not yet settled , star formation is unlikely to occur there .
a comparison of the molecular gas and the radio continuum can be found in fig 2 .    *
ngc 5666 * - this is an isolated , low - luminosity early - type galaxy .
this galaxy is known to have an extended hi disk ( lake , schommer , & van gorkom 1987 ) , and extended diffuse radio emission ( wrobel & heeschen 1988 ) . the radio continuum in ngc 5666 is extended on scales of a few kpc , and matches the co extent along both axes rather well .
the co is centrally peaked whereas the continuum intensity peak ( @xmath02 mjy beam@xmath2 ) is not located at the center , but to the southeast .
a slice through the major and minor axes of the radio emission shows that the radio emission drops to @xmath01.5 mjy beam@xmath2 inside a roughly circular region , @xmath0 7   in diameter , centered on the optical nucleus .
the co and radio continuum data are at matched resolution and so the absence of a central hole in the co emission suggests that the underlying emitting structures are different , and so the mismatch at the center is not the result of a resolution effect .
a comparison of the molecular gas and the radio continuum can be found in fig 3 .
a recent @xmath15-band image of ngc 5666 from donzelli & davoust ( 2003 ) shows a ring of star formation , 10  in diameter and coincident with the ring of radio continuum emission in figure 3 .
this irregular ring contains some fragments which are suggestive of being spiral arm segments , all within the region where co emission is found .
the galaxy s spectrum is dominated by an old stellar population , however ( donzelli & davoust 2003 ) .
thus it appears that some flocculent spiral structure may be developing in the center of this gas rich early type galaxy .    *
ugc 1503 * - this is a an isolated early - type galaxy .
the radio emission in this galaxy is extended on scales of @xmath05 kpc , and is centered on the optical center of the galaxy .
all of the emission is contained in a ring with a central circular depression ( @xmath05@xmath16  in diameter ) where the radio emission nearly drops to zero .
the peak of the radio emission ( @xmath0 0.3 mjy beam@xmath2 ) occurs in several places along the ring whereas the co emission is centrally peaked
. the extent of the co and radio emission are roughly coincident along both the galaxy s major and minor axes . a comparison of the molecular gas and the radio continuum can be found in fig 4 .
* ngc 807 * - this is an isolated giant early - type galaxy .
the co distribution is asymmetric with most of the molecular gas ( @xmath0 70@xmath14 ) located in the southeast section of the galaxy .
young ( 2002 ) suggests that shear from the differential rotation in the disk would destroy the asymmetry on a time scale of 10@xmath17 years .
since the asymmetry is still present , this may be an indication that the gas was acquired in an interaction less than 10@xmath17 years ago .
the radio continuum in ngc 807 is extended on scales of  9 kpc , and matches the co extent along both axes rather well .
there is an asymmetry in the distribution of the continuum much like that of the co distribution ( most of the emission , @xmath0 60@xmath14 , is located in the southeast ) .
in fact the radio emission ( @xmath00.3 mjy beam@xmath2 ) peaks at roughly the same location as the co emission .
however , it may be possible that the perceived asymmetry in the co is in reality due to differing co - to - h@xmath18 conversion factors existing on alternate sides of the galaxy .
if this is the case , then the correlation of the two emission peaks may be nothing more than a coincidence .
however , the high degree of morphological similarity between the co and radio continuum present in several of the other sample galaxies suggest it is more than a coincidence .
a comparison of the molecular gas and the radio continuum can be found in fig 5 .
* ngc 185 * - this is a dwarf elliptical companion of m31 . the molecular gas in this galaxy
is found in large complexes similar to those observed in the milky way ( young 2000 , 2001 ) .
an hi and optical study of the interstellar medium ( ism ) by young & lo ( 1996 ) shows that complex dust regions near the center of the galaxy are coincident with the detected co emission , and the hi emission exhibits a very clumpy distribution .
we detect a small amount of radio emission ( peak intensity is 5@xmath19 ) in our image of ngc 185 .
the small patch of emission is slightly offset from the center of the galaxy and has a flux of @xmath20jy ( see figure 6 ) .
the emission is roughly coincident with a region of extended h@xmath21 emission detected by young & lo ( 1997 ) .
the surface brightness , luminosity , diameter , and morphology of the h@xmath11 emission are within the ranges of known supernova remnants ( e.g. , long _
et al . _  1990
supernova remnants in m31 have 20 cm radio fluxes that range from s@xmath22=150 - 2600 @xmath6jy ( braun & walterbos 1993 ) . since m31 and ngc 185 are at relatively the same distance we conclude that the small patch of radio emission probably originates from a supernova remnant inside the galaxy
. the small spatial offsets between the co detections and the radio continuum emission would not be noticeable in the other sample galaxies at 1860 mpc .
* ngc 4476 * - this galaxy is a small early - type 12@xmath13   in projection from m87 near the center of the virgo cluster .
the stellar kinematics are normal for an elliptical galaxy with a central velocity dispersion of 132  km  s@xmath2 .
dust is clearly seen on optical photographs , showing the presence of a nuclear dust - ring ( tomita _ et al . _  2000 ) .
young ( 2002 ) favors an external origin for the detected co emission based on the fact that the specific angular momentum of the gas is three times the specific angular momentum of stars .
vla 20 cm continuum observations of ngc 4476 in the b , c , and d configurations yield an upper limit of 5@xmath2310@xmath24 w m@xmath25 hz@xmath2 or 0.5 mjy ( f. owen , priv . comm .
no previous radio detections have been reported for this galaxy .
* ngc 205 * - this galaxy is another dwarf elliptical companion to m31 .
similarly to ngc 185 , the molecular gas in this galaxy is found in large complexes similar to those observed in the milky way ( young 2000 , 2001 ) .
an hi and optical study of the interstellar medium ( ism ) by young & lo ( 1996 ) shows that complex dust regions near the center of the galaxies are coincident with the detected co emission , and the hi emission exhibits a very clumpy distribution .
no radio continuum emission was detected in the image for ngc 205 .
an upper limit to the radio continuum brightness in this galaxy is determined from a 3 sigma upper limit of the rms noise ( @xmath0 0.02 mjy beam@xmath2 ) .
the flux for ngc 205 is located in table 4 .
no previous radio detections have been reported for this galaxy .
there are two major components of the radio continuum emission in galaxies . these are synchrotron emission ( non - thermal ) and free - free emission ( thermal bremsstrahlung ) . in star - forming galaxies ,
synchrotron emission comes from cosmic ray ( cr ) electrons spiraling in the magnetic field of a galaxy .
type ii and type ib supernovae and their remnants are thought to be the main electron accelerators ( kennicutt 1989 , condon 1992 ) .
often , the synchrotron spectra can be approximated by a power law in some limited frequency range . when this occurs , the spectral index @xmath11 is defined as s@xmath26 , where s is the integrated flux density and @xmath27 is the frequency .
thermal emission originates from shocked regions around supermassive stars where superheated electrons collide with the ism .
usually the optically thin thermal emission dominates over the optically thick .
optically thin thermal radiation has a characteristic spectrum ( @xmath280.5 ) , and so it should be distinguishable from the steeper - spectrum non - thermal emission .
typically @xmath11 is determined by comparing the total flux densities obtained at two frequencies with closely matching resolutions .
for example : @xmath29 the available data allow us to calculate the spectral index for three of our sample galaxies : ngc 5666 , ngc 7468 , and ngc 3656 ( see table 4 ) .
the spectral indices determined from the total fluxes in the two frequency bands for these three galaxies range from 0.5 to 0.8 .
pixel by pixel maps of the spectral index do not show any evidence of a central point source of significantly different spectral index .
we note here that in order to avoid variability effects , this type of comparison should be done using observations made with identical visibility coverage on the same day .
we have not done this , so the calculated spectral indices are only rough estimates .    the high resolution radio images rule out nuclear monsters both spectroscopically ( @xmath30  1 for luminous agn ; zhang , lin , & fan 2003 ) and morphologically ( only ngc 3656 exhibits a centrally peaked radio distribution ) .
it is possible that our galaxies contain very low radio power agn .
we estimate the limit to the flux that an unresolved agn would contribute to the total radio power by measuring the radio flux contained inside one beam centered on the optical centers of our galaxies .
we find that these limits are only @xmath010 percent of the total flux and so any contribution from a central point source to the total radio power is small .
the correlation between the fir luminosity and the radio continuum emission observed for normal spiral galaxies is extremely tight and roughly linear ( @xmath4q@xmath312.3 at 1.4 ghz ; helou , soifer , @xmath32 rowan - robinson 1985 ) which implies that the process of star formation produces both types of emission .
it is assumed that the fir and radio emission comes from dust heated by massive stars and cr from supernovae remnants of the same massive stars , respectively .
the fir - radio correlation has been observed to be non - linear at low frequencies where nonthermal emission dominates ( @xmath335 ghz ) as well as for optically - selected samples containing low - luminosity ( @xmath34 ) galaxies ( fitt , alexander , & cox 1988 , devereux @xmath32 eales 1989 , condon 1992 ) .
one explanation for this is that low - luminosity galaxies are less massive than normal galaxies , and are more likely to lose their cr to convection ( chi & wolfendale 1990 ) .
subsequently , a small galaxy would appear under - luminous in the radio with respect to the fir emission , and the correlation would become non - linear even if the low luminosity galaxies are forming stars .
it has also been shown that galaxies with low l@xmath35 also drop off the fir - radio correlation ( yun reddy & condon 2001 ) .
deviations from the correlation may also be due to changes in the surrounding environment such as the presence of an intercluster medium ( icm ) .
the magnetic fields in galaxies can become compressed as they move through a dense icm .
this compression may significantly alter the amount of observed synchrotron emission ( reddy & yun 2003 ) and should enhance the radio continuum emission in the galaxy , making it appear over - luminous in the radio .
the fir flux can be approximated from the iras 100 and 60 micron flux densities using the following relation : @xmath36 where s@xmath37 and s@xmath38 are the 100@xmath6 m and 60@xmath6 m fir fluxes in janskys ( condon 1992 ) .
helou , soifer , & rowan - robinson ( 1985 ) defined the ratio between the far - infrared and radio continuum emission as follows : @xmath39 for @xmath27=1.4 ghz .
an evaluation of the ratio q for the galaxies in our sample show that five ( ngc 3656 , ngc 5666 , ngc 7468 , ngc 807 , and ugc1503 ) of the eight galaxies sit within 3@xmath19 of the radio - fir correlation defined for star forming spirals . the median q value for the sample galaxies detected in radio continuum is @xmath4q@xmath312.7@xmath400.1 ( see figure 7 ) .
we would also like to note that our sample galaxies with the lowest optical luminosities ( @xmath41 ) also have the lowest fir luminosities .
condon , cotton , and broderick ( 2002 ) find that nearly all galaxies with spectroscopically or morphologically determined agns have a fir spectral index @xmath42 @xmath4 1.5 .
the fir spectral index can be determined as follows : @xmath43 the galaxies in our sample have fir spectral indices greater than 1.5 ( see table 5 ) , which is another indication that these galaxies do not contain agn or at least the fir emission does not have a significant contribution from a dust enshrouded agn .
this result is consistent with the radio morphologies of the sample galaxies , which also indicates that luminous agn are not present in these galaxies .
the star formation rate ( sfr ) is calculated using the method derived by condon , cotton , and broderick ( 2002 ) for spiral galaxies .
the fir luminosity is proportional to the sfr as : @xmath44\ ] ] this is the star formation rate for stars of m @xmath45 5 m@xmath1 . to calculate the star formation rate for all the stars between 0.1 to 100 m@xmath1 we simply multiply this result by a factor of 5.5 ( for a salpeter initial mass function ) .
the average star formation rate for stars between 0.1 and 100 m@xmath1 for the sample galaxies detected in radio continuum is @xmath0 2 m@xmath1 year@xmath2 ( see table 5 ) .
the gas depletion time scale can be determined by dividing the h@xmath46 masses found by young ( 2002 ) by the star formation rate .
the h@xmath46 masses for the sample galaxies can be found in table 6 .
the average depletion time scale for our sample galaxies detected in co is @xmath0  0.7 gyr ( see table 5 ) which is much smaller compared to typical gas depletion times scales found in spiral galaxies ( @xmath0  2 gyr ; kennicutt 1983 ) .
the smaller depletion timescales may indicate that star formation in early - type galaxies is more efficient than in normal spiral galaxies .
there are several other possibilities for the smaller depletion time scales
. the star formation could be episodic with prolonged periods between star formation events either due to differing ism conditions throughout the galaxies or to the fact that some molecular gas is still in falling from outside of the galaxies ( pipino _ et al . _  2005 ;
ferreras & silk 2003 ; and many others ) or both .
it is also possible that the gas masses are incorrectly estimated ( this will be discussed further in section 10 ) .
the five best candidates for star formation in our sample are ngc 3656 , ngc 5666 , ngc 807 , ngc 7468 , and ugc 1503 .
these galaxies sit within 3@xmath47 of the radio - fir correlation defined for star forming spirals .
ngc 3656 , ngc 5666 , ngc 807 , and ugc 1503 have co emission that is roughly coincident with the observed radio continuum emission .
in fact , the distribution of the continuum emission in these galaxies resembles those of normal spiral galaxies , and probably traces a ring of star formation .
ngc 3656 has a star formation rate a factor of two greater than the average for our sample .
we might expect this since ngc 3656 has undergone some sort of interaction in the past , which may have produced tidally - induced star formation .
optical and h@xmath11 images show that ngc 7468 possess star forming regions much like those in our own galaxy ( cair@xmath12s _ et al . _
the radio emission is coincident with two large patches of h@xmath48 emission .
ngc 7468 sits right on the radio - fir correlation .
however , the star formation rate for ngc 7468 ( 0.8 m@xmath1 year@xmath2 ) is approximately a factor of three less than the average of our sample .
this smaller sfr could indicate that most of the molecular gas has been used up in a recent star formation event , and the sfr has slowed down to its present value .
the lower star formation rate is also consistent with the fact that this galaxy contains two orders of magnitude less molecular gas than the other galaxies in the sample .
our results indicate that star formation is occurring in more than half of our sample galaxies .
we note here that our galaxy sample is fir and co selected , and so it may not be representative of early - type galaxies .
however , if we consider all of the early - type galaxies in the sample of cotton , condon & broderick ( 2002 ) we find that many of them do lie on the radio - fir correlation ( see figure 7 ) .
their galaxies are both radio and fir selected , not co selected .
full detections that sit within 3@xmath19 of @xmath4q@xmath7@xmath492.3 include @xmath015 galaxies , suggesting that our sample galaxies may not be that unusual .
interferometric images of the molecular gas in ngc 185 and ngc 205 find objects much like the giant molecular clouds in our own galaxy , in terms of sizes and masses ( young 2000 , 2001 ) .
strangely , these galaxies show significant deviations from the spiral galaxy radio - fir correlation in the sense that they are radio faint ( they are the two symbols at the bottom left of fig 7 ) .
the radio nondetection found for ngc 205 and the small amount of radio emission detected in ngc 185 are consistent with studies of low surface brightness dwarfs ( hoeppe et al .
ngc 205 and ngc 185 most likely deviate from the radio - fir correlation because they do not posses enough mass to retain all of their cr .
the superbubbles created by supernovae are able to expand out of the galaxy , compressing the magnetic field .
it is easier for the cr to escape if the magnetic fields are located in shocks produced by the expanding bubbles ( chi @xmath32 wolfendale 1990 ) .
also , the sfr may be slow enough so that the cr from defunct snr have already long since escaped the galaxy ( klein _ et al . _  1992 ) .
ngc 4476 is the most puzzling galaxy in our sample .
there is certainly enough molecular gas in this galaxy to form stars .
why then is there no detected radio emission ?
the star formation rate in this galaxy is low , but @xmath50 is similar to the star - forming ngc 807 .
one explanation for the lack of radio continuum may be that most of the atomic hydrogen has been removed due to ram pressure as it travels through the icm surrounding m87 ( lucero , young , & van gorkom 2005 ) .
perhaps the bulk of the magnetic field has been removed along with the hi in the stripping event , and consequentially the radio is not observed .
a similar phenomenon appears to to be occuring in the virgo cluster spiral galaxy ngc 4522 .
this galaxy has a large amount of extraplanar hi , and radio continuum is indeed associated with the extraplanar gas ( vollmer _ et al . _  2007 ) .
okuda _ et al .
_  ( 2005 ) recently have shown that the giant radio galaxy 3c  31 ( ngc 383 ) contains a massive molecular gas disk , though little star formation activity has been found ( owen , odea , & keel 1990 ) .
et al._(2005 ) argue that star formation is not occurring because the galaxy s deep potential stabilizes the gas disk against collapse , an idea which was also put forward by kennicutt ( 1989 ) for early type galaxies in general .
the stability of a rotating gas disk can be assessed quantitatively by toomre s criterion as revised by sakamoto _
et al . _  ( 1999 ) , @xmath51 where @xmath52 and where @xmath19 is the velocity dispersion , g is the gravitational constant , @xmath53 is the gas surface density , @xmath54 is the angular velocity , v is the gas velocity , and @xmath55 is the epicyclic frequency . writing the @xmath3 parameter in this way elucidates the role of the galaxy s dynamical mass and the gas to dynamical mass ratio ; values of @xmath56 indicate gravitational stability .    for the estimation of @xmath3
we assume our galaxies have an axisymmetric potential in which the velocity of the gas either rises linearly ( @xmath57 ) or is constant ( @xmath58 ) out to the edge of the gas disk ( binney & tremaine 1987 , pp .
120 - 122 ) .
gas and dynamical masses are taken from young ( 2002 ) and can be found in table 6 .
molecular hydrogen masses are multiplied up by a factor of 1.36 to account for the presence of helium .
the true velocity dispersion of the gas is not known so a typical value of 10 km s@xmath2 is assumed .
the ratios m@xmath59/m@xmath60 for this sample range from 0.02 to 0.1 , similar to the values found by koda _
et al . _  ( 2005 ) in early type spirals .
the values of the @xmath3 parameter range from 0.8 to 4 under these assumptions .
only one of the estimated @xmath3 values is less than 1 ( indicating instability ) , but the remainder are sufficiently close to 1 that they may also be consistent with gravitationally unstable disks . the simplified analysis described above neglects the destabilizing influence of the stars in the galaxy .
jog & solomon ( 1984 ) show that the interactions between gas and stars can make the two - fluid system gravitationally unstable even if the gas is stable when considered by itself ( as is done above ) .
thus the values of @xmath3 which are somewhat higher than 1 should not necessarily be interpreted as stability .
we must also note that local surface densities are undoubtedly higher than average at some points in the galaxies and lower at others .
a more detailed analysis is precluded by the fact that @xmath3 scales with the velocity dispersion , which might be a factor of two different from what we have assumed ( e.g. okuda _ et al . _
2005 ; jog & solomon 1984 ) .
we have shown that the existence of a toomre - type gravitational instability in these molecular disks can not be ruled out , and the moderate to low @xmath3 values found here may be consistent with the other evidence suggesting star formation activity in these galaxies .
in particular , ngc 3656 has the lowest @xmath3 value in the sample , and it has one of the highest star formation rates and @xmath50 values .
recently it has been suggested that that the star formation process itself may have an effect on a galaxies position on the radio - fir correlation .
several studies comparing the spatial extent of radio continuum and fir emission in spiral galaxies suggest that the diffusion scale length of cr electrons is much larger than the mean free path of the dust heating photons ( bicay & helou 1990 ; murphy _ et al . _
this means that the continuum emission will become more and more smoothed out compared to the fir emission as the cr diffusion timescales increase .
these studies also show that the cr diffusion timescales tend to increase as the disk averaged sfr in the galaxy decreases ( murphy _ et al . _
so if a galaxy has a sfr below some threshold value , its radio continuum emission could become so smoothed out that it may be hard to detect .
this effect would be amplified if a galaxy is not massive enough to retain its cr , has had its magnetic field removed , or experiences episodic star formation .
it is interesting to note that the galaxies in our sample with the smallest sfr ( @xmath61  0.2 m@xmath62 yr@xmath2 ) also appear to contain little to no radio continuum .
it may be that these galaxies deviate from the radio - fir correlation simply because their sfr are too low .
if the star formation proceeds at rates inferred from the far infrared emission , the transformation from gas to stellar disks will occur over a few billion years . according to young ( 2002 ) , these disks will be rotationally supported , on the order of a few kpc , will contain  1@xmath14 of the total stellar mass , and will be very similar to the stellar disks known to exist in many early - types ( scorza & bender 1995 ; emsellem _ et al . _
our results clearly show that stellar disks can grow out of gas disks in at least some cases .
the gas in these disks could be produced by either internal processes ( temi _ et al . _
2007 ) or it could be left over from mergers ( barnes 2002 ; naab , jesseit & burkert 2006 ) . in other words ,
the presence of a star forming gas disks in spheroidal galaxies is consistent with both the early - assembly ( monolithic collapse ) and the late assembly ( hierarchical growth ) pictures .
while our results can neither confirm of rule out either scenario they do emphasize that the morphology of an early - type galaxy , parameterized by its bulge - to - disk ratio , may evolve somewhat through the process of star formation even while the galaxy remains on the red sequence .
it would be very interesting to see if stellar disks are detectable in early - type galaxies at high redshift .
the h@xmath46 mass in galaxies is estimated using co(1 - 0 ) line emission . the standard conversion factor used for molecular clouds in the milky way is given by n(h@xmath46)/i@xmath63 = 3@xmath2310@xmath64 cm@xmath25 k
km s@xmath2)@xmath2 ( sanders , solomon , @xmath32 scoville 1984 ) .
the use of this conversion factor assumes that the mean properties of the molecular gas in the early - type galaxies in the sample are similar to those in the milky way .
the molecular clouds in early - types might not exhibit the same properties , such as temperature and density , as those found in the milky way ( georgakakis _ et al . _
it is expected that such differences will modify the co to h@xmath46 conversion factor .
if the molecular gas densities in early - types are larger than in the milky way , the conversion will underestimate the the h@xmath46 mass . if instead the densities are the same , but the temperatures are lower than those in the milky way , the conversion factor will overestimate the h@xmath46 mass .
a modified h@xmath46 mass will affect the depletion timescales as well as the gas mass to dynamical mass ratios .
if gas masses are overestimated then the depletion timescales become even faster , and the gravitational stability greater . in order to achieve depletion time scales like those in spiral galaxies
the gas mass would have to be underestimated by a factor of 4 .
in this paper we present new high resolution radio continuum observations of eight fir and co selected early - type galaxies , and compare them to high resolution co maps in order to investigate the importance of star formation in these galaxies . in summary : +   + @xmath66 five of the eight galaxies ( ngc 807 , ngc 3656 , ngc 5666 , ngc 7468 and ugc 1503 ) sit within 3@xmath47 of the radio - fir correlation defined for star forming spirals .
+   + @xmath66 the radio continuum is extended on scales of a few to @xmath09 kpc .
good spatial correlations between the co and radio emission exist in four of the sample galaxies ( ngc 807 , ngc 5666 , ugc 1503 , and ngc 3656 ) .
this morphological similarity confirms wrobel & heeschen s ( 1988 ) suggestion that early - type galaxies on the fir - radio correlation may be forming stars .
many of the early - type galaxies in the nvss catalogue also sit on the radio - fir correlation .
these results are in agreement with other recent work suggesting that star forming early - type galaxies are more common than has been thought in the past .
+   + @xmath66 the average star formation rate for the galaxies detected in radio continuum is @xmath0 2 m@xmath1 year@xmath2 .
the average gas depletion time scale for the galaxies detected in co is @xmath0 0.7 gyr , which is much less than the gas depletion time scales observed in star forming spiral galaxies ( @xmath0  2 gyr ) .
+   + @xmath66 a simple gravitational stability estimate suggests that the molecular disks are likely to be unstable , driving the inferred star formation activity .
+   + @xmath66 the molecular gas in our star forming early - type sample galaxies will produce stellar disks containing @xmath0  1@xmath14 of the total stellar population , and radii of a few kpc .
+   + @xmath66 a few of the sample galaxies ( ngc 4476 , ngc 185 , and ngc 205 ) deviate significantly from the fir - radio correlation in that they are radio faint . in the case of the two dwarf early - types
, their small masses may allow a significant number of cosmic rays to escape . in the case of ngc 4476
, the large scale magnetic field may have been removed along with hi gas due to the passage through a hot intercluster medium .
we thank the anoymous referee for a very constructive report , and joan wrobel for her insightful comments which have made this work much stronger .
this research has made use of the nasa / ipac extragalactic database ( ned ) which is operated by the jet propulsion laboratory , california institute of technology , under contract to the national aeronautics and space administration .
last but not least , this work has been supported by nsf grants ast-0074709 and ast-0507432 to new mexico institute of mining and technology .
lisa m. young thanks new mexico tech for a sabbatical leave and is grateful to the university of oxford sub - department of astrophysics for their hospitality .    lcccccl + galaxy & type@xmath67 & ra & dec & l@xmath68 & distance & environment + & & & & ( 10@xmath17  l@xmath69 ) & ( mpc ) & + ngc185 & de3p & 00h38m57.5s & 48d20m12s & 0.057 & 0.63 & group + ngc205 & e5p & 00h40m22.0s & 41d41m07s & 0.30 & 0.85 & group + ugc1503 & e & 02h01m19.8s & 33d19m46s & 17 & 69 & field + ngc807 & e & 02h04m55.7s & 28d59m15s & 32 & 64 & field + ngc3656 & io & 11h23m38.4s & 53d50m31s & 16 & 45 & merger remnant + ngc4476 & e / so & 12h29m59.1s & 12d20m55s & 3.5 & 18 & virgo cluster + ngc5666 & ec & 14h33m09.2s & 10d30m39s & 6.3 & 35 & field + ngc7468 & e3p & 23h02m59.2s & 16s36m19s & 4.0 & 28 & group +