 is at $ k \sim 1.5 ~h {\rm Mpc}^{-1}$. Initial velocities for WDM are drawn from a Fermi-Dirac distribution and added to the proper velocity assigned by linear theory.

In our second prescription, WDM is implemented starting directly from a properly designed CAMB power spectrum, by setting $\Delta N_{\rm eff}$ to express departures from the canonical $N_{\rm eff}$ value according to the relic mass – where $N_{\rm eff} = 3.046 + \Delta N_{\rm eff}$. These runs are simply indicated as WDM in Table \[table\_supporting\_sims\_wdm\], and we adopt this implementation for our [*Grid Suite*]{}, since it is more physically motivated and it is also straightforward (although delicate) to include massive neutrinos within this framework. The inclusion is simply achieved by accounting for the neutrino mass degeneracy, where: $$\Delta N_{\rm eff} = \Big ( {93.14 \cdot \Omega_{\rm WDM} h^2 \over m_{\rm WDM} } \Big )^{4/3}.$$ We eventually use the one-to-one mapping between thermal relics and keV sterile neutrinos to simulate massive sterile neutrino cosmologies in the DW mechanism – as previously explained.

Within the second prescription, we also simulate for the first time cosmologies with both massive active and sterile neutrinos, in a consistent fashion. This is achieved by introducing two neutrino eigenstates (one for the active ones, and one for the sterile counterpart) in a three-component simulation, and then by considering a neutrino mass splitting – so that $\Omega_{\nu} = \Omega_{\rm m} - \Omega_{\rm b}$ now contains two distinct contributions. In essence, three active massive-with-degenerate-mass neutrinos contribute as $N_{\rm eff} =3.046$, while a non-thermalized massive sterile neutrino contributes with an additional $\Delta N_{\rm eff}$ to the number of effective neutrino species, directly proportional to its mass. Such runs represent a unique novelty in the literature, and allow us to study the combined effects of massive and sterile neutrinos, as well as to put constraints on pure WDM scenarios directly from Ly$\alpha$ forest data. The latter simulations are indicated as $\sum m_{\nu}^{+} m_{\rm WDM}^{+}$ in Table \[table\_grid\_sims\_cross\_terms\] ([*Grid*]{} runs), and as [*NU\_WDM*]{} in Table \[table\_
is at $ k \sim 1.5 ~h { \rm Mpc}^{-1}$. Initial velocities for WDM are drawn from a Fermi - Dirac distribution and added to the proper velocity impute by analogue theory. 

 In our second prescription, WDM is follow through begin directly from a properly design CAMB might spectrum, by setting $ \Delta N_{\rm eff}$ to carry departures from the canonical $ N_{\rm eff}$ value accord to the relic mass – where $ N_{\rm eff } = 3.046 + \Delta N_{\rm eff}$. These runs are just argue as WDM in Table \[table\_supporting\_sims\_wdm\ ], and we adopt this execution for our [ * Grid Suite * ] { }, since it is more physically motivated and it is also square (although finespun) to include massive neutrinos within this framework. The inclusion body is simply achieved by accounting for the neutrino mass degeneracy, where: $ $ \Delta N_{\rm eff } = \Big ({ 93.14 \cdot \Omega_{\rm WDM } h^2 \over m_{\rm WDM } } \Big) ^{4/3}.$$ We eventually use the one - to - one function between thermal relics and keV aseptic neutrinos to simulate massive sterile neutrino cosmologies in the DW mechanism – as previously explain. 

 Within the second prescription, we also simulate for the first time cosmology with both massive active and sterile neutrinos, in a consistent fashion. This is achieved by introducing two neutrino eigenstates (one for the active ones, and one for the sterile counterpart) in a three - component simulation, and then by consider a neutrino mass splitting – so that $ \Omega_{\nu } = \Omega_{\rm m } - \Omega_{\rm b}$ now control two clear-cut contributions. In essence, three active massive - with - debauched - mass neutrino contribute as $ N_{\rm eff } = 3.046 $, while a non - thermalized massive sterile neutrino contributes with an extra $ \Delta N_{\rm eff}$ to the number of effective neutrino species, directly proportional to its mass. Such runs represent a singular novelty in the literature, and allow us to study the combined impression of massive and aseptic neutrino, equally well as to put constraints on pure WDM scenarios directly from Ly$\alpha$ forest data. The latter simulations are indicated as $ \sum m_{\nu}^{+ } m_{\rm WDM}^{+}$ in Table \[table\_grid\_sims\_cross\_terms\ ] ([ * Grid * ] { } rivulet), and as [ * NU\_WDM * ] { } in Table \[table\ _
 is at $ k \sim 1.5 ~h {\rm Mpc}^{-1}$. Iniuial velocities for WDM ere drasn from x Fermi-Dirac distribution anv adeed ti the proper velocity xssigned hy lineae thtory.

In our second prescriicion, SFM iv implemented sjarting direwtly from a prmpdrpy designed CAMB power spectrum, by fetting $\Dflta N_{\rm eff}$ tj exkrefs dsiavtures from the canonical $N_{\rm efr}$ value according to the relic mass – where $N_{\rm eff} = 3.046 + \Delta N_{\rm eff}$. Hhese runs qre fumply indicaged as WDM in Table \[tagle\_supporting\_sims\_wdm\], and we adoot thns implemenjcrioj for our [*Grmd Suine*]{}, since it is more pvysicalky motivated akd it is also straightforward (although delicate) to incluge massive neutrinow qithit thhs ffqmedori. Vhe includioi is simply achieved bt accounting for tht nqlyrino mass dsgenerwcr, where: $$\Delta N_{\rm eff} = \Big ( {93.14 \cdot \Omega_{\gm WSM} h^2 \over m_{\rm WDM} } \Big )^{4/3}.$$ We eventually use thg one-to-one mapping between thermal relics and keV sterile nautriior ti fknupate massive sterile neutrino cosmologies in ege DE mechanism – af previouslu fxljained.

Within jhe second prescription, we apso simolate dor the fyrst time cosmologies with both massive actpve qnd sterile neutriuos, in a conristgnt faxhion. This is achieved yy intdoducing twl neutrink eigenstates (one fog tha active ones, and one for ehe sterioe cpunterpxrt) on a trree-componfnt slkulation, and then hy cousidesing a neuhrino mass splitting – so that $\Omega_{\nu} = \Omega_{\rk k} - \Mmega_{\rm y}$ now gontains two diftinct contribotions. In essevce, three zctive kassive-with-qegenerate-masv neutrinos cmntributq as $N_{\rm eff} =3.046$, wfkle a non-thermslized massive sterule neutrino contrlbuter with an additiinao $\Delta N_{\rm eff}$ to thq jukbqs of effectiee ndutfono soecies, diregtlh prpportional to its mavs. Shch runs represent a unique bovelty yn the literayure, and allow us uo stuvy the comboneq effects of massive and sterime neutrijos, as well as eo pmt cjnstraints on pure WDM scenarios directly from Ly$\elpha$ forest data. The lqtter simulations ate lndicated as $\sum i_{\nu}^{+} m_{\rm WGM}^{+}$ in Table \[table\_grie\_sims\_cross\_terms\] ([*Gvid*]{} runs), and as [*NU\_WDM*]{} jn Tabne \[tahle\_
is at $ k \sim 1.5 ~h Initial for WDM drawn from a the velocity assigned by theory. In our prescription, WDM is implemented starting directly a properly designed CAMB power spectrum, by setting $\Delta N_{\rm eff}$ to express from the canonical $N_{\rm eff}$ value according to the relic mass – where eff} 3.046 \Delta eff}$. These runs are simply indicated as WDM in Table \[table\_supporting\_sims\_wdm\], and we adopt this implementation our [*Grid Suite*]{}, since it is more physically and it is also (although delicate) to include massive within framework. The is achieved accounting for the mass degeneracy, where: $$\Delta N_{\rm eff} = \Big ( {93.14 \cdot \Omega_{\rm WDM} h^2 \over m_{\rm WDM} \Big )^{4/3}.$$ use the mapping thermal and keV sterile simulate massive sterile neutrino cosmologies in – as previously explained. Within the second prescription, also simulate the first time cosmologies with both active and sterile neutrinos, in a consistent fashion. is achieved by introducing two neutrino eigenstates (one for the active ones, and one for counterpart) in a three-component and then by a mass – that $\Omega_{\nu} \Omega_{\rm m} - \Omega_{\rm b}$ now contains two distinct contributions. In three active massive-with-degenerate-mass neutrinos contribute as $N_{\rm eff} =3.046$, while massive neutrino contributes with additional $\Delta N_{\rm eff}$ the of effective neutrino species, to mass. a novelty the literature, and allow to study the combined effects massive and sterile neutrinos, constraints on pure WDM scenarios directly from Ly$\alpha$ data. The latter simulations are indicated as m_{\nu}^{+} m_{\rm WDM}^{+}$ in Table \[table\_grid\_sims\_cross\_terms\] ([*Grid*]{} runs), and as [*NU\_WDM*]{} in \[table\_
 is at $ k \sim 1.5 ~h {\rm Mpc}^{-1}$. Initial veloCities for WdM are DraWn fRoM a FeRmi-DIrac distributiON and Added to the proper velociTy assIgNEd by LInEar thEory.

In oUR sECOnd PrEsCriPtIOn, wDM is ImpLementeD starting dIreCtLy from a propeRLy Designed CAmB pOwer spectrum, By sEtting $\deLta n_{\Rm eff}$ To eXpresS deparTUres frOm the canoNiCAl $N_{\rm eFF}$ value aCCOrDing To the relic mass – wheRE $N_{\RM eff} = 3.046 + \Delta N_{\rm efF}$. These RuNS aRE SimPly Indicated aS WdM in TABle \[tablE\_SuPPORtiNG\_sims\_wdm\], and we Adopt this imPLemEntatiOn For OUr [*Grid suite*]{}, SiNCe iT is more physIcalLy motivatEd and iT Is also sTRaightfOrward (AltHouGh deLIcAtE) to InCLudE MaSsiVE neUtrinos wItHiN this FramEWORK. The IncLusiOn is sImply achieved By aCcouNTinG for tHe neuTrinO mAss deGeneraCy, wheRe: $$\delta N_{\rm eff} = \Big ( {93.14 \cDot \OMega_{\rm WDM} H^2 \ovEr M_{\rm wDm} } \Big )^{4/3}.$$ WE EventuAllY usE the one-To-one maPPinG bETWEeN thermal relics and kEV STErIle neutrInos to SImUlATe massivE sTerIle nEUTrino CosmOLoGies in thE DW mecHAnIsM – as prevIoUsly exPlAinEd.

WIthin THe seCond prEscriptiOn, we aLSo simulate for tHE first time cosMOlOGIeS With BotH massive actIve aND steRile NEuTriNOs, in a ConsiStENt FAshion. This is achieveD bY introDucinG two neutrino eIgenstates (ONE For the acTive ONeS, And one for the stErile CounterparT) In a three-CompoNent simuLation, and THEn by consIdeRinG a nEutRINo Mass splitting – SO That $\omEga_{\nu} = \OmEga_{\Rm m} - \OmegA_{\rm B}$ noW coNtaInS two distiNct contrIbUtIoNs. in eSsencE, Three actIvE maSsIve-With-dEGeneraTe-masS neuTrInOS coNtributE As $n_{\RM eff} =3.046$, WhIlE a noN-thErMalizEd maSSivE sterilE neutrino ConTRibuTeS wIth an adDitional $\Delta n_{\rM eff}$ to the nUmBer Of effeCTIve neutrIno species, directly propoRTional tO itS mass. such Runs repreSenT a uniqUe nOVelty iN the liTeratUrE, anD ALlow uS TO sTudY tHe combined EFFecTs of mAsSive And sterIle neutrinos, as well AS to Put constraintS on Pure wdm sCenARiOS diReCTly FROm Ly$\alpha$ forest Data. The latTeR SiMulations aRE inDiCated as $\Sum m_{\nu}^{+} m_{\Rm WDM}^{+}$ IN Table \[tAble\_grid\_sIms\_cross\_tErMs\] ([*GrID*]{} RunS), and as [*NU\_WDm*]{} in Table \[Table\_
 is at $ k \sim 1.5 ~h {\r m Mpc}^{-1 }$. I nit ial  v eloc itie s for WDM ared rawn  from a Fermi-Dirac di strib ut i on a n dadded  to the pr o p erve lo cit ya ss igned  by  linear  theory.
Inou r second pre s cr iption, WD M i s implemente d s tartin gdir e ctlyfro m a p roperl y  desig ned CAMBpo w er spe c trum, b y se ttin g $\Delta N_{\rme ff } $ to express d epartu re s  f r o m t hecanonical$N _{\rm eff}$ v a lu e a cco r ding to the r elic mass – whe re $N_ {\ rme ff} =3.046  + \De lta N_{\rmeff} $. Theseruns a r e simpl y  indica ted as  WD M i n Ta b le  \ [ta bl e \_s u pp ort i ng\ _sims\_w dm \] , and  wea d o p t th isimpl ement ation for our  [* Grid Sui te*]{ }, si nceit  is m ore ph ysica ll y motivated and  itis also s tra ig htf or ward( althou ghdel icate)to incl u dema s s i ve  neutrinos withinth i s  f ramework . Thei nc lu s ion is s im ply  ach i e ved b y ac c ou nting fo r then eu tr ino mas sdegene ra cy,  wh ere:$ $\De lta N_ {\rm eff } = \ B ig ( {93.14 \c d ot \Omega_{\r m  W D M }h ^2 \ ove r m_{\rm WD M} } \Big  )^{ 4 /3 }.$ $  We e ventu al l yu se the one-to-one m ap ping b etwee n thermal rel ics and ke V s terile n eutr i no s  to simulate m assiv e sterilen eutrinocosmo logies i n the DWm e chanism– a s p rev iou s l yexplained.

W i t hinth e secon d p rescrip tio n,weals osimulatefor thefi rs tti mecosmo l ogies wi th  bo th  ma ssive active  andster il en eut rinos,i na cons is te nt f ash io n. Th is i s  ac hievedby introd uci n g tw one utrinoeigenstates ( on e for theac tiv e ones , and onefor the sterile counter p art) in  athree -com ponent si mul ation,  an d  thenby con sider in g a n eutri n o  m ass  s plitting – s o t hat $ \O mega _{\nu}= \Omega_{\rm m} - \Om ega_{\rm b}$now  con t a in s t w od ist in c t c o n tributions. Inessence, t hr e eactive mas s ive -w ith-deg enerate -mass neutrin os contri bute as $ N_ {\rm e ff}  =3.046$,while anon-therm a lized ma ssive  st erilene utr ino c ontrib u tes  with  an ad di tional  $\De lt a N_{\rm  eff}$ to the number of  effec tiveneu trino spe cie s , d irectly p ropo rtional to  it s m ass.Suc h  runs  rep r es ent a uni quen ovelty in th e l i t er ature, anda l l owus to  st u dy the  com bined effects ofm assive and ste rile n eut rin o s, a swell as to put  co ns t r aints on  p ure WDM sce narios d ir e ctlyfrom L y$\alp ha$ for e s td ata. T he l att er simula tio ns are ind ic at e d as $ \sum  m _{\nu} ^{+} m _ {\rm W DM}^{+}$ in Tabl e \[t a b le\_g r id\ _sims \_ cross\_ t erms \] ([*Grid *]{} runs),  and a s [* NU\_W DM*]{}in  Table  \[ ta ble\_
 is_at $_k \sim 1.5 ~h_{\rm Mpc}^{-1}$._Initial_velocities for_WDM_are drawn from_a Fermi-Dirac distribution_and added to the_proper velocity assigned_by_linear theory.

In our second prescription, WDM is implemented starting directly from a properly designed_CAMB_power spectrum,_by_setting_$\Delta N_{\rm eff}$ to express_departures from the canonical $N_{\rm_eff}$ value_according to the relic mass – where $N_{\rm_eff}_= 3.046 +_\Delta N_{\rm eff}$. These runs are simply indicated as_WDM in Table \[table\_supporting\_sims\_wdm\], and we_adopt this implementation_for_our_[*Grid Suite*]{}, since it_is more physically motivated and it_is also straightforward (although delicate) to_include massive neutrinos within this framework. The_inclusion is simply achieved by accounting_for the neutrino mass degeneracy,_where: $$\Delta_N_{\rm eff} = \Big (_{93.14 \cdot \Omega_{\rm_WDM} h^2_\over m_{\rm WDM}_} \Big )^{4/3}.$$ We eventually use_the one-to-one mapping_between thermal relics and keV sterile_neutrinos_to simulate massive_sterile_neutrino_cosmologies in_the DW mechanism_–_as previously_explained.

Within_the second prescription, we also simulate_for_the first time cosmologies with both massive_active and sterile neutrinos,_in_a consistent fashion. This_is achieved by introducing two_neutrino eigenstates (one for the active_ones, and_one for_the sterile counterpart) in a three-component simulation, and then by considering_a neutrino mass splitting – so_that $\Omega_{\nu} = \Omega_{\rm_m} -_\Omega_{\rm_b}$ now contains_two_distinct contributions._In essence, three active massive-with-degenerate-mass neutrinos contribute_as $N_{\rm_eff} =3.046$, while a non-thermalized massive_sterile neutrino contributes with_an_additional $\Delta N_{\rm eff}$ to the_number of effective neutrino species, directly_proportional to its mass. Such_runs_represent_a unique novelty in the_literature, and allow us to study_the combined effects_of massive and sterile neutrinos, as well_as_to put constraints on pure WDM_scenarios_directly from Ly$\alpha$ forest data. The_latter_simulations_are indicated as $\sum m_{\nu}^{+}_m_{\rm WDM}^{+}$ in Table \[table\_grid\_sims\_cross\_terms\] ([*Grid*]{}_runs), and as [*NU\_WDM*]{} in Table \[table\_