rare decays stand among the most powerful probes of physics beyond the standard model ( sm ) .
indeed , most new physics ( np ) scenarios suffer from severe constraints due to their potentially large contributions to flavor observables .
this has motivated an intense experimental search , with special focus on observables which are strongly suppressed in the sm .    in 2013 ,
the lhcb collaboration reported on the measurement of several observables in processes involving @xmath3 transitions .
a significant tension with the sm was found in some cases .
these include angular observables in @xmath4  @xcite , particularly large in case of the popular @xmath5 @xcite , as well as a decrease , with respect to the sm expectation , in several branching ratios  @xcite .
these anomalies received immediate attention in the flavor community , and soon several independent global fits @xcite showed that the tension could be alleviated in the presence of new physics contributions .
interestingly , in 2014 the lhcb collaboration also found an indication of lepton universality violation in the theoretically rather clean ratio  @xcite @xmath6 measured in the low dilepton invariant mass regime , which departs from the sm result @xmath7 by @xmath8  @xcite .
again , this led to some excitement in the community , in particular after it was found that this hint is compatible with the previous anomalies in @xmath3 transitions : they can be explained by the same type of np contributions to muonic operators @xcite .
so far , the discussion has focused on results obtained by lhcb with an integrated luminosity of @xmath9 fb@xmath10 .
the latest chapter of this story is the recent announcement of the lhcb collaboration of new results using the full lhc run i dataset  @xcite , with an integrated luminosity of @xmath11 fb@xmath10 , which confirm the robustness of the lhcb data .
indeed , the new results are compatible with those found with @xmath9 fb@xmath10 , and several theorists have used the new data to update their analyses . according to @xcite , the hypothesis of non - zero np contributions is preferred over the sm by @xmath12 , @xmath13 if the 2014 measurement of @xmath14 is included , whereas the analysis of @xcite finds a slightly larger statistical significance for np ( slightly larger than @xmath15 ) even in the absence of @xmath14 . in any case , assuming that hadronic effects @xcite are not behind these anomalies ( something impossible in case of @xmath14 ) , there is clear evidence of new physics in @xmath16 meson decays . for a complete review of the subject see @xcite and references therein .
several approaches have been considered in order to explain the @xmath17 anomalies .
some papers @xcite consider generic @xmath2 bosons with flavor violating couplings .
these have been shown to provide a simple way to reconcile theory predictions with experimental data .
there are also some works that were built to address the first anomalies but fail to address the lepton universality violating @xmath14 measurement , see for example refs .
@xcite in the context of 331 scenarios or ref .
@xcite , where 4-quark scalar interactions were considered .
finally , a few recent models can also account for the lepton universality violating @xmath14 measurement .
one finds three types of working models : models with scalar or vector leptoquarks ( or similarly , r - parity violating supersymmetry , where the squarks can play the role of the leptoquarks )  @xcite , composite higgs models  @xcite and models with a @xmath2 .
the latter is generally considered the easiest approach , since the anomalies can be solved with new contributions to vectorial operators . however , to the best of our knowledge , the only ( complete ) models that have been put forward in order to account for these anomalies using a @xmath2 are : ( 1 ) the one introduced in ref .
@xcite ( and the two - higgs - doublet version in @xcite ) , which makes use of a @xmath18 gauge symmetry , and ( 2 ) the one proposed in the recent paper @xcite , which extends the horizontal gauge symmetry to the quark sector .
the existence of dark matter ( dm ) is a well established evidence of new physics , supported by a plethora of astrophysical and cosmological observations .
this motivates further extensions of the sm of particle physics that include valid dm candidates . in dm models ,
the current planck @xmath19 limits for the dm relic density , @xmath20 @xcite , typically imply strong constraints on the mass and couplings of the dm candidate .
this , combined with additional constraints ( such as those from flavor physics ) , sometimes leads to very predictive scenarios where the parameter space of the dm model shrinks to small regions compatible with all observations .
we propose a simple model which captures the main ingredients required to explain the @xmath3 anomalies and , simultaneously , provides a dm candidate . in order to do so
, we extend the sm with a dark sector charged under an additional @xmath1 gauge symmetry .
the spontaneous breaking of this symmetry gives rise to a massive @xmath2 boson , which communicates the sm particles with the dark matter particle and solves the @xmath3 anomalies with contributions to the relevant flavor observables .
the interplay between dm and flavor physics leads to a very constrained scenario .
in particular , achieving the measured dm relic density while providing the required np contributions to @xmath16 meson decays turns out to be very restrictive on the model parameters .
the rest of the paper is organized as follows : in sec .
[ sec : general ] we discuss some general aspects about the @xmath0 anomalies and introduce the basic language to be used throughout the paper . in sec .
[ sec : model ] we define our model , show how it addresses the @xmath3 anomalies while providing a dm candidate and discuss extensions to generate non - zero neutrino masses .
[ sec : dm ] is devoted to the dm phenomenology of our model and the @xmath2 portal responsible for its production in the early universe .
[ sec : constraints ] is a review of the most relevant constraints in our model , to be considered in sec .
[ sec : results ] , where we present our numerical results . finally , in sec .
[ sec : conclusions ] we summarize our results , discuss some related aspects and derive some general conclusions .
in this section we discuss some general aspects about the @xmath3 anomalies . for a complete review of the subject ,
we refer the reader to @xcite .      the effective hamiltonian for @xmath3 transitions is usually written as @xmath22 where @xmath23 is the fermi constant , @xmath24 the electric charge and @xmath25 the ckm matrix . @xmath26 and @xmath27 are the effective operators that contribute to @xmath3 transitions , and @xmath28 and @xmath29 their wilson coefficients .
since the most important anomalies have been found in semileptonic @xmath16 meson decays , we will consider the following set of operators , @xmath30 here @xmath31   and @xmath32 , for muons . ] .
it is also customary to split the wilson coefficients in two parts : the sm contributions and the np contributions . since the primed operators , @xmath33 and @xmath34
, do not receive significant sm contributions , this is usually applied only to the unprimed wilson coefficients , which can be written as @xmath35 the sm contributions , @xmath36 and @xmath37 , have been computed by different groups . assuming that these are the only contributions to the wilson coefficients
, several independent global fits have found a sizable tension with experimental data on @xmath3 transitions .
this motivates the addition of np contributions .
when this is done , the global fits are clearly improved . according to @xcite ( we will mainly consider the results of this global fit ) ,
the best improvements are found in two cases :    * * scenario 1 : * np provides a negative contribution to @xmath38 , with @xmath39 , leading to a wilson coefficient @xmath40 significantly smaller than the one in the sm . *
* scenario 2 : * np enters in the @xmath41 invariant direction @xmath42 , with @xmath43 .    in both cases ,
the rest of operators involving muons are perfectly compatible with the sm expectations .
similarly , no np is required for operators involving electrons or tau leptons .      once determined the type of contributions to @xmath44 and @xmath45
a np model has to induce , one can figure out the main ingredients of a simple working model .
arguably , the simplest one contains the following elements :    * a massive @xmath2 boson , responsible for the vectorial operators @xmath46 and @xmath47 * the @xmath2 must have flavor violating couplings to quarks * the @xmath2 must couple differently to electrons and muons    this setup can be easily parameterized by the lagrangian @xcite @xmath48 in order to account for the anomalies , one requires @xmath49 and either ( 1 ) @xmath50 , or ( 2 ) @xmath51 and @xmath52 , depending on the scenario one wants to consider . the rest of the @xmath2 couplings to
sm fermions can be set to zero . as we will see in sec .
[ sec : model ] , the model we are going to consider belongs to scenario ( 2 ) . in this case , the quark and lepton currents contributing to @xmath46 and @xmath47 are both left - handed and one finds at tree - level  @xcite @xmath53 with @xmath54 where @xmath55 is the fine structure constant .
note that @xmath56 and the ckm elements appear in eq . in order to normalize the wilson coefficients as defined in eqs . and .
+ we now introduce a complete renormalizable model with these properties .
we extend the sm gauge group with a new dark @xmath57 factor , under which all the sm particles are assumed to be singlets .
the only particles charged under the @xmath57 group are the following vector - like fermions , @xmath58 as well as the complex scalar fields @xmath59 where we denote the gauge charges under @xmath60 and the @xmath41 doublets can be decomposed as @xmath61 and @xmath62 .    besides canonical kinetic terms
, the new vector - like fermions have dirac mass terms , @xmath63 as well as yukawa couplings with the sm fermions @xmath64 where @xmath65 and @xmath66 are @xmath11 component vectors .
the scalar potential takes the form @xmath67 here @xmath68 is the sm higgs doublet and latexmath:[$\mathcal v_{\text{sm } } = m_h^2    involving the @xmath57 charged scalars are @xmath70 and @xmath71 we will assume that the scalar potential is such that only the standard higgs boson and the @xmath72 field acquire non - zero vacuum expectation values ( vevs ) , @xmath73 therefore , the @xmath72 field will be responsible for the spontaneous breaking of @xmath57 , giving a mass to the @xmath2 , @xmath74 , where @xmath75 is the @xmath57 gauge coupling , and inducing mixings between the vector - like fermions and their sm counterparts thanks to the yukawa interactions in eq . .
furthermore , after spontaneous symmetry breaking , the resulting lagrangian contains a remnant @xmath76 symmetry , under which @xmath77 is odd and all the other fields are even .
therefore , @xmath77 is a stable neutral scalar , and thus a potentially valid dm candidate .
it is worth noting that the mechanism to stabilize the dm particle does not introduce additional ad - hoc symmetries , but simply makes use of the same @xmath57 symmetry that is required in order to give an explanation to the lhcb observations .
this goal has been achieved by breaking the continuous @xmath57 symmetry to a remnant @xmath76 , something that can be easily accomplished with a proper choice of @xmath57 charges  @xcite .    before concluding this section
we must comment on @xmath1 mixing .
it is well known that nothing prevents @xmath1 factors from mixing . in the model under consideration , this would be given by the lagrangian term @xcite @xmath78 where @xmath79 are the usual field strength tensors for the @xmath80 groups .
in the presence of a non - zero @xmath81 coupling , kinetic mixing between the @xmath57 and @xmath82 gauge bosons is induced . as a consequence of this ,
the physical @xmath2 boson would couple to all particles that carry hypercharge , this is , to all the sm fermions .
this would lead to phenomenological problems since couplings to the first generation are strongly constrained .
therefore , we will assume that the tree - level @xmath81 coupling vanishes .
this is easily justified in our model because this term is not induced via renormalization group running if it is zero at some high - energy scale ( where one may speculate about a ultraviolet completion ) .
the reason is our choice of the @xmath57 charges  @xcite .
in addition , we must keep the 1-loop induced @xmath81 coupling , generated in loops including heavy vector - like quarks and leptons , under control .
we find @xmath83 therefore , @xmath84 would ensure small @xmath81 couplings , while allowing large @xmath75 . something required by dm constraints ( see sec . [
sec : dm ] ) .       and
@xmath47 in our model . ]
this model solves the @xmath3 anomalies in a similar fashion as the one in ref .
the @xmath2 couplings to the sm fermions are generated after their mixing with the corresponding vector - like quarks and leptons , as shown in fig .
[ fig : couplings ] . neglecting @xmath85 and @xmath86 ,
the resulting @xmath2 couplings are found to be @xmath87 in case @xmath88 , the @xmath2 coupling to a pair of sm quarks can be further approximated to @xmath89 .
these results can be combined with eq . in order to determine the allowed ranges for the model parameters that explain the @xmath3 anomalies found by lhcb .
is a new scalar field with @xmath57 charge @xmath90 , necessary in order to make the operator gauge invariant . ]
we can further extend the model to get non - zero masses for the sm neutrinos .
this can be done trivially by adding new particles , singlets under @xmath57 , which mediate the standard mechanisms .
for example , the addition of right - handed neutrino singlets , @xmath91 , allows for the usual type - i seesaw mechanism .
it is , however , more interesting to consider a mechanism that involves the @xmath57 sector of the model .
this can be done by means of the effective operator @xmath92 as shown in fig .
[ fig : numass ] . here
@xmath93 is a new scalar field with with @xmath57 charge @xmath94 , necessary in order to make the operator gauge invariant . an example model that can serve as ultraviolet completion of @xmath95
is obtained with the addition of the scalar @xmath96 , together with a vector - like ( dirac ) fermion @xmath97 .
this allows for the yukawa couplings @xmath98 and @xmath99 , which lead to @xmath100 after integrating out @xmath101 and @xmath102 .
although @xmath93 must get a non - zero vev in order to break lepton number and generate neutrino masses , we note that our choice @xmath94 guarantees that the remnant @xmath76 symmetry that stabilizes the dm particle @xmath77 is preserved .
at high temperatures @xmath77 attains thermal equilibrium in the heat bath through its reactions with the different degrees of freedom to which it couples . at decoupling ,
the dominant processes determining the dm yield are the following @xmath103 reactions : , @xmath104 plays a subdominant role . ]
* @xmath105 ( higgs portal ) : processes enabled by the @xmath106 coupling . * @xmath107 , with @xmath101 standing for the sm and the new vector - like quarks and leptons ( @xmath2 portal ) : the former enabled by @xmath108 mixing and @xmath57 coupling , while the latter solely by @xmath75 .    for very heavy dm with @xmath109 , also @xmath110 can take place .
depending on the relative size of these processes ( pure scalar and @xmath2-mediated ) one can then distinguish several scenarios . among them
those solely involving @xmath2-mediated processes , are  arguably  expected to be dominant ( they are driven by a gauge coupling ) .
interestingly enough , in that case a clear correlation with flavor physics must exist .
note that these processes match those in fig .
[ fig : couplings ] if one trades one of the fermion pairs for @xmath111 . therefore , under the fairly reasonable assumption that the @xmath2-mediated processes play a dominant role , an interplay between flavor and dm physics is possible establishing .
thus , in turn , further constraining the flavor - transition parameters through the restrictions imposed by the condition of generating the correct dm relic density ( see the useful ref .
@xcite ) .
although our results rely on micromegas @xcite , a simple analytical discussion is worth doing .
the cross section for the @xmath107 processes can be estimated to be @xmath112 where the @xmath113 and @xmath114 is a kinematic function .
depending on the relative size of @xmath115 and @xmath116 ( @xmath117 ) , one can distinguish two regimes . for @xmath118 ,
dm annihilation processes involve dominantly sm quarks and leptons ( @xmath119 and @xmath120 ) . in that case , however , the corresponding cross sections are suppressed by chiral / vector - like mixing ( see eq .
( [ eq : delta - parameters ] ) ) , annihilation is rather inefficient and therefore leads to an overpopulation of @xmath77 scalars . for @xmath121
, processes involving the vector - like fields are also possible , and since they are not suppressed by mixing parameters they can lead to the appropriate dm relic abundance .
note that if the dm and vector - like fermion mass splitting is small , then the cross section will be phase space suppressed : @xmath122 . in this case
, resonant annihilation is needed to efficiently deplete the scalar @xmath77 population ( see fig .
[ fig : c9 dm ] ) . for large mass splittings , in contrast , annihilation is very efficient , and so the correct relic density can be readily obtained .
however , reconciling this `` scenario '' with flavor constraints turns out to be tricky .
avoiding the resonance requires @xmath123 , which means @xmath124 , and a large @xmath75 close to 1 . in that limit ,
@xmath125 which for typical vector - like quark masses , @xmath126  tev , implies large chiral / vector - like quark mixing ( order one @xmath127 couplings ) in order to be compatible with the lhcb observations .
this , however , is forbidden by quark flavor constraints ( see next section ) .    in summary , the correct dm relic density can be easily produced within our setup provided @xmath128 and dm annihilation proceeds resonantly .
resonant annihilation can be avoided for @xmath129 , but finding spots in parameter space consistent with both , dm and flavor constraints seems challenging . the right relic density might as well be produced when @xmath118 , but probably this would require being sharply at the @xmath2 resonance .
there are constraints on the mass of additional gauge bosons coupled to sm states .
the ballpark of these limits is about 2.53.0  tev for @xmath1 extensions which predict a coupling of the @xmath130 to light quarks of @xmath131 @xcite .
however , in our model these couplings are suppressed by the small mixing between the sm quarks and the vector - like states .
thus , we nearly get any constraint on the mass of the @xmath130 from lhc searches .
in addition , we are going to assume in the lepton sector only a sizable mixing for muons with the new states .
therefore , we get also hardly any limit on @xmath132 from lep searches .    thus , we can safely assume in the following that the @xmath2 couplings to the first generation of sm fermions can be neglected and evade all lep and lhc limits .
similarly , @xmath2 contributions to flavor observables related to the first generation can be safely neglected .
let us now review several relevant constraints in our model .
the masses of the vector - like quarks have strong bounds from the lhc .
being colored particles , vector - like quarks can be efficiently produced in @xmath133 collisions , which typically pushes their masses towards the tev scale , see @xcite and references therein .
since our setup works with vector - like quarks with masses at or above the tev scale , these bounds are easily satisfied .
vector - like leptons can be searched for at the lhc in the standard multi - lepton channels . using the cms analysis in @xcite , based on searches for final states including @xmath11 or more leptons with an integrated luminosity of about @xmath134  fb@xmath10 at @xmath135  tev , ref .
@xcite obtained a lower limit on the mass of doublet vector - like leptons of about @xmath136  gev , in case they decay to electrons or muons , and about @xmath137  gev , in case they decay to tau leptons .
more recently , ref .
@xcite considered the analogous atlas multi - lepton search @xcite , with @xmath138  fb@xmath10 at @xmath135  tev , also looking for final states with @xmath11 or more leptons .
the results were similar to the ones derived from the cms analysis , with lower limits on the mass of doublet vector - like leptons of about @xmath139  gev .
finally , one can also derive limits from @xmath140 searches at the lhc .
however , according to @xcite , these are rather mild , only relevant for @xmath2 masses below @xmath141  gev .
contrary to the model in @xcite , where the muonic current contributing to the @xmath3 observables is purely vectorial ( @xmath142 ) , in our model the current is left - handed and thus contains an axial vector contribution as well ( @xmath143 ) .
this is relevant for the @xmath144 decay , especially sensitive to axial - vector currents .
currently , there is a little tension between the sm prediction for the average time - integrated branching ratio , @xmath145 , and the cms and lhcb measurements .
while the sm prediction is found to be @xmath146  @xcite , the combination of the cms and lhcb measurements leads to @xmath147  @xcite . in view of this little deficit
, the new @xmath2 mediated contribution might be potentially welcome .
the @xmath144 decay is known for its reach in probing new physics scenarios ( see for example the recent papers @xcite ) . in the sm , the amplitude is dominated by the axial vector leptonic current , and thus by the @xmath148 wilson coefficient ( @xmath149 ) . similarly , in the model under consideration , the main np contribution is given by a modification of @xmath148 .
we note that vector contributions vanish for lepton flavor conserving channels .
therefore , our prediction for the average time - integrated branching ratio is simply given by @xmath150 nevertheless , the large experimental error in @xmath145 precludes from obtaining useful bounds . using the sm and the experimental average time - integrated branching ratios given above
, one finds @xmath151 ( at the @xmath152 level ) . by combining @xmath153 with the sm values for @xmath40 and @xmath154
, one finds that in our setup @xmath155 , which is slightly above the @xmath156 limit , and thus perfectly compatible with the experimental measurement of @xmath145 at the @xmath157 level .
another relevant flavor constraint comes from @xmath158 mixing , induced at tree - level by @xmath2 exchange as soon as a non - zero @xmath159 coupling is considered .
allowing for a @xmath160 deviation from the sm expectation in the mixing amplitude , @xmath161 , one finds  @xcite @xmath162      in principle , the mixings with the vector - like leptons can induce lepton flavor violating ( lfv ) processes of the type @xmath163 , mediated by the @xmath2 boson .
however , one can suppress them with a proper parameter choice . since all
@xmath3 anomalies can be simultaneously explained with ( only ) new contributions to @xmath32 , @xmath164 is required , but one can choose @xmath165 .
this would eliminate the coupling of the @xmath2 boson to electrons and tau leptons , and thus all lfv processes mediated by the @xmath2 boson .
this includes lfv in @xmath16 decays , recently suggested in @xcite and further studied in @xcite .    besides the anomalies in @xmath3 transitions
, the lhc might have found additional hints for new physics .
recently , the cms collaboration found a 2.4@xmath166 excess in the @xmath167 channel which translates into br@xmath168%  @xcite .
for this result , the collaboration made use of the 2012 dataset taken at @xmath169  tev with an integrated luminosity of 19.7  fb@xmath10 .
this large higgs lfv branching ratio can not be accommodated in our setup , since vector - like leptons are known to be unable to reach such lfv rates in higgs decays due to limits from the radiative @xmath170 decay  @xcite .
however , as suggested in @xcite and recently confirmed in several works @xcite , a simple extension with a second higgs doublet suffices to explain the cms hint .
finally , we comment on limits from the ( former ) non - observation of lepton universality volation .
indeed , very strong bounds have been derived in many scenarios due to the non - observation of lepton universality violating effects in , for example , pion and kaon decays ( see for example @xcite ) .
however , these are absent in our model due to the suppression of the vector - like and first generation sm quark mixing .
there are several precision measurements in the lepton sector which might be potentially sensitive to the @xmath2 interactions considered in this model .
first of all , there is the so - called _ neutrino trident production _  @xcite .
this is the production of a muon anti - muon pair by scattering of muon neutrinos in the coulomb field of a target nucleus .
the cross - section gets additional contributions in the presence of a new @xmath2 boson that couples to the muons , and these can be used to constrain @xmath171 and @xmath132 .
the impact of this bound on @xmath2 models designed to solve the @xmath3 transitions was studied in @xcite . using the ccrf measurement of the trident cross - section  @xcite
, one can derive bounds on the mass of the @xmath2 boson and its coupling to muons . for scenarios with left - handed @xmath2 coupling to muons , one finds  @xcite @xmath172    finally , there are two other constraints related to precision measurements in leptonic observables : @xmath173 and the @xmath174 boson couplings to the muon ( which are modified at the 1-loop level due to the @xmath2 interactions ) . however , in @xcite these are found to be relatively mild . for this reason we will not consider them in our numerical analysis .
for the numerical analysis we have implemented the model into the mathematica package sarah@xcite and generated fortran code to link this model to spheno@xcite .
the big advantage of this setup is that we can cross - check our tree - level approximations with a full numerical calculation .
moreover , based on the flavorkit interface @xcite , the generated sphenoversion does not only calculate the tree - level contributions to the wilson coefficients but also all 1-loop corrections .
this way , we can see in what parameter range the tree - level approximation works well and where it breaks down .
in general , one can expect loop effects to become important when the mixing between sm quarks and the vector - like states is sizable : this leads to significant changes in the sm - like box and penguin contributions to @xmath44 and @xmath45
. this alone would not be necessarily a problem .
however , the same mixing influences also loop contributions to the other wilson coefficients , in particular @xmath175 and @xmath176 .
these effects have to be small in order to satisfy the results of global fits to @xmath3 data . + as second step , we used the calchep @xcite output of sarahto implement the model in micromegas to calculate the relic density . interfacing the parameter values between sphenoand micromegas
is done via the slha+ functionality of calchep @xcite .     and @xmath177 ( dashed black ) in the @xmath178 plane . for @xmath179
we show the full 1-loop results via the black lines , while the dotted grey lines give the values using the tree - level approximation only . ]     ( full black lines ) and @xmath180 ( dotted grey lines ) in the @xmath178 plane . ]
our main numerical results are summarized in figs .
[ fig : c9 dm ] and [ fig : c7c9 ] . in fig .
[ fig : c9 dm ] we show the dark matter relic density in the @xmath181 plane , together with the ratio @xmath179 .
the other important parameters entering this calculation have been fixed to @xmath182 first of all , we see that with these parameters the model is perfectly compatible with the constraints discussed in section [ sec : constraints ] .
furthermore , there is a region for moderately large @xmath183 where the dm constraint can be satisfied and @xmath184 holds . to reduce the relic density to @xmath185 one has to be rather close the the @xmath130 resonance around 2  tev . in the same figure
we show also @xmath186 when using the tree - level approximation only . in the interesting region
where the flavor anomalies can be explained , this approximation is working quite well . however , for decreasing @xmath75 one would expect smaller values for @xmath186 at tree - level , while we find a large increment .
the reason is that for constant @xmath132 and smaller @xmath75 , @xmath187 is increasing , leading to a larger mixing in the quark sector with the vector - like states .
thus , loop effects become important and quickly dominate the behavior .
these effects are not only expected to show up in @xmath188 but also in the other wilson coefficients as stated above . to demonstrate this , we plot in fig .
[ fig : c7c9 ] the contours for @xmath189 and @xmath190 in the same plane .
obviously , we find a similar enhancement for small @xmath75 for both coefficients .
however , for the most interesting region where tree - level contributions to @xmath188 dominates and can explain the anomalies , the change in @xmath175 is very moderate : @xmath191 .
we have presented a dm model that successfully addresses the anomalies in @xmath3 transitions recently found by the lhcb collaboration . in our setup , @xmath16 mesons decay through the @xmath2 into the dark sector , which mixes with the sm via yukawa couplings .
an interesting connection between dm and flavor physics emerges , leading to a very constrained scenario where the np masses and parameters are restricted to lie in thin regions of the full parameter space .
the model offers several interesting phenomenological possibilities . due to the mixing between the sm muons and the vector - like lepton doublet ,
our model predicts a reduction of the @xmath192 coupling , @xmath193 , with respect to the sm prediction .
the most direct probe of this coupling is , of course , the higgs boson decay into a pair of muons , @xmath194 , experimentally very challenging in a hadron collider due to the expected low branching ratio ( @xmath195 in the sm ) . in 2013 ,
the atlas collaboration presented results based on an integrated luminosity of @xmath196 fb@xmath10 in collisions at @xmath135 tev , without any evidence of a signal  @xcite . for @xmath197 gev
, they obtained an upper limit on the @xmath198 signal strength , @xmath199 ( at @xmath200 cl ) . assuming a sm - like higgs boson production
, this limit translates into @xmath201 .
obviously , our model is well within the limit , since it predicts a reduction of the coupling .
nevertheless , it is worth pointing out that a precise determination of this coupling in a future linear collider would be a strong test of the model   coupling with an uncertainty of about @xmath202 . ] .
let us also comment on the mixings between the vector - like fermions and the sm ones .
these are determined by free yukawa parameters which , in principle , are naturally expected to be of the same order for all sm generations . however , in order to avoid conflict with existing experimental data , one is forced to strongly suppress the mixings with the first generation , implying @xmath203 .
this might seems like an ad - hoc assumption .
however , it just reflects the _ constant _ need for a theory of flavor , also required to understand the yukawa structure of the sm  .
finally , one can envisage several phenomenological directions in which this model can be further explored . by allowing for non - zero couplings to first generation quarks and leptons ,
the impact of the model gets extended to many additional observables of interest .
in particular , the lhc would be able to probe the parameter space of the model via the production of the @xmath2 boson .
furthermore , since the scalar dm carries , at early times , a conserved @xmath57 charge it will develop an asymmetry , provided chemical equilibrium with the thermal bath is guaranteed .
in such a case ( asymmetric dm scenario @xcite ) , a link between the baryon and dm asymmetries is something worth exploring .
an update of the global fit @xcite was recently presented in @xcite .
while this would change slightly the numerical output of our analysis , our conclusions would remain unchanged .
we are grateful to javier virto and jorge martin camalich for enlightening discussions .
das wants to thank diego restrepo and carlos yaguna for very useful conversations on dark matter .
das is supported by a `` charg de recherches '' contract funded by the belgian fnrs agency .
av acknowledges partial support from the expl / fis - nuc/0460/2013 project financed by the portuguese fct .
r.  aaij _ et al . _
[ lhcb collaboration ] , phys .
* 111 * ( 2013 ) 19 , 191801 [ arxiv:1308.1707 [ hep - ex ] ] .
j.  matias , f.  mescia , m.  ramon and j.  virto , jhep * 1204 * ( 2012 ) 104 [ arxiv:1202.4266 [ hep - ph ] ] .
s.  descotes - genon , j.  matias , m.  ramon and j.  virto , jhep * 1301 * ( 2013 ) 048 [ arxiv:1207.2753 [ hep - ph ] ] .
s.  descotes - genon , t.  hurth , j.  matias and j.  virto , jhep * 1305 * ( 2013 ) 137 [ arxiv:1303.5794 [ hep - ph ] ] .
r.  aaij _ et al . _
[ lhcb collaboration ] , jhep * 1406 * ( 2014 ) 133 [ arxiv:1403.8044 [ hep - ex ] ] .
r.  aaij _ et al . _
[ lhcb collaboration ] , jhep * 1307 * ( 2013 ) 084 [ arxiv:1305.2168 [ hep - ex ] ] .
s.  descotes - genon , j.  matias and j.  virto , phys .
d * 88 * ( 2013 ) 7 , 074002 [ arxiv:1307.5683 [ hep - ph ] ] .
w.  altmannshofer and d.  m.  straub , eur .
j.  c * 73 * ( 2013 ) 12 , 2646 [ arxiv:1308.1501 [ hep - ph ] ] .
f.  beaujean , c.  bobeth and d.  van dyk , eur .
j.  c * 74 * ( 2014 ) 6 , 2897 [ erratum - ibid .
c * 74 * ( 2014 ) 12 , 3179 ] [ arxiv:1310.2478 [ hep - ph ] ] .
t.  hurth and f.  mahmoudi , jhep * 1404 * ( 2014 ) 097 [ arxiv:1312.5267 [ hep - ph ] ] .
r.  aaij _ et al . _
[ lhcb collaboration ] , phys .
* 113 * ( 2014 ) 15 , 151601 [ arxiv:1406.6482 [ hep - ex ] ] . g.  hiller and f.  kruger ,
d * 69 * ( 2004 ) 074020 [ hep - ph/0310219 ] .
r.  alonso , b.  grinstein and j.  martin camalich , phys .
* 113 * ( 2014 ) 24 , 241802 [ arxiv:1407.7044 [ hep - ph ] ] . g.  hiller and m.  schmaltz , phys .
d * 90 * ( 2014 ) 5 , 054014 [ arxiv:1408.1627 [ hep - ph ] ] .
d.  ghosh , m.  nardecchia and s.  a.  renner , jhep * 1412 * ( 2014 ) 131 [ arxiv:1408.4097 [ hep - ph ] ] .
t.  hurth , f.  mahmoudi and s.  neshatpour , jhep * 1412 * ( 2014 ) 053 [ arxiv:1410.4545 [ hep - ph ] ] . w.  altmannshofer and d.  m.  straub , arxiv:1411.3161 [ hep - ph ] .
c.  langenbruch , talk presented in the ew session of the _ 50th rencontres de moriond _ , march 20th , 2015 .
s.  jger and j.  martin camalich , jhep * 1305 * ( 2013 ) 043 [ arxiv:1212.2263 [ hep - ph ] ] .
s.  jger and j.  martin camalich , arxiv:1412.3183 [ hep - ph ] .
r.  gauld , f.  goertz and u.  haisch , phys .  rev .
d * 89 * ( 2014 ) 1 , 015005 [ arxiv:1308.1959 [ hep - ph ] ] .
a.  j.  buras and j.  girrbach , jhep * 1312 * ( 2013 ) 009 [ arxiv:1309.2466 [ hep - ph ] ] . a.  j.  buras , j.  girrbach - noe , c.  niehoff and d.  m.  straub , jhep * 1502 * ( 2015 ) 184 [ arxiv:1409.4557 [ hep - ph ] ] .
r.  gauld , f.  goertz and u.  haisch , jhep * 1401 * ( 2014 ) 069 [ arxiv:1310.1082 [ hep - ph ] ] .
a.  j.  buras , f.  de fazio and j.  girrbach , jhep * 1402 * ( 2014 ) 112 [ arxiv:1311.6729 [ hep - ph ] , arxiv:1311.6729 ] .
a.  datta , m.  duraisamy and d.  ghosh , phys .
d * 89 * ( 2014 ) 7 , 071501 [ arxiv:1310.1937 [ hep - ph ] ] .
s.  biswas , d.  chowdhury , s.  han and s.  j.  lee , jhep * 1502 * ( 2015 ) 142 [ arxiv:1409.0882 [ hep - ph ] ] .
b.  gripaios , m.  nardecchia and s.  a.  renner , arxiv:1412.1791 [ hep - ph ] .
i.  d.  m.  varzielas and g.  hiller , arxiv:1503.01084 [ hep - ph ] . c.  niehoff , p.  stangl and d.  m.  straub , arxiv:1503.03865 [ hep - ph ] . w.  altmannshofer , s.  gori , m.  pospelov and i.  yavin , phys .  rev .
d * 89 * ( 2014 ) 9 , 095033 [ arxiv:1403.1269 [ hep - ph ] ] .
a.  crivellin , g.  dambrosio and j.  heeck , arxiv:1501.00993 [ hep - ph ] .
a.  crivellin , g.  dambrosio and j.  heeck , arxiv:1503.03477 [ hep - ph ] .
p.  a.  r.  ade _ et al .
_ [ planck collaboration ] , astron .
astrophys .
* 571 * ( 2014 ) a16 [ arxiv:1303.5076 [ astro-ph.co ] ] .
a.  j.  buras , f.  de fazio and j.  girrbach , jhep * 1302 * ( 2013 ) 116 [ arxiv:1211.1896 [ hep - ph ] ] . l.  m.  krauss and f.  wilczek , phys .
* 62 * ( 1989 ) 1221 . b.  petersen , m.  ratz and r.  schieren , jhep * 0908 * ( 2009 ) 111 [ arxiv:0907.4049 [ hep - ph ] ] .
d.  a.  sierra , m.  dhen , c.  s.  fong and a.  vicente , arxiv:1412.5600 [ hep - ph ] .
b.  holdom , phys .
b * 166 * ( 1986 ) 196 .
r.  m.  fonseca , m.  malinsk and f.  staub , phys .
b * 726 * ( 2013 ) 882 [ arxiv:1308.1674 [ hep - ph ] ] .
x.  chu , y.  mambrini , j.  quevillon and b.  zaldivar , jcap * 1401 * ( 2014 ) 01 , 034 [ arxiv:1306.4677 [ hep - ph ] ] .
a.  alves , s.  profumo and f.  s.  queiroz , jhep * 1404 * ( 2014 ) 063 [ arxiv:1312.5281 [ hep - ph ] ] .
a.  alves , a.  berlin , s.  profumo and f.  s.  queiroz , arxiv:1501.03490 [ hep - ph ] .
g.  belanger , f.  boudjema , a.  pukhov and a.  semenov , comput .
phys .  commun .
* 185 * , 960 ( 2014 ) [ arxiv:1305.0237 [ hep - ph ] ] .
e.  accomando , a.  belyaev , j.  fiaschi , k.  mimasu , s.  moretti and c.  shepherd - themistocleous , arxiv:1503.02672 [ hep - ph ] .
d.  barducci , a.  belyaev , m.  buchkremer , g.  cacciapaglia , a.  deandrea , s.  de curtis , j.  marrouche and s.  moretti _ et al .
_ , jhep * 1412 * ( 2014 ) 080 [ arxiv:1405.0737 [ hep - ph ] ] .
collaboration , _ a search for anomalous production of events with three or more leptons using 19.5/fb of sqrt(s)=8 tev lhc data _ , tech .
cms - pas - sus-13 - 002 , cern , geneva , 2013 .
a.  falkowski , d.  m.  straub and a.  vicente , jhep * 1405 * ( 2014 ) 092 [ arxiv:1312.5329 [ hep - ph ] ] .
r.  dermek , j.  p.  hall , e.  lunghi and s.  shin , jhep * 1412 * ( 2014 ) 013 [ arxiv:1408.3123 [ hep - ph ] ] .
collaboration , _
search for new phenomena in events with three charged leptons at a center - of - mass energy of 8 tev with the atlas detector _ , atlas - conf-2013 - 070 , cern , geneva , 2013 .
c.  bobeth , m.  gorbahn , t.  hermann , m.  misiak , e.  stamou and m.  steinhauser , phys .
* 112 * ( 2014 ) 101801 [ arxiv:1311.0903 [ hep - ph ] ] . and * lhcb * collaborations , eps - hep 2013 european physical society conference on high energy physics , stockholm , sweden , 2013 , conference report no .
cms - pas - bph-13 - 007 , lhcb - conf-2013- 012 , http://cds.cern.ch/record/1564324 .
k.  de bruyn , r.  fleischer , r.  knegjens , p.  koppenburg , m.  merk , a.  pellegrino and n.  tuning , phys .
* 109 * ( 2012 ) 041801 [ arxiv:1204.1737 [ hep - ph ] ] .
a.  j.  buras , r.  fleischer , j.  girrbach and r.  knegjens , jhep * 1307 * ( 2013 ) 77 [ arxiv:1303.3820 [ hep - ph ] ] .
g.  yeghiyan , phys .
d * 88 * ( 2013 ) 3 , 034020 [ arxiv:1305.0852 [ hep - ph ] ] .
d.  becirevic , n.  kosnik , f.  mescia and e.  schneider , phys .
d * 86 * ( 2012 ) 034034 [ arxiv:1205.5811 [ hep - ph ] ] .
h.  k.  dreiner , k.  nickel and f.  staub , phys .
d * 88 * ( 2013 ) 11 , 115001 [ arxiv:1309.1735 [ hep - ph ] ] .
w.  altmannshofer , s.  gori , m.  pospelov and i.  yavin , phys .
* 113 * ( 2014 ) 091801 [ arxiv:1406.2332 [ hep - ph ] ] .
s.  r.  mishra _ et al .
_ [ ccfr collaboration ] , phys .  rev .
* 66 * ( 1991 ) 3117 .
s.  l.  glashow , d.  guadagnoli and k.  lane , phys .
lett .   * 114 * ( 2015 ) 9 , 091801 [ arxiv:1411.0565 [ hep - ph ] ] .
b.  bhattacharya , a.  datta , d.  london and s.  shivashankara , phys .
b * 742 * ( 2015 ) 370 [ arxiv:1412.7164 [ hep - ph ] ] .
s.  sahoo and r.  mohanta , arxiv:1501.05193 [ hep - ph ] .
s.  m.  boucenna , j.  w.  f.  valle and a.  vicente , arxiv:1503.07099 [ hep - ph ] .
v.  khachatryan _ et al . _
[ cms collaboration ] , arxiv:1502.07400 [ hep - ex ] .
s.  davidson and g.  j.  grenier , phys .
d * 81 * ( 2010 ) 095016 [ arxiv:1001.0434 [ hep - ph ] ] .
j.  kopp and m.  nardecchia , jhep * 1410 * ( 2014 ) 156 [ arxiv:1406.5303 [ hep - ph ] ] .
d.  aristizabal sierra and a.  vicente , phys .
d * 90 * ( 2014 ) 11 , 115004 [ arxiv:1409.7690 [ hep - ph ] ] .
j.  heeck , m.  holthausen , w.  rodejohann and y.  shimizu , arxiv:1412.3671 [ hep - ph ] .
i.  dorsner , s.  fajfer , a.  greljo , j.  f.  kamenik , n.  kosnik and i.  nisandzic , arxiv:1502.07784 [ hep - ph ] .
a.  abada , d.  das , a.  m.  teixeira , a.  vicente and c.  weiland , jhep * 1302 * ( 2013 ) 048 [ arxiv:1211.3052 [ hep - ph ] ] .
a.  abada , a.  m.  teixeira , a.  vicente and c.  weiland , jhep * 1402 * ( 2014 ) 091 [ arxiv:1311.2830 [ hep - ph ] ] .
t.  asaka , s.  eijima and k.  takeda , phys .
b * 742 * ( 2015 ) 303 [ arxiv:1410.0432 [ hep - ph ] ] .
[ atlas collaboration ] , atlas - conf-2013 - 010 , atlas - com - conf-2013 - 003 .
f.  staub , arxiv:0806.0538 [ hep - ph ] .
f.  staub , comput .
commun .   * 181 * ( 2010 ) 1077 [ arxiv:0909.2863 [ hep - ph ] ] .
f.  staub , comput .
commun .
* 182 * ( 2011 ) 808 [ arxiv:1002.0840 [ hep - ph ] ] .
f.  staub , comput .
commun .
* 185 * ( 2014 ) 1773 [ arxiv:1309.7223 [ hep - ph ] ] .
w.  porod , comput .
commun .
* 153 * ( 2003 ) 275 [ hep - ph/0301101 ] .
w.  porod and f.  staub , comput .
commun .
* 183 * ( 2012 ) 2458 [ arxiv:1104.1573 [ hep - ph ] ] .
w.  porod , f.  staub and a.  vicente , eur .
j.  c * 74 * ( 2014 ) 8 , 2992 [ arxiv:1405.1434 [ hep - ph ] ]
. e.  e.  boos , m.  n.  dubinin , v.  a.  ilyin , a.  e.  pukhov and v.  i.  savrin , hep - ph/9503280 .
a.  pukhov , hep - ph/0412191 .
g.  belanger , n.  d.  christensen , a.  pukhov and a.  semenov , comput .
commun .
* 182 * ( 2011 ) 763 [ arxiv:1008.0181 [ hep - ph ] ] .
g.  milutinovi - dumbelovi , i.  boovi - jelisavi , c.  grefe , s.  luki , m.  pandurovi and p.  roloff , arxiv:1412.5791 [ hep - ex ] . s.  m.  boucenna and s.  morisi , front .
* 1 * ( 2014 ) 33 [ arxiv:1310.1904 [ hep - ph ] ] . w.  altmannshofer and d.  m.  straub , arxiv:1503.06199 [ hep - ph ] .


observations indicate that the escape of ly@xmath0 photons is facilitated enormously by the presence of outflowing interstellar gas ( e.g. * ? ? ?
* ) . scattering off these outflows effectively doppler shifts ly@xmath0 photons out of resonance into the red wing of the absorption line profile . in models of this process
the outflow is often represented with a geometrically thin shell of gas ( e.g. * ? ? ?
* ; * ? ? ?
* ) . in spite of its simplicity
, this so - called ` shell - model ' has been very successful at reproducing observed spectra line shapes @xcite .
however , the shell - model has recently been shown to have difficulties , especially in reproducing the strength of ` blue bumps ' in a subset of observed spectra @xcite .
@xcite studied what may be viewed as an extreme example of blueshifted ly@xmath0 emission in which the ly@xmath0 line of spatially extended ly@xmath0 emission around a radio galaxy as a whole is blueshifted relative to the observed 21-cm absorption , with emission extending beyond 1500 km s@xmath2 into the blue wing of the line profile .
blueshifted ly@xmath0 radiation arises naturally when the photons scatter through optically thick inflowing gas @xcite . indeed ,
@xcite show that their data can be reproduced if more than @xmath3 of cold inflowing gas is present .
this large mass of cold neutral gas inside a massive dark matter halo ( @xmath4 ) though , is at odds with expectations from theory , which predicts that the gas should be predominantly accreted in the hot - mode ( e.g. @xcite ; fig .  7 of @xcite ) .
an alternative process which gives rise to blueshifted ly@xmath0 emission is described by neufeld & mckee ( 1988 ; henceforth nm88 ) , who show that scattering across a shock front can give rise to such blueshift .
nm88 presented analytical calculations , and therefore were forced to adopt simplifying assumptions , namely that the shock - crossing acceleration process is modelled by ly@xmath0photons bouncing between two partially transparent , frequency - preserving mirrors . in reality , however , the frequency of a ly@xmath0photon is not preserved at each ` reflection ' . instead , the frequency of the photon diffuses through frequency - space as the photon is scattered by the material on either side of the shock front , in turn changing the optical depth of the gas to the photon .
in addition , they assume that photons cross the shock front isotropically .
finally , the analysis presented by nm88 was restricted to a simplified geometry of two adjacent semi - infinite slabs .    in spite of an increasing prevalence of ly@xmath0 radiative transfer monte - carlo codes @xcite , they have so far not been used to investigate this fermi - like acceleration mechanism ( while acknowledging the inaccuracy of doing so , we henceforth refer to this as ` fermi acceleration ' for brevity ) .
the goals of this paper are two - fold : ( _ i _ ) numerically study fermi acceleration of ly@xmath0 photons across a shock front without the simplifying assumptions that were required in the analytical treatment , and ( _ ii _ ) investigate whether this mechanism can help to reproduce observed blue bumps in spectra that have been difficult to explain with conventional shell - models .
the outline of this paper is as follows :
  [ mechanism ] describes the basic acceleration mechanism ,   [ simulations ] details the simulations which we use to investigate the physical mechanism ,   [ simresults ] presents the output of our numerical experiments ,   [ applications ] discusses the applications of our findings , and in   [ conclusions ] we summarise the conclusions .
ly@xmath0photons are resonantly scattered by hi , making the ly@xmath0radiative transfer complex .
the propagation of ly@xmath0photons is affected by both bulk gas motion and , in detail , microscopic motion of individual hydrogen atoms .
the result is that ly@xmath0photons undergo a random - walk like motion in both physical and frequency - space ( see * ? ? ?
* for a detailed review ) .    as ly@xmath0
photons diffuse away from the resonance frequency , their mean free path - and hence their escape probability - increases . as a result ,
the spectrum of ly@xmath0 photons emerging from static optically thick media consists of two peaks , which are distributed symmetrically around the resonance frequency ( as frequency diffusion can occur to lower and higher energies with equal probability k ) the energy deposited in the recoiling scattering atom becomes important , and the red peak is enhanced . ] ) .
when the scattering medium is contracting [ expanding ] however , the frequency diffusion preferentially occurs towards higher [ lower ] frequencies . as a result , the spectrum of ly@xmath0 photons emerging from contracting optically thick media is blueshifted .
an alternative way to see this is that the converging flow of contracting gas is doing work on the ly@xmath0 photons as they scatter outwards , which increases their mean energy @xcite .
the fermi acceleration mechanism described in nm88 invokes similar converging flows of gas .
the mechanism is illustrated with a simplified geometry in figure [ setupsfig ] , which shows two semi - infinite adjacent slabs of neutral gas .
the slab on the left is moving to the right with velocity @xmath5 into a stationary slab . if @xmath6 exceeds the sound speed in the left slab , then the two slabs are separated by a shock front . when a ly@xmath0source is in the vicinity of this shock front
we expect some ly@xmath0photons to diffuse through space and cross the front .
when a photon traverses the shock front , the doppler boost will impart a blueshift in the local gas frame .
partially coherent scattering off atoms with thermal velocity @xmath7 mostly preserves the blueshift of this photon : for a photon with frequency @xmath8 each scattering event pushes back the photon to the line centre by an average amount @xmath9 @xcite . when the photon scatters back across the shock front before this ` restoring force ' has returned the photon to line center , the blueshift of the photon increases with each shell crossing .
this is what nm88 referred to as ` fermi acceleration ' .    in order for fermi acceleration of a ly@xmath0photon to actually occur we require it to scatter after crossing the shock front .
the optical depth to a ly@xmath0 photon that is blueshifted to @xmath10 through a slab of gas with column density @xmath11 equals @xmath12 .
the condition @xmath13 thus translates to @xmath14 @xmath15 .
inversely , a photon can be fermi accelerated to a maximum blueshift of @xmath16 km s@xmath2 , _ in the frame of the scattering medium .
_ after scattering , the ly@xmath0 photon can get an additional doppler shift @xmath17 depending on whether the scattering medium is moving or not .
in this section we describe the setup and results of the numerical simulations which we use to study the fermi acceleration mechanism briefly summarised in   [ mechanism ] .       and [ simulations ] ) .
the ly@xmath0source is denoted by the cross in the left slab .
the left slab has hi column density @xmath18 , and the right slab has hi column density @xmath19 .
the gas in the left slab moves with velocity @xmath20 , photons cross the shock front with velocity @xmath21 , and the angle between @xmath20 and @xmath21 is labelled @xmath22 .
@xmath23 is then the projection of @xmath21 onto the unit direction vector of the shock front , @xmath24 .
, width=302 ]    we use modified versions of the slaf 3d monte - carlo ly@xmath0radiative transfer code ( chung et al . 2015 , in prep ) , as well as the code described in @xcite .
the vanilla version of slaf is a grid - based code , able to handle arbitrary gas and source distributions within a finite volume of 3d space . for this work
we have modified it , removing the dependence on an underlying grid - structure , so that the code can handle the specific setup described in nm88 .
we choose to focus on the geometry discussed in nm88 to facilitate a straightforward comparison .
explicitly , there are two distinct regions of gas , represented by slabs of infinite width ( y ) and height ( z ) . the slab on the left side represents outflowing gas , and has a bulk velocity of @xmath25 km s@xmath2 from left to right .
the hi column density across the outflowing slab is @xmath26@xmath15 .
the slab on the right side represents static gas with an hi column density @xmath27 @xmath15 @xmath28 .
a shock front exists at the interface of the two slabs .
ly@xmath0photons are emitted at line centre in the rest frame of the outflowing slab on the left , which represents ly@xmath0 emission powered by star formation triggered by the passing shockwave ( see nm88 ) .
we ignore dust throughout the analysis in
[ simresults ] , which simplifies the interpretation of our results and does not affect our main results .
we discuss the impact of this assumption separately when we discuss applications of our results in   [ bluewingsrg ] and   [ bluebumpspectra ] .
figure [ setupsfig ] shows the described experimental setup .
figure 2 was produced from a run tracking @xmath29 crossing events , while figures 3 and 4 were produced from runs of @xmath29 photons . in all cases
the results were checked for convergence , and found to already be converged with much lower photon counts .
table [ crossingproportions ] shows the fraction of photons escaping the simulation volume , @xmath30 , as a function of the total number of shock crossings , @xmath31 .
the second and third columns indicate the @xmath32 found in our simulations and the analytic estimate , respectively .
the following paragraphs explain how we obtained these analytic estimates .
photons emitted in the centre of the outflowing slab are equally likely to leave this slab on the left or right hand side .
we therefore expect @xmath33 of all photons not to cross the shock front , which is in agreement with the simulated fraction .
photons that enter the static slab for the first time are doppler boosted by an average amount @xmath34 ( see   3.2.2 ) .
the probability of these photons being transmitted all the way through the slab is given by @xmath35 , where @xmath36 denotes the optical depth of the entire static slab to photons that enter at frequency @xmath37 ( in the frame of the slab , see @xcite ; note that this transmission probability only applies when the optical depth of the slab to incoming photons is @xmath38 ) .
substituting numbers yields @xmath39 , where @xmath40 denotes the hi column density of the slab in units of @xmath41 @xmath15 , and @xmath42 denotes @xmath10 in units of 100 km s@xmath2 ( notation adopted from nm88 ) .
we then get @xmath43 , which compares reasonably to the fraction @xmath44 we found in our simulation ( exact agreement is not expected as photons cross the slab over a range of frequencies ) .
photons that cross the front twice are back in the low column density outflowing slab , and appear doppler shifted by an average amount @xmath45 in the slab frame ( see   3.2.2 ) .
the slab will appear optically thin to these photons , and we estimate the transmission probability from @xmath46 $ ] . here , @xmath47 denotes the optical depth through the outflowing shell for the incoming photons , where @xmath48 is the ly@xmath0 absorption cross - section at @xmath10 .
the factor of @xmath49 in the exponent accounts for the fact that the photons enter the slab under an angle @xmath22 with probability @xmath50 ( see   3.2.2 ) .
the transmission probability is thus the sum of the probability ( _ i _ ) that photons are transmitted directly through the slab ( the @xmath51 term ) , and ( _ ii _ ) that photons that do scatter eventually escape the slab without crossing the shock again . for the latter
, we assume that photons are equally likely to escape on the left and right hand side after scattering , which seems reasonable given that the slab appears optically thin on average . substituting numbers we obtain @xmath52 , which agrees with the simulation result .    with @xmath35 and @xmath53 , we estimate @xmath54 , 48.5% larger than what we obtain from the simulation .
we attribute this discrepancy to the approximation @xmath35 breaking down at larger shifts from line centre .
the previous analysis allows us to understand quantitatively the simulation results .
more photons escape after undergoing an even number of shock crossings than odd .
this is simply a reflection of the fact that in our experimental setup we have specified that @xmath55 .    as we pointed out in
[ introduction ] , it is theoretically possible that a ly@xmath0photon that has been fermi accelerated into the line wing can subsequently diffuse back into the core of the line prior to escaping or crossing the shock front again .
when this happens , any memory of previous fermi acceleration is erased .
however , our numerical simulations indicate that in practice this almost never occurs , and therefore we can identify odd [ even ] numbered contiguous shock crossings with photons which exit on the right [ left ] side of the simulation .
.fraction of photons undergoing _
n _ successive shock crossings before exiting the simulation ( middle column ) .
the right column refers to the same fraction calculated with an analytic method . [ cols="^,^,^",options="header " , ]      we denote the angle at which a photon crosses the shock front with @xmath22 ( see fig  [ setupsfig ] ) . for an isotropic distribution of shock crossing directions we have @xmath56 , where @xmath57 .
the average doppler shift experienced by a photon as it crosses the shock front is @xmath58 .
after @xmath59 crossings the photon experiences an average blueshift of @xmath60 ( nm88 ) .    with our code
we are able to track individual photons as they cross the shock front , which allows us to directly measure @xmath61 .
figure [ crossingdistr ] shows the number count , @xmath62 , for photons that cross the @xmath63 time for various @xmath59 . @xmath61
only differs from @xmath62 by a normalisation factor . here , @xmath59 refers to the current ( rather than final ) contiguous crossing count of a particular photon as it crosses the shock front .
a photon which finally exits the simulation after @xmath64 crossings is represented @xmath64 times ( once for each crossing ) .
figure [ crossingdistr ] shows clearly that the distribution of photon crossing projections is not isotropic @xmath65 [ @xmath66 for photons that cross the shock front from left - to - right [ right - to - left ]
. in figure [ crossingdistr ] only photons on an odd [ even ] shock crossing contribute to @xmath62 for positive [ negative ] values of @xmath67 . ] .
instead , we find that @xmath68 . this distribution has been found in previous analyses ( e.g. * ? ? ? * ; * ? ? ?
* ) , and its origin is discussed in appendix  [ app : pmu ] . for this distribution
we expect @xmath69 ( see also fig .  [ projvsscattering ] ) , and hence that a photon experiences a doppler boost @xmath70 after @xmath59 shock crossings , then the total hydrogen column density to the edge of the slab is @xmath71 , where @xmath72 [ @xmath73 for the slab on the left [ right ] in figure  [ setupsfig ] .
the angle - averaged column to the edge of the slab is given by @xmath74 , as used in
3.2.1 . ] .    .
for each value of @xmath59 the plotted distribution includes all photons that make an @xmath63 crossing regardless of what their final crossing count , @xmath64 , is.,width=302 ]    figure [ velshiftvsscatterings ] shows the velocity shift @xmath10 of exiting photons as a function of the total number of shell crossings @xmath64 , in the frame of the gas into which the photon is crossing ( i.e. for odd [ even ] numbered shock crossings the velocity shift is given in the static [ outflowing ] frame ) .
results from our monte - carlo simulations are represented by the data points .
the red dashed line shows the analytic result under the assumption that photons cross the shock front isotropically ( as in * ? ? ?
* ) , while the green solid line shows the analytic result for photons crossing the shock front according to @xmath68 .
for @xmath75 the simulation results follow the @xmath76-relation ( where here @xmath64 is a good proxy for @xmath59 ) , after which it approaches the isotropic shock - crossing case .
the reason for this transition is that as the photons get blueshifted further into the wing of the line profile their mean free path increases .
as we discuss in detail in appendix  [ app : pmu ] , we expect a transition to isotropic shock crossing when the mean free path becomes comparable to the thickness of the slab .
+     shock crossings as found in our simulations ( points ) , as measured in the frame of the gas .
the error bars show the standard deviation within the distribution .
the dashed red line shows the @xmath77-relation which applies to isotropic shock crossings , whereas the green line refers to the @xmath78-relation which applies to anisotropic shock - crossings of the form @xmath68 ( see text).,width=321 ]    finally , figure [ spectrumsinglesource ] shows the spectrum outcoming from the slab , as measured in the lab frame .
the result can be easily understood in terms of our previous analysis . more specifically , the total spectrum ( black histogram ) , can be explained as follows .
* the peak at @xmath79 km s@xmath2 is composed primarily of the @xmath80 of all photons that did not cross the shock front .
these photons diffused outward of the left slab , and would have emerged with a characteristic double peaked emission line profile centered around the resonance frequency ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) . however , when we doppler boost these photons back into the lab frame this double peak is diluted by the fact that the photons escaping the slab do so at various angles .
the doppler boost is dependent on the exit angle , with an average doppler boost of @xmath81 km s@xmath2 .
the grey crosshatched histogram shows the spectrum of these photons . *
the small fraction of photons that escape after a single shock crossing , @xmath82 , is indicated by the purple hatch - filled histogram and escapes with a mean blueshift of @xmath83 km s@xmath2 ( also see fig .
[ velshiftvsscatterings ] , these photons escape from the static shell , and no additional doppler boost into the lab frame is required ) . *
the photons that escape after two shock crossings ( green dotted line histogram ) account for @xmath84 of the total .
figure  [ velshiftvsscatterings ] shows that these photons have accumulated a blueshift of @xmath85 km s@xmath2 in the frame of the outflowing gas .
a doppler boost back into the lab - frame transforms these photons back to a blueshift of @xmath83 km s@xmath2 .
a similar reasoning applies to photons that escape after a larger number of crossings , and shows that it is easy to understand the shape of the total spectrum : the redshifted peak ( grey crosshatched histogram ) at @xmath86 is composed of photons that did not cross the shock front .
the first peak blueward ( @xmath87 ) of the ly@xmath0resonance is at @xmath88 and consists of photons that crossed the shock front @xmath89 times .
similarly , the second blueward peak at @xmath90 consists of photons that crossed the shock 3 - 4 times , etc .
the prevalence of these peaks depends on a number of factors , including @xmath20 , @xmath18 , @xmath91 , and also the distribution of ly@xmath0 sources relative to the hi gas .
finally we note that our results here are in excellent agreement with the analytic estimate of the maximum blueshift detailed in   [ mechanism ] .
setting @xmath92@xmath15 in the analytic formula gives @xmath93 km s@xmath2 , which lies in the blue tail of figure [ spectrumsinglesource ] .
is the number of shock crossings that the photon undergoes prior to exiting the simulation ( see   [ crossingstats]).,width=302 ]
we now explore , by means of our monte - carlo codes and assuming some simple models , whether the fermi acceleration mechanism can provide a viable explanation for some observed systems .
the analysis by nm88 was motivated by observations of 3c326.1 , a radio source with a ly@xmath0 spectral line profile with emission extending far into the blue wing of the line profile ( @xcite , @xcite ) .
nm88 modelled 3c326.1 with a shock front that was propagating into a collection of dense clumps .
the shock passage triggers star formation and ly@xmath0 emission inside the post - shock gas .
the subsequent radiative processes in each clump proceed just as in the slab models we discussed previously .
the details of the emerging spectrum depend on the distribution and covering factor of the clumps , but for simplicity and to facilitate a direct comparison to nm88 we follow their prescription and assume that the emerging spectrum is a superposition of the spectra of individual clumps .
nm88 further show that dust does not affect the emerging spectra , provided that the gas - to - dust ratio @xmath94 , where @xmath95 denotes the dust - to - gas ratio relative to the local interstellar value .
they argue that this limit is acceptable because the dust grains would be destroyed by the passage of the shock front .
we finally assume all clumps - and therefore the radiative transfer processes inside them - to be identical , then we can simply adopt the spectrum from figure  [ spectrumsinglesource ] . in reality , from our findings in   3.2.2 , we expect that a distribution of clump properties would give rise to a superposition of spectra with peaks in different locations , and therefore that the emission in the blue wing of the line profile would have less prominent features spectrum depends on the placement of ly@xmath0 sources . in particular , changing only the source position leaves the peak positions unchanged , but modulates their relative amplitudes .
we have explicitly verified , however , that our spectrum barely changed if we assumed a uniform distribution of ly@xmath0 sources throughout the outflowing slab . ] .
putting these reservations aside , we observe that qualitatively our simulated spectrum shows similarities to that predicted by nm88 .
we should note though that in figure  [ spectrumsinglesource ] we omit the direct sources at @xmath96 , which nm88 put in by hand as a delta function at @xmath96 .
we note that the two largest peaks immediately redward and blueward of line centre can not be attributed to photons exiting without crossing the shock front as in nm88 ( i.e. the standard double - peak profile produced by a static , optically thick medium ) .
instead , our analysis shows clearly that these peaks are primarily caused by photons which have undergone 0 and \{1 , 2 } shock crossings respectively .
the smaller , secondary blueward peak is primarily composed of photons which have undergone 3 or 4 shock crossings .
the primary differences between the spectra predicted by us and by nm88 are in the asymmetry of the two main peaks , and the ratio of the peaks .
as we have previously mentioned , the source distribution and physical parameters of the clumps affect the outcoming spectrum , and it is likely that tuning of either or both of these will alleviate this discrepancy .
the two blue peaks in our model appear to be shifted redward compared to the model predictions from nm88 .
however , a comparison of the observed spectrum in nm88 and our model shows that the positions of our blue peaks are also a good fit to the data .
in fact , if the amplitude of the peaks is ignored ( this is dependent on the exact source distribution ) , the position of the peaks in our model is arguably a better fit to the data than the nm88 model .
we now consider the observed ly@xmath0 spectra of @xmath97 radio galaxy b2 0902 + 34 @xcite .
there is broad agreement in the shape of the observed and predicted spectra , with the broad wing extending well beyond @xmath98 km s@xmath2 , and some spectra displaying prominent peaks blueshifted by @xmath98 km s@xmath2 .
the spectra from b2 0902 + 34 vary widely depending on where exactly on the galaxy the fiber is placed .
it is however true to say that for the majority of fibers the most prominent peak has been observed to lie either at the systemic velocity , which was determined from the observed 21-cm absorption signature , or slightly blueward of it .
this would appear to be in tension with our model here which , following nm88 , predicts the most prominent peak to appear at @xmath96 ( not shown in figure [ spectrumsinglesource ] ) as a result of direct ` blister ' sources on the edge of the clump .
recalling that the amplitude ratios of the peaks is determined by the source distribution , we speculate that this tension could be relieved if the most prominent peak in the b2 0902 + 34 observations is instead identified as the first blueshifted peak in our model , and the absence of a pronounced peak at systemic velocity attributed to a lower prevalence of ` blister ' sources in this system .      as we mentioned in
[ introduction ] , observed ly@xmath0 spectral line profiles can often be reproduced surprisingly well with shell - models . in these models ,
a central ly@xmath0 source is surrounded by a geometrically thin shell of gas .
the shell - models contain two parameters describing the ly@xmath0 source : ( _ i _ ) the assumed fwhm of the ly@xmath0 line prior to scattering , and ( _ ii _ ) the strength of the ly@xmath0 emission line , which is quantified by the equivalent width ( ew ) .
the shell itself is described by four additional parameters : ( _ i _ ) the hi column density of the shell @xmath99 , ( _ ii _ ) its outflow velocity @xmath100 , ( _ iii _ ) its dust content @xmath101 , and ( _ iv _ ) its velocity dispersion @xmath102 .
some recent analyses have pointed out that shell - models have difficulties explaining blue bumps in observed spectra ( e.g. * ? ? ?
* ; * ? ? ?
nice examples can be seen e.g. in figure  7 of @xcite . more recently @xcite pointed out that the shell - models , when applied to blue - bump objects , consistently require the intrinsic ly@xmath0 line to have a fwhm which is too large ( an excessively high fwhm was only required for blue - bump objects ) .    because fermi acceleration naturally gives rise to blueshifted emission
, we investigate a minor modification of the shell - model in which we embed the outflowing shell within a static gas cloud .
this modification can be interpreted as the situation in which the outflow is still propagating into the static interstellar medium ( possibly prior to breaking out of it ) .
this configuration now includes a shock front as in our previous analysis .
we specifically study a shell - model with parameters based loosely on those inferred by @xcite : @xmath103 @xmath15 , @xmath104 km s@xmath2 , and fwhm=200 km s@xmath2 .
we further assume @xmath105 , ew@xmath106 ( i.e. pure ly@xmath0emission ) , and @xmath107 km s@xmath2 ( corresponding to gas at @xmath108k ) , but note that these assumptions do not affect our results at all .
these results show clearly that dust barely affects the blue bumps .
this result is not surprising : our modification only adds a small amount of additional hydrogen , and this additional gas triggers fermi - acceleration into the ( blue ) wings of the ly@xmath0 line profile , where ly@xmath0photons escape more easily . ] . finally , the key new model ingredient is the addition of a static shell of gas adjacent to the outflowing shell .
this shell is characterised primarily by its column density @xmath109 .
we assume that it has no dust and that it has the same temperature as the outflowing shell .
the results of this analysis for @xmath110 ( i.e. the original shell - model ; solid black histogram ) , @xmath111 @xmath15 ( red dashed line ) , and @xmath112 @xmath15 ( blue dotted line ) are shown in figure  [ fig : shellresults ] ,     spectra emerging from a ly@xmath0 source surrounded by a shell of hi gas outflowing at 200  km  s@xmath2 ( the shell - model , see text ) embedded within a static neutral shell with a column density @xmath109 .
the solid black histogram shows the spectrum for @xmath110 ( i.e. the standard shell - model ) , while the red dashed and blue dotted line shows the spectrum for @xmath111 @xmath15 and @xmath113 @xmath15 , respectively .
the plot demonstrates that adding even a small amount of hydrogen ( 1/10@xmath114 of that in the shell ) triggers the on - set of fermi acceleration , which gives rise to a blue bump.,width=302 ]    where it is clear that a small additional column of hydrogen ( @xmath115 ) dramatically affects the spectrum blueward of the systemic velocity .
this large change of the spectrum can be easily understood , as the outflowing shell with @xmath103 @xmath15 directly transmits a significant fraction of ly@xmath0 photons .
however , the surrounding static shell is optically thick to photons emitted near line centre , because these photons still appear close to the centre of the line in the frame of this gas .
the static shell therefore effectively reflects back photons into the outflowing shell we show that this mechanism produces blue bumps for column densities as low as @xmath116 @xmath15 .
this is because a static shell of gas with @xmath116 @xmath15 remains optically thick to ly@xmath0 photons emitted close to line centre , and can therefore reflect back photons into the expanding shell . ] .
the reflected photons appear blueshifted by @xmath117 in the frame of the outflowing shell , where their newly acquired large blueshift makes them escape efficiently .
for example , photons that are scattered by @xmath118 after being reflected back into the outflowing shell escape with a blueshift of @xmath119 , which is indeed where the new peaks in the spectrum emerge .
this process is depicted schematically in figure  [ fig : shell ] . the analytic estimate for the maximum blueshift in   [ mechanism ] provides a decent estimate of the maximum blueshift @xmath120 km s@xmath2 ( where the @xmath121 km s@xmath2 was the maximum shift for a column density of @xmath122 @xmath15 in the frame of the scattering medium ) .
spectra in the shell - model .
the outflowing shell is surrounded by a static ism .
, width=321 ]
we have presented a detailed analysis of fermi acceleration of ly@xmath0 photons across a shock front , a process that was first studied analytically by nm88 .
our monte - carlo radiative transfer simulations of the slab model ( figure [ setupsfig ] ) confirm the basic results in nm88 .    in particular , we find that each time a photon crosses a shock front that is propagating at @xmath20 , it experiences a doppler boost @xmath123 in the gas frame .
the precise numerical coefficient depends on the opacity of the gas to the ly@xmath0 photons .
we also find the blueshift of the ly@xmath0 photons to increase in proportion with the number of shock crossings @xmath59 , i.e. @xmath124 .
we discussed how our results can help to explain extended blue wings observed in spectra of radio galaxy 3c326.1 ( which was proposed previously by nm88 ) , but we note that our line as a whole appears redshifted by @xmath125 km s@xmath2 .
we reached the same conclusions for the radio galaxy b2 0902 + 34 @xcite .
we nevertheless consider fermi acceleration a plausible alternative to the model proposed by @xcite , which involves the collapse of @xmath126 of neutral gas .
last , we show that fermi acceleration naturally gives rise to blue bumps in ly@xmath0 spectra , which are difficult to reproduce with conventional shell - models ( at least those with reasonable values for the fwhm of the intrinsic ly@xmath0 line , see @xcite ) .
we presented a natural extension to shell - models in which the shell is expanding into static gas that contains a low column density of hi .
these models can give rise to blue bumps in ly@xmath0 spectra without abandoning the simplicity of the conventional shell - model transfer through more complex , simulated gas distributions ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , though results from these calculations have not been compared to observed spectra in detail .
note however , that the sensitivity of the ly@xmath0 transfer scattering process to the distribution and kinematics of neutral gas implies that these predictions depend on sub - grid feedback prescriptions as well as the spatial resolution of the simulation ( which ideally needs to be sub - pc , see dijkstra & kramer 2012 for a discussion ) .
simulations therefore have their own uncertainties , and thus provide a complementary route to addressing ly@xmath0 scattering on interstellar scales . ] . in our model ,
the bump consists of photons that initially streamed through the outflowing shell , but which were reflected back into this shell by the surrounding static interstellar medium .
this suggests that blue bumps in ly@xmath0 spectra are associated with outflows that are still confined to the ism of the galaxy , which may represent an earlier stage in the evolution of the galaxy .
we thank koki kakiichi and michele sasdelli for enlightening discussions , and the referee for critical , helpful comments .
adams , t.  f.  1972 , apj , 174 , 439    j.  j. , hill g.  j. , macqueen p.  j. , 2009 , apj , 694 , 314 ahn , s .- h . , lee , h .- w . , & lee , h.  m.  2001 , apj , 554 , 604 barnes , l.  a. , haehnelt , m.  g. , tescari , e. , & viel , m.  2011 , mnras , 416 , 1723 behrens , c. , & braun , h.  2014 , a&a , 572 , a74    cantalupo , s. , porciani , c. , lilly , s.  j. , miniati , f. , 2005 , apj , 628 , 61 - 75    chonis , t.  s. , blanc , g.  a. , hill , g.  j. , et al .
2013 , apj , 775 , 99 dekel , a. , & birnboim , y.  2006 , mnras , 368 , 2    dijkstra , m. , haiman , z. , & spaans , m.  2006a , apj , 649 , 14    dijkstra , m. , haiman , z. , & spaans , m.  2006b , apj , 649 , 37 dijkstra , m. , & kramer , r.  2012 , mnras , 424 , 1672    dijkstra , m.  2014 , pasa , 31 , e040    s. , 1988 , in kron r.  g. , renzini a. , eds , towards understanding galaxies at large redshift vol .
141 of astrophysics and space science library , lyman - alpha galaxies at large redshifts .
pp 259274 garavito - camargo , j.  n. , forero - romero , j.  e. , & dijkstra , m.  2014 , apj , 795 , 120 harrington , j.  p.  1973
, mnras , 162 , 43    hashimoto , t. , verhamme , a. , ouchi , m. , et al .
2015 , arxiv:1504.03693 kere , d. , katz , n. , weinberg , d.  h. , & dav , r.  2005 , mnras , 363 , 2    kulas , k.  r. , shapley , a.  e. , kollmeier , j.  a. , et al .
2012 , apj , 745 , 33    kunth , d. , mas - hesse , j.  m. , terlevich , e. , et al .
1998 , a&a , 334 , 11 lake , e. , zheng , z. , cen , r. , et al .
2015 , apj , 806 , 46    laursen , p. , razoumov , a.  o. , & sommer - larsen , j.,2009 , apj , 696 , 853    d.  a. , mckee c.  f. , 1988 , apjl , 331 , l87 neufeld , d.  a.  1990 , apj , 350 , 216    osterbrock , d.  e.  1962 , apj , 135 , 195    pierleoni , m. , maselli , a. , ciardi , b. ,  2009 , mnras , 393 , 872    schaerer , d. , hayes , m. , verhamme , a. , & teyssier , r.  2011 , a&a , 531 , a12
m.  a. , mccarthy p.  j. , spinrad h. , djorgovski s. , van breugel w. , liebert j. , dickinson m. , 1987 , in bergeron j. , kunth d. , rocca - volmerange b. , tran thanh van j. , eds , high redshift and primeval galaxies 3c 326.1 - a forming galaxy at a redshift of 1.82 ? .
pp 8592    tasitsiomi , a. , 2006 , apj , 645 , 792    verhamme , a. , schaerer , d. , & maselli , a.  2006 , a&a , 460 , 397    verhamme , a. , schaerer , d. , atek , h. , & tapken , c.  2008 , a&a , 491 , 89 verhamme , a. , dubois , y. , blaizot , j. , et al .
2012 , a&a , 546 , a111    yajima , h. , li , y. , zhu , q. , abel , t.,2012,mnras , 424 , 884    yajima , h. , li , y. , zhu , q. , abel , t.,2015,apj , 801 , 52    zheng , z. , & miralda - escud , j.  2002 , apj , 578 , 33
the goal of this appendix is to clarify why the distribution of angles at which ly@xmath0 photons cross the shock front scales as @xmath68 ( where @xmath57 ) , rather than isotropically .
photons that cross the shock front will on average have traversed an optical depth @xmath127 .
depending on @xmath67 , @xmath127 corresponds to a different ( frequency dependent ) physical depth , @xmath128 , away from the shock front .
it is easy to see that @xmath129 , so that at fixed frequency the volume of gas that the photon is likely to have last scattered in is @xmath130 .
if we further assume that the density of scattering events is homogeneously distributed in the gas , then this leads directly to the relation @xmath131 , where @xmath132 is the intensity of photons .
therefore @xmath133 , which we showed in figure [ crossingdistr ] , and which was found previously by @xcite and @xcite .
figure  [ projvsscattering ] shows the average angle at which photons cross the shock front , @xmath134 , as a function of crossing number @xmath59 . as we mentioned in the paper , odd [ even ] number of shock crossings correspond to crossings from left - to - right [ right - to - left ] .
this figure clearly shows that @xmath135 for even shock - crossings , while it approaches @xmath136 for odd - shock crossings when @xmath137 .    .
the solid lines represent the expected values for isotropic scattering , while the dashed lines show @xmath138 .
red triangular points are for odd-@xmath59 shock crossings , blue square points are for even-@xmath59 shock crossings .
, width=302 ]    an explanation for this is offered in figure [ depthfromshock ] , which shows the mean distance from the shock front ( scattering depth ) of the last scattering before crossing the shock as a function of the shock crossing number @xmath59 .
this distance follows an inverse exponential distribution .
as the photon becomes more blueshifted with increasing @xmath59 , the mean free path of the photon also increases , so that for each successive @xmath59 the photon is coming from deeper within the slab , explaining the increasing trend in figure [ depthfromshock ] .
however , as the mean free path increases , the photons which propagate the farthest actually propagate all the way through the gas slab and exit the simulation . therefore , the distribution of scattering depth is truncated for the next scattering count , resulting in a flattening in the growth of scattering depth as a function of @xmath59 .
this is precisely the behaviour we observe for odd-@xmath59 in figure [ depthfromshock ]
. the scattering depth distribution is truncated , but the photon blueshift is unaffected .
therefore , the gas between the next scattering event and the shock front becomes effectively optically thin to the photon , and we expect photons to cross the shock front isotropically , which is why @xmath139 approaches @xmath140 in figure [ projvsscattering ] .    for even-@xmath59 we do not observe this transition to isotropic crossing .
this is because @xmath141 , so that at a given @xmath59 and blueshift the mean free path of a photon is much lower in the right ( static ) slab and the scattering depth distribution for even-@xmath59 ( i.e. crossings where the previous scattering was in the right , static slab ) has not yet become truncated .
as discussed above it is the truncation of the scattering depth distribution that causes the transition to isotropic shock crossing , and so in its absence for even-@xmath59 the ly@xmath0 photons still cross the shock front following @xmath68 .    .
d is shown in simulation units , i.e. cells ( each slab is 128 cells wide ) .
red triangular points and the red dotted line show the odd - l crossings , while blue square points and the blue dot - dash line show the even - l crossings.,width=302 ]    finally , we note that the highest shock crossing counts that occur in our simulation are rare events , and the resultant poor statistics are responsible for the turnaround in figure [ depthfromshock ] ( recall that we are sampling from an inverse exponential distribution ) .
in dusty media the increase in path length caused by scattering can allow the dust to significantly attenuate the ly@xmath0 flux that escapes .
because photons of different frequencies have different scattering cross - sections the average number of scatterings they undergo will be different .
therefore they will have different changes in path - length , and subsequently undergo different levels of dust attenuation .
we take a direct approach to address this issue and perform numerical experiments where grey dust is added to either the inner , outer , or both shells .
the results are shown in figure [ fig : b1 ] .
it is clear that the differences in the emergent spectrum caused by the different dust prescriptions are very small .
the addition of dust slightly enhances the blue bump .
this occurs because once a photon is fermi accelerated the optical depth of the hi shells to the photon is vastly reduced and the photon easily escapes , avoiding or reducing the number of further scatterings .
thus , a fermi acclerated photon has a shorter total path - length so that the presence of dust affects it less than an unaccelerated photon .
this results in a slightly enhanced blue peak in the normalised spectra shown in figure [ fig : b1 ] .
including three different dust prescriptions .
for each histogram @xmath142 refers to the optical depth of the inner shell , @xmath143 refers to the optical depth of the outer shell .
the reference dust - free spectrum is shown with a heavy black line.,width=302 ]
figure  [ fig : b2 ] shows that the blue bump remains visible even if we further reduce @xmath109 by order(s ) of magnitude .    , but here we vary @xmath109 while keeping the other parameters fixed .
the prominence of the blue peak reduces with hi column density , but a blue bump clearly remains even when @xmath144 @xmath15.,width=302 ]