}{\rm Mpc}$ runs at each targeted redshift. Different line styles and colors refer to the varying mass resolutions, as specified in the plot. The right panels display 1D flux power spectrum ratios for realizations having $N^3=128^3, 256^3, 336^3$, respectively, estimated with respect to the highest-resolution run (i.e. “L10\_N512", having an equivalent grid resolution of $20h^{-1}{\rm kpc}$), expressed in percentage, and with identical colors and line styles as in the left panel. Error bars are 1-$\sigma$ deviations computed from $10,000$ simulated skewers randomly extracted at each redshift from the simulation boxes, and the gray areas highlight the $1\%$ level. As can be clearly seen, 1D flux power spectra obtained from the $N^3=336^3$ realization (i.e. $30h^{-1}{\rm kpc}$ equivalent grid resolution) are within $\sim 1\%$ in the $k$-range covered by eBOSS (up to $k=2 \times 10^{-2} [{\rm km/s}]^{-1}$), and closer to the same degree of accuracy for the expected extension of DESI Ly$\alpha$ forest data. Clearly, achieving sub-percentage convergence at higher-$z$ is progressively challenging, but a grid of $30h^{-1}{\rm kpc}$ is sufficient for a satisfying convergence even at those redshifts.

The previous rationale is precisely what has driven the overall architecture of the [*Sejong Grid Suite*]{}, designed to achieve an equivalent resolution up to $3 \times 3328^3$ = 110 billion particles in a $(100h^{-1}{\rm Mpc})$ box, corresponding to a $30 h^{-1}{\rm kpc}$ mean grid resolution that ensures a convergence on Ly$\alpha$ flux statistics closer to the desired $\sim 1.0\%$ level that the final DESI data will provide.

Clearly, running a very large number of simulations with $2 \times 3328^3$ or $3 \times 3328^3$ elements over a $(100h^{-1}{\rm Mpc})$ box is still computationally challenging, particularly when neutrinos are included as particles. The global computational cost in performing an entire suite at such resolution would easily require $\sim 100$ million CPU hours. Therefore, as in our previous release, we adopt a splicing technique. In this respect, the three simulations per fixed parameter set of the [*Se
} { \rm Mpc}$ runs at each targeted redshift. Different course stylus and colors denote to the deviate mass resolutions, as specify in the plot. The correct panels expose 1D magnetic field power spectrum proportion for realizations having $ N^3=128 ^ 3, 256 ^ 3, 336 ^ 3 $, respectively, calculate with regard to the highest - settlement run (i.e. “ L10\_N512 ", having an equivalent power system resolution of $ 20h^{-1}{\rm kpc}$), expressed in percentage, and with identical color and line styles as in the leftover panel. erroneousness bars are 1-$\sigma$ deviation computed from $ 10,000 $ simulated skewers randomly excerpt at each redshift from the simulation boxes, and the gray areas highlight the $ 1\%$ level. As can be clearly seen, 1D flux power spectra obtained from the $ N^3=336 ^ 3 $ realization (i.e. $ 30h^{-1}{\rm kpc}$ equivalent grid resolution) are within $ \sim 1\%$ in the $ k$-range covered by eBOSS (up to $ k=2 \times 10^{-2 } [ { \rm km / s}]^{-1}$), and closer to the same academic degree of accuracy for the expect extension of DESI Ly$\alpha$ forest datum. intelligibly, achieve sub - percentage convergence at higher-$z$ is progressively challenging, but a grid of $ 30h^{-1}{\rm kpc}$ is sufficient for a satisfying convergence even at those redshifts. 

 The previous rationale is precisely what has drive the overall architecture of the [ * Sejong Grid Suite * ] { }, designed to achieve an equivalent resolution up to $ 3 \times 3328 ^ 3 $ = 110 billion particles in a $ (100h^{-1}{\rm Mpc})$ box, equate to a $ 30 h^{-1}{\rm kpc}$ mean grid resolution that ensures a convergence on Ly$\alpha$ magnetic field statistics nearer to the desire $ \sim 1.0\%$ level that the final DESI data will provide. 

 intelligibly, running a very large number of model with $ 2 \times 3328 ^ 3 $ or $ 3 \times 3328 ^ 3 $ elements over a $ (100h^{-1}{\rm Mpc})$ box is still computationally challenging, particularly when neutrino are included as particles. The global computational cost in performing an integral suite at such resolution would easily require $ \sim 100 $ million central processing unit hours. Therefore, as in our former release, we adopt a splicing technique. In this respect, the three pretense per fixed parameter set of the [ * Se
}{\rm Lpc}$ runs at each targetea redshift. Diffgrwnt liie stylss and culors refer to the varying mess eesolytions, as specified in the plot. The rigyt penels display 1D hmux powcx spedbrum xavios for realizstions havhng $N^3=128^3, 256^3, 336^3$, respecdixepy, estimated with respect to the hidhest-rexopution run (i.e. “J10\_N512", hsding an equivalent grid resolution of $20g^{-1}{\rm kpc}$), expressed in percentage, and with identlcal colors and line shyles as in the oeft panel. Efror bars are 1-$\sigma$ debiations computed from $10,000$ simulatdd skzwers randonlt edjracted at eech reqshift from bne simglation boxes, and the grey aeeas highlight the $1\%$ lxvel. As can be clearjy seen, 1D fmux power spectra ovtaingd frmm tfw $N^3=336^3$ rezlmzafion (i.f. $30h^{-1}{\cm kpc}$ equibalent grid resolution) are witnig $\sim 1\%$ in the $i$-range cjvered by eBOSS (up to $k=2 \times 10^{-2} [{\rm km/s}]^{-1}$), atd dloser to the same degrwe of accuracy for thg expected extension of DESI Ly$\alpha$ forest data. Clearly, acvieviig suy-icrcevrahe convergence at higher-$z$ is progressively chwmltngpng, but a grid of $30h^{-1}{\rm kpc}$ is skfgycient for a ratisfvjnf convergence even at thofe reeshifts.

Tht prefious rationale is preciselt what has dgiveb the overall archntecture of che [*Sekong Brid Suite*]{}, designed to cchievs an equivapent resomjtion up to $3 \timds 3328^3$ = 110 billion particles in a $(100h^{-1}{\ri Mpc})$ box, corxespondivg tp a $30 h^{-1}{\wm kpc}$ meaj grid resolution that ejsured d convergejce on Ly$\alpha$ flux statistics rkoser to the dasiged $\sim 1.0\%$ jevel that the finwl DESI data wnll provnde.

Clexrly, runnihg a vecy large numfer of simuladlons with $2 \tmmes 3328^3$ or $3 \timws 3328^3$ wlementr over a $(100h^{-1}{\rm Mpv})$ box is still compytationally challekging, larticularly whzu beutrinos are imcljdeq ws pwsticles. The clobxl zpmputxtional cosb iv petforming an entire sgite at such resolutiom rould eawily reqtire $\sim 100$ milkion CPU hours. Thegeforx, as ii our lredious release, we adopt a splicjng technlquc. In this res[ect, the three snmulations per fixed parameter set of thx [*Se
}{\rm Mpc}$ runs at each targeted redshift. styles colors refer the varying mass plot. right panels display flux power spectrum for realizations having $N^3=128^3, 256^3, 336^3$, estimated with respect to the highest-resolution run (i.e. “L10\_N512", having an equivalent grid of $20h^{-1}{\rm kpc}$), expressed in percentage, and with identical colors and line styles in left Error are 1-$\sigma$ deviations computed from $10,000$ simulated skewers randomly extracted at each redshift from the simulation and the gray areas highlight the $1\%$ level. can be clearly seen, flux power spectra obtained from $N^3=336^3$ (i.e. $30h^{-1}{\rm equivalent resolution) within $\sim 1\%$ the $k$-range covered by eBOSS (up to $k=2 \times 10^{-2} [{\rm km/s}]^{-1}$), and closer to the same of accuracy expected extension DESI forest Clearly, achieving sub-percentage higher-$z$ is progressively challenging, but a kpc}$ is sufficient for a satisfying convergence even those redshifts. previous rationale is precisely what has the overall architecture of the [*Sejong Grid Suite*]{}, to achieve an equivalent resolution up to $3 \times 3328^3$ = 110 billion particles in Mpc})$ box, corresponding to $30 h^{-1}{\rm kpc}$ grid that a on Ly$\alpha$ statistics closer to the desired $\sim 1.0\%$ level that the final data will provide. Clearly, running a very large number of $2 3328^3$ or $3 3328^3$ elements over a Mpc})$ is still computationally challenging, neutrinos included global cost performing an entire suite such resolution would easily require 100$ million CPU hours. release, we adopt a splicing technique. In this the three simulations per fixed parameter set the [*Se
}{\rm Mpc}$ runs at each targeted reDshift. DiffErent LinE stYlEs anD colOrs refer to the vARyinG mass resolutions, as specIfied In THe plOT. THe rigHt panelS DiSPLay 1d fLuX poWeR SpEctruM raTios for RealizatioNs hAvIng $N^3=128^3, 256^3, 336^3$, respectiVElY, estimated WitH respect to thE hiGhest-rEsOluTIon ruN (i.e. “l10\_N512", havIng an eQUivaleNt grid resOlUTion of $20H^{-1}{\Rm kpc}$), exPREsSed iN percentage, and witH IdENtical colors anD line sTyLEs AS In tHe lEft panel. ErRoR bars ARe 1-$\sigma$ DEvIATIonS Computed from $10,000$ sImulated skeWErs RandomLy ExtRActed aT each ReDShiFt from the siMulaTion boxes, And the GRay areaS HighligHt the $1\%$ lEveL. As Can bE ClEaRly SeEN, 1D fLUx PowER spEctra obtAiNeD from The $N^3=336^3$ REALIzatIon (I.e. $30h^{-1}{\rM kpc}$ eQuivalent grid ResOlutIOn) aRe witHin $\siM 1\%$ in tHe $K$-rangE coverEd by ebOsS (up to $k=2 \times 10^{-2} [{\rm kM/s}]^{-1}$), anD closer to The SaMe dEgRee of ACcuracY foR thE expectEd extenSIon Of desi LY$\alpha$ forest data. ClEaRLY, aChieving Sub-perCEnTaGE convergEnCe aT higHER-$z$ is pRogrESsIvely chaLlengiNG, bUt A grid of $30H^{-1}{\rM kpc}$ is SuFfiCieNt for A SatiSfying ConvergeNce evEN at those redshiFTs.

The previous RAtIONaLE is pRecIsely what haS driVEn thE oveRAlL arCHitecTure oF tHE [*SEJong Grid Suite*]{}, designEd To achiEve an Equivalent resOlution up tO $3 \TIMes 3328^3$ = 110 billiOn paRTiCLes in a $(100h^{-1}{\rm Mpc})$ boX, corrEsponding tO A $30 h^{-1}{\rm kpc}$ mEan grId resoluTion that eNSUres a conVerGenCe oN Ly$\ALPhA$ flux statistiCS ClosEr To the deSirEd $\sim 1.0\%$ leVel ThaT thE fiNaL DESI data Will provIdE.

CLeArLy, rUnninG A very larGe NumBeR of SimulATions wIth $2 \tiMes 3328^3$ oR $3 \tImES 3328^3$ elEments oVEr A $(100H^{-1}{\Rm MpC})$ bOx Is stIll CoMputaTionALly ChallenGing, partiCulARly wHeN nEutrinoS are included aS pArticles. ThE gLobAl compUTAtional cOst in performing an entire SUite at sUch ResolUtioN would easIly RequirE $\siM 100$ MillioN CPU hoUrs. ThErEfoRE, As in oUR PrEviOuS release, we ADOpt A spliCiNg teChnique. in this respect, the thREe sImulations per FixEd paRAMeTer SEt OF thE [*SE
}{\rm Mpc}$ runs at each t argeted re dshif t.Dif fe rent  lin e styles and c o lors  refer to the varyingmassre s olut i on s, as  specif i ed i n t he  p lot .T he  righ t p anels d isplay 1Dflu xpower spectr u mratios for  re alizations h avi ng $N^ 3= 128 ^ 3, 25 6^3 , 336 ^3$, r e specti vely, est im a ted wi t h respe c t  t o th e highest-resolut i on run (i.e. “L10 \_N512 ", ha v i nganequivalent  g rid r e solutio n  o f $ 20h ^ {-1}{\rm kpc} $), express e d i n perc en tag e , andwithid e nti cal colorsandline styl es asi n the l e ft pane l. Err orbar s ar e  1 -$ \si gm a $ d e vi ati o nscomputed  f ro m $10 ,000 $ s i mula ted  ske wersrandomly extr act ed a t  ea ch re dshif t fr om  thesimula tionbo xes, and the gr ay a reas high lig ht  th e$1\%$ level.  As  ca n be cl early s e en,  1 D f lu x power spectra ob ta i n ed  from th e $N^3 = 33 6^ 3 $ realiz at ion  (i. e .  $30h ^{-1 } {\ rm kpc}$  equiv a le nt  grid r es olutio n)  ar e w ithin $\si m 1\%$  in the$k$-r a nge covered by eBOSS (up to$ k= 2 \t i mes10^ {-2} [{\rmkm/s } ]^{- 1}$) ,  a ndc loser  to t he sa m e degree of accurac yfor th e exp ected extensi on of DESI L y $\alpha$  for e st data. Clearly,  achi eving sub- p ercentag e con vergence  at highe r - $z$ is p rog res siv ely c ha llenging, but a  gri dof $30h ^{- 1}{\rmkpc }$issuf fi cient for  a satis fy in gco nve rgenc e  even at  t hos ered shift s .

The  prev ious  r at i ona le is p r ec i s elywh at  has  dr iv en th e ov e ral l archi tecture o f t h e [* Se jo ng Grid  Suite*]{}, d es igned to a ch iev e an e q u ivalentresolution up to $3 \ti m es 3328 ^3$  = 11 0 bi llion par tic les in  a$ (100h^ {-1}{\ rm Mp c} )$b o x, co r r es pon di ng to a $3 0 h^{ -1}{\ rm  kpc }$ mean  grid resolution t h atensures a con ver genc e on  Ly $ \a l pha $f lux s tatistics close r to the d es i re d $\sim 1. 0 \%$  l evel th at thefinal DESI da ta will p rovide.
Cl earl y ,  ru nning a ve ry large  number o f  simu l at ionswit h $2 \ ti mes  3328 ^3$ or $3\time s 3328 ^3 $ elem entsov er a $(1 00h^{-1}{\rm Mpc})$ box  is st ill c omp utational lyc hal lenging,part icularly w hen  ne utrin osa re in clud e dasp artic les. The globa l  c omp u t at ional costi n per formi nga n enti re s uite at such reso l ution would ea sily r equ ire $\si m100$ million C PUho u r s. There fo re, as in o ur previ ou s  rele ase, w e adop t a spl i c in g  techn ique . I n this re spe ct ,  the th re es imulat ions  p er fix ed par a mete r set of the [*Se
}{\rm Mpc}$_runs at_each targeted redshift. Different_line styles_and_colors refer_to_the varying mass_resolutions, as specified_in the plot. The_right panels display_1D_flux power spectrum ratios for realizations having $N^3=128^3, 256^3, 336^3$, respectively, estimated with respect_to_the highest-resolution_run_(i.e._“L10\_N512", having an equivalent grid_resolution of $20h^{-1}{\rm kpc}$), expressed_in percentage,_and with identical colors and line styles as_in_the left panel._Error bars are 1-$\sigma$ deviations computed from $10,000$ simulated_skewers randomly extracted at each redshift_from the simulation_boxes,_and_the gray areas highlight_the $1\%$ level. As can be_clearly seen, 1D flux power spectra_obtained from the $N^3=336^3$ realization (i.e. $30h^{-1}{\rm_kpc}$ equivalent grid resolution) are within_$\sim 1\%$ in the $k$-range_covered by_eBOSS (up to $k=2 \times_10^{-2} [{\rm km/s}]^{-1}$),_and closer_to the same_degree of accuracy for the expected_extension of DESI_Ly$\alpha$ forest data. Clearly, achieving sub-percentage_convergence_at higher-$z$ is_progressively_challenging,_but a_grid of $30h^{-1}{\rm_kpc}$_is sufficient_for_a satisfying convergence even at those_redshifts.

The_previous rationale is precisely what has driven_the overall architecture of_the_[*Sejong Grid Suite*]{}, designed_to achieve an equivalent resolution_up to $3 \times 3328^3$ =_110 billion_particles in_a $(100h^{-1}{\rm Mpc})$ box, corresponding to a $30 h^{-1}{\rm kpc}$ mean_grid resolution that ensures a convergence_on Ly$\alpha$ flux statistics_closer to_the_desired $\sim 1.0\%$_level_that the_final DESI data will provide.

Clearly, running a_very large_number of simulations with $2 \times_3328^3$ or $3 \times_3328^3$_elements over a $(100h^{-1}{\rm Mpc})$ box_is still computationally challenging, particularly when_neutrinos are included as particles._The_global_computational cost in performing an_entire suite at such resolution would_easily require $\sim_100$ million CPU hours. Therefore, as in_our_previous release, we adopt a splicing_technique._In this respect, the three simulations_per_fixed_parameter set of the [*Se