it is now well - established that early - type galaxies possess substantial ism , with neutral hydrogen masses on the order of @xmath0 @xmath5  @xmath6 @xcite . accompanying this neutral gas
are emission - line regions extending out to several kiloparsecs , with line ratios characteristic of gas in a relatively low state of ionization .
they appear to be powered by the diffuse galactic background , rather than being a superposition of individual h ii regions around bright compact sources . in the past
, much effort has been devoted to understanding the sources of ionizing radiation powering these nebulae . in passively - evolving ( non - star - forming ) galaxies , it is the old stellar population itself which appears likely to provide the dominant contribution to the local extreme - uv background ( e.g * ? ? ?
* ) , at least outside of their innermost nuclei . there
, low - luminosity agn had until recently been thought to provide at least a significant contribution , at least in many cases ( e.g. * ? ? ? * ) . however
, increasing evidence suggests that the old stellar population may still provide a significant , or even dominant role in powering liner emission ( e.g. * ? ? ?
* ) . within the stellar population ,
post - asymptotic giant branch stars ( pagbs ) remain the favoured candidate for powering the observed emission lines .
in particular , pagb stars are expected to provide sufficient ionizing photons in order to account for the observed h@xmath7 flux @xcite , and the radial surface brightness profiles of the observed emission lines are consistent with ionization by a diffuse population , rather than a central source ( e.g. * ? ? ? * ; * ? ? ? * ) .    at the same time ,
similar arguments have proven capable of ruling out other possible ionizing sources .
recently , @xcite proposed that understanding the mechanism powering the diffuse line - emission in passively - evolving galaxies may shed light on the origin of type ia supernovae ( sne ia ) .
there exist two families of models for how these tremendous explosions come about : in the double - degenerate scenario @xcite , these events result from the merger of a binary pair of white dwarfs ( wds ) . in the standard model of the single - degenerate channel @xcite
, a carbon - oxygen wd accretes matter from a main sequence or red giant companion until triggering a thermonuclear explosion .
if the latter is correct , then accreting , nuclear - burning wds should be an extremely luminous component of any stellar population . with temperatures on the order of @xmath8
@xmath4 , they should provide the dominant ionizing background in relatively young ( @xmath3 4gyr ) passively - evolving stellar populations .
if this is the case , it should be readily apparent from spectroscopic observations of those ellipticals which host low - ionization emission - line regions @xcite .
this provides a unique opportunity to constrain plausible progenitor models for sne ia .
in particular , the observed line flux in recombination lines of ionized helium can provide an excellent test for the presence of relatively high temperature sources ( @xmath9 2  6 @xmath10 @xmath2 ) , and their total ionizing luminosity .
this covers much of the temperature range we expect from accreting , nuclear - burning white dwarfs , such as supersoft sources ( ssss ) , or in particular , the so - called ` ultra - soft sources ' ( those accreting white dwarfs with greatly inflated photospheres , see * ? ? ?
* ) . however , we expect higher temperatures ( up to @xmath0 @xmath4 ) for the most massive ( @xmath11 1.0 @xmath12 ) nuclear - burning white dwarfs . in the past ,
any possible population of such objects has been severely constrained based on the observed paucity of soft x - ray emission in early type galaxies , and the lack of individual soft x - ray sources ( e.g. * ? ? ?
* ; * ? ? ?
however , independent diagnostic tools remain extremely important .    here
we extend the work of @xcite to include a number of emission - line diagnostic tests for the presence of a population of very hot sources ( @xmath13 @xmath11 @xmath14 ) , as well as those too cool to strongly ionize he ii ( @xmath13 @xmath9 @xmath15 ) . in @xmath162 , we review our model assumptions for the ism in early - type galaxies , as well as the possible sources of ionizing radiation in passively - evolving stellar populations . in @xmath163 , we discuss those lines which are most sensitive to ionization by high temperature sources , in particular those of [ o i ] , [ n i ] , and [ c ii ] .
together with observations of recombination lines of he ii , this extends the range over which emission lines can meaningfully constrain the hardness of the ionizing background in passively - evolving galaxies to cover all temperatures expected from steadily nuclear - burning white dwarfs , including those with inflated photospheres ( @xmath17 @xmath18 @xmath13 @xmath18 @xmath4 ) . in @xmath164 , we outline the even more pronounced effect of any putative sd progenitor population ( with @xmath19 @xmath11 @xmath2 ) in young ( t @xmath3 1gyr ) , post - starburst populations ( such as e+a galaxies ) . in this case
the ionizing luminosity provided by sd progenitors would exceed that available from any other stellar source by up to two orders of magnitude , indicating the specific luminosity in any of the classical strong optical lines ( h@xmath20 , [ o ii ] 3727 , etc . )
should be similarly enhanced .
we follow this in @xmath16 5 with estimates of the observational prospects for producing limits using a number of instruments and surveys , as well as how the presence of a significant population of accreting wds may compromise measures of the sfr in galaxies .
finally , we discuss implications for the evolution of binary stellar populations and the progenitors of type ia supernovae .
we model the emission - line regions in early - type galaxies using the 1-d photoionization code mappings iii ( see e.g. * ? ? ?
our procedure follows closely that of @xcite ; here for convenience we summarize some of the essential ingredients .
unless otherwise stated , we assume a solar @xcite metallicity for the gas phase , and 2.5 solar for the stellar population . the former is consistent with the mean oxygen abundance measured in the warm ism of nearby retired galaxies @xcite . for the ionizing continuum from the ` normally - evolving ' ( single star ) stellar population , we use the population synthesis calculations of @xcite .
the most common morphology of the cold ism in such galaxies is a disky distribution extending from the nucleus out to several tens of kpc @xcite , with emission - line nebulae typically ionized by the diffuse background of the stellar population rather than nearby individual sources .
therefore , we assume plane parallel geometry throughout . for our standard case ,
we assume a constant hydrogen number density of 100@xmath21 , consistent with the observed [ s ii ] 6717/6731  ratio observed in stacked passively - evolving galaxies from the sdss @xcite , although this may only be an approximate upper bound @xcite . in principle , the majority of the optical emission lines considered in this work are relatively insensitive to variations in the density .    in our photoionization calculations ,
we assume that the modeled nebulae are ionization - bounded ( see * ? ? ? * and references therein ) , and terminate our calculations when the ionized hydrogen fraction falls below 1% .
the ionization parameter ( u = @xmath22 ) is similarly well constrained by the observed ratio of [ o iii ] 5007 / h@xmath23 , generally lying in the range -4 @xmath3 log(u ) @xmath3 -3.5 ( e.g. * ? ? ?
* ; * ? ? ?
* ) .    where appropriate
, we normalize the line flux per unit stellar mass of the host population , and assume a covering fraction of unity .
this can of course be scaled as needed , since the line luminosity is directly proportional to the stellar ionizing luminosity absorbed by the ism in ellipticals . in practice , it is often more convenient to provide line predictions normalized to the nearest recombination line of hydrogen , in order to minimize the importance of reddening in comparison with observed values .
this has the added benefit of removing any dependence on the covering fraction of the gas , therefore we also provide predictions for each optical line luminosity relative to the nearest hydrogen recombination line .    ]      in early - type galaxies , pagbs remain favoured to be the predominant ionizing source ( at least within the stellar population ) .
empirically , this has been strongly supported by energetic considerations ( e.g. * ? ? ?
* ) , as well as more detailed emission - line studies @xcite .
the spectral emission from these sources has been studied in detail by @xcite , whose models are included in the spectral synthesis calculations of @xcite used here . despite mounting evidence ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , some doubt remains about the longevity of post - agb stars , and therefore their expected numbers and contribution to the ionizing background @xcite .
therefore , where appropriate , we will consider in our calculations the ionizing background with and without the contribution normally expected from pagbs .
non - stellar sources , such as low - luminosity agn , likely also play a role in many galaxies @xcite , at least in their central regions .
fast shocks have been thought to play a role in some cases , as well as ionization by the hot , x - ray emitting diffuse ism @xcite .
however , such contaminating cases can be excluded with the careful use of line diagnostics ( e.g * ? ? ?
* ) , allowing one to confine their attention to so - called retired galaxies where the old stellar population is the principle source of ionizing emission .
other stellar sources , such as low - mass x - ray binaries , have been suggested , but do not appear to be significant in this regard @xcite .
accreting , nuclear - burning wds have , however , been found to likely provide a significant contribution , and would be expected to dominate the ionizing contribution for relatively young , passively - evolving stellar populations , at least if they exist in sufficient numbers to account for the observed sn ia rate @xcite .
even if they are not the primary evolutionary channel for sne ia , accreting wds may still account for a substantial contribution to the ionizing background ( see chen et al . in prep ) .
the number of nuclear - burning wds expected in any stellar population , assuming the sd channel accounts for all sne ia , can be estimated from the sn ia rate , the mass accretion rate , and the characteristic mass accreted per sn ia ( as in * ? ? ?
* ; * ? ? ?
* ) :    @xmath24    where @xmath25 is sn ia rate of the galaxy , and
@xmath26 is the total mass accreted prior to explosion ( @xmath27  @xmath28 ) .
this depends sensitively on our assumptions regarding the typical accretion rate , which in turn depends on as - yet poorly understood details of binary evolution .
fortunately , we can remove any dependence on the accretion rate if we consider instead the total bolometric luminosity of all accreting , nuclear - burning wds , @xmath29 .
the luminosity of any individual source is simply @xmath30 , with @xmath31erg / s the energy release from nuclear - burning of hydrogen , and @xmath32 the hydrogen fraction of solar metallicity gas .
therefore , the total luminosity of any putative sd progenitor population can easily be estimated as :    @xmath33    independent of the assumed accretion rate @xcite . for the snia rate we will use the delay time distribution derived for passively evolving galaxies by @xcite
: @xmath34snia / century/@xmath35 ( assuming a chabrier imf and z = 0.05 ) .    the spectral shape of the ionizing background in ellipticals
can then be found from the sum of two terms :    @xmath36    .predicted specific line luminosities ( in units of erg / s/@xmath12 ) for those lines most strongly effected by the introduction of a high temperature source population .
we assume ionization by a 3 gyr - old ` normal ' stellar population only ( column 3 ) , and with the addition of a @xmath13 = @xmath4 sd sn ia progenitor population ( column 4 ) . in both cases
, we take log(u ) = -3.5 .
[ line_list ] [ cols="^,^,^,^,^ " , ]     [ fits ]    in figure [ abs_lums ] , we plot the total luminosity in ly@xmath20 , h@xmath20 , [ n ii ] 6583 , [ o iii ] 5007  and the [ o ii ] 3727  doublet , as a function of source temperature for ionization by sd progenitors alone .
note that the predicted ly@xmath20 line luminosity should be only be considered an indicative estimate , as discussed in @xmath16 3 .
we take a 500 myr - old ssp as our example of a post - starburst stellar population and assume our standard characteristic accreted mass @xmath37 . from fig.5
one can see that the presence of a high - temperature sd progenitor population leads to an enhancement for all lines by well over an order of magnitude ( except for the [ o iii ] 5007  line for high source temperatures ) .
the [ n ii ] 6583   line luminosity exceeds that expected in the case of photoionization by the `` normally - evolving '' stellar population ( @xmath9 3.6 @xmath10 @xmath38 @xmath39 for @xmath40 = @xmath41 ) by a factor of @xmath42 across the considered temperature range .
the most luminous of the considered lines is the [ o ii ] 3727  doublet , which lies in the near - uv on the boundary of the wavelength range typically considered the `` optical '' .
however , a number of ground - based surveys have already been carried out specifically investigating this line ( typically as a proxy for the sfr , see @xmath16 [ prospects ] ) .    the total luminosity in any line is determined by the incident ionizing flux , which in the sd scenario is proportional to the accreted ( and nuclear - processed ) mass ( recall eq .
[ eq : ltot ] ) .
therefore any measurement ( or upper limit ) on their luminosities associated with a given stellar population will constrain the total ionizing luminosity of any putative sd progenitors .
this is particularly useful in the lower range of temperatures predicted for nuclear - burning wds ( @xmath43 ) , where both the he ii diagnostic and tracers of the str " omgren boundary fail . the significant enhancement in the line luminosity expected given photoionization by sd progenitors should allow the total mass processed by nuclear - burning accreting wds to be constrained to @xmath3 0.01
0.001 @xmath44/sn ia ( for example , using the [ o ii ] 3727  line , see fig . [ oii_ev_lums ] ) , should line fluxes be found to be consistent with ionization by pagbs alone .
6300  line and the [ o ii ] 3727  doublet as a function of stellar population age .
we assume @xmath45 and an ionization parameter of log(u ) = -3.5 , with line fluxes given ionization by the `` normally - evolving '' stellar population only ( red , solid ) , and with the inclusion of our example sd progenitor population ( blue , dashed ) . ]
a number of optical spectroscopic surveys of extended low - ionization emission line regions in early - type galaxies have already been conducted ( e.g. * ? ? ? * ; * ? ? ?
* ) , or are presently ongoing @xcite
. therefore , a large sample of such galaxies already exists which can be compared with our predictions . in estimating the feasibility of constraining the sd channel using optical forbidden lines ,
the most useful quantity is not the absolute luminosity per se , but rather the ew , as it characterizes the observability of the line . in fig .
[ ew_hight ] , we show for a 3 gyr - old ssp the ew of [ n i ] 5200 , [ o i ] 6300 , and [ c ii ] 1335 ) , taking the continuum emission as predicted from @xcite . for significantly older stellar populations ( @xmath46 5gyr )
our diagnostics become signigicantly less discriminating , at least for the presence of sd sn ia progenitors , given their much lower expected luminosity at late times ( recall figs .
[ 0.1__1gyr ] , [ oii_ev_lums ] ) . with the ew of the line and its broadening @xmath47
, we can estimate a minimum s / n of the spectrum needed to detect any line at a desired amplitude - to - noise ratio a / n using the formula given by @xcite :    @xmath48{2\pi}\sigma}{\rm{ew}}\label{s_n}\ ] ]    where we will take their characteristic broadening and minimum a / n ( @xmath47 = 120km / s , a / n = 4 ) .
for sd source temperatures above @xmath49 , we predict for the [ o i ] line an ew of @xmath0 1  for the realistic case @xmath40 = 0.5 ( see fig . [ ew_hight ] ) .
therefore , in order to barely detect [ o i ] 6300  emission , we would require a spectrum with s / n @xmath9 25 .
such a s / n is only just out of reach of typical sdss spectra ( s / n @xmath0 10  15 ) , and a very robust detection should esily be possible with stacking .
note that the [ o i ] line is also frequently detected in individual galaxies in the sample of @xcite .
given the similar continuum level at 6563  and the much greater intrinsic strength of the h@xmath20 line ( see fig .
[ abs_lums ] ) , the latter should similarly be quite easily detectable .    as discussed in @xmath16 3 , normalizing our predictions for the flux in any line to that of the nearest hydrogen recombination line has the advantage of removing any dependence on the covering fraction or total mass of the stellar population . in fig . [ oi_dm ]
, we plot the predicted [ o i ] 6300 / h@xmath20 ratio as a function of the total mass accreted by sd progenitors at an effective temperature of @xmath4 .
note that , whereas we would expect a linear relation for the total line flux as a function of @xmath26 , this is not the case for diagnostic line ratios .
this is primarily due to the contribution of the ionizing luminosity from sd progenitors to the ionization of hydrogen .    as an illustrative example
, we can compare our models with the observed [ o i ] 6300 / h@xmath20 ratio in the nucleus of the post - starburst galaxy ngc 3489 @xcite . as a conservative estimate ( see e.g. * ? ? ?
* ) , we assume a 2 gyr - old ssp , with the available photoionizing continuum originating from both pagbs and sd progenitors . we see that , for a wide range of gas phase metallicities ( 0.5@xmath50 @xmath18 @xmath51 @xmath18 1.5@xmath50 ) , the [ o i ] 6300 / h@xmath20 diagnostic can in practice distinguish ( at @xmath9 1@xmath47 confidence ) between ionization by pagbs alone and the inclusion of a sd progenitor population which accretes more than @xmath90.04@xmath12 per sn ia at @xmath19 = @xmath4 .
[ ni_dm ]     which may be accreted by the typical wd in the sd scenario at a given photospheric temperature .
upper limits are obtained assuming an upper limit of 10% on the he ii 4686 / h@xmath7 ratio ( red line ) , and for the measured [ o i ] 6300 / h@xmath20 ratio in the nucleus of ngc 3489 @xcite ( black , solid ) .
for the latter , we also show 1@xmath47 upper limit ( black , dashed line ) .
note that any upper limit on the accreted mass is equivalent to an upper limit on the total luminosity of the population of high temperature sources ( right y - axis ) .
we assume 2 gyr - old stellar population , log(u ) = -3.5 , and solar metallicity for the gas phase . ]
a similar case can be made for the [ n i ] 5200 / h@xmath7 line diagnostic ( see fig .
[ ni_dm ] ) . though the [ n i ] 5200  line strength is always weaker than [ o i ] 6300  ( and , unfortunately , not detected in the sample of * ? ? ?
* ) , it is still observed with confidence in a number of liers ( e.g. * ? ? ?
[ ew_hight ] , we see that for sd source temperatures @xmath52 , the [ n i ] 5200  ew is roughly @xmath0 0.5  ( for @xmath53 and log(u)=-3.5 ) . using eq .
[ s_n ] above , we find that this line could then be detected in a 3 gyr - old stellar population in spectra with s / n @xmath9 40 .
this is well within what is typically achieved in individual sauron spectra ; with stacking , very robust detections ( or upper limits ) should be feasible .
stacking of sdss spectra should allow for a confident detection with only @xmath0 20 galaxies . from fig [ ni_dm ]
, we see that similar to the [ o i ] line , a detection or upper limit on the [ n i ] 5200 / h@xmath7 ratio of @xmath9 10% would constrain the average mass proccessed by nuclear - burning wds in the sd channel to a few @xmath54 per supernova .    turning again to the [ o i ] 6300 / h@xmath20 ratio in ngc 3489
, we can extend the exercise in fig .
[ ni_dm ] to investigate its dependence on wd temperature . in fig .
[ dm_teff ] , we plot the maximum mass which can be accreted by nuclear - burning wds in the sd scenario for effective temperatures @xmath8 @xmath3 @xmath13 @xmath3 @xmath4 , should our models , we include the expected contribution from pagbs to the ionizing background . ]
remain consistent with the observed ratio [ o i ] 6300 / h@xmath20 .
note that the upper limit for any temperature in fig , [ dm_teff ] excludes any matter accreted at any other temperature .
for example , accretion of @xmath55 of matter at @xmath56 k will produce an [ o i ] line luminosity equal to the observed value , therefore no more material at any other temperature could additionally be accreted without conflicting with observations . in practice ,
any realistic sd scenario would imply a population of wds in varying accretion states , with a range in temperatures . therefore any detailed sd scenario may be constrained considering the integrated luminosity and the spectrum of the ionizing radiation produced by the population of sd progenitors .    for ngc 3489 ,
given the observed value of the [ o i ] 6300 / h@xmath20 ratio , no more than @xmath9 @xmath57 per supernova can be accreted and nuclear - processed by wds with photospheric temperatures
@xmath13 @xmath11 @xmath58 . taking into account the [ o i ] line flux error
, this translates to an upper limit of @xmath59 at @xmath60 .
note that the flux errors in @xcite are not dominated by systematics , and in principle dedicated observations should be able to reduce the uncertainty considerably .
if we compare this figure with the sensitivity of the he ii diagnostic in @xcite ( see their fig . 8 , or the red curve in our fig . [ dm_teff ] ) , we see that the two complement each other remarkably well ; upper limits on the he ii line luminosity ( here we assume he ii 4686 / h@xmath7 @xmath61 10% ) reach their peak effectiveness at @xmath9 2.5@xmath62 , just as the [ o i ] 6300 / h@xmath20 diagnostic becomes ineffective .
an important point to consider is that , while we typically quote limits here in terms of the average mass accreted by nuclear - burning wds ( @xmath26 ) , in principle it is rather the total luminosity of any high - temperature sources which we constrain ( see right y - axis in fig . [ dm_teff ] ) .
even if sd progenitors are not primarily responsible for the sn ia rate , accreting wds may still provide a significant contribution to the ionizing background in passively - evolving galaxies .
therefore , their consideration would be of vital importance in any exercise in understanding the role of pagbs in low - ionization emission line regions , or their relative importance with respect to llagn ( as in , e.g. * ? ? ?
therefore , careful modelling of the evolution of wd binaries and any phases of stable mass transfer are essential ( see chen et al , in prep ) .
5007 / h@xmath7 and [ o i ] 6300 / h@xmath20 for varying ionization parameter , assuming ionization by a 3 gyr - old ssp alone ( red ) , and with the include of a sd progenitor population which accretes @xmath26 = @xmath63 ( green ) and @xmath26 = @xmath64 ( blue ) with an effective temperature @xmath13 = @xmath4 .
the ionization parameter increases along these lines from log(u ) @xmath9 -4.5 to -3 , in the direction indicated by the black arrow . shown in black
are a subset of the results from @xcite . ]     3869/[o ii ] 3727and the ew of [ o ii ] 3727 / h for varying ionization parameter , assuming ionization by a 3 gyr - old ssp alone ( black ) , and with the inclusion of a sd progenitor population which accretes varying @xmath26 with photospheric temperature @xmath13 = @xmath14 and log(u ) = -3.5 ( red , solid ) and @xmath13 = @xmath4 with log(u ) = -3.5 ( blue , solid ) and log(u ) = -3.0 ( blue , dashed ) . for the tracks which include sd progenitors , @xmath26 increases from 0 to 0.3@xmath12 from left to right in the figure . ]
finally , we note that in the presence of a hot sd progenitor population , the higher values of [ o i ] 6300 / h@xmath20 and [
n i ] 5200 / h@xmath7 would place most passively - evolving galaxies in the region of the bpt diagram associated strictly with shock ionization ( -0.5@xmath65([o iii ] 5007 / h@xmath7 ) @xmath66 , and @xmath67([o i ] 6300 / h@xmath20 ) @xmath46 -0.5 ) .
however , this has not been observed ( e.g. * ? ? ?
* ; * ? ? ?
* ) . in fig .
[ oiii_oi ] , we plot the predicted values of the [ o iii ] 5007 / h@xmath7 and [ o i ] 6300 / h@xmath20 ratios with varying ionization parameter assuming ionization by a 3gyr - old ssp alone , and with the inclusion of our standard sd progenitor population .
we also include a subsample of the results from @xcite ; namely , those galaxies with stellar population ages @xmath61 4 gyr ( from * ? ? ?
* ) , with ew(h@xmath20 ) < 1.5  @xcite .
this restricts us to relatively young early - type galaxies for which the `` normally - evolving '' stellar population may plausibly power the observed nebular emission . the observed sample in fig .
[ oiii_oi ] appears to be largely reproduceable with ionization by the ssp alone , although any interpretation is complicated by the considerable uncertainty in the majority of the emission line measurements .    for a number of reasons , emission - line diagnostics in
the near - uv provides an attractive alternative .
as discussed in the previous section , the [ o ii ] 3727  doublet provides a much higher luminosity than h@xmath20 , in a waveband where the continuum is instrisically lower .
the very nearby [ ne iii ] 3869  line is also available as a gauge of the ionization parameter ( through the ratio , e.g. * ? ? ? * ) .
this allows one to construct a bpt - style diagnostic diagram using a remarkably narrow wavelength range ( as shown in fig .
[ neiii_oii ] ) , using the ew of the [ o ii ] doublet to measure the specific ionizing luminosity of the observed stellar population , and [ ne iii ] 3869/[o ii ] 3727  to fix the ionization parameter . in fig , [ neiii_oii ] , we see that the ssp - only and sd progenitor ionizing populations occupy clearly distict regions at t = 3gyr . for a young , post - starburst galaxy ( t @xmath9 500myr ) ,
the ew of [ o ii ] 3727  would be enhanced by a further order of magnitude for the same ionization parameter ( see previous section ) .
the nearby 5@xmath683 transition of he ii ( 3203 ) 3 times weaker than the 4 @xmath683 transtion at 4686 , the continuum is predicted to be @xmath04 times weaker at @xmath93200 .
therefore the observability and utility of this transition is comparable ( see eq .
[ s_n ] ) ] . would also allow one to constrain the hardness of the ionizing source population .
together , this presents an excellent opportunity for near - uv spectroscopy to constrain the progenitors of sne ia .
the prospect of calorimetry of the ionizing background in post - starburst galaxies offers a highly promising avenue for constraining the sd channel at the earliest delay - times . in fig .
[ ew_stronglines ] , we plot the ew of the strongest optical lines for a 500 myr - old starburst . as discussed in @xmath16 4 , the most intuitive ( and least model - dependent ) line for this purpose may be the h@xmath20 recombination line . in fig .
[ ew_stronglines ] , we see that in our standard sd case , we predict an h@xmath20 ew ranging from @xmath9 3.5  at @xmath69 to @xmath9 12  at @xmath70 , assuming @xmath40 = 0.5 .
this may be compared with an h@xmath20 ew of @xmath9 0.2  expected from the ssp alone , giving a factor of @xmath71 enhancement .
thus , if observations of post - starburst galaxies gave h@xmath20 line intensities consistent with ssp - only case , the contribution of the sd channel could be strongly constrained to a few per cent level , or , equivalently , to @xmath72 per supernova .
the [ n ii ] 6583  line can provide a similar constraint ( see fig . [ ew_stronglines ] ) .
however , the nitrogen abundance in the ism of external galaxies is known to be rather uncertain @xcite .
therefore , there is little benefit in focussing on this line for calorimetry , unless it is unresolved from the h@xmath20 line .
the [ o iii ] 5007  specific line luminosity drops off dramatically at the highest sd temperatures ; this is because with harder spectra , oxygen is increasingly ionized to @xmath73 near the illuminated face of the gas , reducing the emitting volume for [ o iii ] . at the same time
, high - energy photons penetrate deeper into the nebula , preferentially ionizing neutral oxygen .
note , however , that this itself may pose a problem for any high - temperature sd channel , as the ism in post - starburst ellipticals are observed to have relatively high values of the [ o iii ] 5007 / h@xmath7 ratio @xcite .
the [ o ii ] 3727  cooling line , originating in the more low - ionized regions of nebulae ( the interior , in the geometry discussed here ) , is expected to have an ew @xmath11 15  ( for @xmath2 @xmath3 @xmath13 @xmath3 @xmath4 ) .
this is a @xmath74-fold increase over that expected from the ssp alone ( @xmath75 , detectable with a minimum s / n @xmath9 40 ) .
the [ o ii ] 3727  line ( as well as h@xmath20 ) has been invoked as a proxy for the star - formation rate ( more on this in the following section ) , and already the volumetric [ o ii ] 3727  luminosity function has been mapped at least to z @xmath18 5.4 ( see * ? ? ? * and references therein ) .
individual nearby galaxies which have undergone recent , possibly quenched starbursts would make ideal targets for follow - up studies . as an example , ngc 3156 is a nearby ( @xmath9 22mpc ) early - type galaxy which the sauron survey found to host a post - starburst subpopulation ( with an estimated mass of @xmath9 @xmath6 ) .
the proposed cubes nuv spectrograph would be able to detect the [ o ii ] 3727  line ( given our standard sd case from above , and assuming @xmath40 = 0.5 ) with a s / n @xmath9 10 @xmath76 in just 1800 seconds .    although its specific luminosity is greatly enhanced with the presence of a high temperature source population ( see fig .
[ ew_hight ] ) , the [ c ii ] 1335  line remains intrinsically weak relative to the strong uv recombination lines ( e.g. ly@xmath20 , he ii 1640 ) .
the difficulty in detecting this line can best be illustrated with an example . at a distance of @xmath0 21 mpc ,
ngc 3607 is a nearby , relatively young ( mean stellar age @xmath0 3gyr , e.g. * ? ? ?
* ) elliptical galaxy which hosts an extended low - ionization emission - line region . assuming our standard sd population from above and a covering fraction of @xmath77
, we can expect a total integrated [ c ii ] 1335  line flux of @xmath0 4.5 @xmath10 @xmath78 erg / s/@xmath79 , using the k - band luminosity to estimate the mass as in @xcite .
the galex fuv spectrograph ( now no longer operational ) would have required a prohibitive @xmath0 51 ksec in order to reach s / n @xmath0 5 , with 1335  lying outside the optimal response of the instrument . however
, this remains a promising target for future uv missions , such as the planned wso . until such time , the aforementioned he ii 1640  line remains an extremely promising target for observations ( see * ? ? ?
* ) , as does the ly@xmath20 transition of hydrogen at 1216 , the latter being easily detected in any emission - line elliptical galaxy . in the limit where sd progenitors provide the only contribution to the ionizing continuum
, the total luminosity expected in ly@xmath20 exceeds that of the he ii 1640  line by 1  2 orders of magnitude for sd source temperatures @xmath13
@xmath9 @xmath2  @xmath4 .      in young stellar populations ,
the presence of strong emission lines indicates ongoing star formation ( in particular , for low values of n ii 6583 / h@xmath20 , e.g. * ? ? ?
* ) , where interstellar matter is ionized by highly luminous o stars . in this case
, the total luminosity in any given line may be used as a tracer of the total ionizing photon flux , and therefore a measure of the number of o stars , indirectly providing the sfr ( for an assumed imf ) . in practice ,
the h@xmath20 and [ o ii ] 3727  fluxes are the most practical in the optical , with the standard calibrations being @xcite :    @xmath80\ ] ]    @xmath81}[erg / s]\ ] ]    however , this does not account for the possibility of an additional ionizing source within the stellar population , such as possible sd progenitors of sne ia .
this would produce an offset in the calibration of any such sfr indicator , which would be dependent on age and mass of the stellar population . for example
, for @xmath8 @xmath18 @xmath13 @xmath18 @xmath4 sd progenitor population , a 3 gyr - old @xmath82 ssp would still appear to have an ongoing star - formation rate of @xmath83 @xmath10 4@xmath12/yr from the [ o ii ] 3727  indicator .
similarly , for a @xmath13 = @xmath2 sd progenitor population , a 3 gyr - old @xmath82 ssp would still appear to have an ongoing star - formation rate of @xmath84 @xmath10 1@xmath12/yr from the h@xmath20 indicator .
because of their insensitivity to interstellar reddening , ir lines and the ir continuum itself are extremely useful sfr indicators . as an important coolant in the neutral ism , the [ c ii ]
157@xmath85 m line has been proposed as a useful sfr indicator , with the present best calibration given by @xcite    @xmath86}[erg / s])^{0.983}}{1.028 \times 10^{40}}\ ] ]    yet there remains considerable scatter in the observed relation .
if there is indeed a significant contribution from populations of accreting , nuclear - burning wds , this should also provide an age - dependent offset .
for example , for our standard case sd scenario , a 3 gyr - old @xmath82 ssp would still appear to have an ongoing star - formation rate of @xmath84 @xmath10 1@xmath12/yr .
the [ o i ] 63@xmath85 m and 146@xmath85 m ir lines would provide only a small contribution to any sfr - indicator based on the continuum flux .
however , as individual lines , both may in principle be a useful measure of the population of accreting wds in galaxies , in particular the 63@xmath85 m line , which is @xmath0 5 times stronger than the [ c ii ] 157@xmath85 m line for sd progenitor effective temperatures of @xmath8 @xmath18 @xmath13 @xmath18 @xmath4 . with the end of the herschel space observatory s operational lifespan , there is at present no instrument which could provide new observations in this waveband , however ample archival data is available .
the field - imaging far - infrared line spectrometer ( fifi - ls ) , will also soon be available , observing in the 42 - 210 micron band .    as a final word of caution , we recall that our results for the ir and fir forbidden lines discussed above are sensitive to the assumed density of the gas , which is observed to vary radially ( e.g. * ? ? ? * ) .
therefore , any quantitative assessment of the contribution of a possible sd channel to the emission of these lines should consider carefully the environment under study .
furthermore , these lines do not originate solely from photoionized regions , but arise also in the neutral medium @xcite .
so far , we have considered only single - age stellar populations , and assumed the same metallicity throughout .
although we have previously demonstrated that varying the metallicity of the stellar population will not significantly alter the ionizing continuum expected from sd progenitors @xcite , we have not tested the effect of varying star - formation histories . in a real galaxy , a much more complex superposition of populations of differing ages
should be expected .
this is particularly concerning , should the presence of a significantly older subpopulation be able to strongly dilute the observed ew in a galaxy with an otherwise young mean stellar age .
this could in turn significantly reduce the detectability of our diagnostic lines .    in order to evaluate this possibility , we perform the following experiment . we mix spectra of varying fractions of 1 gyr - old , z = 2.5@xmath87 and 10 gyr - old , z = @xmath87 stellar populations from @xcite . using the method of @xcite based on the lick indices @xcite
, we find the luminosity - weighted age depending on the relative fractions of young and old populations ( johansson , private communication , see also johansson et al in prep ) . at the same time , we compute the expected flux in any of emission lines discussed above for the composite stellar population , using the combined spectra in our photoionization calculations . here
we assume the same gas density @xmath88 with a gas phase metallicity of z = @xmath50 , as above , and with ionization parameter log(u ) = -3.5 .
6300  emission line as a function of time for single - age ssps ( blue ) , and for composite populations ( 1gyr and 10gyr , in red ) as a function of their  measured  age ( johansson et al , in prep ) . ]
the result of this experiment is shown in fig .
[ ew_mixed_pop ] , where we plot the ew of the [ o i ] 6300  emission line for such composite stellar populations , shown as a function of the derived `` mean age . '' for comparison , we also plot the [ o i ] ew for a single - aged ssp with z = 0.05 .
note that the two tracks do not agree at 10 gyrs because of the different assumptions for the metallicity of the stellar population .
as can be clearly seen , the ew is not significantly affected in the case of composite stellar populations .
we have demonstrated how observations of forbidden lines of the most abundant metals , as well as recombination lines of hydrogen and ionized helium , can be used to constrain any population of sources with effective temperatures @xmath13 @xmath11 @xmath2 in early - type galaxies .
these sources may not be directly detectable due to local interstellar absorption , as well as along the line of sight in the milky way , but may be revealed from the dramatic effect on the emission line fluxes predicted for their host galaxies .    a confident detection , or upper limits on
, any recombination line of he ii in the extended , low - ionization emission line regions of passively - evolving galaxies can constrain the total luminosity of sd progenitors with effective temperatures in the range @xmath89 @xmath3 @xmath13 @xmath3 @xmath90 . while this comfortably covers the range expected for symbiotics , many supersoft sourcess , and the so - called ultrasoft sources , other diagnostics
must be invoked in order to constrain all possible sd progenitors .
we have demonstrated that , in relatively young early - type galaxies which host significant neutral and weakly ionized ism , forbidden lines of [ o i ] , [ n i ] , and [ cii ] can be enhanced by up to @xmath91 times in the presence of a significant ( @xmath92(sd ) @xmath9 @xmath92(total ) , very hot ( @xmath13 @xmath11 @xmath14 ) population of accreting , nuclear - burning wds consistent with the sd progenitor scenario . in particular , past and presently ongoing optical spectroscopic surveys
should already be able to either detect evidence of the sd channel or constrain its contribution to the sn ia rate down to a few per cent of the observed total for inferred delay - times 1 gyr @xmath3 t @xmath3 4 gyr .
future observations with the planned wso ( the world space observatory * ? ? ?
* ) should also be capable of making similar progress with observations of the [ c ii ] 1335  ( and he ii 1640 ) lines .
the possible contribution of the sd channel at the earliest delay - times ( @xmath93 gyr ) can also be remarkably well - constrained with measurements of the specific ( per unit stellar mass ) luminosities of the `` classical '' strong emission lines , such as h@xmath94 , [ o ii ] 3727  , [ n ii ] 6583 .
for a 500 myr - old stellar population , the h@xmath20 recombination line luminosity is predicted to be enhanced over the ssp - only case by up to @xmath0 two orders of magnitude for sd progenitor temperatures t @xmath3 2@xmath62 .
a similar enhancement is expected for the [ o ii ] 3727  line luminosity for sd temperatures spanning the @xmath95 k range , suggesting calorimetry of post - starburst galaxies may be able to constrain the sd contribution to less than @xmath0 1% .
alternatively , observations of emission - line nebulae consistent with ionization by the old ssp alone may be interpreted as a constraint on the total mass accreted in the steady nuclear - burning regime . assuming a hypothetical sd scenrio accounted for all sne ia
, this would restrict the mass accreted at the @xmath95 k range of photospheric temperatures to @xmath3 @xmath96/yr .    finally ,
if the sd scenario did indeed account for a significant fraction of the sn ia rate , than one might expect the resulting enhancement in the luminosities of emission - lines h@xmath20 , [ o ii ] 3727 , and [ c ii ] 157@xmath85 m to introduce an age - dependent offset in their calibration as indicators of the sfr in galaxies . at the same time
, one would expect nuclear - burning wds to provide a vital ingredient in calculating the ionization budget in low - ionization emission line regions , necessitating a careful understanding of the evolution of accreting wd populations .
we thank the referee and jonas johansson for their helpful comments and discussion , and brent groves for making his current version of mappings iii available at http://www.mpia-hd.mpg.de/@xmath0brent/mapiii.html .


networks of semiflexible filaments and fibers are common in biological systems , where dynamic structures of tunable strength , elasticity , and adjustable response are required @xcite . from the microscopic scale of cytoskeleton to the macroscopic scale of fibrous tissue ,
such networks utilize stiff or semiflexible filaments , allowing for rich mechanical response @xcite , dynamic remodelling @xcite , and controlled structural failure @xcite , while remaining an open structure allowing easy through access for solvents and solutes @xcite .
inspired by these observations , filament networks are produced in industry , and now stand as a promising area in manufacturing functional materials  @xcite .
the physics of a single semiflexible chain or filament is well understood with the aid of worm - like chain model , in various implementations and approximations  @xcite .
the full theory is capable of accurately describing the force - extension relationship , as well as the magnitude of transverse fluctuations , for a range of filament stiffness spanning from very high ( a rigid elastic rod ) to very flexible ( a classical polymer coil ) and a range of end - to - end extensions from zero up to the full stretch where the entropic force has a characteristic divergence .
when such filaments form a macroscopic crosslinked network , the main scaling features of its nonlinear mechanical response are still determined by the single filament properties , as reviewed in  @xcite .
there have been many key contributions to the elastic theory of stiff and semiflexible networks , including ones highlighting the issues of strain non - affinity  @xcite .    in order to obtain the constitutive relationship of a semiflexible network ,
a widely used approach follows the following procedure ( explained in detail in chapters 3,4 of ref .
@xcite ) : the shear stress on @xmath1 plane is calculated by summing the contributions of the tension along the @xmath2 axis of all the chains crossing the @xmath3 plane .
this approach has been very successful in deriving viscoelastic properties of polymer solutions and melts under shear flow  @xcite ; storm et al .
@xcite have applied the same stress - strain relation to a crosslinked network by assuming the affine elastic strain acting on each filament and the probability distribution of semiflexible filament length of wilhelm and frey  @xcite .
this model qualitatively describes how the network behaves when being deformed ; however , it can only be used numerically ( since no closed expressions are possible ) and it omits the pre - stress acting on the chains , as in this model the tension is assumed as zero when there is no deformation .
wilhelm and frey also produced a numerical simulation of filament network elasticity on their own  @xcite , using the mikado model of connectivity and athermal rigid rods as elastic elements .
although important issues of percolation rigidity threshold are exposed , this work can not be used to describe most experimentally relevant ( and most biological ) filaments .
more recently , palmer and boyce  @xcite proposed a closed analytical form of elastic free energy , with the corresponding constitutive relation , using the same force - extension filament relationship as storm et al .
they applied the ` 8-chain model ' originally introduced in the context of ordinary rubber elasticity  @xcite with an approximate expression for individual filament elasticity applied for each strand .
this might be the best attempt in formulating the nonlinear elasticity of semiflexible network to date . on the other hand , unterberger et.al .
@xcite developed a ` 1-chain ' model for a network , by assuming filaments have a homogenous orientational distribution in the equilibrium system , and the average ( imposed ) stretch of the network @xmath4 is a @xmath5-root average of deformations of individual filaments : @xmath6^{1/p } / a^{1/p}$ ] , where @xmath7 is the deformation of an individual filament , @xmath5 is the averaging parameter and @xmath8 the area of the unit sphere .
this procedure allows the local stretch to fluctuate around its average value , and may partly account for the nonaffinity .
this model could only be applied numerically to reproduce several key mechanical features ( shear and normal stress ) .
most good models discussed here make a successful fitting of shear stress data of different filament networks , which captures the generic stress - stiffening effect originating from the characteristic divergence of the force - extension curve of an individual filament .
usually , the authors employ the versions of celebrated marco - siggia interpolation formula  @xcite , which gives the correct response for filaments near full - extension ( but much less so for more flexible filaments ) .
in fact , the well - known mackintosh scaling of differential shear modulus with stress , @xmath9 , in the stress - stiffening regime is entirely based on the mentioned divergence of an individual filament near full extension  @xcite ( it holds even for an athermal network of undulating filaments , as long as the tensile force is proportional to the inverse - square of the compression , as is the case near full extension  @xcite ) . therefore , a mere agreement ( good fitting ) of shear stress - strain curves is not a sufficient test of different theories .
in particular , they must simultaneously descibe the effect of negative normal stress in the network of stiffer filaments , with a positive normal stress ( also known as the weissenberg effect ) for networks made of more flexible chains  @xcite .
it turns out that none of the mentioned theories , including the 8-chain model of palmer and boyce , are able to produce the correct normal stress , or address the network stability increasing with filament pre - tension@xcite . in this paper
we develop and discuss a continuum theory of semiflexible network ( in closed analytical form ) that addresses these issues , while retaining accurate fitting of a wide range of shear stress data .
before introducing free energy of a semiflexible network , we need to discuss the properties of a single semiflexible filament . a filament connecting two neighboring crosslinks in a network is sketched in fig .
[ fig_1](a ) , with coordinates @xmath10 , where @xmath11 is the arc - length coordinate along the chain .
different approaches to the chain ( in)extensibility have been tested over the years  @xcite , from the strict constraint to the requirement that the length of filament remains constant on average ( while small local fluctuations are allowed )  @xcite , to the models that explicitly include filament stretching  @xcite .
it turns out that in the regime of high extension , when there are no foldbacks ( hairpins ) on the chain , all length - constraining models give the same divergence of the force - extension , @xmath12 , with the end - to - end ratio @xmath13 , where @xmath14 is the end - to - end length of the filament , fig .
[ fig_1](a ) .
this was used by the famous interpolation formula of marko and siggia  @xcite .
more recently , a complete theory of filament entropic elasticity has been developed , which spans the full range of extensions and the full range of bending modulus  @xcite .    in the worm - like chain model  @xcite
, the bending energy can be expressed by @xmath15 , where @xmath16 is the bending rigidity .
the key physical quantity for describing the stiffness of a polymer chain is the persistence length @xmath17 .
when the contour length of the filament , @xmath18 , is comparable with its persistence length @xmath19 , the chain is regarded as  semiflexible " . combining the effects of enthalpy arising from bending and entropy of conformation fluctuations , the closed form of the single chain
free energy  @xcite can be expressed as a function of its end - to - end factor , @xmath13 : @xmath20 where @xmath21 is a dimensionless stiffness parameter reflecting the competition between bending and thermal energy ; in our notation @xmath22 . in ref.@xcite
one can find the comparison with several notable models of semiflexible filament ( marko - siggia and ha - thirumalai ) and where they deviate from the accurate expression ( [ achain ] ) .    , and end - to - end length @xmath14 being the mesh size .
( b ) the relationship @xcite between the tensile force , @xmath23 , and end - to - end factor , @xmath24 , plotted for semiflexible chains with different stiffness @xmath25 , and marking the position @xmath26 for a chain to stay at the force - free state . ]
as shown in fig  [ fig_1](b ) , if the value of @xmath25 is smaller than a critical value @xmath27 , the minimum of the free energy ( or the force - free natural length ) will be at @xmath28 ; such a chain can be regarded as flexible . when @xmath29 , one recovers the gaussian entropic spring form : @xmath30 .
the marko - siggia ( in fact , fixman - kovac ) limit commonly used for quick fitting of force - extension curves is reached for flexible chains at high extension @xmath31 : @xmath32 .
when @xmath33 , the chain can be regarded as a stiff rod with a natural length @xmath34 , and its energy is dominated by elastic bending : @xmath35 $ ] . under compression
@xmath36 such a filament undergoes euler buckling instability .
note that on forming a crosslinked network , @xmath37 does not need to be equal to @xmath26 for each strand in equilibrium , meaning there can be pre - tension in the network .
we now construct the continuum elastic free energy of a network of such filaments using the methodology that was successfully developed in rubber elasticity  @xcite .
if the sample shape is changed with a deformation tensor @xmath38 , then the corresponding cauchy - green tensor  @xcite is @xmath39 , with eigenvalues : @xmath40 , @xmath41 and @xmath42 .
the values @xmath43 , @xmath44 , @xmath45 can be interpreted as stretching ratios along the principal directions of deformation .
a class of simple , yet powerful theories of rubber elasticity is formulated in terms of invariants of the c - g tensor , @xmath46 , @xmath47 and @xmath48 .
famous examples of this class are neo - hookean ( gaussian ) , mooney - rivlin , and gent models ( see @xcite for review ) . note that polymer networks are frequently treated as incompressible , so @xmath49 .        for a polymer network , the chains are distributed and crosslinked randomly throughout the material , which is the main difficulty in obtaining an analytical expression of the total free energy of the network ( unless the individual chain is gaussian ) . to incorporate a more complicated chain , such as eq .
( [ achain ] ) , into a constitutive framework , it is necessary to have a model that relates the chain deformation to the applied affine strain .
this can be accomplished by representing an element of volume in the average network as a `` mesh cell '' , which is then symmetrically multiplied to fill the volume . here
we need to distinguish  affinity " and  non - affinity " in a filament network constructed by the unit - cell method .
a  mesh cell " is deformed differently from the material , as these cells are constructed in orientation aligned with the principal stretching direction of the material , which is referred to as  non - affinity " by palmer and boyce  @xcite .
all mesh cells deform in the same way , with repeated deformed structure .
however , the real non - affinity can arise from the different responses among different  mesh cells " , due to local force relaxation . in this work
we refer to the concepts of  affinity " or  non - affinity " in the local context of individual filament junctions and mesh cells , rather than the global one between the orientation of a mesh cell and the macroscopic deformation of the material . in assuming a symmetry of repeated mesh cells
, we automatically discard the effects of the local non - affine deformations  @xcite , i.e. take all mesh cells deforming uniformly , which we know must be the case at least for stiff filaments @xcite . since there is no quantitative way of assessing the degree of the error introduced by the affine approximation ( in our interpretation ) , we will have to look at the fits to the experimental data for validation . within a mesh cell , chains or crosslinks
have several possible arrangements , reflecting what one assumes about the topology of the network mesh , see fig .
[ models ] .
most acceptable structures include the homogeneous sphere ( hs ) in 1-chain network model , the primitive cubic ( pc ) in 3-chain model , tetrahedral ( th ) in 4-chain model , and body - centered cubic ( bcc ) in 8-chain model ( see @xcite and a review @xcite for detail ) .    in 1-chain model  @xcite , one end of a polymer chain
is fixed at the center of a sphere , while the other end is on the sphere surface at an arbitrary orientation @xmath50 , distributed isotropically .
when deformed by @xmath38 , with stretching ratios along principal directions @xmath43 , @xmath44 , and @xmath51 , the lengths of the three semi - axes will change from @xmath14 , to @xmath52 , @xmath53 and @xmath54 , respectively .
the elastic energy density of the network can be expressed as an orientational average of a deformed filament : @xmath55,\end{aligned}\ ] ] where @xmath56 is the density of crosslinked chains and @xmath57 .
if one takes the ` entropic spring ' limit of the gaussian chain , the average @xmath58 reduces to the classical neo - hookean rubber - elastic expression @xmath59 $ ] .
however , in the general case the free energy in 1-chain model can not be expressed in analytical form as a function of strain invariants , which renders it less convenient . in the following ,
we compare the 3- and 8-chain models , and decide on the 3-chain model preference , in particular due to the failure of 8-chain model in reproducing the normal stress .    in the 3-chain model , a primitive cubic
is constructed with lattice points representing the crosslinking sites , and the edges are aligned along the principle directions of deformation tensor @xmath38 .
three chains are linked with their end - to - end vectors along the edges and the equilibrium mesh size @xmath14 .
on deformation , the lengths of three perpendicular edges at one lattice point become @xmath60 , @xmath61 and @xmath62 , respectively
. then the free energy density of a semiflexible network can be expressed as @xmath63 equation  ( [ 3chain1 ] ) can be rearranged as a function of the strain invariants : @xmath64\end{aligned}\ ] ] where @xmath22 and @xmath13 as used in eq .
( [ achain ] ) .
the body - centered cubic cell of 8-chain model is constructed with eight filaments connected from the center point to all eight lattice points @xcite .
the edges of the cell are @xmath14 as we define the mesh size , while the chains are shorter by factor @xmath65 .
the important feature of the high - symmetry 8-chain model is that on deformation all chains change their distance by exactly the same amount , from @xmath66 to simply @xmath67 , since the mesh cell is aligned along the principal axes of deformation tensor @xmath38 .
the free energy density of the network in this case is given by the single - chain expression directly : @xmath68    when the elastic energy is expressed as a function of strain invariants , the stress tensor of an incompressible material ( with @xmath49 ) can be obtained as  @xcite : @xmath69 - p\delta_{ij } ,   \label{strain}\end{aligned}\ ] ] where @xmath70 is the cauchy strain , and @xmath71 the lagrangian multiplier for incompressibility , the value of which determined by the boundary conditions .
let us consider the simple shear deformation such that in cartesian coordinates a point @xmath72 in an original material will change to @xmath73 after being sheared ; @xmath74 is the shear strain .
the incompressibility is satisfied automatically , and the remaining strain invariants are : @xmath75 .
the shear stress in 3-chain model can be obtained from the general constitutive relation ( [ strain ] ) as : @xmath76_{\mathrm{3c } } = \frac{2}{3}n k_\mathrm{b } t\gamma x^2 \left [   \frac{\left(1-x^4\right)}{c \pi   \left[1-\left(2+\gamma^2\right ) x^2+x^4\right]^2}-    c \pi ^2\right],\end{aligned}\ ] ] while the corresponding expression in 8-chain model is : @xmath77_{\mathrm{8c}}= \frac{2}{3}nk_\mathrm{b}t\gamma x^2\left [   \frac{9}{c \pi   \left[3-\left(3+\gamma^2\right ) x^2\right]^2}-    c \pi ^2\right].\end{aligned}\ ] ] one must distinguish the nominal shear modulus @xmath78 , from the differential shear modulus @xmath79  @xcite , while the linear shear modulus of the network @xmath80 is equal to the differential modulus @xmath81 at @xmath82 .    , which has a limit of 1 at @xmath82 ; the fitting parameters for 3-chain model are given in table  [ thenetwork ] . ]
take actin and fibrin network under simple shear deformation as an example .
both 3-chain and 8-chain models fit the shear - experiment data from refs .
@xcite equally well , in fact  perfectly , as shown in fig .
[ constitutiverelation](a ) .
the stiffness , @xmath25 , and the initial end - to - end factor , @xmath83 , obtained by fitting with 3-chain model are a bit smaller than those in 8-chain model , but both in the semiflexible regime .
this is because the eight chains in a bcc are stretched equally ( in principal axes ) and share the deformation , while the three chains in a pc are stretched differently and the most stretched one contributes the most to the nonlinear elastic energy .
we believe the intrinsic heterogeneity in the 3-chain model , and the fact that the chain lying along the maximum principle stretch direction dominates the response of the whole cubic in semiflexible networks , is closer to the realistic case of filament network .
in addition , we shall see below that the 8-chain model can not explain negative normal stress .
hence we apply the 3-chain model in the following .
figure  [ constitutiverelation](b ) shows fits to experimental data for a wide variety of semiflexible filaments . in this figure
we plot the shear modulus @xmath84 instead of stress , because the fitting convergence is much better when the initial data section is close to 1 .
fitted parameters @xmath85 are listed in table  [ thenetwork ] , along with other parameters for each material that we list from the literature .
.fitting parameters @xmath85 for collagen , actin , vimentin , fibrin and neurofilament data obtained from @xcite . also shown
are the linear modulus @xmath0 extracted from the original data and used for scaling in fig .
[ constitutiverelation ] , the literature values of @xmath86 , and the calculated mesh size @xmath87 . [ cols="^,^,^,^,^,^",options="header " , ]     first of all , one may be surprised that the persistence length of collagen fibers is quoted as @xmath88 m , when there is a large body of literature that would claim that collagen fibers are stiff athermal rods with persistence length of centimeters .
this is all to do with the way a sample is prepared , and we use / quote the data  @xcite where the collagen was apparently less aggregated than in a typical extra - cellular matrix .
the same ambiguity will apply to actin networks as well , below .
another point to note about the fitted values in table  [ thenetwork ] is about neurofilament , which has @xmath89 , that is , rather flexible chains . taking the fitted value for @xmath90
, this leans that @xmath91 in this network was @xmath92 m , i.e. about 3.6 times longer than @xmath86 . by calculating @xmath34 for different filaments listed in table  [ thenetwork ]
, we see that the fitted @xmath37 is smaller than @xmath26 , indicating all of the filaments in fig .
[ constitutiverelation](b ) are pre - compressed , rather than pre - stretched in the equilibrium state of the network .
as table  [ thenetwork ] shows , no matter what the effective stiffness of the examined biofilaments , the ratio of fitted parameters @xmath93 remains close to 1 . later , when we discuss the network stability ( fig .
6 ) , it will become clear that all the networks we examined here lie very near the stability boundary .
it is not clear to us whether the fact that the mesh size is close to the filament persistence length is an unintended result of different crosslinking density in experiments  @xcite , or is a relevant and universal biological feature .
it is clear that networks of biological filaments have a great variety even within the same substance .
figure [ actin32 ] shows the published data for in - vitro crosslinked actin networks reported by different groups , all performing the simple - shear experiment ( the data is digitized from a reference given in the plot ) .
the stress - strain plotted in log - log format allows a clear identification of the linear regime @xmath94 ( the modulus varying between 95pa and 1pa for different sets ) , and the subsequent stiffening at higher shear .
all curves in fig .
[ actin32](a ) are fitted by the same eq .
( [ const3 ] ) with the linear modulus @xmath0 and the two parameters @xmath85 taking values : @xmath95 pa , @xmath96 ; @xmath97 pa , @xmath98 ; @xmath99 pa , @xmath100 ; @xmath101 pa , @xmath102 ; @xmath103 pa , @xmath104 and @xmath105 pa , @xmath106 for the six sets from top to bottom .
the difference between two data sets from @xcite is the density of heavy meromyosin ( hmm ) crosslinker ( labelled on the plot ) ; the difference between two data sets from @xcite is the degree of actin filament bundling : an initial network of f - actin filaments turned into a more sparse collection of bundles under repeated shear cycles @xcite . as a result
the network at ` cycle 7 ' has lower @xmath0 but higher stiffening .    , before the stiffening sets in at higher shear .
( b ) the stress - stiffening of actin networks represented by the @xmath107 scaling relation .
all data sets are the same as labelled in the plot ( a ) , and solid lines are theoretical curves plotted with eq .
( [ const3 ] ) using the same fitted parameters . ]    to find the onset of non - linearity in our theory , the differential shear modulus @xmath108 can be expanded in powers of shear strain @xmath74 : @xmath109     \label{ktaylor } \\ & & + 4 nk_\mathrm{b}t x^4 \frac{1 + x^2}{c \pi ( 1 - x^2)^5 } \ , \gamma^2 .
\nonumber\end{aligned}\ ] ] the crossover point when these two terms are comparable with each other , in a more stiff network with @xmath110 , can be approximated as @xmath111 .
the stress at this point is @xmath112 $ ] .
this crossover stress @xmath113 increases with the temperature as @xmath114 , closely matching the observations  @xcite . at the end of range ,
when the shear strain approaches the point of divergence in stress @xmath115 , ( @xmath116 in eq .
[ const3 ] ) , the stress and the differential modulus can be approximated as : @xmath117 ^ 2 } , \label{limit }   \\
\frac{k(\gamma)}{n k_\mathrm{b } t } & \simeq & \frac{2 ( 1 - x^4 ) \left[x^2 - ( 2 - 3 \gamma^2 ) x^4 + x^6\right ] } {   3 \pi c \left[1 - ( 2 + \gamma^2 ) x^2 + x^4\right]^3}. \label{limit2}\end{aligned}\ ] ] both expression diverge due to the same vanishing denominator , while maintaining the obvious scaling relation @xmath107 , which has been reported by different theories and experiments  @xcite .
we see this limit exposed clearly in fig .
[ actin32](b ) for very different actin networks . in this plot , the data points are the same as in fig .
[ actin32](a ) , and we also plot the predicted curves of @xmath118 obtained by differentiating eq .
( [ const3 ] ) , using parameters @xmath119 and @xmath37 from the fitting in fig .
[ actin32](a ) , for each data set .
one has to make a comment here , in the context of @xmath107 scaling and different experiments . in many cases , such as in unterberger et.al .
@xcite , the actin network was crosslinked by hmm and apparently retains some transient activity , producing the stress - softening and plasticity at higher stress ( this is also the case in  @xcite before the actin filaments formed bundles , or with f - actin crosslinked by filamin ) .
of course , our theory is not intended to deal with network plasicity ( we assumed all crosslinkes permanent ) and therefore we only retained the experimental data points in the early stress - stiffening regime to see the 3/2 scaling .    on the other hand , just by examining the actual data for the actin network of storm et.al .
@xcite , one might conclude that it strongly and systematically deviates from the 3/2 scaling .
in fact , the authors of  @xcite develop their own theory invoking various additional factors ( e.g. filament extensibility ) to account for this data .
however , our basic theory , assuming permanent crosslinks , bulk incompressibility and inextensible filaments , evidently fits both @xmath120 and @xmath118 data very well .
the fact that the data ( and the predicted curve ) do not appear to follow the 3/2 scaling is be due to the fact that , for these values of @xmath25 and @xmath37 , the final crossover to this characteristic scaling regime would occur at an even higher stress ( at which point the network would probably not survive in practice ) .
though simple shear deformation was helpful in this analysis , the important issue of normal stress remains controversial .
this is mainly because of the uncertain boundary conditions , see ref.@xcite for detail .
since the experiments reporting normal stress measurements are most commonly conducted in rotating cylindrical geometry of a standard rheometer   @xcite , we will consider this geometry and realistic boundary conditions to describe the response of the material in its normal direction when a shear is applied , as shown in fig .
[ poynting](a ) .
suppose the height and the radius of the undeformed cylinder are @xmath121 and @xmath122 , respectively ; on deformation they may become @xmath123 and @xmath124 .
incompressibility maintains @xmath125 . in the @xmath126
coordinate system , after rotating the top plate by the angle @xmath127 , the coordinates change as : @xmath128 : @xmath129 where @xmath130 denotes the shear strain , which is a function of the radial position in this parallel - plate geometry .
the strain invariants become : @xmath131/\lambda_h$ ] and @xmath132/\lambda_h^{2}$ ] . given the shear strain at the outermost surface , @xmath133 ,
the total free energy then becomes a function of the stretching ratio @xmath134 along the @xmath135 axis , after integration over radius : @xmath136 the equilibrium @xmath134 can be obtained by minimizing this free energy . when the cylinder radius becomes smaller ( or larger ) upon shear , with its height becoming correspondingly larger ( or smaller )  this phenomenon is called the positive ( or negative ) poynting effect @xcite . this geometric effect with stress - free top plate corresponds to the positive ( negative ) normal stress required to maintain the fixed plate separation in a more common rheometry experiment @xcite .
the 8-chain model fails in obtaining the correct normal stress in a sheared network : due to its core assumption that eight chains in a cell are identically deformed ( stretched ) upon shear , the elongation ratio of the network along the stretching direction is always larger than 1 , making the normal stress always positive .     at the outer surface of the cylinder ( black curve ) and the normal strain @xmath137 for a flexible network ( red curve ) with @xmath138 , and a more stiff network ( blue curve ) with @xmath139 . ]
figure  [ poynting](b ) shows the results for equilibrium @xmath134 , presented as a response to an oscillating imposed shear @xmath140 , for two model materials with different filament stiffness and pre - tension . a flexible network ( @xmath138 ) has a positive poynting effect , or positive normal stress is required to counter the expansion along the height ( in other setting , this is called the weissenberg effect ) . in a flexible network ( @xmath141 ) , when the mesh size is much smaller than the contour length of the subchains connecting the neighboring crosslinks , _
i.e. _ , @xmath142 , the entropic energy plays a main role . though chains are stretched along the principal extension direction in the shear geometry illustrated in fig .
[ poynting](a ) , the chains in other two principle directions are more likely being compressed , leading to the material contraction along the radial and circumferential directions . however , if the mesh size is comparable with the contour length of the subchains , @xmath143 , the material can behave in a negative poynting effect manner : the height of the cylindrical sample contracts , or negative normal stress is required to counter that and maintain the fixed height .
this is because the force acting on the chains directed along the principal extension in fig .  [
poynting](a ) is close to a divergence if @xmath144 , and it causes less energy when the material is compressed in the longitudinal direction , rather than stretched .
similar reason works also for a network of more stiff filaments ( see red curve for @xmath145 ) .
we did not specifically calculate the normal stress here ( this would be a cumbersome process involving the full tensor form of eq .
( [ strain ] ) , first fixing the pressure @xmath71 from the condition of zero radial stress on the free outer surface ) .
however , the magnitude of normal stress can be easily estimated from the linear relationship @xmath146 , which uses the fact that the normal strain is quite small ( i.e. the linear regime is justified ) and the young modulus is @xmath147 . taking the fibrin values in table  [ thenetwork ] as an example
, the normal stress is about @xmath148pa under a shear strain of @xmath149 .
these are very close to the observations in experiments @xcite .
the magnitude of @xmath150 pa for the same shear also accurately matches the actin results obtained in @xcite .
figure  [ map ] gives the full ` phase diagram ' of the stiffness - tension parameter space ( @xmath151 ) with phase boundaries separating positive / negative normal stress regions . in order to generate this diagram , for each parameter set ( @xmath152 )
we have calculated the value of @xmath153 by minimizing @xmath154 , under the given shear . in this way the boundary separating the positive ( @xmath155 , @xmath156 ) and the negative ( @xmath157 , @xmath158 ) normal stress regions in fig .  [ map ]
is calculated .
there is also a weak dependence of this boundary on the magnitude of the applied shear strain , which is due to the inherent non - linearity of stress - strain response ; however , we are not showing this in the figure to avoid clutter .
figure  [ map ] shows that a loose flexible network usually has a positive poynting effect , while a network of more stiff filaments has a negative poynting effect , especially when the filaments are crosslinked with increasing pre - tension .    , and filament pre - tension , @xmath37 , showing the boundaries of positive / negative poynting effect and the boundary of network stability ( the dashed line of the neutral filament , @xmath26 ,
was defined in fig .
[ fig_1 ] ) .
three ` softer ' filament networks from table [ thenetwork ] are shown in this map : @xmath159 - vimentin , @xmath160 - fibrin , and @xmath161 - neurofilament networks .
three actin networks from fig .
[ actin32 ] also fit on this map : @xmath162 - @xcite , @xmath163 - @xcite , and @xmath164 - @xcite . ]
one can see from figs .  [
constitutiverelation](b ) and [ actin32](b ) that the linear regime when @xmath94 persists for a different range of strain in different filaments ; the onset of hyperelastic regime is determined by how much pre - tension is in the filaments , i.e. how @xmath37 compares with @xmath165 from eq .
( [ achain ] ) .
the linear shear modulus @xmath80 can be easily obtained from our theory : in 3-chain model it is given by the first term in eq .
( [ ktaylor ] ) .
the marginal rigidity condition @xmath166 determines the strand pre - tension that is required for achieving a stable network .
the stability criterion here is : @xmath167 note that the expression in the right hand side is always greater than @xmath168 as defined in section 2 .
flexible chains with @xmath169 are therefore always stable in the rubbery network , that is , @xmath170 always .
on the other hand , stiff filaments have to be crosslinked with @xmath171 exceeding the pre - tension threshold given by eq .
( [ stabil ] ) , which turns out to be slightly lower than @xmath165 defined in fig .  [ fig_1](b ) . in other words
, there have to be tensile forces acting on the crosslinked filaments in the network in order for it to be mechanically stable with a non - zero shear modulus .
this notion is familiar from the `` tensegrity '' concept in biology and engineering  @xcite . for stiff athermal filament network ,
the window of pre - tension between the linear modulus @xmath172 at @xmath173 , and @xmath174 at @xmath110 ( we assume inextensible chains ) is very narrow .
the stability boundary of the network is plotted in fig .
the condition for @xmath172 gives the equilibrium case labelled as @xmath175 in the phase diagram .
however , the full analysis shows that the magnitude of the shear strain modifies the stability condition , as represented by the coloured curves for @xmath176 and @xmath177 .
this shift means that the region of mechanical stability of filament network expands on increasing deformation , which matches exactly what the recent paper by sharma _
states  @xcite .
to summarize , in this work we develop a continuum elastic theory of a network of semiflexible filaments , by implementing the general free energy of one semiflexible chain into that of a disordered network . on reflection
, we choose the 3-chain model as most closely matching the realistic system  and achieve various quite stringent fits a number of different experimental data sets over the full range of nonlinear stress - strain range .
we demonstrate that the general theory produces the conditions for positive / negative poynting effect ( normal stress ) in a network under imposed shear .
the greatest weakness of this model is omitting the effects of local strain non - affinity , which might play an important role when crosslink density of the network is small and filaments stiff .
the effect of tensegrity , or linking of filament pre - tension to the network stability and the magnitude of the shear modulus @xmath178 is an unexpected result of this theory .
the stringent match of experimental data is reassuring , and we believe that the presented analytical continuum model can provide as an efficient and portable tool for studying the mechanical properties of semiflexible networks . +
* acknowledgments : *   this work has been funded by the tcm critical mass grant from epsrc ( ep / j017639 ) .
we are grateful for informative discussions with masao doi and zhongcan ouyang , and thank cornelius storm for providing their raw experimental data @xcite .
b.  wickstead and k.  gull , _ j. cell .
biol . _ , 2011 , * 194 * , 513525 g.  faury , _ pathol .
_ , 2001 , * 49 * , 310325 s.  zbek , p.
g. balasubramanian , r.  chiquet - ehrismann , r.  p. tucker and j.  c. adams , _ mol .
biol . cell _ , 2010 , * 21 * , 43004305 k.  e. kasza , a.  c. rowat , j.  liu , t.  e. angelini , c.  p. brangwynne , g.  h. koenderink and d.  a. weitz , _ curr .
cell biol .
_ , 2007 , * 19 * , 101107 p.
bursac , g.  lenormand , b.  fabry , m.  oliver , d.  a. weitz , v.  viasnoff , j.  p. butler and j.  j. fredberg , _ nature mater .
_ , 2005 , * 4 * , 557561 d.  a. fletcher and r.  d. mullins , _ nature _ , 2010 , * 463 * , 485492 x.  trepat , g.  lenormand and j.  j. fredberg , _ soft matter _ , 2008 , * 4 * , 17501759 t.  lecuit , p .- f .
lenne and e.  munro , _ annu .
biol . _ , 2011 , * 27 * , 157184 e.  w. wong , p.  e. sheehan and c.  m. lieber , _ science _ , 1997 , * 277 * , 19711975 m.  r. falvo , g.  j. clary , r.  m. taylor , v.  chi , f.  p. brooks , s.  washburn and r.  superfine , _ nature _ , 1997 , * 389 * , 582584 m.  fixman and j.  kovac , _ j. chem . phys . _ , 1973 , * 58 * , 15641568 j.  f. marko and e.  d. siggia , _ macromolecules _ , 1995 , * 28 * , 87598770 b.  y. ha and d.  thirumalai , _ j. chem . phys . _ , 1996 , * 106 * , 14 c.  bouchiat , m.  d. wang , j .- f .
allemand , t.  strick , s.  m. block and v.  croquette , _ biophys .
j. _ , 1999 ,
* 76 * , 409413 j.  r. blundell and e.  m. terentjev , _ macromolecules _ , 2009 , * 42 * , 53885394 f.  c. mackintosh , _ cytoskeletal mechanics : models and measurements _ , cambridge university press , cambridge , 2006 c.  p. broedersz and f.  c. mackintosh , _ rev .
_ , 2014 , * 86 * , 9951036 d.  a. head , a.  j. levine and f.  c. mackintosh , _ phys . rev
. lett . _ , 2003 , * 91 * , 108102 c.  heussinger , b.  schaefer and e.  frey , _ phys .
e _ , 2007 , * 76 * , 031906 m.  doi and s.  f. edwards , _ the theory of polymer dynamics _ , oxford university press , oxford , 1986 d.  c. morse , _
macromolecules _ , 1999 , * 32 * , 59345943 c.  storm , j.  j. pastore , f.  mackintosh , t.  lubensky and p.  a. jamney , _ nature _ , 2005 , * 435 * , 191194 j.  wilhelm and e.  frey , _ phys .
_ , 1996 , * 77 * , 25812584 j.  wilhelm and e.  frey , _ phys . rev
_ , 2003 , * 91 * , 108103 j.  s. palmer and m.  c. boyce , _ acta biomater . _ , 2008 ,
* 4 * , 597612 m.  j. unterberger , k.  m. schmoller , a.  r. bausch and g.  a. holzapfel , _ j. mech
_ , 2013 , * 22 * , 95114 m.  l. gardel , j.  h. shin , f.  c. mackintosh , l.  mahadevan , p.  matsudaira and d.  a. weitz , _ science _ , 2004 ,
* 304 * , 13011305 t.  van  dillen , p.  onck and e.  van  der giessen , _ j. mech
solids _ , 2008 , * 56 * , 22402264 g.  agar , p.
r. onck and e.  van der giessen , _ macromolecules _ , 2011 , * 44 * , 70267033 p.  a. janmey , m.  e. mccormick , s.  rammensee , j.  l. leight , p.  c. georges and f.  c. mackintosh , _ nature mater . _ , 2007 , * 6 * , 4851 h.  kang , q.  wen , p.  a. janmey , j.  x. tang , e.  conti and f.  c. mackintosh , _ j. phys . chem .
b _ , 2009 , * 113 * , 37993805 a.  sharma , a.  j. licup , k.  a. jansen , r.  rens , m.  sheinman , g.  h. koenderink and f.  c. mackintosh , _ nature phys . _ , 2016 ,
* 12 * , 584587 m.  otto , j.  eckert and t.  a. vilgis , _ macromol . theory sim . _ , 1994 , * 3 * , 543555 d.  c. morse , _ macromolecules _ , 1998 , * 31 * , 70307043 o.  kratky and g.  porod , _ recueil des travaux chimiques des pays - bas _ , 1949 , * 68 * , 11061122 m.  c. boyce and e.  m. arruda , _ rubber chem .
_ , 2000 , * 73 * , 504523 l.  r.  g. treloar , _ the physics of rubber elasticity _
, oxford university press , oxford , 1975 a.  f. bower , _ applied mechanics of solids _ , crc press , boca raton , 2009 c.  o. horgan , _ intl . j. non - lin .
mech . _ , 2015 , * 68 * , 916 a.  zaccone , j.  r. blundell and e.  m. terentjev , _ phys . rev .
b _ , 2011 , * 84 * , 174119 c.  p. goodrich , a.  j. liu and s.  r. nagel , _ nature phys . _ , 2014 , * 10 * , 578581 a.  m. stein , d.  a. vader , l.  m. jawerth , d.  a. weitz and l.  m. sander , _ j. microscopy _ , 2008 , * 232 * , 463475 f.  gittes , b.  mickey , j.  nettleton and j.  howard , _ j. cell biol .
_ , 1993 , * 120 * , 923934 n.  mcke , l.  kreplak , r.  kirmse , t.  wedig , h.  herrmann , u.  aebi and j.  langowski , _ j. mol
. biol . _ , 2004 , * 335 * , 12411250 o.  i. wagner , s.  rammensee , n.  korde , q.  wen , j.  f. leterrier and p.  a. janmey , _ exp .
_ , 2007 , * 313 * , 22282235 k.  schmoller , p.  fernandez , r.  arevalo , d.  blair and a.  bausch , _ nature comm . _ , 2010 ,
* 1 * , 134 p.  h.  j. kouwer , m.  koepf , v.  le sage , m.  jaspers , a.  m. van buul , z.  h. eksteen - akeroyd , t.  woltinge , e.  schwartz , h.  j. kitto , r.  hoogenboom , s.  j. picken , r.  j.  m. nolte , e.  mendes and a.  e. rowan , _ nature _ , 2013 , * 493 * , 651655 m.  l. gardel , j.  h. shin , f.  c. mackintosh , l.  mahadevan , p.  a. matsudaira and d.  a. weitz , _ phys
lett . _ , 2004 , * 93 * , 188102 y.  c. lin , n.  y. yao , c.  p. broedersz , h.  herrmann , f.  c. mackintosh and d.  a. weitz , _ phys .
lett . _ , 2010 , * 104 * , 058101 y.  c. lin , c.  p.
broedersz , a.  c. rowat , t.  wedig , h.  herrmann , f.  c. mackintosh and d.  a. weitz , _ j. mol .
biol . _ , 2010 , * 399 * , 637644 c.  o. horgan and j.  g. murphy , _ j. elast .
_ , 2010 , * 98 * , 205221 l.  a. mihai and a.  goriely , _ proc .
a _ , 2011 ,
* 467 * , 36333646 j.  h. poynting , _ proc .
a _ , 1909 , * 82 * , 546559 d.  e. ingber , _ annu .
physiol . _ , 1997 , * 59 * , 579599