 redshift evolution of baryonic features in the galaxy power spectrum; iii) weak-lensing tomography derived from measurement of galaxy shear patterns and iv) the luminosity distance–redshift relation of type Ia SNe. We focus our attention on theoretical issues related to the first three of these tests.

The power spectrum of fluctuations at recombination is calculated to high accuracy from linear theory$^{\ino}$, so the problem of realizing, through direct simulation, the evolution of a finite patch of a particular world model from this epoch forward is well posed. To support DES-like observations, one would like to evolve multiple regions of Hubble Length dimension with the principal clustered matter components — dark matter and multiple phases of baryons (stars and cold gas in galaxies, warm/hot gas surrounding galaxies and in groups/clusters) — represented by multiple fluids. Mapping observable signatures of the theoretical solution along the past light-cone of synthetic observers in the computational volume allows ‘clean’ mock surveys to be created, which can further be ‘dirtied’ by real-world observational effects in support of survey data analysis.

Two fundamental barriers stand in the way of achieving a complete solution to this problem. One is the wide dynamic range of non-linear structures (sub-parsecs to gigaparsecs in length, for example), and the other is the complexity of astrophysical processes controlling the baryonic phases. Since DES-like surveys probe only the higher mass portion of the full spectrum of cosmic structures, the first issue is not strongly limiting. The DES will, however, require understanding galaxy and cluster observable signatures, so uncertainties in the treatment of baryonic physics will play a central role. In a companion paper$^{\ino}$, we outline theoretical uncertainties associated with the large-scale structure channels DES will use to test dark energy. Here, we offer a critique of the computational methods that provide theoretical support for DES and similar surveys.

Challenges for Computational Cosmology 
=======================================

Given the wide scale and scope of large-scale structure, a number of approaches have evolved to address restricted elements of the full problem. Since cold dark matter dominates the matter density, N-body methods that follow collisionless clustering have played an important role in defining the overall evolution of the cosmic web and the structure of collapsed halos formed within it. Combined N-body and gas dynamics techniques explore gravitationally coupled evolution of baryons and dark matter. Knowledge gained from these ‘direct’ approaches led to the development of ‘semi-analytic’
redshift evolution of baryonic features in the galaxy power spectrum; iii) weak - lensing imaging derive from measurement of galaxy shear patterns and iv) the luminosity distance – red shift relation of character Ia SNe. We focus our attention on theoretical issue related to the first three of these trial. 

 The power spectrum of fluctuations at recombination is forecast to high accuracy from linear theory$^{\ino}$, so the problem of realize, through lineal simulation, the evolution of a finite while of a particular world exemplar from this epoch forward is well posed. To support DES - comparable observations, one would like to evolve multiple regions of Hubble Length property with the principal clustered matter components — dark topic and multiple phases of baryons (stars and cold gas in galaxies, warm / blistering gasoline surrounding galaxies and in groups / clusters) — represented by multiple fluids. Mapping observable signatures of the theoretical solution along the past luminosity - cone of synthetic observers in the computational book allow ‘ clean ’ mock surveys to be created, which can far be ‘ dirty ’ by real - world observational effect in support of survey data analysis. 

 Two fundamental barrier stand in the way of achieving a complete solution to this trouble. One is the wide active range of non - linear structures (sub - parsecs to gigaparsecs in length, for case), and the other is the complexity of astrophysical processes controlling the baryonic phase. Since DES - like surveys probe only the higher multitude portion of the full spectrum of cosmic structure, the first topic is not strongly restrict. The DES will, however, require reason galaxy and cluster discernible signatures, so uncertainties in the treatment of baryonic physics will play a central role. In a companion paper$^{\ino}$, we outline theoretical uncertainties associated with the large - scale social organization channels DES will use to screen dark energy. Here, we offer a critique of the computational methods that provide theoretical support for DES and exchangeable surveys. 

 Challenges for Computational Cosmology 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 Given the wide scale and scope of large - plate structure, a number of approaches have evolved to address restricted elements of the full trouble. Since cold dark matter dominates the matter density, N - body methods that follow collisionless clustering have played an significant role in defining the overall evolution of the cosmic web and the structure of collapsed aura formed within it. Combined N - body and gas dynamics techniques explore gravitationally copulate evolution of baryons and dark matter. Knowledge gained from these ‘ direct ’ approaches led to the development of ‘ semi - analytic ’
 refshift evolution of baryunic features iu the gelaxy pkwer speztrum; iii) weak-lensing tomogrephy dericed from measurement ow galaxy dhear parteris and iv) the luminosity distahge–redviift relation on type Ia STe. We focus ous xtcention on theoretical issues relateq to thr virst three of thexq teans.

Uhe power spectrum of fluctuationa at rebombination is cakculated to high accuracy vrom linear theory$^{\ino}$, do the proboem jd realizing, ghrough digzct simulatjon, the evolution of a finite pxtch pf a partixuoar forld model from this epoch nprward is welk posed. To supiort VES-luke observations, one xould like to evolve multiple rzgions of Hubble Lengrh dimetsiot wigy tfe kriicilal clkstxred matter components — dark matter and mileplle phases or baryjnf (stars and cold gas in galaxies, warm/hou gas surrounding galaxies abd in groups/clusters) — represeneed by multiple fluids. Mapping observable signatuses oh ghe tmeordrifal solution along the past light-cone of syntrstoc observers in bhe computational fopukg allows ‘clean’ mock surbeys to be created, which san fyrther be ‘diryied’ by real-world observatiinal effects un support of survzy data analvsis.

Twp funcamental barriers stand in fhe way of wchieving x complete solutkon tm this problem. One is the ride dynanic xange of non-kinear structured (sub-idrsecs to gigaparsfcs iu lencth, for exwmple), and the other is the complexity of astrppvyspcal proczsses gontrolling the baryonic phasgs. Since BES-likd surveys irobe onlb the higher mass portion lf the full vpectrum of xosmuc struzgures, the firsy issue is not strobgly limiting. The VES djll, however, reqbnrt ynderstanding balxxy ajd cjgster observdble siematurds, so uncerbaivtiex in the treatment ox badyonic physics wilk ilay a cebtral roje. In a compamion paper$^{\ino}$, we oktlinx theoceticak ugcertainties associated with tge large-sfalc structure cranncls QES will uxe to test dark energy. Here, we offer a rritique of the computarional methods that krovide theorxtical support xor DES and similar wurveys.

Challenges for Computational Coamologf 
=======================================

Givfn the wide scale and scope of large-scale structure, a number of approachew have evolved to zddrrss rvstxicced eleientd of the full pronlem. Since cold dark matter dommnates the mctter density, N-body methods tnag follow collksionless clustering havs played an important role in defining yhe overall evolution of the cosnic wsb and the structure of collspsed kalos dormed qithln it. Combined N-bpdn and gas fynakies techniques explore gravitationqlly ciuoled evolution oh baryons anf dark mattet. Knowledge gaines frkl there ‘direct’ apptoaches led to the development of ‘veii-analytic’
redshift evolution of baryonic features in the spectrum; weak-lensing tomography from measurement of the distance–redshift relation of Ia SNe. We our attention on theoretical issues related the first three of these tests. The power spectrum of fluctuations at recombination calculated to high accuracy from linear theory$^{\ino}$, so the problem of realizing, through simulation, evolution a patch of a particular world model from this epoch forward is well posed. To support DES-like one would like to evolve multiple regions of Length dimension with the clustered matter components — dark and phases of (stars cold in galaxies, warm/hot surrounding galaxies and in groups/clusters) — represented by multiple fluids. Mapping observable signatures of the theoretical solution the past synthetic observers the volume ‘clean’ mock surveys created, which can further be ‘dirtied’ effects in support of survey data analysis. Two barriers stand the way of achieving a complete to this problem. One is the wide dynamic of non-linear structures (sub-parsecs to gigaparsecs in length, for example), and the other is the astrophysical processes controlling the phases. Since DES-like probe the mass of the spectrum of cosmic structures, the first issue is not strongly limiting. DES will, however, require understanding galaxy and cluster observable signatures, in treatment of baryonic will play a central In companion paper$^{\ino}$, we outline associated the DES use test dark energy. Here, offer a critique of the methods that provide theoretical surveys. Challenges for Computational Cosmology ======================================= Given the scale and scope of large-scale structure, a of approaches have evolved to address restricted elements of the full problem. cold dark the matter density, N-body methods that follow collisionless have played an important in defining the overall evolution of the cosmic web the of collapsed formed within it. N-body and gas techniques explore gravitationally of baryons dark Knowledge ‘direct’ approaches led to the development ‘semi-analytic’
 redshift evolution of baryonIc features In the GalAxy PoWer sPectRum; iii) weak-lensINg toMography derived from meaSuremEnT Of gaLAxY sheaR patterNS aND Iv) tHe LuMinOsITy DistaNce–RedshifT relation oF tyPe ia SNe. We focus OUr Attention oN thEoretical issUes RelateD tO thE First ThrEe of tHese teSTs.

The pOwer spectRuM Of flucTUations AT ReCombInation is calculatED tO High accuracy frOm lineAr THeORY$^{\inO}$, so The problem Of RealiZIng, throUGh DIREct SImulation, the eVolution of a FIniTe patcH oF a pARticulAr worLd MOdeL from this epOch fOrward is wEll posED. To suppORt DES-liKe obseRvaTioNs, onE WoUlD liKe TO evOLvE muLTipLe regionS oF HUbble lengTH DIMensIon With The prIncipal clusteRed MattER coMponeNts — daRk maTtEr and MultipLe phaSeS of baryons (stars And cOld gas in gAlaXiEs, wArM/hot gAS surroUndIng GalaxieS and in gROupS/cLUSTeRs) — represented by mulTiPLE fLuids. MapPing obSErVaBLe signatUrEs oF the THEoretIcal SOlUtion aloNg the pASt LiGht-cone Of SyntheTiC obSerVers iN The cOmputaTional voLume aLLows ‘clean’ mock sURveys to be creaTEd, WHIcH Can fUrtHer be ‘dirtieD’ by rEAl-woRld oBSeRvaTIonal EffecTs IN sUPport of survey data anAlYsis.

TwO fundAmental barrieRs stand in tHE WAy of achiEvinG A cOMplete solution To thiS problem. OnE Is the widE dynaMic range Of non-lineAR StructurEs (sUb-pArsEcs TO GiGaparsecs in leNGTh, foR eXample), aNd tHe other Is tHe cOmpLexItY of astropHysical pRoCeSsEs ConTrollINg the barYoNic PhAseS. SincE dES-likE survEys pRoBe ONly The highER mASS porTiOn Of thE fuLl SpectRum oF CosMic struCtures, the FirST issUe Is Not stroNgly limiting. THe dES will, howEvEr, rEquire UNDerstandIng galaxy and cluster obseRVable siGnaTures, So unCertaintiEs iN the trEatMEnt of bAryoniC physIcS wiLL Play a CENtRal RoLe. In a compaNIOn pAper$^{\iNo}$, We ouTline thEoretical uncertainTIes Associated witH thE larGE-ScAle STrUCtuRe CHanNELs DES will use to tEst dark eneRgY. heRe, we offer a CRitIqUe of the ComputaTionaL Methods That proviDe theoretIcAl suPPOrt For DES and sImilar suRveys.

ChalLEnges FOr compuTatIonal COsMolOgy 
=======================================

GiVen the WIde Scale And scoPe Of largE-scalE sTructure, A number of approaches have EvolveD to adDreSs restricTed ELemEnts of the Full Problem. SinCe cOld Dark mAttER domiNateS ThE maTTer deNsitY, n-body methODs ThaT FOlLow collisioNLESs cLusteRinG Have plAyed An important role in DEfining the overAll eVOLutIon OF the CoSmic web and the sTruCtURE of collaPsEd halos formEd within It. cOmbinEd N-bodY and gaS dynamiCS TeCHniqueS expLorE gravitatIonAlLY coupleD eVoLUtion oF barYoNs and dArk matTEr. KnOWLedge gained from tHese ‘dIREct’ apPRoaChes lEd To the deVElopMent of ‘semi-Analytic’
 redshift evolution of bar yonic feat uresinthe  g alax y po wer spectrum;i ii)weak-lensing tomograph y der iv e d fr o mmeasu remento fg a lax ysh ear  p a tt ernsand  iv) th e luminosi tydi stance–redsh i ft  relationoftype Ia SNe.  We  focus  o ura ttent ion  on t heoret i cal is sues rela te d  to th e  firstt h re e of  these tests.

Th e  p o wer spectrum o f fluc tu a ti o n s a t r ecombinati on  is c a lculate d  t o h igh accuracy from  linear the o ry$ ^{\ino }$ , s o  the p roble mo f r ealizing, t hrou gh direct  simul a tion, t h e evolu tion o f a  fi nite pa tc h o fa  pa r ti cul a r w orld mod el  f rom t hise p o c h fo rwa rd i s wel l posed. To s upp ortD ES- likeobser vati on s, on e woul d lik eto evolve multi pleregions o f H ub ble  L ength dimens ion  wi th theprincip a l c lu s t e re d matter component s— da rk matte r andm ul ti p le phase sofbary o n s (st arsa nd  cold ga s in g a la xi es, war m/ hot ga ssur rou nding gala xies a nd in gr oups/ c lusters) — rep r esented by mu l ti p l ef luid s.Mapping obs erva b le s igna t ur eso f the  theo re t ic a l solution along th epast l ight- cone of synth etic obser v e r s in the  com p ut a tional volumeallow s ‘clean’m ock surv eys t o be cre ated, whi c h  can fur the r b e ‘ dir t i ed ’ by real-wor l d  obs er vationa l e ffectsinsup por t o fsurvey da ta analy si s. 
Tw o f undam e ntal bar ri ers  s tan d int he way  of a chie vi ng a c omplete so l u tion  t othis  pr ob lem.Onei s t he wide  dynamicran g e of  n on -linear  structures ( su b-parsecsto  gi gapars e c s in len gth, for example), andt he othe r i s the  com plexity o f a stroph ysi c al pro cesses  cont ro lli n g  theb a ry oni cphases. Si n c e D ES-li ke  sur veys pr obe only the highe r  ma ss portion of  th e fu l l  s pec t ru m  of  c o smi c structures, the  first iss ue is  not stron g lyli miting.  The DE S wil l , howev er, requi re unders ta ndin g gal axy and cl uster ob servables ignat u re s, so  un certai nt ies  in t he tre a tme nt of  baryo ni c phys ics w il l play a  central role. In a com panion  pape r$^ {\ino}$,weo utl ine theor etic al uncerta int ies  asso cia t ed wi th t h elar g e-sca le s t ructure c h an nel s DE S will uset o tes t dar k e n ergy.Here , we offer a crit i que of the com puta t i ona l m e thod sthat provide t heo re t i cal supp or t for DES a nd simil ar surve ys.

C hallen ges for C om p utatio nalCos mology 
= === == = ======= == == = ====== ==== == ====== ==

Gi v en t h e  wide scale andscope o f lar g e-s calest ructure ,  a n umber of a pproaches h ave ev olve d toaddress  r estric ted  e lements of the fullprobl em. Sin ce  col d d ark ma tter d omina testh e m atter den s i ty ,  N -b o dymeth ods t ha t fo llow coll i sionless  cl u stering  h ave p layeda ni m portant ro lein de f i ning the o v eral l  e v oluti on ofthe co smic we b  an dthe str uct u r e of coll apsed hal o s f or medwithin i t.  Com bi ned  N -body and gas d ynamicstechniqu e s  explo reg ravi t ationa ll y coupled ev oluti on of bary o n s a nd dark mat ter . K nowl edg e ga i ned from these ‘direct’appr oache s  le d tothe d ev elopment of ‘ semi-ana l ytic’
 redshift_evolution of_baryonic features in the_galaxy power_spectrum;_iii) weak-lensing_tomography_derived from measurement_of galaxy shear_patterns and iv) the_luminosity distance–redshift relation_of_type Ia SNe. We focus our attention on theoretical issues related to the first_three_of these_tests.

The_power_spectrum of fluctuations at recombination_is calculated to high accuracy_from linear_theory$^{\ino}$, so the problem of realizing, through direct_simulation,_the evolution of_a finite patch of a particular world model from_this epoch forward is well posed._To support DES-like_observations,_one_would like to evolve_multiple regions of Hubble Length dimension_with the principal clustered matter components_— dark matter and multiple phases of_baryons (stars and cold gas in_galaxies, warm/hot gas surrounding galaxies_and in_groups/clusters) — represented by multiple_fluids. Mapping observable_signatures of_the theoretical solution_along the past light-cone of synthetic_observers in the_computational volume allows ‘clean’ mock surveys_to_be created, which_can_further_be ‘dirtied’_by real-world observational_effects_in support_of_survey data analysis.

Two fundamental barriers stand_in_the way of achieving a complete solution_to this problem. One_is_the wide dynamic range_of non-linear structures (sub-parsecs to_gigaparsecs in length, for example), and_the other_is the_complexity of astrophysical processes controlling the baryonic phases. Since DES-like surveys_probe only the higher mass portion_of the full spectrum_of cosmic_structures,_the first issue_is_not strongly_limiting. The DES will, however, require understanding_galaxy and_cluster observable signatures, so uncertainties in_the treatment of baryonic_physics_will play a central role. In_a companion paper$^{\ino}$, we outline theoretical_uncertainties associated with the large-scale_structure_channels_DES will use to test_dark energy. Here, we offer a_critique of the_computational methods that provide theoretical support for_DES_and similar surveys.

Challenges for Computational Cosmology_
=======================================

Given_the wide scale and scope of_large-scale_structure,_a number of approaches have_evolved to address restricted elements of_the full problem. Since cold dark matter dominates the_matter density, N-body_methods that follow collisionless clustering_have_played_an important role in defining the overall evolution of the_cosmic web_and the structure_of collapsed halos formed within it. Combined N-body and gas_dynamics techniques explore gravitationally coupled evolution of_baryons and dark matter. Knowledge gained from these ‘direct’ approaches led_to the development of ‘semi-analytic’