strong ovi absorptions observed around starforming galaxies at @xmath8 with impact parameter as high as 150 kpc reveal huge gaseous galactic haloes @xcite .
even conservative estimates lead to amount of gas in them far exceeding the gas reservoir in galaxies themselves .
such a conclusion is based on the assumption that gas in the haloes has solar metallicity .
the arguments underlying this assumption stem from the standard estimate of oxygen mass from the observed column density @xcite [ mo ] m_o = 5r^2n__ovim__o f_hit ( 0.2 f__ovi ) = 1.210 ^ 7 ( 0.2f__ovi)m _ , where a typical column density @xmath9 @xmath3 , the halo radius @xmath10 kpc , the hit rate correction factor ( covering factor ) @xmath11 are assumed following by @xcite , @xmath12 is atomic mass of oxygen .
the fractional abundance of ovi under thermal collisional ionization equilibrium never exceeds @xmath13 ( e.g. * ? ? ?
* ; * ? ? ?
* ) , such that within this assumption eq .
( [ mo ] ) provides a lower estimate of oxygen mass in the halo @xcite .
it is readily seen from here that the mass @xmath14 , i.e. around 10 to 70% of oxygen mass in the ism , is indeed a conservative estimate .    in principle , this conclusion for galactic haloes to bear such a large gas mass would might solve the problem of missing baryons and metals @xcite , though requiring enormously high oxygen production and mass ejection rates .
moreover , the fraction @xmath13 under collisional ionization equilibrium is kept only in a very narrow temperature range : @xmath15
k ( e.g. * ? ? ?
* ; * ? ? ?
it is therefore totally unrealistic to assume that all observed haloes keep their temperatures within such a narrow range independent on distance from the host galaxy .    in their estimates
@xcite used the ionic fractions calculated under collisional and/or photoionization equilibrium conditions , i.e. independent of time @xcite .
however , ionization state of gas situated in time - dependent ( nonequilibrium ) environment can differ qualitatively from the one settled on to equilibrium , particularly , for solar metallicity @xcite .
this difference can be smaller for low - density gas exposed to a strong ionizing field .
however , in the process of galactic evolution both the magnitude and the shape of the radiation spectrum can change . in addition the extragalactic spectrum , which can be important on the periphery of galactic haloes ,
does also evolve @xcite . under such conditions
it is natural to expect the ionic composition to experience time variations .    in this paper
we therefore concentrate on the question of whether evolution of the ionizing radiation field can result in considerable changes of fractional ionization of oxygen such to make estimates of gas mass in galactic haloes more reliable . in next sections
we will demonstrate that indeed under time - dependent conditions the observed column densities of ovi correspond to at least half order of magnitude less massive haloes , which can stay in a much wider range of physical parameters .
the paper is organized as follows . in section 2
we describe the details of the model . in section 3
we present our results .
section 4 summarizes the results .
thermal and ionization state in our model is fully time - dependent : the model involves ionization and thermal evolution of gas located at radii @xmath16  kpc in the galactic halo , exposed to extragalactic and galactic time - dependent ionizing radiation field .      in this paper
we only briefly touch main features of calculation of the ionization and thermal evolution of gas immersed into external time - dependent ionizing radiation .
the details can be found in @xcite .
we study ionization and thermal evolution of a lagrangian gas element : a gas parcel is assumed to be optically thin to external ionizing radiation . in the calculations we include all ionization states of the elements h , he , c , n , o , ne , mg , si and fe . we take into account the following major processes : photoionization , multi - electron auger ionization process , collisional ionization , radiative and dielectronic recombination , collisional excitation as well as charge transfer in collisions with hydrogen and helium atoms and ions .    the total cooling and heating rates are calculated using as a subroutine the photoionization code cloudy ( ver . 10.00 ,
more specifically , we input into cloudy code a given set of all ionic fractions @xmath17 calculated at temperature @xmath18 , gas density @xmath19 and external ionization flux @xmath20 and obtain the corresponding cooling and heating rates .
the latter also includes compton heating from x - rays .
for the solar metallicity we adopt the abundances reported by @xcite , except ne for which the enhanced abundance is adopted @xcite . in all calculations
we assume the helium mass fraction @xmath21 , which corresponds to [ he / h]=0.081 , and is close to the observed one @xcite .
we solve a set of 96 coupled equations ( 95 for ionization states and one for temperature ) using a variable - coefficient ordinary differential equation solver @xcite .
we considered the two regimes of gas evolution : isochoric and isobaric .
the isochoric regime suggests that gas density in a cloud is kept constant , while in isobaric model gas pressure is assumed constant .
the two regimes correspond to two opposite limits of the ratio between the cooling and the sound crossing times : @xmath22 and @xmath23 , correspondingly . in the external heating radiation field
isobaric models show essentially similar thermal and ionization evolution , though on longer time scales due to decreasing density coupled to increasing temperature .
it results in an increase of the cross - section of gas clouds such that their covering factor increases as well .      in the process of galaxy evolution
the stellar content changes : massive stars produce enormous number of uv photons and form ultimately compact objects which emit hard ionizing photons
. spectrum of the ionizing radiation escaping galactic interstellar medium ( ism ) and exposing further the halo depends on amount of metals in the ism disk absorbing ionizing photons , and thus on chemical evolution of the galaxy . in order to follow evolution of stellar mass , metallicity and galaxy spectrum , we use the spectro - photometric code pegase @xcite .
we assume a schmidt - like power - law starformation rate ( sfr ) : @xmath24 , typical sfr for massive starforming galaxies , where @xmath25 is the normalized mass of gas in @xmath26 . in some regions of the galaxy sfr
can be inhibited by gas outflows from the disk , however , when averaged over the whole disk the sfr remains sufficiently high over the whole period of active star formation . in our models
we assumed a closed - box regime .
in general , though this can not be applied to galaxies with active star formation . however , many parameters related to mass and energy exchange between galaxies and the intergalactic medium , such as the rates of mass ejection from and mass accretion from the ambient medium , gas metallicity and corresponding cooling rate , clustering of sne explosions , etc are highly uncertain and hard to be coherently described phenomenologically .
gas in galactic haloes is exposed by a cumulative ionizing background consisting of the extragalactic and the galactic components .
the extragalactic component is uniform on galactic halo scales , and is nearly constant on timescales @xmath27  myr , while changes significantly on longer times . for the extragalactic background
we accept the spectrum described by @xcite .
its evolution covers redshifts @xmath28 to @xmath29 divided by 49 equally spaced log bins .
contribution from possible flux sources varies over cosmological time , and irregular changes in different bands of the spectrum with redshift can be met .
for this reason simple linear approximation between neighbour redshift bins was used .    the galactic component instead ,
may change on much shorter timescales . in general , the corresponding time scale is close to the lifetime of massive stars , i.e. @xmath30  myr .
this galactic radiation component originates in the central star forming region , and is seen from outer parts of the halo at distances @xmath31  kpc as a nearly spherical domain of size of @xmath32  kpc .    the ionizing radiation from stellar population is partly modified by absorption in the interstellar medium of underlying bulge and disk .
in order to account this absorption we assume that ionizing photons pass through a layer of neutral gas in galactic disc with the optical depth @xmath33 , throughout the paper @xmath34  @xmath3 , @xmath35  @xmath3 are considered as fiducial .
the corresponding optical depths at the hi and hei lyman limits are as high as @xmath36 and @xmath37 , respectively . as a result only photons with @xmath38  ev escape the galaxy and penetrate into halo ; the ionizing flux with photons of @xmath39  ev decreases as @xmath40 . in what follows
we will discuss the dependence of our results on the hi column density .    in our model
the galactic uv spectrum is calculated with making use the pegase code @xcite , which gives spectral luminosity in the range from 91  to 160  @xmath41 m . in order to extend the spectrum to higher energies ( up to @xmath42 ev , responsible for ionization of highly charged ionic species ) we use the empirical relation between the x - ray luminosity and the star formation rate ` @xmath43 ' @xcite .
this relation is well established for massive starforming galaxies considered here .
overall the cumulative spectrum varies on time scales from several to hundreds of millions years .
we consider gas in outer haloes of massive ( milky way type ) starforming galaxies with stellar mass of several @xmath44 .
recent simulations of the milky way halo show that it extends up to the virial radius of the milky way ( i.e. , @xmath16  kpc ) with densities ranging within @xmath45  @xmath7 @xcite .
observational estimates of the circumgalactic gas density around the milky way and other local group galaxies give similar numbers @xmath46
@xmath7 at @xmath47  kpc @xcite . in our calculations
we follow these numbers and set @xmath48  @xmath7 in the circumgalactic volume .
we consider both isochoric and isobaric regimes .
we start the calculations at @xmath49 ( the lookback time is around 10 gyrs ) .
this timescale is nearly the cooling time for hot gas with @xmath50  k and @xmath45  @xmath7 @xcite .
the last major merging for milky way - type galaxies is thought to occur earlier than @xmath51 ( e.g. , * ? ? ?
the initial ionic composition and the temperature are set equal to the ones corresponding to photoequilibrium in gas exposed to the extragalactic haardt & madau spectrum at @xmath49 .
such radiation field is sufficient to force low density gas with @xmath19 in the accepted range to settle quickly onto photoequilibrium @xcite .
calculations cover physical time scales much longer than the relaxation time scale of ionization and thermal processes in gas exposed to the time - dependent spectrum adopted here .
gas metallicity in our models is assumed to range within @xmath52 to @xmath53 .
higher metallicities are rare in dla qso and dla grb absorbers @xcite . the lower limit lies around ten times above the upper limit of the igm metallicity at @xmath54 ( e.g. , * ? ? ?
* ; * ? ? ?
( solid black ) , and the stellar mass , @xmath55 ( solid grey ) ; the initial time corresponds to @xmath49 .
the four models with the rate @xmath56 are shown : two models with fixed @xmath57  myr  @xmath58 and the initial gaseous mass @xmath59 ( model a  solid line ) , and @xmath60 ( model c  dashed line ) .
the other two models are with fixed @xmath60 and different @xmath61  myr  @xmath58 ( model d  dotted line ) , and @xmath62  myr  @xmath58 ( model e  dash - dotted line ) . , width=302 ]    from the chemical evolution models calculated with the pegase code
we obtain the time - dependent starformation rate , stellar mass , gas metallicity and spectral luminosity .
we simulated several models corresponding to massive galaxies and chose two of them , whose parameters , namely the sfr and the stellar mass , are close to those of the starforming galaxies described by @xcite . with the accepted sfr law @xmath63 , where @xmath25 is the normalized mass of gas in @xmath26 and through the paper @xmath64
is assumed , the notations are as in @xcite . in total
we consider four models .
the first two suggest a fixed normalization factor @xmath57  myr  @xmath58 , and two different values of the initial gaseous mass : @xmath59  model a , and @xmath65  model c. the other two models assume a fixed initial gaseous mass @xmath60 , and two different normalization constants @xmath66 : one with @xmath61  myr  @xmath58  model d , and the other with @xmath62  myr  @xmath58  model e.    figure  [ figssfr - time ] presents the star formation rate , @xmath67 , the specific starformation rate , @xmath68 , and the stellar mass , @xmath55 : model a ( solid line ) , model c ( dashed line ) , model d ( dotted line ) and model e ( dash - dotted line ) .
two features are to note : first , the ssfr reveals differences between the models only after several hundreds of myrs , and second , models a and c show practically equal ssfr , and as a result gas in these models is equally converted into stars by @xmath69  gyr .
instead , model d ( @xmath61  myr  @xmath58 ) leaves the galaxies gas - rich  only @xmath70  % of gaseous mass exhausts by this time .
the characteristic time for the gas to exhaust is @xmath71 , which gives @xmath72 myr for all considered models , resulting in a
 nearly constant sfr ( upper panel in figure  [ figssfr - time ] ) over first 200  myr with the sfr@xmath73 , and decreasing on later times as sfr@xmath74 .
figure  [ figsfr - mstar ] shows the dependence of specific starformation rate on the stellar mass for the models a ( pentagons ) , c ( circles ) , d ( up - triangles ) and e ( down - triangles ) . filled large symbols
mark time moments shown nearby .
data for the starforming and passive galaxies studied by @xcite are depicted by small squares and rhombs , correspondingly .
the grey - scale map shows sdss+galex galaxies @xcite .
it is readily seen that almost all points for the starforming galaxies are locked between tracks of the models considered here .
one can expect therefore that spectral properties of starforming galaxies from @xcite are similar to those in the models a and c during the latest @xmath753 - 4  gyrs of their evolution .
passive galaxies have an order of magnitude lower ssfr .    . , width=302 ]     and a distance from the galaxy @xmath76  kpc , which consists of the uv and x - ray galactic spectrum , attenuated by galactic neutral gas ( dash line ) , the extragalactic ionizing background ( dotted line ) .
the galactic spectral luminosity is shown by dash - dotted line ( the right axis ) .
, width=302 ]    dash - doted line in figure  [ figspec ] shows the galactic spectral luminosity at @xmath77  gyr , which is equal time elapsed from @xmath49 to 0.2 .
a significant break is clearly seen at the wavelength 91   , which corresponds to the minimum wavelength in the spectrum reached in the pegase code .
the break is due to an exponential decrease of the number of such hard photons emitted by stellar population ( see e.g. * ? ? ?
* ) , as they are only produced by very massive stars , whose number is very small .
as soon as we consider galactic evolution on timescales longer than 1 - 3  gyr x - ray binaries are expected to have already formed , and we extend the spectrum to higher energies assuming the empirical relation between x - ray luminosity and star - formation rate `` @xmath43 '' @xcite .    by the thick grey line figure  [ figspec ] shows also an example of the total spectral radiation flux exposing a given gas parcel located at distance @xmath76  kpc from the galaxy evolved till @xmath78 .
the total spectrum consists of the galactic ( dash line ) and extragalactic ( dotted line ) ionizing photons . at low energies , @xmath79  ev
, the extragalactic contribution dominates at large distances , @xmath80  kpc , while the stellar population turns into play in uv range at smaller distances .
strong absorption of the galactic ionizing photons ( @xmath81  ev ) in the galactic disk with the fiducial @xmath82 and @xmath83 values results in the total dependence of the ionic composition of halo gas on the extragalactic background . the significance of absorption by the galactic disk can be understood from comparison of the galactic spectral luminosity shown by dash - dotted line and the cumulative flux ( thick grey line ) .    a narrow bump at @xmath84  ev is due to the galactic photons survived against absorption in thick disk .
its magnitude is obviously determined by our choice of the neutral column densities @xmath82 and @xmath83 in the disk . for the fiducial values of @xmath82 and @xmath83
the optical depth is about 3 for photons with @xmath85  ev .
higher column densities can erase the bump and the extragalactic flux becomes dominating in the whole galactic halo .
decrease of the galactic contribution is also seen at larger distances from the galaxy due to dilution @xmath40 . note that the ionization potential of ov , @xmath86  ev falls exactly in this the range @xmath84  ev .
this means that the fraction of ovi ions can increase in galaxies with lower column densities ( @xmath82 and @xmath83 ) in their disks .
( dash lines ) and @xmath87 ( dash - dotted lines ) , where the superscipts g and e refer to galactic and extragalactic backgrounds . thick and thin lines correspond to redshifts 1.8 and 0.2 , respectively .
, width=302 ]    in general , the excess of galactic photons with energies higher than 113.9  ev contributes crucially into ionization kinetics of ovi ion .
figure  [ figflux - ratio ] presents the ratio of monochromatic galactic and extragalactic fluxes at 113.9  ev ( solid lines ) in model a for our fiducial column densities @xmath82 and @xmath83 .
it is clearly seen that the galactic flux dominates within @xmath88  kpc at redshifts @xmath89 . at lower redshifts
the region of predominance of galactic flux widens mostly due to a steep decrease of the extragalactic flux in these epochs .
the monochromatic radiation at the hi lyman limit is fully absorbed in the disk for the fiducial @xmath82 and @xmath83 column densities .    in order to estimate both the escape of hi ionizing photons and the efficiency of ov ionization we calculate radial dependences of the ratios of ionizing fluxes @xmath90 and @xmath87 , where the superscipts _g_ and _e_ refer to the galactic and extragalactic contributions .
figure  [ figflux - ratio ] shows these ratios at @xmath91 and 0.2 for model a. in the energy range @xmath92  ev the galactic contribution dominates up to distance @xmath93  kpc at @xmath91 , and extends even to @xmath88  kpc at @xmath78 due to strong decrease of the extragalactic background at low redshifts .
it is obvious , that increase of ionizing flux enlarges radius of the zone of predominance of galactic ionizing photons : for instance , this zone increases from @xmath94  kpc in model a to @xmath95  kpc in model c at @xmath91 .
decrease of absorption in the disk does also promote the zone of the galactic predominance to grow : it reaches @xmath96  kpc for @xmath97  @xmath3 even in model  a.       and @xmath83 .
, width=302 ]    figure  [ figt - ovi - evol ] shows evolution of temperature ( lower panel ) and ovi fraction ( upper panel ) in gas located at different distances from the galactic center ; here we assume the sfr as in model a and our fiducial column densities @xmath82 and @xmath83 .    in the beginning , @xmath51 , oxygen is mainly locked in the ovii state due to high ionizing flux .
its fraction establishes around @xmath98 under the influence of the extragalactic ionizing background .
after several hundreds million years the extragalactic background starts decreasing following the cosmic star formation rate .
it should result in a quick transition from ovii to ovi . in a cooling plasma ovi recombines rapidly into lower ionization states and practically disappears shortly  a well - known ovi `` fragility '' .
however , in our case the excess of photons with @xmath99  ev emitted by starforming galaxies does not allow ovi to recombine .
consequently , its fraction remains almost frozen at the level @xmath98 in the region @xmath100 kpc over the range @xmath101 till 0 . in our models
@xmath102 reaches @xmath0 in a low - metallicity ( @xmath53 ) gas within @xmath103  kpc from @xmath104 to @xmath105 .
this value is several times higher than the maximum ovi fraction , @xmath106 , reached in gas exposed only to the extragalactic background ( e.g. * ? ? ?
* ; * ? ? ?
* ) , and in nonequilibrium collisional gas evolved from @xmath107  k ( see e.g. * ? ? ? * ) .
note that even at large distances @xmath108 kpc @xmath102 remains higher than @xmath109 .
the temperature of gas with such high ovi fraction is within @xmath110  k ( see lower panel in figure  [ figt - ovi - evol ] ) .
the dependence of @xmath102 on metallicity is weak : for instance , ovi fraction in gas with @xmath111 is @xmath112 within distances @xmath113  kpc in model  a.    cooling of gas with @xmath114 exposed only to the extragalactic background is mainly due to hydrogen and helium , whereas metals ( oxygen and carbon ) play a minor role in radiation losses ( e.g. * ? ? ?
* ; * ? ? ?
* ) . a considerable increase of ovi fraction in the zones where galactic ionization dominates enhances the contribution of oxygen into cooling . as a result ,
gas temperature is lower in these regions as seen on the lower panel in figure  [ figt - ovi - evol ] .
deviation of column densities from the fiducial values accepted above may result in a considerable change of the overall picture due to changes of the interrelation between the fractions of galactic and extragalactic ionizing photons . in order to understand how sensitive is the oxygen ionization state to surface density of the underlying galactic gaseous discs we calculate several models with different @xmath115 and @xmath116 . for simplicity
we assume that @xmath117 .
figure  [ fig - ovi - hi ] presents the galactic part of the spectrum at distance @xmath76  kpc and redshift @xmath118 for several values @xmath115 ( upper panel ) , and the evolution of ovi fraction in gas exposed to this radiation field ( lower panel ) . for column density
@xmath119  @xmath3 the galactic spectrum in model a is fully absorbed in the range 13.6 - 80  ev . only for a hundredth of this column density
the attenuation becomes less depressing .
a lower column density @xmath120  @xmath3 allows a higher fraction of ionizing photons with the energy above the ionization potential of ov ( 113  ev ) to penetrate into the halo and support a higher fraction of ovi ( lower panel ) : the ovi fraction reaches @xmath121 and slightly increases for lower column densities .
this value is considerably higher than in collisionally ionized gas and in gas photoionized only by the extragalactic background .
the reason for such a high fraction of ovi is in the break at the energy slightly lower than the ovi ionization potential ( 136  ev ) .
it is worth noting that the uncertainty of the spectral energy distribution in this energy range within stellar population codes is rather high .
this partly stems from the stellar atmosphere models and the rarity of extremely massive stars @xcite . on one hand
this is a shortcoming of our model , though on the other , the rarity of massive stars and their short lifetime allows us to think that stellar contribution into circumgalactic ionizing field falls steeply down at energies @xmath122  ev ( see the spectra of the hottest stars in * ? ? ?
* ) , which is lower than the ovi ionization potential ( 136  ev ) .
however it is important to note that that the predominance of galactic flux over the extragalactic background in the range around 110 - 130  ev results in ovi fraction as high as 0.6 - 0.8 .
higher hi column density in underlying galactic discs , @xmath123  @xmath3 , heavily erases galactic flux in the range @xmath124  ev , which , for instance in model a becomes lower than the extragalactic background and the effect of an enhanced ovi fraction vanishes .
kpc and redshift @xmath118 for several values @xmath115 ( upper panel ) , and the evolution of ovi fraction in gas exposed this spectrum ( lower panel ) .
the thick black line on the upper panel shows the extragalactic part of the spectrum .
redshift dependence of the fraction of ovi ions on the lower panel for different column densities on lower panel shows a catastrophic change at @xmath125 @xmath3 with a heavily depressed @xmath102 at lower @xmath126 . the ratio @xmath127 is assumed .
, title="fig:",width=302 ]  kpc and redshift @xmath118 for several values @xmath115 ( upper panel ) , and the evolution of ovi fraction in gas exposed this spectrum ( lower panel ) .
the thick black line on the upper panel shows the extragalactic part of the spectrum .
redshift dependence of the fraction of ovi ions on the lower panel for different column densities on lower panel shows a catastrophic change at @xmath125 @xmath3 with a heavily depressed @xmath102 at lower @xmath126 . the ratio @xmath127 is assumed . , title="fig:",width=302 ]      k ( thick grey ) and @xmath129  k ( thin grey lines ) for surface gravity @xmath130 , @xmath131 and @xmath132  cm  s@xmath133 ( dark grey , grey and light grey lines ) at solar metallicity taken from the stellar library calculated by @xcite are shown by grey lines ( the left axis ) .
the galactic spectral luminosity accepted in model a is depicted by thick black line ( the right axis ) .
several power - law functions @xmath134 are shown for comparison . _
middle panel : _ the spectral luminosity in models a ( black line ) and a+1 .. a+5 at @xmath78 , corresponding to different power - law index @xmath135 as shown by grey lines . _ bottom panel : _ the cumulative ionizing background flux at a distance from the galaxy @xmath76  kpc and redshift @xmath78 in models a , a+1 .. a+5 ( thick lines ) , the extragalactic ionizing background ( thin dotted line ) .
models c+1 , c+2 , ... are defined as the model c with the spectra accepted for models a+1 , a+2 , ... , correspondingly .
, title="fig:",width=291 ]  k ( thick grey ) and @xmath129
k ( thin grey lines ) for surface gravity @xmath130 , @xmath131 and @xmath132  cm  s@xmath133 ( dark grey , grey and light grey lines ) at solar metallicity taken from the stellar library calculated by @xcite are shown by grey lines ( the left axis ) .
the galactic spectral luminosity accepted in model a is depicted by thick black line ( the right axis ) .
several power - law functions @xmath134 are shown for comparison . _
middle panel : _ the spectral luminosity in models a ( black line ) and a+1 .. a+5 at @xmath78 , corresponding to different power - law index @xmath135 as shown by grey lines . _ bottom panel : _ the cumulative ionizing background flux at a distance from the galaxy @xmath76  kpc and redshift @xmath78 in models a , a+1 .. a+5 ( thick lines ) , the extragalactic ionizing background ( thin dotted line ) .
models c+1 , c+2 , ... are defined as the model c with the spectra accepted for models a+1 , a+2 , ... , correspondingly . , title="fig:",width=283 ]  k ( thick grey ) and @xmath129  k ( thin grey lines ) for surface gravity @xmath130 , @xmath131 and @xmath132  cm  s@xmath133 ( dark grey , grey and light grey lines ) at solar metallicity taken from the stellar library calculated by @xcite are shown by grey lines ( the left axis ) .
the galactic spectral luminosity accepted in model a is depicted by thick black line ( the right axis ) .
several power - law functions @xmath134 are shown for comparison . _
middle panel : _ the spectral luminosity in models a ( black line ) and a+1 .. a+5 at @xmath78 , corresponding to different power - law index @xmath135 as shown by grey lines . _ bottom panel : _ the cumulative ionizing background flux at a distance from the galaxy @xmath76  kpc and redshift @xmath78 in models a , a+1 .. a+5 ( thick lines ) , the extragalactic ionizing background ( thin dotted line ) .
models c+1 , c+2 , ... are defined as the model c with the spectra accepted for models a+1 , a+2 , ... , correspondingly . , title="fig:",width=283 ]    kpc from the galaxy exposed to the spectra corresponding to models a , a+1 .. a+5 .
, width=302 ]     it is obviously clear from discussion in section 3.2 that the enhanced ovi fraction is mainly due to the excess of galactic photons with @xmath136  ev in the bump at @xmath84  ev .
this bump maintains the ovi fraction at high level and keeps them against recombination to lower ionic states .
the left border of this spectral feature , formed by to strong absorption of the galactic ionizing radiation in the disc , is irrelevant from the point of view of the maintainance of ovi  even in the absence of any absorption it remains enhanced ( see figure  [ fig - ovi - hi ] ) .
the right edge of the bump is , however , artificial in origin : as a matter of fact , the cut at 136  ev ( 91   ) is a choice of the developers of stellar population codes ( e.g. pegase , starburst99 , galaxev ) .
therefore , in the following we discuss convergence of our results against possible galactic spectral variations at energy higher than 136  ev .
extreme ultraviolet and soft x - ray photons
( @xmath137  ev ) can be produced by post - agb stars before entering the white dwarf cooling phase ( e.g. * ? ? ?
during this very bright evolution phase lasting @xmath138  yrs the effective temperature may reach more than @xmath139  k , while the surface gravity varies as @xmath140 to @xmath141  cm  s@xmath133 .
figure  [ fig - spec - z02 ] ( upper panel ) presents the two spectra for @xmath142  k and @xmath129  k for solar metallicity taken from the stellar library @xcite .
it is seen that the amount of photons with @xmath143  ev increases when the surface gravity @xmath144 decreases from @xmath131 to @xmath132  cm  s@xmath133 for @xmath145  k. on the other hand even a small decrease of effective temperature results in a shortage of photons with @xmath146  ev independent on @xmath144 .
the majority of observed post - agb stars is known to have masses below @xmath147 and temperatures @xmath139 to @xmath129  k ( e.g. * ? ? ?
* ) , thus falling into the range shown here . for several
spectra a break around @xmath148  ev is clearly seen .
one can think that in order to conservatively account this break in stellar population codes the value 136  ev ( 91   ) was choosen as an upper limit of energy in their stellar library .
it is clearly seen though that even the spectra with a break at @xmath148  ev continue up to @xmath149  ev and do not show such a strong decrease of flux at @xmath143  ev as used in our calculations .
one can think that this gradually decreasing flux above @xmath148  ev produces qualitatively similar effects as does the flux with the break at @xmath143  ev .     in order to check this tentative conclusion we performed calculations of ionic composition for several models of the ionizing flux above 136  ev
we construct several spectra and test the appearence of enhanced ovi fraction .
we smooth the spectra calculated in pegase code at energies @xmath150  ev and match them at @xmath151  ev to the spectrum obtained from the `` @xmath43 '' relation as described in section 2.3 ( see * ? ? ?
we proceed in the manner that the resulting spectrum would be qualitatively close to the spectra depicted in the upper panel .
figure  [ fig - spec - z02 ] ( middle panel ) shows these spectra normalized on the spectral luminosity in model a. we calculate the ionization and thermal evolution of gas exposed to a cumulative spectrum which includes the extragalactic and the galactic radiation field .
figure  [ fig - spec - z02 ] ( bottom panel ) presents an example of the cumulative ionizing background flux at a distance @xmath76  kpc from the galaxy in models a , a+1 .. a+5 located at a redshift @xmath78 ( thick lines ) . figure  [ fig - ovi - spec ] shows evolution of ovi fraction in gas exposed to the ionizing flux as in models a , a+1 .. a+5 .
it is readily seen that in all a+ models the ovi fraction is greater than 0.1 at redshifts @xmath152 .
the ovi fraction is several times lower than in model a , but this is still obviously higher than @xmath153 , which establishes in gas exposed only to the extragalactic radiation ( compare to the dot - dot line for log  n(hi)=20.5 in lower panel of figure  [ fig - ovi - hi ] ) . in models a+1 and a+2 oxygen at @xmath154 is mainly locked in ovii ions .
starting from @xmath104 the ionization production rate of ovii ions decreases such that recombinations become efficient to replenish lower ionization states , and as a result the fraction of ovi increases at @xmath155 substantially : from @xmath156 at @xmath157 to @xmath158 in a+1 and @xmath159 in a+2 at @xmath152 .
this balance between ovii and ovi is kept then quasi - steady down to @xmath160 . at large distances
the contribution from galactic photons decreases resulting in a higher fraction of ovi in models a+1 and a+2 .
note that contrary to models a+1 and a+2 the flux in model a+5 is low and oxygen accumulates mostly in oiv - ov states .
higher galactic luminosities presented in a set of models c and c+3 .. c+5 , result in a higher ovi fraction : @xmath161 in gas located at @xmath76  kpc at @xmath152 . however , in models c+1 and c+2 the ionization flux at @xmath146 ev is excessive and the fraction of ovi decreases to 0.03 and 0.07 at @xmath78 , and reaches only 0.05 and 0.1 , at @xmath28 . at larger distances from a galaxy where the flux drops models c+1 and c+2 show an increase of ovi fraction : e.g. , in model c+1 it increases to @xmath162 ( 0.11 ) at 150  kpc and @xmath163 ( 0.17 ) at 200  kpc at @xmath78 ( @xmath28 ) .
it is readily seen thus that in the majority of these models ovi fraction remains high and reaches @xmath164 at distance @xmath80  kpc from a galaxy located at @xmath152 .
this is more than an order of magnitude higher than can be reached in gas exposed only to the extragalactic radiation at @xmath152 . in model
c+1 the ovi fraction reaches 0.1 only for @xmath165  kpc ( @xmath152 ) , while at closer distances most of oxygen is confined into ovii .
it is worth stressing though that the model c+1 can be considered as a model with an extreme possible contribution from post - agb stars , in the sense that models with harder spectrum are unlikely ( e.g. * ? ? ?
* ) , thus one can expect that within conservative models of galactic x - ray spectrum the effect of enhanced ovi remains stable .      in the previous section we have found that ovi fraction can reach high values of @xmath166 under the action of ionizing radiation from the underlying galactic stellar population along with the extragalactic radiation field .
such ovi fraction is at least an order higher than the maximum @xmath163 reached in commonly used models with gas ionized collisionally and/or by the extragalactic radiation . assuming the spherical symmetry and using the time - dependent ovi radial distribution around a host galaxy one can find the ovi column densities along the line of sight crossing the galactic halo at impact parameter  @xmath167 .
we integrate from @xmath168  kpc to 300  kpc .
@xmath7 , the middle and lower panels show the column density for the isobaric ( ib ) regime with the initial density @xmath169  @xmath7 and @xmath170  @xmath7 , respectively ; the metallicity is @xmath53 .
open rhombs and squares depict the observed column densities in the passive and starforming galaxies , correspondingly @xcite ; integration is over @xmath168  kpc to 300  kpc , different symbols correspond to different redshifts . , width=302 ]    , but for gas metallicity @xmath111 .
, width=302 ]     figure  [ figcd_ovi ] ( upper panel ) presents the dependence of @xmath171 on impact parameter @xmath167 in a @xmath172 isochoric gas within the models a and c. the ovi column density ranges from @xmath173 to @xmath174  @xmath3 at the impact parameter @xmath175  kpc  a factor of 3 - 10 higher than observed by @xcite .
therefore circumgalactic gas of an order of magnitude lower metallicity @xmath176 would provide good agreement with observations as shown in figure  [ figcd_ovi2 ] ( upper panel ) .
it is worth noting though that spatial distribution of metals in circumgalactic environment is highly nonhomogeneous ( e.g. * ? ? ?
* ; * ? ? ?
* ) , and it is most likely for metals to be locked in clumps of smaller covering factor resulting in a proportional decrease of the column density .
figure  [ figcd_ovi ] ( middle and lower panels ) show the column density of ovi in isobaric gas with the initial density @xmath169 and @xmath170  @xmath7 , respectively .
isobaric regime is supported by radiation losses and heating from ionizing radiation .
photo - heating grows at large distances from the galaxy , @xmath177  kpc where contribution from ovi ions to cooling decreases . at lower distances
ovi fraction increases and it becomes a dominant coolant resulting in a rather efficient cooling .
consequently , column densities grow under isobaric compression : @xmath171 can reach @xmath178  @xmath3 for the initial density @xmath169  @xmath7 , and @xmath179  @xmath3 for @xmath170  @xmath7 . similarly to the isochoric case gas with lower metallicity
( 0.01  @xmath180 ) fits observational data better .
such high column densities are reached for galactic spectra with a cut at 91    ( see footnote in section 3.1 ) . in reality though galactic spectra do not show such sharp breaks . in section 3.3
we considered galactic spectra approximated smoothly around @xmath181  ev to account possible variations of contribution from post - agb stars .
we have found that in the low redshift range , @xmath152 , ovi fraction in gas exposed to such approximated spectra remains sufficiently high in the majority of models considered here . in all cases
we have a considerable ( an order of magnitude ) excess of the column densities caluculated in models a , c and with respect to the observed ones as seen in figure  [ figcd_ovi ] ) .
@xcite have argued that the circumgalactic medium can be a significant reservoir of the ejected material from galaxies . in order to explain large column densities of ovi observed in the haloes they have had to assume that the circumgalactic gas has nearly solar metallicity .
this assumption meets though difficulties because the minimum oxygen mass in the halo obtained within such an assumption reaches around 10 - 70% of the total oxygen mass in the ism @xcite , which requires in turn unprecedentedly high mass exchange between the galactic ism disk and a huge circumgalactic reservoir extending up to 150  kpc .
higher ovi fraction obtained in our model reduces the estimates to a more reasonable level , and as a consequence weakens constrains on the sources of oxygen in underlying galaxies .
contrary , an order of magnitude lower metallicity results in a more consistent scenario of metal production and mass exchange between galactic and circumgalactic gas .
indeed , the total mass of the circumgalactic gas in ovi haloes within @xmath182  kpc is about @xmath183 independent on redshift as estimated from hi data @xcite .
if such gas would be enriched to the metallicity as high as @xmath53 , then the oxygen mass in it equals @xmath184 .
this amount can be produced during a 0.2  gyr period in a galaxy with sfr @xmath185 ( as in our model a ) , with 0.015  @xmath26 of oxygen returned to the ism per each 1  @xmath26 star formed ( e.g. , * ? ? ? * ) .
this timescale is significantly shorter than the whole period of the evolution
several gyrs , and comparable to the initial period of evolution in our models with the sfr kept nearly constant ( upper panel on figure  [ figssfr - time ] ) .
because of the high sfr during initial several hundreds myrs , gas can be ejected from the disk with velocities @xmath186  km / s sufficient to reach @xmath75100  kpc in approximately 1  gyr .
we assumed throughout that metals are honogeneously distributed in a spherical layer between 50 and 300  kpc .
this can be an overestimate , because in general the volume occupied by ejected metals depends on the galaxy mass and for dwarf galaxies can be smaller ( e.g. * ? ? ?
moreover , metals are locked most likely into small - size dense fragments ( e.g. * ? ? ?
* ) , where radiative cooling is efficient .
it can result in a rapid cooling and compression and make such fragments to escape detection .
dwarf satellites may add metals in the halo , though their contribution is apparently small because of low metallicity in their ism : @xmath187 ( e.g. , * ? ? ?
* ; * ? ? ?
* ) and massive dwarfs are rare ( e.g. , * ? ? ?
* ) .
in extended galactic haloes gas with a low density , @xmath188  @xmath7 , is exposed to rather a high ionizing background consisting of galactic and extragalactic components . under such conditions the ionic state may reach photoequilibrium .
ionic composition under photoequilibrium is calculated with making use the cloudy .
we assume the time - dependent ionizing background described in section 2.3 .
the difference between our calculations and the ones performed by @xcite , is that besides the haardt - madau extragalactic radiation field used in @xcite we add the galactic contribution , which as we showed above may have crucial consequences : the excess of photons with @xmath189  ev from the underlying stellar population competes the extragalactic radiation , while x - ray photons in the excess with @xmath136  ev substantially enhance ovi fraction ( section 3.2 ) .     for the fiducial column densities
@xmath190  @xmath3 and @xmath191  @xmath3 ovi fraction in photoequilibrium is several magnitudes smaller than in our time - dependent models . in photoequilibrium oxygen
is mainly locked in less ionized states , oiv - ov , and ovi column density drops to @xmath192  @xmath3 .
when @xmath82 and @xmath83 decrease the fraction of ovi instead increases . even for @xmath193  @xmath3 ( here @xmath194 is assumed ) the ovi column density reaches @xmath195  @xmath3 at the impact distances @xmath196  kpc and @xmath118 .
further decrease of @xmath82 results in ovi fraction to grow up to @xmath197 , and consequently in higher ovi column densities .
figure  [ figcd_ovi - eq ] presents the dependence of the ovi column density on the impact parameter for photoequilibrium model with @xmath198  @xmath3 and time - dependent model with @xmath199  @xmath3 .
the clearly seen considerable difference between photoequilibrium and time - dependent models originates from a well - known fact that gas under nonequilibrium conditions is overionized with respect to what occurs under equilibrium ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
when @xmath82 decreases , the absorption of galactic ionizing radiation falls down and the peak in the spectrum around @xmath200  ev grows . at such conditions photoionization time shortens , @xmath201 , and the ionic composition approaches photoequilibrium . at low @xmath126 the extragalactic
flux decreases : for example , for @xmath118 @xmath202  erg  s@xmath203  @xmath3  hz@xmath203 at @xmath204  ev ( figure  [ fig - ovi - hi ] ) . as a result
the phoionization timescale @xmath205 of ovi ions increases up to several hundreds myr , here @xmath206  @xmath3 and @xmath207 are assumed ; @xmath208  ev is the ovi ionization potential and @xmath209 is the width of the peak between 113 - 136  ev .
this is about several times shorter than the recombination time for gas with @xmath210  @xmath7 making nonequilibrium effects important .
smaller absorption in the galactic disc provides higher galactic ionizing flux and shorter photoionization timescale such that ionic composition shifts to photoequilibrium .
thus , nonequilibrium effects are important for ovi ionization kinetics in the circumgalactic gas with @xmath211  @xmath7 at @xmath5 .
t is therefore can be concluded that the main factor which provides high ovi column densities in massive starforming galaxies at low and moderate metallicities is the excess of photons with energies 113 - 130  ev from the stellar population of a host galaxy .
@xmath3 , whereas in time - dependent model @xmath212
@xmath3 is adopted ( @xmath213 is assumed ) .
thin and thick lines correspond to nonequilibrium and photoequilibrium , respectively ; density and metallicity are @xmath210 @xmath7 and @xmath53 , the other notations are as in figure  [ figcd_ovi ] . , width=302 ]
in this paper we have presented the photoequilibrium and nonequilibrium ( time - dependent ) ionization and thermal state of circumgalactic gas located at distances up to @xmath16  kpc around starforming galaxies , and exposed to both extragalactic and galactic time - dependent ionizing background . for the extragalactic background we considered the spectra obtained by @xcite . using the pegase code @xcite
we have calculated chemical and spectro - photometric evolution of galaxies , and have chosen the two models , whose specific star formation rate ( @xmath214 ) and stellar masses are close to the starforming galaxies with large ovi column densities observed in @xcite .
we have found that    * the maximum ovi fraction can reach @xmath0 under physical conditions ( gas density and metallicity , and the spectrum shape ) , which are typical in haloes of starforming galaxies ; such a high ovi fraction is due to galactic photons with @xmath1  ev ; the effect of enhanced ovi remains stable within conservative models of galactic x - ray spectrum fluctuations at @xmath215  ev ; * due to such high fraction of ovi its column density ranges in @xmath216  @xmath3 even for a low metallicity @xmath217 , and @xmath218  @xmath3 for @xmath219 at impact parameters up to @xmath220  kpc ; this results in several times more conservative estimate of the oxygen mass in haloes compared to @xmath221 @xcite .
* we have shown therefore that the large ovi column densities observed in haloes of starforming galaxies @xcite can be found in circumgalactic conditions with nearly 0.01 - 0.1 of solar metallicity , and correspondingly the requirements to the sources of oxygen in the extended haloes become reasonably conservative .
high ovi column densities in haloes of starforming galaxies can emerge under photoequilibrium as well as under nonequilibrium conditions .
the main source is a high radiation flux of photons with @xmath1  ev from the stellar population of starforming galaxies .
nonequilibrium effects for ovi ionization kinetics are important in the circumgalactic environment with @xmath211  @xmath7 at @xmath5 .
very recently @xcite have reported about high ovi column densities in the circumgalactic medium at @xmath54 with @xmath171 reaching @xmath222  @xmath3 ; the observational sample includes absorbers of different type : from lyman limit to damped ly@xmath135 systems .
from our point of view such a high ovi column density can be due to the excess of galactic ionizing photons with @xmath223  ev over the extragalactic background , while the metallicity might be rather low .
we thank jason tumlinson for providing data and discussion , and the anonymous referee for valuable comments and pointing to a mistake . this work is supported by the rfbr through the grants 12 - 02 - 00365 , 12 - 02 - 00917 , 12 - 02 - 92704 , and by the russian federal task program `` research and operations on priority directions of development of the science and technology complex of russia for 2009 - 2013 '' ( state contracts 14.a18.21.1304 , 2.5641.2011 and 14.a18.21.1179 ) .
ev is grateful for support from the `` dynasty '' foundation .
asplund , m. , grevesse , n. , & sauval , a. j. 2005 , in asp conf .
336 , cosmic abundances as records of stellar evolution and nucleosynthesis , ed .
t. g. barnes iii & f. n. bash ( san francisco : asp ) , 25