hii galaxies can be considered as the strong emission line subset of the blue compact dwarf galaxies ( bcd ) .
their optical emission is dominated by strong and narrow emission lines produced by the interstellar gas ionized by a young and luminous cluster .
their integrated optical spectra are indistinguishable from a normal hii region .
hii galaxies are gas rich and metal poor .
these facts had led to the early proposal that these systems may be very young perhaps undergoing their first burst of star formation ( sf ) @xcite . however , many authors have found clear evidences of the existence of a low surface brightness , old , non - ionizing stellar population in the majority of these galaxies , as explained by ( * ? ? ?
* here in after mmdt ) and references therein .
the existence of an old population plus the relative paucity of hii galaxies with extreme low metal content are consistent with a scenario in which hii galaxies are suffering at present a starburst , instead of being old systems with a slow evolution .
therefore we use the term _ starburst _ to denote a violent episode of star formation where a large number of massive stars have been formed in a small volume of space and over a time scale of a few million years .
galaxies called _ starbursts galaxies _ are a very heterogeneous category , including from bcds to ulirgs .
we will handle here only hii galaxies , those dwarf galaxies with active star formation whose effects actually dominates their uv - optical emission . in these galaxies
the starburst is considered a phase , a process of strong star formation in terms of intensity and duration .
nevertheless , the fundamental details of how the star formation history proceeds in hii galaxies is still an unresolved problem .    until now
, three basic evolutionary scenarios have been postulated :    * * bursting star formation * : the stars form in short but intense episodes separated by long quiescent periods of very low or null activity @xcite * * gasping star formation * : the star formation takes places as long episodes of sf of moderate intensity separated by short quiescent periods @xcite . * * continuous star formation * : the process of stellar formation is continuous and of low intensity during the galaxy life , with superimposed sporadic bursts @xcite .
the low metallicities , lack of dust and the optical colors argue in favour of a bursting star formation with long quiescent periods ( e.g. * ? ? ?
however , inactive star forming periods seem to rarely occur , and inter - burst stages characterized by low - levels of star formation are perhaps more plausible than its complete cessation @xcite . if this is the case , a continuous low level state of star formation would dominate during the vast majority of the time during which many or perhaps most of the stars in bcd galaxies would be formed .
chemical evolution models with gasping or bursting star formation modes seem to be appropriate to describe local dwarf irregulars ( dirr ) @xcite , but the chemical evidence alone can not differentiate between gasping and bursting scenarios . the gasping scenario can be viewed as a transition between the continuous star formation and the bursting scenario @xcite .
the most remarkable difference between gasping and bursting scenarios is that short episodes of star formation enrich the ism in a time scale of few tens of myr while long - lasting episodes enrich gradually the ism in a longer time scale , and any further episode of sf does not leave an appreciable imprint on the chemical evolution .
the abundances , which increase rapidly in bursting models , change more slowly and during a longer time in gasping models , but with older ages in their stellar populations @xcite .
this picture , however , does not entirely exclude the possibility that star formation histories appear to be composed by burst cycles when examined with a high time resolution , but looks constant when averaged over longer time scales .
it does not exclude either the possibility that the sf propagates through the galaxy , taking place in independent luminous short lived hii regions : `` the short timescales associated to starbursts in dwarf galaxies may be understood as _
flickering _ events , small components of a larger starburst in the galaxy '' @xcite    these observations of stellar populations suggest that it is most likely that the star formation in dwarf galaxies is sporadic , separated by millions to billions of years , even in isolated systems , as most of hii are .
if the sf of hii galaxies is dominated by intermittent starburst episodes separated by long inactive time spans , quiescent blue dwarf galaxies without the dominant starburst and showing similar properties to those with such intense star formation events , should be relatively common .
support to this framework comes from the fact that blue compact dwarfs ( bcd ) , with emission lines on average much weaker than these of hii galaxies , are more common than those ones .
following this idea @xcite obtain a relationship between the duration time of the starburst phase and the quiescent periods , based on the number of objects of each type of their sample .
they find that , if the duration of the burst phase is 10myr , the time in quiescence must be at least 0.27gyr , implying several mayor starburst episodes along the life of a bcd galaxy .
therefore , the diversity in properties exhibited by dwarfs may be at least partially the consequence of observing them at different times in their star formation cycles , and the frequency distributions of galaxies in the various phases could correspond to the time spent in each of them
.    it would be interesting to see which is the most probable star formation scenario given by the _ ab initio _ cosmological hydrodynamical simulations .
however this is difficult since the dwarf galaxies produced in the existing works are usually poor gas , dwarf elliptical ( de ) or dwarf spheroidals , ( dsp ) objects @xcite .
other simulations , such as @xcite used fixed initial conditions for stellar and gas masses , and this way the mass to halo mass ratio is not a result of the simulation . in @xcite simulations , the ultraviolet ( uv ) background is not included , which probably contributes to a high star formation efficiency and to large final stellar masses .
@xcite follow the evolution of a galaxy with a halo mass of 10@xmath0m@xmath1 but only until z=5 .
@xcite have calculated some complete hydrodynamical simulations , in which baryonic processes , as gas cooling , heating from the cosmic uv field , star formation and supernova - driven gas heating , are included with sufficient spatial resolution , with clumps as low as 10@xmath2m@xmath1 resolved .
the created dwarf galaxies without bulges are similar to the observed dwarf irregular galaxies . however , besides the authors do not show the star formation histories , these galaxies show a maximum star formation rate ( sfr ) of 0.25m@xmath3 and a present sfr of 0.01m@xmath3 , much lower than observational values for hii galaxies , and not comparable to these objects that we will try to model here .
only @xcite has recently presented some star formation histories , obtained from cosmological hydrodynamical simulations for dwarf galaxies , with enough spatial resolution for gravitational masses of  10@xmath0m@xmath1 .
their star formation histories seem to be sporadic , and the stars continue forming sporadically even at late times . however , the simulations are still having some problems to reproduce the adequate number of galaxies and the author claims that more work is still necessary . as @xcite explains , all dwarf galaxies (  10@xmath4m@xmath1 ) formed in the current hydrodynamical simulations are more than an order of magnitude more luminous than expected for these masses . in any case
they are much more massive than bcds and hii galaxies that we want modeling .    for what refers to other models existing in the literature ,
there are not models for the study of bcd or hii galaxies in the way we want to use them .
a number of models have computed purely chemical evolution models for bcds or dwarf irregular galaxies ( * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* and many others ) .
most of them assume that the star formation occurs in bursts and include the effects of galactic winds and/or gas infall .
however , they limit the study to the evolution of n and o abundances and/or to the luminosity - metallicity relation . @xcite and @xcite used the information coming from the chemical evolution models from @xcite to perform the next step and combine chemical and spectral evolution for irregular galaxies .
these models , however , exclude the early stages of evolution , i.e. during the nebular phase when most massive stars dominate the energy output budget . @xcite or @xcite
compute chemical and photometric evolution models in a way more or less consistent for starburst , spiral , early types or low surface brightness galaxies , respectively , any of them applied to hii galaxies .
other works , in turn , are focused on the ionized gas properties and make models using single stellar populations ( ssps ) spectral energy distributions ( seds ) of a given metallicity for applying a photo - ionization code to obtain the emission lines and studying diagnostic diagrams and/or abundances obtained by empirical calibrations @xcite .
most of them ignore the star formation history of the galaxy and not take into account the underlying stellar population , so these studies are only valid for the study of the current stellar generations of the galaxy . summarizing a code combining chemical ,
evolutionary synthesis and photo - ionization models has not been applied for the analysis of bcds and hii galaxies , with the exception of the ours ( mmdt ) , but only to elliptical massive galaxies @xcite .
the aim of this work is to critically analyze the possibility of dwarf galaxies undergo recurrent phases of star formation by comparing the predictions of our self - consistent evolutionary models with the most evolution sensitive observed parameters , i.e. stellar continuum colors , emission line equivalent widths and chemical composition in a large sample of bcd galaxies .
we have used an updated version of our models from mmdt , developed to predict the main characteristics of hii galaxies , such as the emission lines , continuum color , continuum plus the contribution of emission lines colors , equivalent widths and chemical abundances .
our code uses the chemical evolution model results for the computation of the seds , which are , in turn , used as the ionizing sources for a photo - ionization code . in mmdt
we studied the viability of this tool , showing its adequacy in the reproduction of the data defined by the stellar populations ( colors , equivalent width of absorption lines or spectral indices , spectral energy distributions .. ) , which define the time evolution , as well as the characteristics of the gas phase ( emission lines , equivalent widths , elemental abundances , gas densities ) , which defines the present time state of the galaxy .
this is done in a self - consistent way , that is , using the same assumptions regarding stellar evolution , model stellar atmospheres and nucleosynthesis , and using a realistic age - metallicity relation .
one of the most important results obtained from our previous work was that observational data are reproduced only if the mass involved in the last burst is much smaller than the mass of the old underlying stellar population .
now , once checked the possibilities of our tool , we try to wide the number of models , changing the star formation scenarios to see if data are reproduced with more or less success for some of them .
it is a widely assumed idea that starburst galaxies are strongly affected by gas infall and outflow and that chemical abundances and star formation can not modeled without these issues .
however , the existence of outflows in dwarf starburst galaxies is by no means a settled issue ; while winds able to escape the galaxy have been found in two prototypical starburst galaxies like ngc1569 and ngc1705 , recent works about mass loss , for example @xcite , find that ionized and neutral hydrogen expansion velocities measured are , in all cases , too low to allow the gas to escape from the gravitational potential of their studied galaxies some of which , like ngc2363 or ngc4861 , can also be considered as _
prototypical_. in fact , according to @xcite:while the observational support for the presence of galactic winds in massive galaxies and gas - rich mergers is quite strong , the case for galactic winds in dwarf galaxies is much weaker  .
the need of winds is also related with the low oxygen abundance found for dwarf galaxies for their gas fractions , and , in this case , both infall and outflow can help .
however , galaxies loosing a large fraction of their gas should show red colors , contrary to what is observed .
in fact , a good fraction of the galaxies in the sample of @xcite can be reproduced with closed box models .
furthermore , it is necessary to take into account the n / o vs o / h relationship for dwarf galaxies .
selective winds models as @xcite have recognized problems in reproducing the observed n / o ratios , and the work by @xcite demonstrates that the observed n / o ratios , as well as their dispersion , can be reproduced with closed box models while selective winds can be ruled out .
actually , @xcite showed that it is mainly the combination of the different time scales in nitrogen production by stars of different masses and the galaxy star formation histories as proceeding from different star formation efficiencies , what establishes the n / o ratio .
therefore , since there is enough proof to conclude that infall and/or outflows are not necessary to reproduce the general trends of hii galaxies , galactic supernova - driven winds are not included in our models .
this scenario is also in agreement with models from @xcite , who find that scaling relations of dwarf galaxies may be reproduced by simulated galaxies without supernova - driven outflows .    for the present work
, new and updated theoretical codes have been applied for the computation of the models .
the main changes correspond to the evolutionary synthesis models , now from ( * ? ? ?
* hereinafter mgvb09 ) instead @xcite .
as explained before , to guarantee the self - consistency of the approach we have used the same assumptions in the stellar evolution , model stellar atmospheres and nucleosynthesis parts of the code .
the use of these new models allows to be more consistent than in our previous work mmdt , since there the imf used in the chemical evolution models @xcite and in evolutionary synthesis models @xcite was not the same .
the new models mgvb allow to us use the same imf for chemical and evolutionary synthesis ( photometrical ) calculations .
moreover , with the previous models from @xcite we had also a problem for the lowest metallicities stellar populations since models for the youngest ages of low metallicities were not available .
now we have models for the same number of ages for all metallicities . to include the lowest metallicities of the youngest ages change appreciably the colors of the stellar populations and
allows to calculate emission lines for very low metallicity regions , as it is shown in @xcite .
the contribution of the emission lines is quite different than expected from the simple extrapolation from other metallicities , also modifying the final results of colors .
the data shown in the graphs have been extracted mainly from two main sources .
first , the compilation from @xcite which provides emission line measurements , corrected for extinction , published for local hii galaxies .
the sample comprises 450 objects and constitutes a large sample of local hii galaxies with good - quality spectroscopic data .
the sample is rather inhomogeneous in nature , since the data proceed from different instrumental setups , observing conditions and reduction procedures , but have been analysed in a uniform way .
data for these sample objects include the emission line intensities of : [ oii]@xmath5 3727,29  [ oiii]@xmath5 4959,5007  and [ nii]@xmath5 6548,84  all of them relative to h@xmath6 , and the equivalent widths of the [ oii ] and [ oiii ] emission lines  ew([oii ] ) ew([oiii ] )  , and the h@xmath6 line , ew(h@xmath6 ) . as a second source
, we have used the metal poor galaxy data from the data release 3 of sloan digital sky survey , taken from @xcite .
the sloan digital sky survey @xcite constitutes a large data base of galaxies with well defined selection criteria and observed in a homogeneous way .
the sdss dr3 @xcite provides spectra in the wavelength range from 3800 to 9300   for @xmath7 530000 galaxies , quasars and stars . @xcite
extracted @xmath7 2700 spectra of non - active galaxies with the [ oiii]@xmath84363   emission detected above 1@xmath9 level .
this initial sample was further restricted to the objects with an observed flux in the h@xmath6 emission line larger than 10@xmath10 erg s@xmath11 @xmath12 and for which accurate abundances could be derived .
they have also excluded all galaxies with both [ oiii]@xmath84959/h@xmath6 @xmath13 0.7 and [ oii]@xmath83727/h@xmath6 @xmath14 1.0 . applying all these selection criteria
, they obtain a sample of @xmath7 310 sdss objects .
data for these sample objects include the emission line intensities of : [ oiii]@xmath84959,5007   and [ nii]@xmath86584 relative to h@xmath6 and the equivalent width of h@xmath6 .
they also include the intensity of the [ oii ] @xmath5 3727,29 emission line for the lowest redshift objects .    in the next section an explanation of the theoretical star - bursting models
the section 3 presents the results , the meaning of each input parameter and the implication of their variation over the results of the models .
we will discuss the influence of these parameters in the reproduction of the observable characteristics of star - forming galaxies , their impact over possible star - formation scenarios and the connection among different type of dwarf galaxies in section 4 . finally ,
a summary and the conclusions of this work are presented in section 5 .
all tables of the models will be available in electronic format .
the viability of our model technique to reproduce the observable characteristics of hii galaxies was discussed in mmdt : in a first stage the chemical evolution model is computed by using certain input parameters to obtain the star formation history and the gas and stars chemical abundances . secondly , the evolutionary population synthesis code is applied to obtain the seds corresponding to each time step of the chemical evolution . finally , the ionizing part of the sed and the resulting abundances are used as an input to the photo - ionization code to compute the time evolution of the emission lines of the ionized gas .
the chemical evolution is computed with a simplified version of the classical _ multiphase chemical evolution model _ from @xcite . in our version of the code
there is only one region ( without the two zones halo and disk ) with a given mass of gas which form stars ( in this case we do not consider the phase of formation of molecular clouds as in the classical multiphase models ) .
we have run models considering the star formation as a set of successive bursts followed by quiescent periods in a region with a total mass of gas of @xmath15 . in each
burst a given amount of gas is consumed to form stars , and this process comes defined by the star formation efficiency .
the code solves the chemical evolution equations to obtain , in each time step , the abundances of 15 elements : h , d , @xmath16 , @xmath17 , c , @xmath18 , o , n , ne , mg , si , s , ca , fe , and _ nr _ ( where _ nr _ are the isotopes of the neutron rich elements , synthesized from @xmath19 , @xmath18 , @xmath20n and @xmath21 inside the co core ) .
the stellar yields are those from @xcite for massive stars , @xmath22 , and those from @xcite for low and intermediate mass stars .
the supernova ia yields used proceed from @xcite . the initial mass function ( imf ) is taken from @xcite with a range between 0.15 and 100 @xmath23 . more recent imfs such as @xcite or @xcite there exist but , for chemical evolution models , it is necessary to use a combination imf+stellar yields sets calibrated with milky way galaxy data , that is , able to reproduce the elemental abundances and gas and star densities as observed .
this calibration is not still done for those imfs , although an update of models taking these imfs and some new stellar yield sets will be shown in a next future ( molla et al . in preparation ) .
we have taken time steps of @xmath24t@xmath7 0.7myr from the initial time , @xmath25 , up to the final one , @xmath26gyr . at each time step , the star formation rate and the mass in each phase
 low mass , massive stars and remnants , total mass in stars created , and mass of gas are also obtained .
the seds are taken from the popstar evolutionary synthesis models by mgvb09 .
isochrones are an update from those from @xcite for 6 different metallicities : @xmath27 , 0.0004 , 0.004 , 0.008 , 0.02 and 0.05 .
the very low metallicity model of @xmath28 had not been included before in similar works .
the age coverage is from @xmath295.00 to 10.30 with a variable time resolution of @xmath30 in the youngest stellar ages .
we have used the ferrini imf results @xcite with mass limits between 0.15 and 100m@xmath1 in order to avoid any inconsistency with the chemical evolution code , which also uses this imf .    to calculate the sed integrated over the whole history of the galaxy , seds of ssps with the corresponding metallicity and age must be convolved with the star formation history ( sfh ) : @xmath31 where @xmath32 is the age of the stellar population created in a time @xmath33 and @xmath34 is the sed for each ssp of age @xmath35 and metallicity z reached in that time @xmath33 .
a sed from the ssp library , @xmath34 , must be assigned to each time step according to its corresponding age and metallicity .
taking into account the sfh , @xmath36 , and the age - metallicity relation , @xmath37 , obtained from the chemical evolution model , we know the age and metallicity assigned to each time step or stellar generation .
for each one of them we have chosen the sed closest in age among those available in the grid of popstar . however ,
in our models the metallicity changes continuously while the available seds of the library are computed only for 6 possible values .
therefore , we have interpolated logarithmically between the two ssp of the same age @xmath35 and closest in metallicities to @xmath37 to obtain the corresponding @xmath38 .
the final result is the total luminosity at each wavelength @xmath8 , l@xmath39 ) , corresponding to the whole stellar population , including the ionizing continuum proceeding from the last formed stellar population .      to compute the photo - ionization models ,
we have used cloudy ( version c06.02 * ? ? ?
the gas is assumed to be spherically symmetric around a point source of radiation and the pressure or density in the gas is imposed by external conditions .
this allows a plane - parallel geometry treatment in which the gas may be regarded as a thin shell .
a closed geometry has been taken for the calculations .
all the photons which escape from the illuminated face of the cloud towards the star , go on to strike the other side of the nebula , ensuring the case b of recombination and the approximation on the spot .
the number of ionizing photons , q(h ) , striking the illuminated face of the cloud , have been calculated directly from the total resulting sed of the models .
we derive the radius of the modelled region from the mechanical energy from massive stars with strong winds ( taken from mgvb09 ) , instead of using the radius required to maintain a constant density of stellar mass through the successive bursts , as we did in mmdt .
@xcite demonstrated that an early - type star with a strong stellar wind can blow out a large cavity or _
bubble _ in the surrounding gas , if it is assumed to be compressed into a thin spherical shell .
the wind - driven shell begins to evolve with an initial phase of free expansion followed by an adiabatic expansion phase , and then the material collapses into a thin , cold shell as a result of radiative cooling . at this stage
the gas traps the ionization front and the radiative phase begins . in this phase
the ionizing photons are absorbed and the region cools via emission in the balmer lines . in this process
, the radius of the outer shock , r@xmath40 , evolves as : @xmath41 where @xmath42 is the total injected mechanical energy ( sn and stellar winds ) per unit time in units of 10@xmath43ergs s@xmath11 , @xmath44 is the interstellar medium density in units of @xmath45 , and @xmath46 is the age of the shell in units of 10@xmath47yr . since we use the value of the energy at each time - step , this radius represents the instantaneous size of the region obtained by adding the winds and sne from the previous age to this age . therefore l@xmath48 is this energy divided by the time step .
then , the ionized gas is assumed to be located in a thin spherical shell at that distance r@xmath40 from the ionizing source .
this approach has the advantage of eliminating the ionization parameter as a free variable in the models since now it is computed from the physical parameters of the evolving young cluster .
the ionized gas abundances are assumed to be those reached at the end of the starburst previous to the current one .
fifteen element abundances have been introduced in the photo - ionization code : he , c , n , o , ne , na , mg , al , si , s , ar , ca , fe and ni , obtained from the chemical evolution model , except for na , ar and ni which are not computed in the model and are scaled to the solar ratio @xcite .
the models assume that the nebula is ionization bounded and no dust has been included in the chemical evolution models neither in the photo - ionization calculations .
however dust grains mixed with the ionized gas have been partially taken into account , since we have included the depletion in refractory elements ( si , fe , ca , si , mg ) taken from @xcite .
the grains can affect the absorption of the _ uv _ photons and decrease the electronic temperature .
the density has been assumed constant for simplicity and equal to 100@xmath45 , which is appropriate for modeling hii galaxies @xcite and large circumnuclear hii regions @xcite , frequently found around the nuclei of starbursts and agns .
although the constant density hypothesis is probably not realistic , it can be considered representative when the integrated spectrum of the nebula is analyzed .
the shape of the ionizing continuum is defined by the pair of values ( @xmath49(ryd ) , log@xmath50 ) ( eq .
[ sp ] ) obtained from the ionizing spectrum given directly by the evolutionary synthesis code .
taking into account the simplified version used for our chemical evolution models , the free input parameters related with the infall time - scale and the molecular cloud formation disappear .
the total mass of the initial gas is within the region from the initial time .
the star formation rate is zero except during the bursts .
we need to define the intensity of these bursts , and the time elapsed between them .
each model is therefore characterized by three input parameters :    * * 1 - the initial efficiency ( _ @xmath51 _ ) * : it is the amount of gas consumed to form stars in the first burst of star formation , that is , @xmath52 + we present here models computed with 2 values of @xmath51 : 1 .
high efficiency models : the first burst of star formation involves 33 per cent of the total initial mass of gas ( 33@xmath5310@xmath54m@xmath1 ) , 2 .
low efficiency models : the first burst of star formation involves 10 per cent of the total initial mass of gas ( 10@xmath5310@xmath54m@xmath1 ) . * * 2 - attenuation * : the star formation efficiency of the following bursts is attenuated in two different ways : 1 .   by a factor which changes with the number of the burst , _
n _ , according to the expression : + @xmath55 + this attenuation type corresponds to a soft attenuation model ( hereinafter sam ) .
2 .   by a constant factor , _ k@xmath56 _ , according to the expression : + @xmath57 + corresponding to a high attenuation model ( hereinafter ham ) . in this last expression _
k _ is the attenuation factor which , in this case , takes values from 0 to 1 .
the lower the _ k _ value , the stronger the attenuation , the burst being less efficient each time ; the higher the _ k _ value , the weaker the attenuation , and then the stronger the bursts .
we only show models with @xmath58 * * 3 - time between bursts ( * _ @xmath59 t _ * ) * : every burst takes place instantaneously and it is followed by quiet periods , whose duration can change . for this work we have taken @xmath59t= 1.3gyr for the inter - burst time , that is , one burst every 1.3gyr as the generic case , although we will also show some models with @xmath60gyr and @xmath61gyr in order to compare the different results .
the models we show here are a selection of the ones calculated and described more widely in martn - manjn , moll & lpez - snchez , 2012 ( in preparation ) . in that work ,
the results of 20 models are compared with data for generic and particular low mass dwarf galaxies .
we have selected the most relevant of them in order to show the effects of the input parameter variation on the resulting evolutionary history of each galaxy .
the input parameters of these six models are shown in table  [ inputs ] .
.input parameters for the theoretical models . in column 1
it is defined the type of model according to the attenuation : soft ( sam ) or high attenuation ( ham ) ; the second column gives the parameter @xmath62 which defines the attenuation ; column 3 is the star formation initial efficiency @xmath51 , and column 4 shows the time between burst @xmath63 ingyr . [ cols="^,^,^,^,^ " , ]
we will show here the main results obtained with the selected models .
the * efficiency * determines the initial star formation rate and the initial metallicity of the gas .
we have plotted in fig .
[ fig1 ] the sfr for the 6 models .
the first burst is strong in all panels , while the subsequent ones are less intense due to the decrease of the available gas to form stars and the attenuation . in a ) and
b ) the last time computed is @xmath64gyr while in c ) and d ) the final time is around 1.3 and 0.6gyr respectively . in all panels
we show the observational limits as dashed ( green ) lines . taking these limits into account
, the dwarf galaxies may suffer 11 bursts for sam , or 8 - 9 for ham since later bursts show rates lower than observed .
the oxygen abundances for these same models are shown in fig .
[ fig2 ] . in this plot
we can see that the two efficiencies chosen , 33 per cent and 10 per cent , give the upper and lower limits respectively for the observed oxygen abundance range in this type of galaxies .
the models with the same initial efficiency and different attenuation type , are very similar , and if @xmath63 is shorter , the low oxygen abundance limit is reached very quickly .    figs .
[ fig3 ] and [ fig4 ] show diagnostic diagrams involving the ratios of intense emission lines .
the number of the model ( see table 1 ) is given in each panel .
the higher the efficiency , the more ionizing photons produced and the higher ionization parameter , leading to a higher excitation of the gas .
the low efficiency models cover the region of the diagram occupied by the young and less metallic galaxies ( high [ oiii]@xmath55007,4959/h@xmath65 , low [ oii]@xmath83727/h@xmath65 and low [ nii]@xmath86584/h@xmath66 ) , as expected .
the high initial efficiency models reproduce high excitation and high abundance galaxies , with high [ oiii]@xmath55007,4959/h@xmath65 and high [ nii]@xmath86584/h@xmath66 .
the ionization degree is driven by the efficiency of the bursts , not by the attenuation mode or attenuation factor.thus , the models with different attenuation modes and the same initial efficiency do not show noteworthy differences , since it is the current burst of star formation which produces the observed emission lines , and the underlying population does not affect the ionization parameter .
the attenuation sets the sfr of the successive bursts and determines the contribution of the underlying population .
a higher attenuation implies a larger contribution of previous bursts to the total sed .
therefore , it is important to see the attenuation factor effect on the the evolution of the broad - band continuum colors .    in fig .
[ fig5 ] we show the evolution of the modeled continuum colors for each burst of star formation .
the inclusion of nebular emission continuum in the computation of the ssps reddens the colors of very young populations significantly , mainly at low metallicity as explained in mgvb09 .
the ham cases , compared with sam , have a higher contribution from the non ionizing underlying continuum , which makes the colors redder for subsequent bursts and the young and blue population features disappear .
meanwhile , sam maintains blue colors , characteristic of the current burst of star formation .    in fig .
[ fig6 ] the evolution of the colors is shown including the emission lines contribution to the wide band filters .
we have done so taking into account the strongest emission lines that contribute to the color in each broad - band spectral interval at redshift zero . these are mainly [ oii]@xmath53727 in _ u _ , h@xmath65 in _ b _ , [ oiii]@xmath55007,4959 in _ v _ , h@xmath67 in _ r _ and [ siii]@xmath59069,9532 in _ i_. in that case the more intense the last burst , the larger the change in the colors , so the dotted lines are more similar to colors calculated only with the stellar populations than those shown by solid lines for higher efficiencies . on other hand ,
not all bands are equally modified at the same time
.    a fine tuning of the attenuation factor can be obtained from the evolution of the broad - band continuum colors versus the equivalent width of h@xmath65 , ew(h@xmath65 ) .
these are important observations which give also an effective method to uncover the presence of old underlying stellar populations @xcite .
the comparison of models and observations shows that in most star - forming dwarf galaxies , ssps are unable to reproduce the reddest colors shown by the data with the lowest ew(h@xmath65 ) values .
in fact the observed h@xmath65 equivalent width values and colors require a large contribution from previous stellar generations @xcite .
the evolution of ew(h@xmath65 ) _ vs. _ a pseudo - color of the continuum , the intensities of the adjacent continua of [ oii]@xmath83727 and [ oiii]@xmath85007 lines ( similar to u - v ) of our models compared with the data is shown in figure [ fig7 ] . the sam and ham cases for @xmath68gyr are plotted in the two first panels . it can be seen that to reproduce the trend of observational data , shifted to red colors at low values of ew(h@xmath65 ) with respect to the ssp predictions , the contribution to the total continuum of the non - ionizing red population must be much higher than the contribution of the current burst which creates the emission line spectrum .
this trend can not be reproduced just with ssps neither by increasing their metallicity or age separately nor simultaneously ( mmdt ) .
a strong attenuation in the sfr is needed , as can be seen in the center panel , to reproduce the whole range in ew(h@xmath65 ) and continuum colors simultaneously .
the time between bursts ( @xmath59 t ) is a parameter which may also have an effect on the models similar to the attenuation .
the reduction of the time between bursts offsets the effect of increasing the attenuation : the underlying population is less evolved and produces less reddening . in figure
[ fig7 ] , right panel , ham with @xmath69gyr and @xmath70gyr are plotted .
the ew(h@xmath65 ) decreases rapidly while resulting colors are not shifted to the red sufficiently to cover the range shown by the observational sample .
the colors of the models with shorter inter - burst time are more similar to those with soft attenuation ( strong bursts ) , and require an extra reddening to reproduce the needed effects of the underlying non ionizing populations .
however , the ew(h@xmath65 ) decreases more from burst to burst than in the case of a soft attenuation , where it maintains a high value at the beginning of every bursts . in order to reproduce the observed trend
, the inter - burst time must be longer than 100myr .    in brief , to reproduce the observable characteristics of star - forming galaxies we have to adjust the three input parameters :
* the initial efficiency must be between 0.10 and 0.60 in order to produce a first burst which provides oxygen abundances within the observed values .
in fact , efficiencies lower than 10% seems more probable to reproduce most of observations . * most ( perhaps all ) hii galaxies require the contribution of previous stellar generations to explain the observed trends shown by their continuum and emission line properties .
the chemical and spectro - photometrical parameters ( equivalent widths and colors ) obtained by our models reproduce the observed relations if the contribution of the underlying population from previous bursts to the total continuum is higher than the contribution of the current burst of star formation which dominates the observed emission line spectrum .
this implies a history of star formation higher in the past than at present , and , even in that case , attenuated burst along the time .
the fine tuning of observations is obtained by the adjusting of the attenuation factor , once the initial efficiency is fixed .
( see martn - manjn et al .
2012 for details ) * the inter - burst time must be shorter than 1.3gyr and longer than 100myr . with shorter periods than 100myr
the underlying continuum is too blue and its contribution to the color do not reproduce the trend shown by hii galaxies .
three star formation scenario are usually discussed in relation with the evolution of dwarf galaxies , ( a ) burst : short star - formation episodes with large quiescent periods , ( b ) gasp : long moderate star - formation episodes with short quiescent periods and ( c ) continuous or almost continuous star formation with few over - imposed sporadic bursts .
the three cases can be simulated by our models simply with a change of parameters :    * the star - bursting scenario ( a ) is this one shown in the previous sections .
the effects of an instantaneous burst can be seen in the emission lines during 10myr after the burst takes place .
long quiescent periods or , in our case , null star formation periods , have been considered to last more than 1gyr , which would be the minimum age of the underlying population belonging to the previous burst .
bursts occurred before the immediately previous one have very poor contribution to the current total continuum luminosity , and the use of this long time for the inter - burst period produces a similar result to assuming a two - burst model for some parameters , as ew(h@xmath65 ) .
however , differences appear between both scenarios since the ism is enriched by every single burst occurring during the galaxy life time , stellar populations older than 1gyr should be present , and in consequence colors would be shifted to the red . * the gasping star formation scenario ( b ) can be simulated by substantially increasing the attenuation of the burst and reducing the inter - burst time . with the change in these two parameters
we make a smoother and more moderate the star formation .
since the subsequent bursts are closer in time , the quiescent periods are shorter . models with extreme attenuation and inter - burst times shorter than 100myr can be considered a good approximation to gasping models .
this particular choice of parameters has important consequences when compared with the observations . increasing the attenuation a major contribution to the continuum by the previous bursts
is obtained . however ,
if , in addition , the inter - burst time is reduced , the underlying population will become younger and bluer .
this approach would predict the presence of intermediate age stellar populations , younger than 1gyr , even after several star bursts . with larger inter - burst times
only the immediately previous stellar generation contributes substantially to the present continuum , but even older stellar generations would be noticeable if we further reduce the time between bursts .
although complicated , it is also possible to simulate more extended star - burst phases by concatenation of several star - bursts .
however , population synthesis modeling of bcd spectra does not favor extended periods of star formation @xcite .
then , if the starburst phase is short , there should be many of these episodes during the galaxy lifetime . according to @xcite , there should be one bcd phase each 0.3gyr , which agrees with our models .
* under a continuous star formation scenario ( c ) , the sfr must be very low and extended .
if we reduce the inter - burst time to a minimum value and , simultaneously , we reduce the intensity of each burst , we obtain a very low star formation rate , not comparable to a starburst , but similar to a continuous star formation history .
in fact , many of the observable features of star - forming galaxies can be modeled with a young stellar population 3 - 5myr old and most massive stars being essentially coeval . however , as @xcite demonstrated , it is also possible to obtain similar results with a continuous star formation , lasting at least 20myr since the beginning of the star formation .
a low continuous star formation rate can not be neglected , especially in low metallicity galaxies .
otherwise , if the sfr is very low , ( log ( sfr ) @xmath13 -3 approximately , according to * ? ? ?
* models ) , there are not enough ionizing photons to produce emission lines . in this case
we are not reproducing the so - called hii or bcd galaxies , but another type of galaxy or another stage of their evolution .
actually , we may consider that , instead of a different scenario , the low and continuous star formation could be the phase occurring between bursts . in this way we could change the perspective of the problem and study the characteristics of the galaxies for the two phases : the burst and the following @xmath71myr , and the inter - burst , low activity periods , that is , after the first 10myr from the last star formation episode .
in fact this is related with the two different time scales involved in the evolution , and related with different kind of data : on the one hand , a short time - scale in which the emission lines and othe effects of the ionizing stellar populations are observed . on the other hand
there is long time - scale which defines the red color evolution mainly due to the the stellar populations older than 10myr .
since a hii galaxy can be considered just a star - forming phase or a stage in the evolution of a gas rich dwarf galaxy , we can consider the inter - burst time in our models as a quiescent star - forming dwarf ( qbcd ) phase @xcite .
then the hii galaxy may be in the phase 10@xmath72 yr after the star burst , in which its effects are still visible ( emission lines or changes in metallicity due to massive stars ejections ) and the qbcd phase would correspond to the inter - burst periods .
the contribution of the young and old stellar populations in each galaxy may be obtained from the comparison of data coming from the two different time - scales , in particular , from observations as emission lines , proceeding from the ionizing stellar populations , with colors or abundances given by the older than 10myr stars .    as we have already seen in our previous works , mmdt and @xcite , and we also show in fig .
[ fig6 ] , when both types of stellar populations are contributing to the light , as occurs in bcd galaxies , the colors are contaminated by the emission lines , showing in most cases different trends than those standard sequences defined by the stellar populations , either young or old . taking into account this contamination ,
the colors observed are those corresponding to the ionizing continuum of a given galaxy , always during the bcd phase , that is , the colors of the young population which dominates the emitted light . in order to study the colors in the inter - burst periods and to know the properties of the underlying population
, we should further eliminate the continuum arising from the starburst and isolate the properties of the underlying host galaxy .
figures [ fig8 ] , [ fig9 ] and [ fig10 ] show the b - v _ vs. _ v - i , v - i _ vs. _ r - i and g - r _ vs. _ r - z color - color diagrams for our 6 models .
black circles show the stellar continuum colors in all panels , while the colors computed including the contribution by the strongest emission lines are shown as green and red for models 1 and 2 respectively , brown and orange for models 3 and 4 , and blue and magenta for models 5 and 6 . in order to represent these older phases
we have also represented , as colored open squares , the colors corresponding to ssps , as given in mgvb09 , for ages older than 10myr that is , without taking into account the youngest stellar populations phases , and metallicities from z=0.0001 to 0.02 .
we have not computed exactly the colors of our inter - burst periods , but they must fall in this same region of the diagrams .    in the b - v _ vs. _ v - i diagram ( fig .
[ fig8 ] ) the inclusion of the emission lines contribution to the continuum colors shifts the position of the model points almost perpendicularly to the originally computed ones . the location of the points is mainly determined by the contribution of strong [ oiii]@xmath55007,4959 emission lines to the _ v _ band , and , logically , they are closer to the normal stellar populations locus when efficiency is low .
this means that the excursion from locus defined by the main sequence for stellar populations is stronger when the star formation efficiency is higher , leading the points far .
we also show in this figure the data by @xcite as cyan dots with error bars .
this set of observations corresponds to readily observed colors , including both continuum and line emission , and no reddening correction has been applied .
it can be seen that some of the data are impossible to be reproduced by the models which do not take into account the contribution by emission lines , even if some amount of reddening is invoked . on the other hand
the data by @xcite of resolved locations in mrk  370 , shown in the figure as open green triangles are very well reproduced by our continuum colors .
the data of mrk  370 are based on _ ubvri _ broadband and h@xmath66 narrow band observations . in this case
, the authors subtracted the contribution of the underlying continuum from the old stellar population , removing the contribution from emission lines and correcting for extinction , measuring in this way true colors of the young star - forming knots .
this way while open green triangles correspond to the observed mrk  370 colors , including the emission lines contribution and the underlying stellar population , the full green triangles are the stellar population colors after correction of extinction , emission lines and old stellar populations . in this case
we see that some points fall out of the expected region of models which we interpret as probably due to a over - correction with ssps models .
the underlying component has colors redder than the bcd ones , and they are also well reproduced by our models .    in fig .
[ fig9 ] we plot v - i _ vs. _ v - r , with the same data and code from @xcite and @xcite than in the previous figure .
moreover , we have also included the data of bcds and its hosts from @xcite .
this sample consists of 15 bcds for which they observed the total colors .
the data corrected of extinction are the green stars .
when they correct of the emission lines effect , the obtain the colors shown by magenta stars .
these points corresponding to the actual young stellar populations are closer to the region where our models are represented by the same coloured dots as in fig .
they subtracted the contribution of the underlying galaxy by assuming the mean surface brightness of the extensions , or hosts , to be constant and representative of the underlying galaxy within the starburst regions .
summarizing , they give separately the colors of this underlying stellar continuum ( blue stars ) within the starburst , the contributions including the emission lines ( green stars ) and the colors of the bcd stellar continuum ( magenta stars ) , corrected from emission lines , which are close to black dots , corresponding to our continuum bcd phase .
the colors of the isolated host galaxy ( blue stars ) show redder colors than the bcd component , in agreement with qbcds , mrk  370 host and model colors with ages older than 10myr .
again we note some differences with our models that we assign to an over - correction of the emission lines .
finally we represent in fig.[fig10 ] the g - r vs. r - z diagram with the same symbols that before .
the selected qbcds @xcite are assumed to be like bcd host galaxies @xcite .
the colors of the ssps with ages older than 10myr ( colored squares ) can reproduce these observations , indicating that they correspond to the colors of a more evolved population than the one corresponding to the current burst in bcds . however , some of the observed dots lie in the bcd zone , indicating that this data sample also contains a certain proportion of bcds , following a continuous evolutionary sequence . in this color - color diagram
there is a larger difference than in the previous ones between models and data for the bluest points , showing these data higher g - r values than predicted for similar low r - z colors .
probably this is due to the transformation used to calculate sdss colors with our johnson colors , which would be not valid for this type of objects with emission lines contaminating the wide band filter magnitudes .
we see in the three figures that models start to separate off the main sequence of colors at different points following the type of attenuation . in the top panel , there are points with large contribution of lines starting even at the lowest and left corner , where the youngest stellar populations lie .
the same points in medium panel begin to separate at redder colors , this implies that the underlying continuum corresponds to an older age .
these models cover all observational range .
however in the bottom panel the points with the contribution of the lines fall in a region much bluer where there are no observations .
the underlying stellar populations are much younger than the populations of the data .    in any case
it seems that bcds and qbcds colors are well reproduced by the same basic models , during the first 10myr after the beginning of star formation and after 10myr up to more than 1gyr respectively , and therefore we can say that both types of galaxy colors overlap and could be considered as different phases in the evolution of the same object .
historically , the most common scenario assumed for the star formation history of the star - forming galaxies is that stars form in a bursting mode .
this scenario consists in successive instantaneous star formation bursts interposed between long periods of null or very low star formation activity .
if these periods are shorter and the bursts are more extended and moderate , we are under a gasping star formation scenario . on the contrary , if the star formation is low and moderate during the whole life of the galaxy , and only some sporadic bursts take place , a continuous star formation scenario is appealed .
we have analyzed the possibilities of these scenarios by means of theoretical models based on the successive bursts star formation hypothesis .
each galaxy is modeled assuming an initial amount of unprocessed gas of 10@xmath73 m@xmath1 .
the evolution is computed along a total duration of 13.2gyr during which successive star - bursts take place . since any of the possible scenarios of the star formation histories described show evidences of being able to explain the observed data , we have tried to check the hypotheses suggested .
the grid of theoretical models are computed by the combination of three tools : a chemical evolution code , an evolutionary synthesis code and a photo - ionization code , all of them previously calibrated .
they have been used in a self - consistent way , _ i. e. _ taking the same assumptions about stellar evolution and nucleosynthesis , and the resulting metallicity in every time step .
these models have three free input parameters which can be changed to obtain different model results : ( a ) the initial star formation efficiency which determines the initial star formation rate ( sfr ) and the initial metallicity of the gas , ( b ) the attenuation of the successive bursts , which determines the evolution of the gas , keeping metallicity between the observational data limits , and ( c ) the inter - burst time , which sets the age of the non ionizing underlying population .
we have shown that ( a ) also leads the behavior of the ionizing gas since emission lines are determined by the number of massive stars and the metallicity of the emitting gas . with this parameter we can control the number of old pristine stars that we are going to find in the next star formation burst .
the strength of the bursts is controlled by ( b ) .
the more attenuated , the less the star formation rate , which favors the previous stellar generations to have more weight on the total current spectrum , hence modifying the equivalent width and the colors of the continuum .
finally , the contribution of the underlying population can be determined by ( c ) .
these periods are the quiescent phases of the bcd .
the shorter the inter - burst time , the younger the underlying stellar populations .
these stars will be part of the host bcd galaxy , that is , all those generations of stars which have not ionizing stars , showing ages of more than 10@xmath72 yr .    the models can reproduce simultaneously the data relative to the current ionizing population and the data which give us information about the evolutionary history of the star formation . in order to reproduce the observational data ,
the parameters are constrained :    * the initial star formation efficiency must be lower than 50% and more probably is around 20% .
* the subsequent star - forming burst must be attenuated , in order to have the adequate proportion of underlying to young / ionizing stellar populations ; and * the inter - burst time must be longer than 100myr but shorter than 1.3gyr .    with our models we are able to compute the stellar continuum colors as well as the broad band colors corrected from the contribution of the more intense emission lines .
this effect shifts the colors almost perpendicularly to the stellar continuum ones , specially during the first bursts of star formation , when the emission lines are more intense .
the contamination of the continuum by the emission lines can make the observed galaxies appear younger than they really are , hiding evidence of underlying old populations .    on the other hand , the different characteristics shown by the star formation histories of our models ,
imply that we may be observing the same objects in different evolutionary stages , that is , just during the burst , in the post - burst phase or in the inter - burst periods .
if we consider that the different star formation scenarios correspond to different phases in the life of a single galaxy , our models are able to reproduce the characteristics of each one of these phases . considering the bcd phase as the first 10myr after a star formation burst , and the quiescent phase as the inter - burst periods , we can reproduce the galaxies both in a high activity and in a post - starburst ( or pre - starburst ) phase .
the main difference between these phases are the existence of emission lines produced by the ionized gas . their intensity , and their contribution to the total continuum color ,
can be modified by the initial efficiency . with the attenuation factor
we can produce a more or less reddened stellar continuum due to the contribution of the host galaxy , the age of these old stars being determined by the inter - burst time .
this work has been partially supported by dgicyt grants aya200767965-c03 - 03 and aya201021887c0402 . also , by the comunidad de madrid under grant s-0505/esp/000237 ( astrocam ) and by the spanish micinn under the consolider - ingenio 2010 program grant csd2006 - 00070
: first science with the gtc ( http://www.iac.es/consolider-ingenio-gtc ) which are acknowledged .
rt is grateful to the mexican research council ( conacyt ) for supporting this research under grants cb-2006 - 49847 , cb-2007 - 01 - 84746 and cb-2008 - 103365-f                                                        gallagher , j.  s. , mould , j.  r. , de feijter , e. , holtzman , j. , stappers , b. , watson , a. , trauger , j. , ballester , g.  e. , burrows , c.  j. , casertano , s. , clarke , j.  t. , crisp , d. , griffiths , r.  e. , hester , j.  j. , hoessel , j. , krist , j. , matthews , l.  d. , scowen , p.  a. , stapelfeldt , k.  r. and westphal , j.  a. , 1996 apj , 466 , 732 g